Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros.

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability.

BACKGROUND: The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA). RESULTS: We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12-10 ka ago. CONCLUSIONS: Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds.

Levels and variability of metals in soils of the province of Golestan (Iran).

Studying the concentration distribution of metals is necessary for soil pollution monitoring and maintaining environmental quality. To date, very little large-scale research has been performed to investigate metal contamination in developing countries. In this study, the content and spatial distributions of five metals (cadmium [Cd], copper [Cu], nickel [Ni], lead [Pb], and zinc [Zn]) were quantified based on 346 topsoil samples from 12 districts in the province of Golestan (northeast [NE] Iran). The concentration levels (mg/kg) of Cd, Cu, Pb, Ni and Zn varied from 0.02 to 0.36, 9.3 to 93.7, 6.8 to 44, 9.5 to 85.35, and 25 to 417.4, respectively. The average concentrations (mg/kg) obtained were as follows: Cd 0.12 ± 0.07, Cu 23.9 ± 9.07, Ni 34.88 ± 11.59, Pb 15.42 ± 5.81 and Zn 82.08 ± 30.87. Significant differences in the distribution of trace elements among districts were detected. The AzadShahr and BandarGaz districts displayed the highest metal concentrations. Greater metal values were obtained in the central, south, west, and NE areas, although Zn concentration was also high to the north of the province. Values of contamination factor and contamination degree indicated that the metal pollution level was on the order of Zn > Cd > Pb > Cu > Ni and that all of the metals belonged to a low or moderate pollution category. Ni and Cu levels derived from natural sources, whereas Cd, Pb, and Zn derived from anthropogenic activities with greater mean concentrations than reference concentrations.

Detecting hybridization between Iranian wild wolf (Canis lupus pallipes) and free-ranging domestic dog (Canis familiaris) by analysis of microsatellite markers.

The genetic threat due to hybridization with free-ranging dogs is one major concern in wolf conservation. The identification of hybrids and extent of hybridization is important in the conservation and management of wolf populations. Genetic variation was analyzed at 15 unlinked loci in 28 dogs, 28 wolves, four known hybrids, two black wolves, and one dog with abnormal traits in Iran. Pritchard's model, multivariate ordination by principal component analysis and neighbor joining clustering were used for population clustering and individual assignment. Analysis of genetic variation showed that genetic variability is high in both wolf and dog populations in Iran. Values of H(E) in dog and wolf samples ranged from 0.75-0.92 and 0.77-0.92, respectively. The results of AMOVA showed that the two groups of dog and wolf were significantly different (F(ST) = 0.05 and R(ST) = 0.36; P < 0.001). In each of the three methods, wolf and dog samples were separated into two distinct clusters. Two dark wolves were assigned to the wolf cluster. Also these models detected D32 (dog with abnormal traits) and some other samples, which were assigned to more than one cluster and could be a hybrid. This study is the beginning of a genetic study in wolf populations in Iran, and our results reveal that as in other countries, hybridization between wolves and dogs is sporadic in Iran and can be a threat to wolf populations if human perturbations increase.

Phylogeny and evolution of the genus Cervus (Cervidae, Mammalia) as revealed by complete mitochondrial genomes.

Mitochondrial DNA (mtDNA) lineages are recognized as important components of intra- and interspecific biodiversity, and allow to reveal colonization routes and phylogeographic structure of many taxa. Among these is the genus Cervus that is widely distributed across the Holarctic. We obtained sequences of complete mitochondrial genomes from 13 Cervus taxa and included them in global phylogenetic analyses of 71 Cervinae mitogenomes. The well-resolved phylogenetic trees confirmed Cervus to be monophyletic. Molecular dating based on several fossil calibration points revealed that ca. 2.6 Mya two main mitochondrial lineages of Cervus separated in Central Asia, the Western (including C. hanglu and C. elaphus) and the Eastern (comprising C. albirostris, C. canadensis and C. nippon). We also observed convergent changes in the composition of some mitochondrial genes in C. hanglu of the Western lineage and representatives of the Eastern lineage. Several subspecies of C. nippon and C. hanglu have accumulated a large portion of deleterious substitutions in their mitochondrial protein-coding genes, probably due to drift in the wake of decreasing population size. In contrast to previous studies, we found that the relic haplogroup B of C. elaphus was sister to all other red deer lineages and that the Middle-Eastern haplogroup E shared a common ancestor with the Balkan haplogroup C. Comparison of the mtDNA phylogenetic tree with a published nuclear genome tree may imply ancient introgressions of mtDNA between different Cervus species as well as from the common ancestor of South Asian deer, Rusa timorensis and R. unicolor, to the Cervus clade.

