Characteristics of the Genetic Structure of Snow Sheep (<Emphasis Type="Italic">Ovis nivicola lydekkeri</Emphasis>) of the Verkhoyansk Mountain Chain

Genetic characteristics of the allele pool of four groups of the Yakut snow sheep subspecies (Ovis nivicola lydekkeri) inhabiting various parts of the Verkhoyansk Mountain Range such as Kharaulakh Ridge, Orulgan Ridge, ridges of the Central Verkhoyansk, and Suntar-Khayata Ridge is presented. Fragment analysis using 17 microsatellite loci was carried out using the ABI 3131xl genetic analyzer. Significant heterozygote deficiency was detected in all investigated snow sheep populations. Differentiation of the studied groups in accordance to their geographical origin was revealed.     

Genome‐wide SNP analysis unveils genetic structure and phylogeographic history of snow sheep (Ovis nivicola) populations inhabiting the Verkhoyansk Mountains and Momsky Ridge (northeastern Siberia)

Insights into the genetic characteristics of a species provide important information for wildlife conservation programs. Here, we used the OvineSNP50 BeadChip developed for domestic sheep to examine population structure and evaluate genetic diversity of snow sheep (Ovis nivicola) inhabiting Verkhoyansk Range and Momsky Ridge. A total of 1,121 polymorphic SNPs were used to test 80 specimens representing five populations, including four populations of the Verkhoyansk Mountain chain: Kharaulakh Ridge–Tiksi Bay (TIK, n = 22), Orulgan Ridge (ORU, n = 22), the central part of Verkhoyansk Range (VER, n = 15), Suntar‐Khayata Ridge (SKH, n = 13), and Momsky Ridge (MOM, n = 8). We showed that the studied populations were genetically structured according to a geographic pattern. Pairwise F_ST values ranged from 0.044 to 0.205. Admixture analysis identified K = 2 as the most likely number of ancestral populations. A Neighbor‐Net tree showed that TIK was an isolated group related to the main network through ORU. TreeMix analysis revealed that TIK and MOM originated from two different ancestral populations and detected gene flow from MOM to ORU. This was supported by the f3 statistic, which showed that ORU is an admixed population with TIK and MOM/SKH heritage. Genetic diversity in the studied groups was increasing southward. Minimum values of observed (Ho) and expected (He) heterozygosity and allelic richness (Ar) were observed in the most northern population—TIK, and maximum values were observed in the most southern population—SKH. Thus, our results revealed clear genetic structure in the studied populations of snow sheep and showed that TIK has a different origin from MOM, SKH, and VER even though they are conventionally considered a single subspecies known as Yakut snow sheep (Ovis nivicola lydekkeri). Most likely, TIK was an isolated group during the Late Pleistocene glaciations of Verkhoyansk Range.     

Phylogenetic analysis of snow sheep (Ovis nivicola) and closely related taxa

Based on mitochondrial cytochrome b gene sequence analysis, the history of true sheep (Ovis) began approximately 3.12 million years ago (MYA). The evolution of Ovis resulted in three generally accepted genetic groups: Argaliforms, Moufloniforms, and Pachyceriforms. The Pachyceriforms of the subgenus Pachyceros comprise the thin-horn sheep Ovis nivicola (snow sheep), Ovis dalli (Dall and Stone sheep), and Ovis canadensis (Rocky Mountain and desert bighorn). North America wild sheep (O. canadensis and O. dalli) evolved separately from Eurasian wild sheep and diverged from each other about 1.41 MYA. Ancestral stock that gave rise to snow sheep, Moufloniforms, and Argaliforms occurred 2.3 MYA, which then gave rise to two different extant lines of snow sheep that diverged from each other about 1.96 MYA. The more recent nivicola line is genetically closer to the North American wild sheep and may represent a close association during the refugium when Alaska and Siberia were connected by the Bering land bridge. The earlier period of evolution of the Pachyceriforms suggests they may have first evolved in Eurasia, the oldest ancestor then giving rise to North American wild sheep, and that a canadensis-like ancestor most likely gave rise to nivicola. Cytogenetic analysis further validates that the standard diploid number for modern nivicola is 52.     

