Research Article | Open Access
Lorraine Pariset, Marco Mariotti, Maria Gargani, Stephane Joost, Riccardo Negrini, Trinidad Perez, Michael Bruford, Paolo Ajmone Marsan, Alessio Valentini, "Genetic Diversity of Sheep Breeds from Albania, Greece, and Italy Assessed by Mitochondrial DNA and Nuclear Polymorphisms (SNPs)", The Scientific World Journal, vol. 11, Article ID 186342, 19 pages, 2011. https://doi.org/10.1100/2011/186342
Genetic Diversity of Sheep Breeds from Albania, Greece, and Italy Assessed by Mitochondrial DNA and Nuclear Polymorphisms (SNPs)
ABSTRACT
We employed mtDNA and nuclear SNPs to investigate the genetic diversity of sheep breeds of three countries of the Mediterranean basin: Albania, Greece, and Italy. In total, 154 unique mtDNA haplotypes were detected by means of D-loop sequence analysis. The major nucleotide diversity was observed in Albania. We identified haplogroups, A, B, and C in Albanian and Greek samples, while Italian individuals clustered in groups A and B. In general, the data show a pattern reflecting old migrations that occurred in postneolithic and historical times. PCA analysis on SNP data differentiated breeds with good correspondence to geographical locations. This could reflect geographical isolation, selection operated by local sheep farmers, and different flock management and breed admixture that occurred in the last centuries.
1. INTRODUCTION
The earliest archaeozoological evidence of domestic sheep comes from a restricted area of south-western Asia: modern Iran, Turkey, and Cyprus [1]. A pioneering genetic study examining the karyotypes of the various species of extant wild sheep [2, 3] showed that domestic sheep derive from the Asiatic mouflon (Ovis orientalis) of Anatolia, western Iran, and southwest Iran. Afterwards, a probable migration of the Neolithic farmers occurred out of the Near East and across Europe following two main routes, through the continental heartland up the Danube valley or along the Mediterranean coast [4, 5] crossing the sea to the major islands. Archaeological data and radiocarbon dates on seeds or bones provide support for an earlier arrival in Western Europe via the Mediterranean route rather than the âDanubianâ route [6].
Both archaeozoological evidence and genetic evidence indicate that the domestication of wild sheep occurred 8000â9000 years ago. The first appearance of the remains of domestic sheep in the western part of Mediterranean Europe, dating approximately 5400 BC, is believed to reflect a rapid spread by sea [7, 8]. The Mediterranean Sea also had a key role in the history of livestock in postneolithic times, when peoples like Phoenicians, Greeks, Romans, and Berbers probably introduced new species of animals and new breeds of livestock into southwest Europe by sea. Some settlers may have improved local livestock by importing stock from overseas [8], which explained the unexpectedly high within-breed diversity in domestic goats [9, 10], the differential cattle migration along the Mediterranean coast [11], and the close genetic relationship between Tuscan and Near Eastern cattle breeds [12]. The role of the Mediterranean Sea as a natural corridor connecting the South West Europe to the Near East and North Africa is particularly plausible for domestic sheep and goats that were adaptable to various environments and easy to transport due to their size [8]. Subsequently, sheep breeds developed after selective breeding for desirable traits (wool, milk, and meat production) and environmental tolerance. Since domestication, sheep have established a wide geographic range due to their adaptability to poor diets and extreme climatic conditions as well as their manageable size.
The genetic history of sheep has been investigated using three major sources of genomic variation: autosomes, Y chromosome, and mitochondrial genome. Analysis of the nonrecombining region of the Y chromosome has revealed patterns of male-mediated introgression during breed development [13, 14]. Recent surveys have tested collections of animals from southern and northern Europe [15] or Europe and the Middle East [16] using microsatellites and enabled the analysis of genetic partitioning at a continental scale. Interestingly, southern European breeds displayed increased genetic diversity and decreased genetic differentiation compared to their northern European counterparts. This is consistent with the expectation that genetic diversity will be maintained high in populations close to the centre of domestication but decreases with increasing geographic distances. Kijas et al. [17] used a SNP panel to analyse sheep nuclear genome, providing the indication that breeds cluster into large groups based on geographic origin and that SNPs can successfully identify population substructures within individual breeds.
A recent study on retrovirus integrations [18] has provided additional information on the introduction of sheep into Europe, indicating an early arrival of the primitive sheep populations (European mouflons, North-Atlantic Island breeds) and a subsequent advent of wool producing sheep.
However, most of the information about history and domestication of the species have been gathered using mtDNA. The existence of multiple mtDNA lineages and their admixture within breeds [8, 19â22] could be due to multiple domestication events and subsequent human selection or introgression by domestic and wild species.
Mitochondrial DNA analyses in sheep identified an increasing number of maternal lineages: two [23â25], three [20, 26], and then five [22]. The main haplogroups A and B are both found in Asia, while B dominates in Europe. Haplogroup C has been found in Portugal, Turkey, the Caucasus, and China [7]. Haplogroup D, present in Rumanian Karachai and Caucasian animals, is possibly related to the haplogroup A. In contrast to taurine cattle, the sheep haplogroups hardly correlate with geography.
Because of their mode of inheritance, mitochondrial markers are more likely to lead to biased estimates of species phylogeny [27]. Combining nuclear and mitochondrial markers may help in avoiding this problem. Nuclear genome evolves five-to-ten times slower than mtDNA; it is contributed by both parents and its variability is less affected by demographic forces such as bottleneck. Therefore, nuclear markers can detect more recent genetic events that influence the extant divergence of domestic breeds. Several studies have demonstrated that the combination of nuclear and mtDNA markers can increase the information obtained [27â30]. The use of both markers might provide a more accurate and comprehensive understanding of a speciesâ history [31]. SNP markers could help in understanding the recent evolutionary history of domestic animals [10, 32].
We aimed at investigating the geographic distribution of the genetic diversity of sheep breeds in Albania, Greece, and Italy and to gather information on the migration history of the species. To accomplish that, we employed sequence data from the mitochondrial D-loop and 27 nuclear loci (SNPs).
2. MATERIALS AND METHODS
2.1. Sampling and DNA Extraction
We focused on sheep breeds of Albania, Greece, and Italy. Samples of the European mouflon were also included. About twenty unrelated samples per breed were selected. Three animals per flock from 11 farms spread over the traditional rearing area were sampled. A total of 313 animals from 18 sheep breeds were analyzed. Breeds, acronyms used, and country of origin of each breed are reported in Table 1. Part of the samples were obtained from a previous project (Econogene, http://www.econogene.eu/). Blood samples were collected in EDTA tubes and frozen at â20°C until extraction. Genomic DNA was isolated using standard procedures, checked for DNA quality on agarose gel and quantified using a DTX microplate reader (Beckman Coulter) after staining with PicoGreen (Invitrogen).
|
2.2. Amplification and Sequencing of the Mitochondrial D-Loop
To amplify the partial D-loop of 721âbp, primers described by Tapio et al. [7] were used from 15,541 to 16,261 of the complete sequence described by Hiendleder et al. [33] available in GenBank (NC_001941.1). Polymerase chain reaction (PCR) was performed in a total volume of 50âÎŒL containing 20âng of genomic DNA, 40âpMol of each primer (Sigma-Aldrich), 200âÎŒM dNTPs, 5X PCR buffer, and 5 units of Taq DNA polymerase (Promega) on a PCR Thermo Cycler (MJ Research). A 5 minutes denaturation step at 95°C was followed by 14 cycles of denaturation at 95°C for 30âsec, annealing for 30âsec starting at 62°C and decreasing 0.5°C per cycle and extension at 72°C for 120âsec, then by 20 cycles of denaturation at 94°C for 30âsec, annealing at 55°C for 30âsec and extension at 72°C for 120âsec; the final extension step was carried out at 72°C for 5 minutes. PCR products were purified through ExoSap-IT (USB Corporation) to remove residual primers and dNTPs and used as templates for forward and reverse sequencing reactions.
Sequencing was performed using the primers described by Tapio et al. [7] with a CEQ 8800 sequencer using DTCS QuickStart Kit and purifying with Agencourt CleanSEQ 96 (Beckman Coulter), according to the manufacturerâs instructions. After the optimization of the sequencing protocol, sequencing was outsourced to Macrogen (http://www.macrogen.com/). The sequences of D-loop were submitted to GenBank (accession numbers: JN184789âJN184999).
2.3. Mitochondrial Sequence Analysis
A fragment of 435âbp, running from 15,541 to 16,261âbp (NC_0019041.1), was selected excluding a central region rich in tandem repeats (from 15,644 to 15,932âbp). mtDNA variations were identified on a total of 313 sequences of 18 breeds analyzed and aligned with BioEdit software [34]. DnaSP 5.00 software [35] was used to calculate haplotype, sequence variation, average number of nucleotide differences (D), and average number of nucleotide substitutions (Dxy) per site between breeds. Neighbour-joining tree for all haplotypes was constructed using Mega version 5 [36]. Analysis of molecular variance (AMOVA) was performed with Arlequin version 3.11 [37]. Sequences of the same D-loop fragment in wild sheep, published by Hiendleder et al. [33], were obtained from GenBank, Ovis vignei arkal (AY091489.1), Ovis vignei bochariensis (AY091490.1, AY091491.1, and AF039580.1), Ovis ammon collium (AY091492.1), Ovis ammon nigrimontana (AY091493.1 and AY091494.1), and used as outgroups in phylogenetic analysis.
Geographic distribution of eigenvectors was performed to investigate population genetic differences on the basis of their geographic distances. This approach permitted the generation of a synthetic configuration of locations based on the pairwise genetic distances that matched the real geographic configuration. Principal component analysis (PCA) scores for the first two components, obtained using Neiâs 1973 genetic distance, were plotted on a geographic map. As breeds are scattered among several farms, a virtual geographic entity representing the centroid of each breed on geographic maps was created using WGS84 geographical coordinates [38]. For a given component, it is a measure of the variance accounted for by that component. On thematic maps produced with the geographic information system (GIS) Manifold software package (Manifold System, Version 7, Manifold Net Ltd., Carson City, USA, http://www.manifold.net/), all breeds are thus represented according to a geometric distribution (see Figures 3(a) and 3(b)). Breeds showing high eigenvectors contribute sensibly to the explanation of the variance related to the component displayed. Classes were elaborated on the basis of the criterion of the natural breaks (Jenks optimization method). This algorithm reduces the variance within classes and maximizes the variance between classes. Colour classes were chosen in order to support the distinction between the different categories of behaviours observed: green: positive contribution; yellow: intermediary values; red: negative contribution to the component displayed.
2.4. Nuclear Polymorphism Analysis
The same 313 sheep belonging to 18 breeds sequenced at D-loop were genotyped with 37 previously described SNPs [39]. SNP ascertainment bias was minimised by sequencing target DNA in at least 8 individuals from different populations. Large-scale genotyping of all animals was performed by outsourcing to a commercial genotyping company (http://www.Kbioscience.co.uk/).
Allele frequencies, Neiâs estimation of observed and expected heterozygosities (Ho and He, resp.), were calculated using Fstat 2.93 [40]. Weir and Cockerhamâs [41] estimates of per population, per locus, and population pairs were calculated for each locus using Genalex 4.0 [42]. The same software was used to test deviations from Hardy-Weinberg equilibrium (HWE) for each locus and population and for locus over all populations; test for conformity with HWE expectations was assessed by calculating the Chi-squared value.
Correlation between geographic and Neiâs 1973 pairwise genetic distances was tested using Mantel tests (999 permutation) implemented in Genalex 4.0 software [42].
A PCA was performed on the covariance matrix of SNP frequency data to investigate spatial patterns of genetic variation using GENETIX software [43].
Nei [44] and Reynolds [45] genetic distances between population pairs were calculated using PowerMarker v3.25 [46].
Geographic distribution of eigenvectors was performed as described above using pairwise genetic distances [47] calculated on the basis of the selected SNP markers.
3. RESULTS
3.1. Mitochondrial Haplotypes
Ninety-three polymorphic sites and 154 haplotypes were identified from 313 sequences. Relatively high haplotype diversity was found in all three sampled geographic regions; the largest nucleotide diversity is present in Albania (0.02107) while the highest number of haplotypes observed is recorded in Greece (83) (Table 2).
|
The average number of nucleotide differences and the average number of nucleotide substitutions per site were used to calculate the genetic distance between breeds. The lowest distance was observed between Laticauda and Anogeiano (D: 2.357âDxy: 0.006), while the highest distance was observed between Bardhoka and Kymi (D: 12.450âDxy: 0.03) (Table 3).
|
AMOVA revealed that mitochondrial diversity is mainly distributed within breeds (95.04%) and only in part among regions (0.90%); low variability was also found among breeds/within regions (4.06%) (Table 4).
|
3.2. Phylogenetic Analysis and Haplogroups
The NJ tree obtained from mtDNA haplotypes and wild sheep sequences, used as out-group, revealed three of the five haplogroups described in the literature: A, B, and C (Figure 1). Haplogroup B is the most frequent among the analyzed samples (89%), while A and C are less common (8% and 3%, resp.). Greek and Albanian breeds are present in all three haplogroups, while Italian breeds are present only in haplogroups B and A (Table 5). This is shown also in Figure 2, representing the percentage of occurrence of each haplogroup in Albania, Greece, and Italy.
|



(a)

(b)

(c)

(d)
3.3. SNP Analysis
A total of 37 SNPs identified as polymorphic on eighteen sheep breeds selected throughout Europe [39] were applied to genotype 313 individuals of three Albanian, ten Greek, and five Italian sheep breeds. After removing those found monomorphic in the selected breeds, 27 SNPs were used for the analysis.
The frequencies of the major alleles ranged from 0.99 for the locus LEP1 to 0.538 for the locus IL2_1. Except for CAST_1, LEP1, LEP2, GDF8, and PRNP_1, which show frequencies of the rare alleles of 0.035, 0.010, 0.020, 0.021, and 0.036, respectively, the remaining SNPs have rare allele frequency of greater than 5%. Observed heterozygosity of all loci ranged from 0.019 (LEP1) to 0.489 (IL2_1), with a mean of 0.250. Expected heterozygosity of the loci ranged from 0.018 (LEP1) to 0.474 (IL2_1), with a mean of 0.259.
Significant deviations from Hardy-Weinberg equilibrium (HWE) for each locus and population (P value < 0.01) were observed for nine loci. Locus BMPR was not in HWE in Skopelos breed. Locus ACVR2B_2 was not in HWE in Bardhoka population, and locus SERPINA3 was not in HWE in Karagouniko population. Loci ZP2 and KRTAP6 were not in HWE in Altamurana. PRNP_1 was not in HWE in four Italian breeds: Bergamasca, Delle Langhe, Gentile di Puglia, and Laticauda. CSN3 was not in HWE in Ruda, Skopelos, and Bergamasca populations; KRT1 in Delle Langhe, Laticauda, and Karagouniko; GHR in Karagouniko and Laticauda.
per population ranged from â0.102 (Shkordane) to 0.276 (Altamurana) with a mean of 0.033 ( value on 9720 randomisations).
The analysis of population differentiation revealed an overall of 0.048, that is, 4.8% of allelic variation accounted across breeds and 95.2% within breeds. Weir and Cockerhamâs [41] estimate of per locus ranged from 0.0132 (IL2_2) to 0.098 (DESMIN_2), with a mean of 0.047.
Both Neiâs [44] and Reynoldsâ [45] genetic distances were calculated (Table 6). The distances range from 0.022 (Ruda-Karagouniko) to 0.253 (Sfakia-Gentile di Puglia) using Reynoldsâ distances and from 0.008 (Ruda-Karagouniko) to 0.096 (Gentile di Puglia-Anogeiano) using Neiâs distances. Both indices indicate Ruda and Karagouniko as the breeds with the minimum pairwise distance and show the maximum distances between the Italian breed Gentile di Puglia and Anogeiano, Karagouniko, Kefalleneas, Kymi, Lesvos, Pilioritiko, Sfakia, Skopelos, and Bergamasca breeds. The Mantel test showed correspondence between geographic and genetic distances with a .
|
3.4. PCA on SNPs
Genetic relationships were also explored by means of PCA. To examine the overall pattern of population differentiation, we considered the first three axes, which cumulatively explained 48.87% of the total inertia contained in the data set (Figure 4). From PCA, it can be seen that the some breeds are quite differentiated, with good correspondence to geographical locations, even if SNPs were few (4 of the 5 Italian breeds are well separated from the main cluster). Particularly, a differentiation between northern Italian (Bergamasca and Delle Langhe) and southern Italian (Gentile di Puglia, Altamurana, and Laticauda) breeds can be seen. The projection of loci in the space formed by the first three principal components (data not shown) shows that the differentiation of outlying breeds is caused by a small number of SNPs: the Delle Langhe population is mainly affected by alleles in the MSTN gene along the second component.; the Skopelos breed position in the graph is affected by alleles in the PRP gene and the Laticauda breed (lesser) differentiation is mainly due to alleles in the CALPA and LEP genes.

PCA scores calculated on mtDNA marker for the first 2 components were plotted on a geographic map, using the centroids of the sampling area of each breed (Figures 3(a) and 3(b)). The highest eigenvector contribution (coloured in green) was observed for Albanian and Greek breeds, as expected, even if four breeds show unexpectedly low diversity.
As for mtDNA analyses, the results of the PCA based on SNPs data were also used for making inferences about population genetic differences on the basis of their geographic distances. PCA scores for the first two principal components were plotted on a geographic map, using the centroids of the sampling area of each breed (Figures 3(c) and 3(d)). The line separating the map in two regions shows the isoline for an eigenvalue of 0. RUDA breed is an exception showing an eigenvalue above 0 but located in an area where all other breeds show lower eigenvalues. In accordance with the domestication history of the species, genetic diversity was higher in south-eastern populations than in north-western populations. The first component of Figure 3(c), in fact, shows a regular loss of genetic diversity towards North West.
4. DISCUSSION
Agriculture arose mainly within the distribution range of the wild ancestors of the most valuable livestock, such as the Fertile Crescent of southwest Asia, where early farmers were able to outcompete local hunter-gatherers. Once livestock slowly spread northwest across Europe, farming also shifted northwest from the Fertile Crescent to areas where farming had never arisen independentlyâfirst to Greece, then to Italy, and finally to northwest Europe [50]. Therefore, today the most productive farming zones do not correspond to those most productive in the past. Then, founder effect, genetic drift, and natural or artificial selection led to the formation of distinct breeds or varieties [51].
We therefore focused on the analysis of sheep of three countries aligned on this route, to evidence signs of migration. Geographical isolation, natural and artificial selection for physical or productive characters, genetic drift, mutations, and interpopulation gene flows have altered gene frequencies over many generations. The genetic diversity within and across breeds and species forms the basis of our current animal genetic resources for food production and other purposes.
The nature of the markers used for the analysis can affect the detection of geographical structuring, as suggested by Naderi et al. [52]. In fact, mtDNA informativeness is limited because it does not detect male-mediated gene flow and does not predict the nuclear genomic diversity [53]. Moreover, results may be affected by phenomena such as homoplasy, incomplete lineage sorting, effective population sizes, and sex-biased dispersal [27]. By combining markers with different modes of inheritance and rates of evolution this bias can be minimized [54].
Our mtDNA analysis shows higher levels of sheep nucleotide diversity in the South-East, which is congruent with data reported in the literature [22] and congruent with the proximity to the domestication centre. This is confirmed by eigenvector analysis, which showed high contribution to variance by Albanian and Greek breeds, even if four breeds show unexpectedly low diversity. However, this behaviour can be explained by recent isolation or selection for some traits that reduced the overall genetic diversity through bottleneck.
Very high haplotypes diversity was found in all three regions analysed (greater than 0.9), in agreement with previously published works on Portuguese breeds [9], Indian breeds [55], and Balkan sheep [56]. The major mitochondrial variation is distributed within breeds (95.04%), while it is lower among regions (0.90%) and among breeds within regions (4.06%). Phylogenetic methods were employed to examine the evolutionary history of the 18 breeds. Neighbour-joining and median-joining network revealed three of the described haplogroups, A, B, and C. The mouflon shares a haplotype with domestic sheep, as previously reported by Hiendleder et al. [33].
The SNP analysis revealed a rare allele frequency <5% for LEP1 and LEP2 loci, in agreement with those observed on a different European breed panel by Pariset et al. [39]. Observed and expected mean heterozygosity also showed similar values to those reported in the same paper. Expected heterozygosity values, which can indicate response to selection, are higher than observed heterozygosity values (Hs 0.063, 0.07, and 0.042; Ho 0.052, 0.06, and 0.038 in CALPA, PRNP-1, and GDF8, resp.).
Among the breeds tested, Altamurana showed the highest value suggesting the inbreeding in this population.
Regarding the phylogeographic structure we found that the 95.2% of variation occurred within breeds indicating the weak phylogeographic structure in sheep. These data are consistent with those previously published by Kijas et al. using a different SNP panel [17]. Sheep generally do not have a strong geographic structure and show a high genetic variability within breeds.
Anyway, Mantel test analysis using SNPs revealed a correlation between genetic and geographic distance. The possibility to assess the presence of a geographic component in genetic diversity using SNPs was already reported in previous studies on sheep [39] and goats [10, 57, 58].
In the PCA, the breeds appear differentiated with 48.9% of the variance explained by the first three principal components. Also this analysis shows a good correspondence to geographical locations: the breeds remaining separated by the main group are all Italian. PCA supports therefore a westward route to Italy that could indicate that transport of animals made by sea as already proposed for cattle [11, 12] and goats [8, 59]. This is particularly plausible because small sized species as sheep are easy to transport during migration and commercial trade and can adapt easily to various environments [53, 60, 61]. A similar decrease in genetic diversity as well as an increase in the level of differentiation at the breed level from South-East to North-West in European sheep breeds, supporting the hypothesis of livestock migration from the Middle East towards western and northern Europe, was found by Peter et al. [16] and Lawson Handley et al. [15], using other nuclear markers.
The formulation of the modern breed concept during mid-1800s has caused remarkable changes in the livestock sector: large-scale production expanded and its application to breeding and husbandry practices led to the formation of well-defined breeds that were exposed to intense anthropogenic selection. The differentiation of three breeds observed using PCA analysis could be related to a recent selection, which appears to be linked to CALPA (Laticauda), PRNP_1 (Skopelos), and GDF8 (Delle Langhe) SNP markers. Gentile di Puglia breed seems influenced by both CALPA and GDF8. In particular, these two genes have an effect on conformation and therefore are an easy target of selection. Other SNPs related to meat traits were found potentially under the effect of selection and apparently not associated with production attitude of the breeds [58].
5. CONCLUSIONS
We employed mtDNA and nuclear SNPs to investigate the genetic diversity of sheep breeds of three countries of the Mediterranean basin: Albania, Greece, and Italy. Our results showed significant genetic differentiation among the sheep breeds, supported by mtDNA and by SNP. The differentiation identified by nuclear markers could indicate a reduced gene flow due to geographical isolation, associated with different flocks management, or an effect of the introduction of different stocks centuries ago (cf. Figure 3(d), showing the 2nd dimension geographic distribution of eigenvectors). In general, D-loop sequence analysis shows a pattern reflecting migrations that occurred in postneolithic and historical times, with the most divergent mtDNA lineage occurring in the southern breeds, as shown in Figure 2 and Table 5. PCA on SNP data differentiated breeds with good correspondence to geographical locations. It is interesting that the correlation between genetic and geographic distances revealed using nuclear SNPs was not confirmed by mtDNA, for which Mantel test was not significant. Our results seem to indicate a better correlation between geographic distances and autosomal markers.
APPENDIX
Members of the Econogene Consortium
Abo-Shehada Mahamoud, Parasitology Resesarch Laboratory, Faculty of Veterinary Medicine, Jordan University of Science and Technology, Ibid, Jordan; Ajmone Marsan Paolo, Istituto di Zootecnica, UniversitĂ Cattolica del S. Cuore, Piacenza, Italy; Al Tarrayrah Jamil, Department of Animal Production, Aristotle Faculty of Agriculture, University of Thessaloniki, Thessaloniki, Greece; Angiolillo Antonella, Dipartimento di Scienze Animali Vegetali e dellâAmbiente, UniversitĂ del Molise, Campobasso, Italy; Baret Philip, Department of Animal Genetics, Catholic University of Louvain, Louvain-la-Neuve, Belgium; Baumung Roswitha, University of Natural Resources and Applied Life Sciences, Vienna, Vienna, Austria; Beja-Pereira Albano, Laboratoire de Biologie des Populations dâAltitude, UniversitĂ© Joseph Fourier et Centre National de la Recherche Scientifique, Grenoble, France; Bertaglia Marco, Department of Food Economics and Consumption Studies, University of Kiel, Kiel, Germany; Bordonaro Salvatore, Dipartimento di Scienze Agronomiche e delle Produzioni Animali, UniversitĂ degli Studi di Catania, Catania, Italy; Bruford Mike, School of Biosciences, Cardiff University, Cardiff, UK; Caloz RĂ©gis, Laboratoire de SystĂ©mes dâInformation GĂ©ographique (LASIG), Ecole Polytechnique FĂ©dĂ©rale de Lausanne (EPFL), Lausanne, Switzerland; Canali Gabriele, Istituto di Economia Agro-Alimentare, UniversitĂ Cattolica del S. Cuore, Piacenza, Italy; Canon Javier, Universidad Computense de Madrid, Madrid, Spain; Cappuccio Irene, Dipartimento di Produzioni Animali, UniversitĂ della Tuscia, Viterbo, Italy; Carta Antonello, Istituto Zootecnico e Caseario per la Sardegna, Sassari, Italy; Cicogna Mario, Istituto di Zootecnia Generale, UniversitĂ degli Studi di Milano, Milano; Italy, Crepaldi Paola, Istituto di Zootecnia Generale, UniversitĂ degli Studi di Milano, Milano, Italy; Dalamitra Stella, Biological Sciences, University of East Anglia, Norwich, UK; Daniela Krugmann, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Dobi Petrit, Department of Animal Production, Faculty of Agriculture, Tirana, Albania; Dominik Popielarczyk, Department of Sheep and Goat Breeding, Warsaw Agricultural University, Warsaw, Poland; Dunner Susana, Universidad Computense de Madrid, Madrid, Spain; DâUrso Giuseppe, Dipartimento di Scienze Agronomiche e delle Produzioni Animali, UniversitĂ degli Studi di Catania, Catania, Italy; El-Barody M. A. A., Animal Production Department, Faculty of Agriculture, Minia University, Minia, Egypt; England Phillip, Laboratoire de Biologie des Populations dâAltitude, UniversitĂ© Joseph Fourier et CNRS, Grenoble, France; Erhardt Georg, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Ertugrul Okan, Veterinary Faculty, Ankara University, Ankara, Turkey; Eva-Maria Prinzenberg, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Eveline Ibeagha-Awemu, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Ewa Strzelec, Department of Sheep and Goat Breeding, Warsaw Agricultural University, Warsaw, Poland; Fadlaoui Aziz, Department of Animal Genetics, Catholic University of Louvain, Louvain-la-Neuve, Belgium; Fornarelli Francesca, Istituto di Zootecnia Generale, UniversitĂ degli Studi di Milano, Milano, Italy; Garcia David, Universidad Computense de Madrid, Madrid, Spain; Georgoudis Andreas, Department of Animal Production, Faculty of Agriculture, Aristotle University of Thessaloniki, Thessaloniki, Greece; Gesine LĂŒhken, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Giovenzana Stefano, Istituto di Zootecnia Generale, UniversitĂ degli Studi di Milano, Milano, Italy; Gutscher Katja, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Hewitt Godfrey, Biological Sciences, University of East Anglia, Norwich, UK; Hoda Anila, Department of Animal Production, Faculty of Agriculture, Tirana, Albania; Horst Brandt, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Istvan Anton, Research Institute for Animal Breeding and Nutrition, Department of Genetics, Herceghalom, Hungary; Juma Gabriela, School of Biosciences, Cardiff University, Cardiff, UK; Joost StĂ©phane, Laboratoire de SystĂšmes dâInformation GĂ©ographique (LASIG), Ecole Polytechnique FĂ©dĂ©rale de Lausanne (EPFL), Lausanne, Switzerland; Jones Sam, The Sheep Trust, Biology Area-8, University of York, York, England; Karetsou Katerina, Department of Animal Production, Faculty of Agriculture, Aristotle University of Thessaloniki, Thessaloniki, Greece; Kliambas Georgios, Department of Animal Production, Faculty of Agriculture, Aristotle University of Thessaloniki, Thessaloniki, Greece; Koban Evren, Department of Biology, Middle East Technical University, Ankara, Turkey; Kutita Olga, Department of Animal Production, Faculty of Agriculture, Aristotle University of Thessaloniki, Thessaloniki, Greece; Lazlo Fesus, Research Institute for Animal Breeding and Nutrition, Department of Genetics, Herceghalom, Hungary; Lenstra Johannes A., Faculty of Veterinary Medicine, Utrecht University, Utrecht, The Netherlands; Ligda Christina, National Agricultural Research Foundation, Thessaloniki, Greece; Lipsky Shirin, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Luikart Gordon, Laboratoire de Biologie des Populations dâAltitude, UniversitĂ© Joseph Fourier et CNRS, Grenoble, France; Marie-Louise Glowatzki, Institut fĂŒr Genetik ErnĂ€hrung und Haltung von Haustieren, UniversitĂ€t Bern, Bern, Switzerland; Marilli Marta, Istituto di Zootecnia Generale, UniversitĂ degli Studi di Milano, Milano, Italy; Marletta Donata, Dipartimento di Scienze Agronomiche e delle Produzioni Animali, UniversitĂ degli Studi di Catania, Catania, Italy; Milanesi Elisabetta, Istituto di Zootecnica, UniversitĂ Cattolica del S. Cuore, Piacenza, Italy; Negrini Riccardo, Istituto di Zootecnica, UniversitĂ Cattolica del S. Cuore, Piacenza, Italy; Nijman, IsaĂ€c J., Faculty of Veterinary Medicine, Utrecht University, Utrecht, The Netherlands; Obexer-Ruff Gabriela, Institut fĂŒr Genetik ErnĂ€hrung und Haltung von Haustieren, UniversitĂ€t Bern, Bern, Switzerland; Papachristoforou Christos, Agricultural Research Institute, Nicosia, Cyprus; Pariset Lorraine, Dipartimento di Produzioni Animali, UniversitĂ della Tuscia, Viterbo, Italy; Pellecchia Marco, Istituto di Zootecnica, UniversitĂ Cattolica del S. Cuore, Piacenza, Italy; Peter Christina, Institut fĂŒr Tierzucht und Haustiergenetik, Justus Liebig-UniversitĂ€t of GieĂen, GieĂen, Germany; Perez Trinidad, School of Biosciences, Cardiff University, Cardiff, UK; PietrolĂ Emilio, Dipartimento di Scienze Animali, Vegetali e dellâAmbiente, UniversitĂ del Molise, Campobasso, Italy; Pilla Fabio, Dipartimento di Scienze Animali, Vegetali e dellâAmbiente, UniversitĂ del Molise, Campobasso, Italy; Roman Niznikowski, Department of Sheep and Goat Breeding, Warsaw Agricultural University, Warsaw, Poland; Roosen Jutta, Department of Food Economics and Consumption Studies, University of Kiel, Kiel, Germany; Scarpa Riccardo, Environment Department, University of York, York, UK; Sechi Tiziana, Istituto Zootecnico e Caseario per la Sardegna, Sassari, Italy; Taberlet Pierre, Laboratoire dâEcologie Alpine (LECA), UniversitĂ© Joseph Fourier, Grenoble, France; Taylor Martin, Biological Sciences, University of East Anglia, Norwich, UK; Togan Inci, Department of Biology, Middle East Technical University, Ankara, Turkey; Trommetter Michel, Institut de recherchĂ© Agronomique, UnitĂ© dâĂconomie et Sociologie Rurales, Grenoble, France; Valentini Alessio, Dipartimento di Produzioni Animali, UniversitĂ della Tuscia, Viterbo, Italy; Van Cann, Lisette M., Faculty of Veterinary Medicine, Utrecht University, Utrecht, The Netherlands; Vlaic Augustin, Department of Animal Genetics, Faculty of Zootechnics, University of Cluj-Napoca, Cluj-Napoca, Romania; Wiskin Louise, Biological Sciences, University of East Anglia, Norwich, UK; Zundel Stephanie, Laboratoire dâEcologie Alpine (LECA), UniversitĂ© Joseph Fourier, Grenoble, France.
ACKNOWLEDGMENTS
This work has been supported by a PRIN 2007 (prot. 2007XRBAEN_003) of Italian Ministry of University and Research. Part of the sampling was operated in the frame of the EU ECONOGENE project (QLK5-CT2001-02461). The content of the publication does not necessarily represent the views of the commission or its services. The authors wish to thank Gabriella Porcai and Paolo Ciorba for technical assistance.
References
- J. D. Vigne, I. CarrĂšre, and J. Guilaine, âUnstable status of early domestic ungulates in the Near East: the example of Shillourokambos (Cyprus, IX-VIIIth millennia Cal. B.C.),â Bulletin de Correspondance Hellenique, vol. 43, supplement, pp. 239â251, 2003. View at: Google Scholar
- C. F. Nadler, K. V. Korobitsina, R. S. Hoffmann, and N. N. Vorontsov, âCytogenetic differentiation, geographic distribution, and domestication in Palaearctic sheep (Ovis),â Zeitschrift fĂŒr SĂ€ugetierkunde, vol. 38, pp. 109â125, 1973. View at: Google Scholar
- T. D. Bunch, W. C. Foote, and J. J. Spillett, âTranslocations of acrocentric chromosomes and their implications in the evolution of sheep (Ovis),â Cytogenetics and Cell Genetics, vol. 17, no. 3, pp. 122â136, 1976. View at: Google Scholar
- K. Dobney and G. Larson, âGenetics and animal domestication: new windows on an elusive process,â Journal of Zoology, vol. 269, no. 2, pp. 261â271, 2006. View at: Publisher Site | Google Scholar
- M. A. Zeder, âDomestication and early agriculture in the Mediterranean Basin: origins, diffusion, and impact,â Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 33, pp. 11597â11604, 2008. View at: Publisher Site | Google Scholar
- J. ZilhĂŁo, âRadiocarbon evidence for maritime pioneer colonization at the origins of farming in west Mediterranean Europe,â Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 24, pp. 14180â14185, 2001. View at: Publisher Site | Google Scholar
- M. Tapio, N. Marzanov, M. Ozerov et al., âSheep mitochondrial DNA variation in European, Caucasian, and Central Asian areas,â Molecular Biology and Evolution, vol. 23, no. 9, pp. 1776â1783, 2006. View at: Publisher Site | Google Scholar
- F. Pereira, S. J. M. Davis, L. Pereira, B. McEvoy, D. G. Bradley, and A. Amorim, âGenetic signatures of a Mediterranean influence in Iberian Peninsula sheep husbandry,â Molecular Biology and Evolution, vol. 23, no. 7, pp. 1420â1426, 2006. View at: Publisher Site | Google Scholar
- F. Pereira, L. Pereira, B. Van Asch, D. G. Bradley, and A. Amorim, âThe mtDNA catalogue of all Portuguese autochthonous goat (Capra hircus) breeds: high diversity of female lineages at the western fringe of European distribution,â Molecular Ecology, vol. 14, no. 8, pp. 2313â2318, 2005. View at: Publisher Site | Google Scholar
- L. Pariset, A. Cuteri, C. Ligda, P. Ajmone-Marsan, and A. Valentini, âGeographical patterning of sixteen goat breeds from Italy, Albania and Greece assessed by Single Nucleotide Polymorphisms,â BMC Ecology, vol. 9, article 20, 2009. View at: Publisher Site | Google Scholar
- T. Cymbron, A. R. Freeman, M. I. Malheiro, J. D. Vigne, and D. G. Bradley, âMicrosatellite diversity suggests different histories for Mediterranean and Northern European cattle populations,â Proceedings of the Royal Society B, vol. 272, no. 1574, pp. 1837â1843, 2005. View at: Publisher Site | Google Scholar
- M. Pellecchia, R. Negrini, L. Colli et al., âThe mystery of Etruscan origins: novel clues from Bos taurus mitochondrial DNA,â Proceedings of the Royal Society B, vol. 274, no. 1614, pp. 1175â1179, 2007. View at: Publisher Site | Google Scholar
- J. R. S. Meadows, O. Hanotte, C. DrögemĂŒller et al., âGlobally dispersed Y chromosomal haplotypes in wild and domestic sheep,â Animal Genetics, vol. 37, no. 5, pp. 444â453, 2006. View at: Publisher Site | Google Scholar
- J. R. S. Meadows and J. W. Kijas, âRe-sequencing regions of the ovine Y chromosome in domestic and wild sheep reveals novel paternal haplotypes,â Animal Genetics, vol. 40, no. 1, pp. 119â123, 2009. View at: Publisher Site | Google Scholar
- L. J. Lawson Handley, K. Byrne, F. Santucci et al., âGenetic structure of European sheep breeds,â Heredity, vol. 99, no. 6, pp. 620â631, 2007. View at: Publisher Site | Google Scholar
- C. Peter, M. Bruford, T. Perez, S. Dalamitra, G. Hewitt, and G. Erhardt, âGenetic diversity and subdivision of 57 European and Middle-Eastern sheep breeds,â Animal Genetics, vol. 38, no. 1, pp. 37â44, 2007. View at: Publisher Site | Google Scholar
- J. W. Kijas, D. Townley, B. P. Dalrymple et al., âA genome wide survey of SNP variation reveals the genetic structure of sheep breeds,â PLoS ONE, vol. 4, no. 3, Article ID e4668, 2009. View at: Publisher Site | Google Scholar
- B. Chessa, F. Pereira, F. Arnaud et al., âRevealing the history of sheep domestication using retrovirus integrations,â Science, vol. 324, no. 5926, pp. 532â536, 2009. View at: Publisher Site | Google Scholar
- S. Y. Chen, Z. Y. Duan, T. Sha, J. Xiangyu, S. F. Wu, and Y. P. Zhang, âOrigin, genetic diversity, and population structure of Chinese domestic sheep,â Gene, vol. 376, no. 2, pp. 216â223, 2006. View at: Publisher Site | Google Scholar
- S. Pedrosa, M. Uzun, J. J. Arranz, B. GutiĂ©rrez-Gil, F. San Primitivo, and Y. BayĂłn, âEvidence of three maternal lineages in near eastern sheep supporting multiple domestication events,â Proceedings of the Royal Society B, vol. 272, no. 1577, pp. 2211â2217, 2005. View at: Publisher Site | Google Scholar
- J. R. S. Meadows, K. Li, J. Kantanen et al., âMitochondrial sequence reveals high levels of gene flow between breeds of domestic sheep from Asia and Europe,â Journal of Heredity, vol. 96, no. 5, pp. 494â501, 2005. View at: Publisher Site | Google Scholar
- J. R. S. Meadows, I. Cemal, O. Karaca, E. Gootwine, and J. W. Kijas, âFive ovine mitochondrial lineages identified from sheep breeds of the near east,â Genetics, vol. 175, no. 3, pp. 1371â1379, 2007. View at: Publisher Site | Google Scholar
- N. J. Wood and S. H. Phua, âVariation in the control region sequence of the sheep mitochondrial genome,â Animal Genetics, vol. 27, no. 1, pp. 25â33, 1996. View at: Google Scholar
- S. Hiendleder, H. Lewalski, R. Wassmuth, and A. Janke, âThe complete mitochondrial DNA sequence of the domestic sheep (Ovis aries) and comparison with the other major ovine haplotype,â Journal of Molecular Evolution, vol. 47, no. 4, pp. 441â448, 1998. View at: Publisher Site | Google Scholar
- S. Hiendleder, B. Kaupe, R. Wassmuth, and A. Janke, âMolecular analysis of wild and domestic sheep questions current nomenclature and provides evidence for domestication from two different subspecies,â Proceedings of the Royal Society B, vol. 269, no. 1494, pp. 893â904, 2002. View at: Publisher Site | Google Scholar
- J. Guo, L. X. Du, Y. H. Ma et al., âA novel maternal lineage revealed in sheep (Ovis aries),â Animal Genetics, vol. 36, no. 4, pp. 331â336, 2005. View at: Publisher Site | Google Scholar
- J. Flanders, G. Jones, P. Benda et al., âPhylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum: contrasting results from mitochondrial and microsatellite data,â Molecular Ecology, vol. 18, no. 2, pp. 306â318, 2009. View at: Publisher Site | Google Scholar
- G. Heckel, R. Burri, S. Fink, J. F. Desmet, and L. Excoffier, âGenetic structure and colonization processes in European populations of the common vole, Microtus arvalis,â Evolution, vol. 59, no. 10, pp. 2231â2242, 2005. View at: Publisher Site | Google Scholar
- J. Bryja, C. Smith, A. KoneÄnĂœ, and M. Reichard, âRange-wide population genetic structure of the European bitterling (Rhodeus amarus) based on microsatellite and mitochondrial DNA analysis,â Molecular Ecology, vol. 19, no. 21, pp. 4708â4722, 2010. View at: Publisher Site | Google Scholar
- M. M. Tollefsrud, J. H. SĂžnstebĂž, C. Brochmann, O. Johnsen, T. SkrĂžppa, and G. G. Vendramin, âCombined analysis of nuclear and mitochondrial markers provide new insight into the genetic structure of North European Picea abies,â Heredity, vol. 102, no. 6, pp. 549â562, 2009. View at: Publisher Site | Google Scholar
- J. W. O. Ballard and M. C. Whitlock, âThe incomplete natural history of mitochondria,â Molecular Ecology, vol. 13, no. 4, pp. 729â744, 2004. View at: Publisher Site | Google Scholar
- G. L. Gonçalves, G. R. P. Moreira, T. R. O. Freitas, D. Hepp, D. T. Passos, and T. A. Weimer, âMitochondrial and nuclear DNA analyses reveal population differentiation in Brazilian Creole sheep,â Animal Genetics, vol. 41, no. 3, pp. 308â310, 2010. View at: Publisher Site | Google Scholar
- S. Hiendleder, K. Mainz, Y. Plante, and H. Lewalski, â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,â Journal of Heredity, vol. 89, no. 2, pp. 113â120, 1998. View at: Publisher Site | Google Scholar
- T. A. Hall, âBioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT,â Nucleic Acids Symposium Series, vol. 41, pp. 95â98, 1999. View at: Google Scholar
- P. Librado and J. Rozas, âDnaSP v5: a software for comprehensive analysis of DNA polymorphism data,â Bioinformatics, vol. 25, no. 11, pp. 1451â1452, 2009. View at: Publisher Site | Google Scholar
- K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, âMEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,â Molecular Biology and Evolution. In press. View at: Publisher Site | Google Scholar
- L. Excoffier, G. Laval, and S. Schneider, âArlequin ver. 3.0: an integrated software package for population genetics data analysis,â Evolutionary Bioinformatics Online, vol. 1, pp. 47â50, 2005. View at: Google Scholar
- S. Joost, L. Colli, P. V. Baret et al., âIntegrating geo-referenced multiscale and multidisciplinary data for the management of biodiversity in livestock genetic resources,â Animal Genetics, vol. 41, no. 1, pp. 47â63, 2010. View at: Publisher Site | Google Scholar
- L. Pariset, I. Cappuccio, P. Ajmone-Marsan et al., âCharacterization of 37 breed-specific single-nucleotide polymorphisms in sheep,â Journal of Heredity, vol. 97, no. 5, pp. 531â534, 2006. View at: Publisher Site | Google Scholar
- J. Goudet, âFSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3),â 2001, http://www2.unil.ch/popgen/softwares/fstat.htm. View at: Google Scholar
- B. S. Weir and C. C. Cockerham, âEstimating F-statistics for the analysis of population structure,â Evolution, vol. 38, no. 6, pp. 1358â1370, 1984. View at: Google Scholar
- R. Peakall and P. E. Smouse, âGENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research,â Molecular Ecology Notes, vol. 6, no. 1, pp. 288â295, 2006. View at: Publisher Site | Google Scholar
- K. Belkhir, P. Borsa, L. Chikhi, N. Raufaste, and F. Bonhomme, (1996â2004) GENETIX 4.05, logiciel sous Windows TM pour la gĂ©nĂ©tique des populations. Laboratoire GĂ©nome, Populations, Interactions, CNRS UMR 5171, UniversitĂ© de Montpellier II, Montpellier, France.
- M. Nei, âGenetic distance between populations,â American Naturalist, vol. 106, no. 949, p. 283, 1972. View at: Google Scholar
- J. Reynolds, B. S. Weir, and C. C. Cockerham, âEstimation of the coancestry coefficient: basis for a short-term genetic distance,â Genetics, vol. 105, no. 3, pp. 767â779, 1983. View at: Google Scholar
- K. Liu and S. V. Muse, âPowerMaker: an integrated analysis environment for genetic maker analysis,â Bioinformatics, vol. 21, no. 9, pp. 2128â2129, 2005. View at: Publisher Site | Google Scholar
- M. Nei, âAnalysis of gene diversity in subdivided populations,â Proceedings of the National Academy of Sciences of the United States of America, vol. 70, no. 12, pp. 3321â3323, 1973. View at: Google Scholar
- M. Nei, Molecular Evolutionary Genetics, Columbia University Press, New York, NY, USA, 1987.
- M. Slatkin and R. R. Hudson, âPairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations,â Genetics, vol. 129, no. 2, pp. 555â562, 1991. View at: Google Scholar
- J. Diamond, âEvolution, consequences and future of plant and animal domestication,â Nature, vol. 418, no. 6898, pp. 700â707, 2002. View at: Publisher Site | Google Scholar
- C. Berthouly, B. BedâHom, M. Tixier-Boichard et al., âUsing molecular and multivariate methods to study the genetic diversity on local European and Asian chickens breeds,â Animal Genetics, vol. 39, pp. 121â129, 2008. View at: Google Scholar
- S. Naderi, H. R. Rezaei, P. Taberlet et al., âLarge-scale mitochondrial DNA analysis of the domestic goat reveals six haplogroups with high diversity,â PLoS ONE, vol. 2, no. 10, Article ID e1012, 2007. View at: Publisher Site | Google Scholar
- M. W. Bruford, D. G. Bradley, and G. Luikart, âDNA markers reveal the complexity of livestock domestication,â Nature Reviews Genetics, vol. 4, no. 11, pp. 900â910, 2003. View at: Publisher Site | Google Scholar
- G. M. Hewitt, âGenetic consequences of climatic oscillations in the Quaternary,â Philosophical Transactions of the Royal Society B, vol. 359, no. 1442, pp. 183â195, 2004. View at: Publisher Site | Google Scholar
- V. C. Pardeshi, N. Y. Kadoo, M. N. Sainani, J. R. S. Meadows, J. W. Kijas, and V. S. Gupta, âMitochondrial haplotypes reveal a strong genetic structure for three Indian sheep breeds,â Animal Genetics, vol. 38, no. 5, pp. 460â466, 2007. View at: Publisher Site | Google Scholar
- M. Äinkulov, Z. Popovski, K. Porcu et al., âGenetic diversity and structure of the West Balkan Pramenka sheep types as revealed by microsatellite and mitochondrial DNA analysis,â Journal of Animal Breeding and Genetics, vol. 125, no. 6, pp. 417â426, 2008. View at: Publisher Site | Google Scholar
- I. Cappuccio, L. Pariset, P. Ajmone-Marsan et al., âAllele frequencies and diversity parameters of 27 single nucleotide polymorphisms within and across goat breeds,â Molecular Ecology Notes, vol. 6, no. 4, pp. 992â997, 2006. View at: Publisher Site | Google Scholar
- L. Pariset, I. Cappuccio, S. Joost et al., âCharacterization of single nucleotide polymorphisms (SNPs) in sheep and their variation as an evidence of selection,â Animal Genetics, vol. 37, pp. 290â292, 2006. View at: Google Scholar
- J. Cañón, D. GarcĂa, M. A. GarcĂa-Atance et al., âGeographical partitioning of goat diversity in Europe and the Middle East,â Animal Genetics, vol. 37, no. 4, pp. 327â334, 2006. View at: Publisher Site | Google Scholar
- G. Luikart, L. Gielly, L. Excoffier, J. D. Vigne, J. Bouvet, and P. Taberlet, âMultiple maternal origins and weak phylogeographic structure in domestic goats,â Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 10, pp. 5927â5932, 2001. View at: Publisher Site | Google Scholar
- H. FernĂĄndez, S. Hughes, J. D. Vigne et al., âDivergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas,â Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 42, pp. 15375â15379, 2006. View at: Publisher Site | Google Scholar
Copyright
Copyright © 2011 Lorraine Pariset et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.