The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals.

The emergence of farming during the Neolithic transition, including the domestication of livestock, was a critical point in the evolution of human kind. The goat (Capra hircus) was one of the first domesticated ungulates. In this study, we compared the genetic diversity of domestic goats to that of the modern representatives of their wild ancestor, the bezoar, by analyzing 473 samples collected over the whole distribution range of the latter species. This partly confirms and significantly clarifies the goat domestication scenario already proposed by archaeological evidence. All of the mitochondrial DNA haplogroups found in current domestic goats have also been found in the bezoar. The geographic distribution of these haplogroups in the wild ancestor allowed the localization of the main domestication centers. We found no haplotype that could have been domesticated in the eastern half of the Iranian Plateau, nor further to the east. A signature of population expansion in bezoars of the C haplogroup suggests an early domestication center on the Central Iranian Plateau (Yazd and Kerman Provinces) and in the Southern Zagros (Fars Province), possibly corresponding to the management of wild flocks. However, the contribution of this center to the current domestic goat population is rather low (1.4%). We also found a second domestication center covering a large area in Eastern Anatolia, and possibly in Northern and Central Zagros. This last domestication center is the likely origin of almost all domestic goats today. This finding is consistent with archaeological data identifying Eastern Anatolia as an important domestication center.

Gazella arabica dareshurii: a remarkable relict population on Farur Island, Iran.

BACKGROUND: The islands in the Persian Gulf are home to several species of gazelles, i.e., Gazella bennettii, G. subgutturosa, and a new subspecies of Mountain gazelles which was discovered on Farur Island and described for the first time in 1993 as Gazella gazella dareshurii. Later, phylogenetic analyses showed that the Mountain gazelles consist of two species: G. gazella and G. arabica. As the Farur gazelles are more closely related to the Arabian forms of the Mountain gazelles, this subspecies is regarded to be G. arabica dareshurii. Until now, the origin of this subspecies has been an enigma. RESULTS: Here, we used mitochondrial cyt b, two nuclear introns (CHD2 and ZNF618), and morphological data to address this question by investigating the taxonomic position of the Farur gazelles. The results show that this population is monophyletic and split from other G. arabica populations probably 10,000 BP. CONCLUSIONS: It is a natural relict population that was trapped on the island due to the rising sea levels of the Persian Gulf after the Last Glacial Maximum. Intermittent drought and flooding are suggested to be the main factors balancing population growth in the absence of natural predators on this monsoon-influenced island. Conservation actions should focus on preserving the natural situation of the island (cease introducing mesquite tree and other invasive species, stop building new construction and roads, and caution in providing water sources and forage), and possibly introducing individuals to other islands (not inhabited by gazelles) or to fenced areas on the Iranian mainland (strictly isolated from other gazelle populations) when the population reaches the carrying capacity of the island.

Evolution and taxonomy of the wild species of the genus Ovis (Mammalia, Artiodactyla, Bovidae).

New insights for the systematic and evolution of the wild sheep are provided by molecular phylogenies inferred from Maximum parsimony, Bayesian, Maximum likelihood, and Neighbor-Joining methods. The phylogeny of the wild sheep was based on cytochrome b sequences of 290 samples representative of most of the sub-species described in the genus Ovis. The result was confirmed by a combined tree based on cytochrome b and nuclear sequences for 79 Ovis samples representative of the robust clades established with mitochondrial data. Urial and mouflon, which are either considered as a single or two separate species, form two monophyletic groups (O. orientalis and O. vignei). Their hybrids appear in one or the other group, independently from their geographic origin. The European mouflon O. musimon is clearly in the O. orientalis clade. The others species, O. dalli, O. canadensis, O. nivicola, and O. ammon are monophyletic. The results support an Asiatic origin of the genus Ovis, followed by a migration to North America through North-Eastern Asia and the Bering Strait and a diversification of the genus in Eurasia less than 3 million years ago. Our results show that the evolution of the genus Ovis is a striking example of successive speciation events occurring along the migration routes propagating from the ancestral area.