Whole genome SNP scanning of snow sheep (<Emphasis Type="Italic">Ovis nivicola</Emphasis>)

This is the first report performing the whole genome SNP scanning of snow sheep (Ovis nivicola). Samples of snow sheep (n = 18) collected in six different regions of the Republic of Sakha (Yakutia) from 64° to 71° N. For SNP genotyping, we applied Ovine 50K SNP BeadChip (Illumina, United States), designed for domestic sheep. The total number of genotyped SNPs (call rate 90%) was 47796 (88.1% of total SNPs), wherein 1006 SNPs were polymorphic (2.1%). Principal component analysis (PCA) showed the clear differentiation within the species O. nivicola: studied individuals were distributed among five distinct arrays corresponding to the geographical locations of sampling points. Our results demonstrate that the DNA chip designed for domestic sheep can be successfully used to study the allele pool and the genetic structure of snow sheep populations.     

Genetic distinctiveness of endangered dwarf blue sheep (Pseudois nayaur schaeferi): evidence from mitochondrial control region and Y-linked ZFY intron sequences

To elucidate the controversial systematic relationship of blue sheep (Pseudois nayaur) and endangered dwarf blue sheep, we sequenced part of the mtDNA control region and Y-linked ZFY intron, and carried out phylogenetic analyses. Mitochondrial results revealed that the dwarf blue sheep is a strongly supported monophyletic group, with an average of 12.21% sequence divergence from the blue sheep. This is the first genetic evidence for the distinctness of the dwarf blue sheep. ZFY intron results showed an average of 0.51% sequence divergence, and one shared haplotype between the dwarf blue sheep and blue sheep. By analyzing an expanded data set that incorporated ZFY intron sequences of two additional Ovis (sheep) species-O. nivicola and O. ammon--we demonstrated that ZFY intron provides good resolution at the species and genus levels. The ZFY intron sequence divergence between dwarf blue sheep and blue sheep was comparable to that within the two Ovis species. Moreover, we found intraspecific sequence variation in ZFY intron for all three species examined. We propose that dwarf blue sheep be designated as a subspecies of blue sheep, P. n. schaeferi.     

A genome-wide set of SNPs detects population substructure and long range linkage disequilibrium in wild sheep

The development of genomic resources for wild species is still in its infancy. However, cross-species utilization of technologies developed for their domestic counterparts has the potential to unlock the genomes of organisms that currently lack genomic resources. Here, we apply the OvineSNP50 BeadChip, developed for domestic sheep, to two related wild ungulate species: the bighorn sheep (Ovis canadensis) and the thinhorn sheep (Ovis dalli). Over 95% of the domestic sheep markers were successfully genotyped in a sample of fifty-two bighorn sheep while over 90% were genotyped in two thinhorn sheep. Pooling the results from both species identified 868 single-nucleotide polymorphisms (SNPs), 570 were detected in bighorn sheep, while 330 SNPs were identified in thinhorn sheep. The total panel of SNPs was able to discriminate between the two species, assign population of origin for bighorn sheep and detect known relationship classes within one population of bighorn sheep. Using an informative subset of these SNPs (n=308), we examined the extent of genome-wide linkage disequilibrium (LD) within one population of bighorn sheep and found that high levels of LD persist over 4 Mb.     

Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale,Merino and Creole sheep