Diversification and subspecies patterning of the goitered gazelle (Gazella subgutturosa) in Iran.

Goitered gazelles, Gazella subgutturosa, exist in arid and semiarid regions of Asia from the Middle to the Far East. Although large populations were present over a vast area until recently, a decline of the population as a result of hunting, poaching, and habitat loss led to the IUCN classification of G. subgutturosa as "vulnerable." We examined genetic diversity, structure, and phylogeny of G. subgutturosa using mitochondrial cytochrome b sequences from 18 geographically distant populations in Iran. The median-joining network of cyt b haplotypes indicated that three clades of goitered gazelles can be distinguished: a Middle Eastern clade west of the Zagros Mountains (and connected to populations in Turkey and Iraq), a Central Iranian clade (with connection to Azerbaijan), and an Asiatic clade in northeastern Iran (with connection to Turkmenistan, Uzbekistan, and other Asian countries as far as northeastern China and Mongolia). Based on our results, we argue that Iran is the center of diversification of goitered gazelles, due to the presence of large mountain ranges and deserts that lead to the separation of populations. In accordance with previous morphological studies, we identified the Asiatic clade as the subspecies G. s. yarkandensis, and the other two clades as the nominate form G. s. subgutturosa. The new genetic information for goitered gazelles in Iran provides the basis for future national conservation programs of this species.

Sheep genome functional annotation reveals proximal regulatory elements contributed to the evolution of modern breeds.

Domestication fundamentally reshaped animal morphology, physiology and behaviour, offering the opportunity to investigate the molecular processes driving evolutionary change. Here we assess sheep domestication and artificial selection by comparing genome sequence from 43 modern breeds (Ovis aries) and their Asian mouflon ancestor (O. orientalis) to identify selection sweeps. Next, we provide a comparative functional annotation of the sheep genome, validated using experimental ChIP-Seq of sheep tissue. Using these annotations, we evaluate the impact of selection and domestication on regulatory sequences and find that sweeps are significantly enriched for protein coding genes, proximal regulatory elements of genes and genome features associated with active transcription. Finally, we find individual sites displaying strong allele frequency divergence are enriched for the same regulatory features. Our data demonstrate that remodelling of gene expression is likely to have been one of the evolutionary forces that drove phenotypic diversification of this common livestock species.

Convergent genomic signatures of domestication in sheep and goats.

The evolutionary basis of domestication has been a longstanding question and its genetic architecture is becoming more tractable as more domestic species become genome-enabled. Before becoming established worldwide, sheep and goats were domesticated in the fertile crescent 10,500 years before present (YBP) where their wild relatives remain. Here we sequence the genomes of wild Asiatic mouflon and Bezoar ibex in the sheep and goat domestication center and compare their genomes with that of domestics from local, traditional, and improved breeds. Among the genomic regions carrying selective sweeps differentiating domestic breeds from wild populations, which are associated among others to genes involved in nervous system, immunity and productivity traits, 20 are common to Capra and Ovis. The patterns of selection vary between species, suggesting that while common targets of selection related to domestication and improvement exist, different solutions have arisen to achieve similar phenotypic end-points within these closely related livestock species.

Contact Zone of Asian and European Wild Boar at North West of Iran.

Wild boar (Sus scrofa) are widely distributed throughout the Old World. Most studies have focused on Europe and East Asia with the genetic diversity of West Asia being less well studied. In particular, the genetic variability and genetic structure of the Iranian populations are not yet known; gaps which prevent scientists from resolving the genetic relationships of the Eurasian wild boar. This paper is the first attempt to provide information about genetic relationships among modern Iranian populations of the Eurasian wild boar (S. scrofa) by sequencing 572 bp of the mitochondrial (mt) DNA control region. As a result of this investigation, it was discovered that Iran contains not only Middle Eastern haplotypes, but also shares haplotypes with Europe and East Asia. The Italian clade, which is endemic in Italy, is not identified in Iran, while all other clades, including Asiatic, European, Near East 1, and Near East 2 are found based on the phylogenetic tree and median-joining network. The results of this study illustrate that north west of Iran (specifically Southwest Caspian Sea) is the contact zone between the Asian (Near Eastern and Far Eastern), and the European clades. In light of the fact that the domestication of pigs occurs in Anatolia, this finding is important.

Is black coat color in wolves of Iran an evidence of admixed ancestry with dogs?

Melanism is not considered a typical characteristic in wolves of Iran and dark wolves are believed to have originated from crossbreeding with dogs. Such hybrid individuals can be identified with the combined use of genetic and morphological markers. We analyzed two black wolves using a 544 base pairs (bp) fragment of the mtDNA control region and 15 microsatellite loci in comparison with 28 dogs, 28 wolves, and four known hybrids. The artificial neural networks (ANNs) method was applied to microsatellite data to separate genetically differentiated samples of wolves, dogs, and hybrids, and to determine the correct class for the black specimens. Individual assignments based on ANNs showed that black samples were genetically closer to wolves. Also, in the neighbor-joining network of mtDNA haplotypes, wolves and dogs were separated, with the dark specimens located in the wolf branch as two separate haplotypes. Furthermore, we compared 20 craniometrical characters of the two black individuals with 14 other wolves. The results showed that craniometrical measures of the two black wolves fall within the range of wolf skulls. We found no trace of recent hybridization with free-ranging dogs in the two black wolves. Dark coat color might be the result of a natural combination of alleles in the coat-color-determining gene, mutation in the K locus due to past hybridization with free-ranging dogs, or the effect of ecological factors and adaption to habitat conditions.

Reassessment of the taxonomic position of Iranocypris typhlops Bruun & Kaiser, 1944 (Actinopterygii, Cyprinidae).

The Iranian cave barb (Iranocypris typhlops Bruun & Kaiser, 1944) is a rare and endemic species of the family Cyprinidae known from a single locality in the Zagros Mountains, western Iran. This species is "Vulnerable" according to the IUCN Red List and is one of the top four threatened freshwater fish species in Iran. Yet, the taxonomic position of I. typhlops is uncertain. We examined phylogenetic relationships of this species with other species of the family Cyprinidae based on the mitochondrial cytochrome b gene. Our results show that I. typhlops is monophyletic and is sister taxon of a cluster formed by Garra rufa (Heckel, 1843) and Garra barreimiae (Fowler & Steinitz, 1956) within a clade that includes other species of the genus Garra. Based on previous molecular and morphological studies, as well as our new results, we recommend that I. typhlops should be transferred to the genus Garra Hamilton, 1822.

Large-scale mitochondrial DNA analysis of the domestic goat reveals six haplogroups with high diversity.

BACKGROUND: From the beginning of domestication, the transportation of domestic animals resulted in genetic and demographic processes that explain their present distribution and genetic structure. Thus studying the present genetic diversity helps to better understand the history of domestic species. METHODOLOGY/PRINCIPAL FINDINGS: The genetic diversity of domestic goats has been characterized with 2430 individuals from all over the old world, including 946 new individuals from regions poorly studied until now (mainly the Fertile Crescent). These individuals represented 1540 haplotypes for the HVI segment of the mitochondrial DNA (mtDNA) control region. This large-scale study allowed the establishment of a clear nomenclature of the goat maternal haplogroups. Only five of the six previously defined groups of haplotypes were divergent enough to be considered as different haplogroups. Moreover a new mitochondrial group has been localized around the Fertile Crescent. All groups showed very high haplotype diversity. Most of this diversity was distributed among groups and within geographic regions. The weak geographic structure may result from the worldwide distribution of the dominant A haplogroup (more than 90% of the individuals). The large-scale distribution of other haplogroups (except one), may be related to human migration. The recent fragmentation of local goat populations into discrete breeds is not detectable with mitochondrial markers. The estimation of demographic parameters from mismatch analyses showed that all groups had a recent demographic expansion corresponding roughly to the period when domestication took place. But even with a large data set it remains difficult to give relative dates of expansion for different haplogroups because of large confidence intervals. CONCLUSIONS/SIGNIFICANCE: We propose standard criteria for the definition of the different haplogroups based on the result of mismatch analysis and on the use of sequences of reference. Such a method could be also applied for clarifying the nomenclature of mitochondrial haplogroups in other domestic species.