The aim of this study was to investigate the genetic diversity within and among three breeds of sheep: Corriedale, Merino and Creole. Sheep from the three breeds (Merino n = 110, Corriedale n = 108 and Creole n = 10) were genotyped using the Illumina Ovine SNP50 beadchip^®. Genetic diversity was evaluated by comparing the minor allele frequency (MAF) among breeds. Population structure and genetic differentiation were assessed using STRUCTURE software, principal component analysis (PCA) and fixation index (F_ST). Fixed markers (MAF = 0) that were different among breeds were identified as specific breed markers. Using a subset of 18,181 single nucleotide polymorphisms (SNPs), PCA and STUCTURE analysis were able to explain population stratification within breeds. Merino and Corriedale divergent lines showed high levels of polymorphism (89.4% and 86% of polymorphic SNPs, respectively) and moderate genetic differentiation (F_ST = 0.08) between them. In contrast, Creole had only 69% polymorphic SNPs and showed greater genetic differentiation from the other two breeds (F_ST = 0.17 for both breeds). Hence, a subset of molecular markers present in the OvineSNP50 is informative enough for breed assignment and population structure analysis of commercial and Creole breeds.     

Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ～36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.     

Genetic diversity of three Chinese native sheep breeds

The sheep (Ovis aries L.) has been an important farm animal species since its domestication. A wide array of indigenous sheep breeds with abundant phenotypic diversity exists for domestication and selection. Therefore, assessing the genetic diversity of a local sheep resource using a multi-molecular system is helpful for maintaining and conserving those breeds. This study aimed to investigate the genetic diversity of three native Chinese sheep breeds (Tibetan sheep, Sishui Fur sheep, and Small-tailed Han sheep) using 15 microsatellite markers and the second exon of the DRA gene. In regards to the microsatellites, on average, 19 alleles per loci were observed among all individuals. Across loci, the HO within the population was 0.652 ± 0.022 in Tibetan sheep, 0.603 ± 0.023 in Small-tailed Han sheep and 0.635 ± 0.022 in SFS, and for most populations, the H _E and H _O were inconsistent. In addition, affluent private alleles within the breed indicated that the breeds have different domestication histories or sites. In regards to the 2 exon of the DRA gene, three haplotypes were constructed by seven single-nucleotide polymorphisms (SNPs), which were identified in the second DRA exon and inferred the potential for phenotypic variety in these Chinese native sheep. In summary, the current study reveals the importance of implementing effective conservation strategies for these three native Chinese sheep.     

Population genetic structure of North American thinhorn sheep (Ovis dalli)

The thinhorn sheep (Ovis dalli ssp.) provides a rare example of a North American large mammal that occupies most of its native range and maintains close to ancestral population size. There are currently two recognized subspecies, Dall's sheep (O. d. dalli) and Stone's sheep (O. d. stonei), the validity of which remains uncertain. We investigated the spatial genetic structure of thinhorn sheep populations representing both subspecies by genotyping individuals (n = 919) from across the species range at 12 variable microsatellite loci. We found high levels of genetic diversity within (HE = 0.722) and significant genetic structure among the 24 sampled areas (FST = 0.160). Genetic distance measures and Bayesian clustering analyses revealed the presence of at least eight subpopulations that are delineated by mountain range topology. A strong overall pattern of isolation-by-distance is evident across the sampling range (r = 0.75, P < 0.001) suggesting limited dispersal and extensive philopatry. Partial Mantel tests of this relationship showed mountain range distinctions represent significant barriers to gene flow (P = 0.0001), supporting the Bayesian analyses. Genetic structure was more strongly pronounced in southern Yukon and Alaska than elsewhere. We also show evidence for genetic differences between the two currently recognized thinhorn subspecies.     

利用AIMs分析亚洲9个绵羊群体的群体结构

祖先信息标记(ancestry informative makers,AIMs)可用来分析群体的遗传结构.本研究利用Illumina Ovine SNP50芯片上的SNP位点,在云南乌骨绵羊和其他8个亚洲绵羊群体中以Rosenberg等定义的Informativeness统计量为筛选方法,选取Informativeness值最高的前20、50、100、500个SNP位点与相应数目的随机SNP位点分别用来推断群体的遗传结构.通过主成分分析和用fast STRUCURE推断祖先成分的方法,评价AIMs在推断亚洲绵羊群体遗传结构中的作用.研究显示,利用筛选到的高信息含量标记AIMs,可减少群体结构研究中需要的SNP位点数目.前50个AIMs可有效地将绵羊群体分为4个大类,这与利用全基因组SNPs分析得到的群体结构是一致的,即乌骨绵羊群体blackbone单独为一类、西藏群体changthangi和tibetan归为一类,孟加拉的banglandeshi、banglandeshi Garole群体和印度的Indian Garole群体归为一个类群,其余三个群体(印度尼西亚sumatran、garut群体和印度的deccani群体)近似归为一类.这4大类群在AIMs上存在显著的分化,利用这些位点信息可以为研究群体特征和进化关系提供线索.     

Evaluating wildlife translocations using genomics: A bighorn sheep case study

Wildlife restoration often involves translocation efforts to reintroduce species and supplement small, fragmented populations. We examined the genomic consequences of bighorn sheep (Ovis canadensis) translocations and population isolation to enhance understanding of evolutionary processes that affect population genetics and inform future restoration strategies. We conducted a population genomic analysis of 511 bighorn sheep from 17 areas, including native and reintroduced populations that received 0–10 translocations. Using the Illumina High Density Ovine array, we generated datasets of 6,155 to 33,289 single nucleotide polymorphisms and completed clustering, population tree, and kinship analyses. Our analyses determined that natural gene flow did not occur between most populations, including two pairs of native herds that had past connectivity. We synthesized genomic evidence across analyses to evaluate 24 different translocation events and detected eight successful reintroductions (i.e., lack of signal for recolonization from nearby populations) and five successful augmentations (i.e., reproductive success of translocated individuals) based on genetic similarity with the source populations. A single native population founded six of the reintroduced herds, suggesting that environmental conditions did not need to match for populations to persist following reintroduction. Augmentations consisting of 18–57 animals including males and females succeeded, whereas augmentations of two males did not result in a detectable genetic signature. Our results provide insight on genomic distinctiveness of native and reintroduced herds, information on the relative success of reintroduction and augmentation efforts and their associated attributes, and guidance to enhance genetic contribution of augmentations and reintroductions to aid in bighorn sheep restoration.     

MLH1 Region Polymorphisms Show a Significant Association with CpG Island Shore Methylation in a Large Cohort of Healthy Individuals

Single nucleotide polymorphisms (SNPs) are the most common form of genetic variation. We previously demonstrated that SNPs (rs1800734, rs749072, and rs13098279) in the MLH1 gene region are associated with MLH1 promoter island methylation, loss of MLH1 protein expression, and microsatellite instability (MSI) in colorectal cancer (CRC) patients. Recent studies have identified less CpG-dense “shore” regions flanking many CpG islands. These shores often exhibit distinct methylation profiles between different tissues and matched normal versus tumor cells of patients. To date, most epigenetic studies have focused on somatic methylation events occurring within solid tumors; less is known of the contributions of peripheral blood cell (PBC) methylation to processes such as aging and tumorigenesis. To address whether MLH1 methylation in PBCs is correlated with tumorigenesis we utilized the Illumina 450 K microarrays to measure methylation in PBC DNA of 846 healthy controls and 252 CRC patients from Ontario, Canada. Analysis of a region of chromosome 3p21 spanning the MLH1 locus in healthy controls revealed that a CpG island shore 1 kb upstream of the MLH1 gene exhibits different methylation profiles when stratified by SNP genotypes (rs1800734, rs749072, and rs13098279). Individuals with wild-type genotypes incur significantly higher PBC shore methylation than heterozygous or homozygous variant carriers (p<1.1×10⁻⁶; ANOVA). This trend is also seen in CRC cases (p<0.096; ANOVA). Shore methylation also decreases significantly with increasing age in cases and controls. This is the first study of its kind to integrate PBC methylation at a CpG island shore with SNP genotype status in CRC cases and controls. These results indicate that CpG island shore methylation in PBCs may be influenced by genotype as well as the normal aging process.     

SNP discovery and genotyping using restriction‐site‐associated DNA sequencing in chickens

Single nucleotide polymorphisms (SNPs) are essential to the understanding of population genetic variation and diversity. Here, we performed restriction‐site‐associated DNA sequencing (RAD‐seq) on 72 individuals from 13 Chinese indigenous and three introduced chicken breeds. A total of 620 million reads were obtained using an Illumina Hiseq2000 sequencer. An average of 75 587 SNPs were identified from each individual. Further filtering strictly validated 28 895 SNPs candidates for all populations. When compared with the NCBI dbSNP (chicken_9031), 15 404 SNPs were new discoveries. In this study, RAD‐seq was performed for the first time on chickens, implicating the remarkable effectiveness and potential applications on genetic analysis and breeding technique for whole‐genome selection in chicken and other agricultural animals.     

Genome-Wide Genetic Diversity and Differentially Selected Regions among Suffolk,Rambouillet, Columbia,Polypay, and Targhee Sheep

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F _ST, significant genetic differentiation appeared between Suffolk and Rambouillet (F _ST = 0.1621), while Rambouillet and Targhee had the closest relationship (F _ST = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.     

Evaluating sample size to estimate genetic management metrics in the genomics era

Inbreeding and relationship metrics among and within populations are useful measures for genetic management of wild populations, but accuracy and precision of estimates can be influenced by the number of individual genotypes analysed. Biologists are confronted with varied advice regarding the sample size necessary for reliable estimates when using genomic tools. We developed a simulation framework to identify the optimal sample size for three widely used metrics to enable quantification of expected variance and relative bias of estimates and a comparison of results among populations. We applied this approach to analyse empirical genomic data for 30 individuals from each of four different free‐ranging Rocky Mountain bighorn sheep (Ovis canadensis canadensis) populations in Montana and Wyoming, USA, through cross‐species application of an Ovine array and analysis of approximately 14,000 single nucleotide polymorphisms (SNPs) after filtering. We examined intra‐ and interpopulation relationships using kinship and identity by state metrics, as well as F_ST between populations. By evaluating our simulation results, we concluded that a sample size of 25 was adequate for assessing these metrics using the Ovine array to genotype Rocky Mountain bighorn sheep herds. However, we conclude that a universal sample size rule may not be able to sufficiently address the complexities that impact genomic kinship and inbreeding estimates. Thus, we recommend that a pilot study and sample size simulation using R code we developed that includes empirical genotypes from a subset of populations of interest would be an effective approach to ensure rigour in estimating genomic kinship and population differentiation.     

Molecular analysis of wild and domestic sheep questions current nomenclature and provides evidence for domestication from two different subspecies

Complete mitochondrial DNA (mtDNA) control regions (CR) were sequenced and analysed in order to investigate wild sheep taxonomy and the origin of domestic sheep (Ovis aries). The dataset for phylogenetic analyses includes 63 unique CR sequences from wild sheep of the mouflon (O. musimon, O. orientalis), urial (O. vignei), argali (O. ammon) and bighorn (O. canadensis) groups, and from domestic sheep of Asia, Europe and New Zealand. Domestic sheep occurred in two clearly separated branches with mouflon (O. musimon) mixed into one of the domestic sheep clusters. Genetic distances and molecular datings based on O. canadensis CR and mtDNA protein-coding sequences provide strong evidence for domestications from two mouflon subspecies. Other wild sheep sequences are in two additional well-separated branches. Ovis ammon collium and O. ammon nigrimontana are joined with a specimen from the transkaspian Ust-Urt plateau currently named O. vignei arkal. Ovis ammon ammon, O. ammon darwini and O. vignei bochariensis represent a separate clade and the earliest divergence from the mouflon group. Therefore, O. musimon, O. vignei bochariensis and Ust-Urt sheep are not members of a 'moufloniform' or O. orientalis species, but belong to different clades. Furthermore, Ust-Urt sheep could be a hybrid population or an O. ammon subspecies closely related to O. ammon nigrimontana.     

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.     

Genome‐wide set of SNPs reveals evidence for two glacial refugia and admixture from postglacial recolonization in an alpine ungulate

Past glaciation events have played a major role in shaping the genetic diversity and distribution of wild sheep in North America. The advancement of glaciers can isolate populations in ice‐free refugia, where they can survive until the recession of ice sheets. The major Beringian refugium is thought to have held thinhorn sheep (Ovis dalli) populations during times of glacial advance. While isolation in the major refugium can account for much of the genetic and morphological diversity seen in extant thinhorn sheep populations, mounting evidence suggests the persistence of populations in smaller minor refugia. We investigated the refugial origins of thinhorn sheep using ~10 000 SNPs obtained via a cross‐species application of the domestic sheep ovine HD BeadChip to genotype 52 thinhorn sheep and five bighorn sheep (O. canadensis) samples. Phylogenetic inference revealed a distinct lineage of thinhorn sheep inhabiting British Columbia, which is consistent with the survival of a group of thinhorn sheep in a minor refugium separate from the Beringian refugium. Isolation in separate glacial refugia probably mediated the evolution of the two thinhorn sheep subspecies, the white Dall's sheep (O. d. dalli), which persisted in Beringia, and the dark Stone's sheep (O. d. stonei), which utilized the minor refugium. We also found the first genetic evidence for admixture between sheep from different glacial refugia in south‐central Yukon as a consequence of post glacial expansion and recolonization. These results show that glaciation events can have a major role in the evolution of species inhabiting previously glaciated habitats and the need to look beyond established refugia when examining the evolutionary history of such species.     

Mitochondrial DNA Diversity of Modern,Ancient and Wild Sheep (Ovis gmelinii anatolica) from Turkey: New Insights on the Evolutionary History of Sheep

In the present study, to contribute to the understanding of the evolutionary history of sheep, the mitochondrial (mt) DNA polymorphisms occurring in modern Turkish native domestic (n = 628), modern wild (Ovis gmelinii anatolica) (n = 30) and ancient domestic sheep from Oylum Höyük in Kilis (n = 33) were examined comparatively with the accumulated data in the literature. The lengths (75 bp/76 bp) of the second and subsequent repeat units of the mtDNA control region (CR) sequences differentiated the five haplogroups (HPGs) observed in the domestic sheep into two genetic clusters as was already implied by other mtDNA markers: the first cluster being composed of HPGs A, B, D and the second cluster harboring HPGs C, E.To manifest genetic relatedness between wild Ovis gmelinii and domestic sheep haplogroups, their partial cytochrome B sequences were examined together on a median-joining network. The two parallel but wider aforementioned clusters were observed also on the network of Ovis gmelenii individuals, within which domestic haplogroups were embedded. The Ovis gmelinii wilds of the present day appeared to be distributed on two partially overlapping geographic areas parallel to the genetic clusters that they belong to (the first cluster being in the western part of the overall distribution). Thus, the analyses suggested that the domestic sheep may be the products of two maternally distinct ancestral Ovis gmelinii populations.Furthermore, Ovis gmelinii anatolica individuals exhibited a haplotype of HPG A (n = 22) and another haplotype (n = 8) from the second cluster which was not observed among the modern domestic sheep. HPG E, with the newly observed members (n = 11), showed signs of expansion. Studies of ancient and modern mtDNA suggest that HPG C frequency increased in the Southeast Anatolia from 6% to 22% some time after the beginning of the Hellenistic period, 500 years Before Common Era (BCE).