Large-Scale Mitochondrial DNA Analysis of the Domestic Goat Reveals Six Haplogroups with High Diversity

Abstract

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.

About the Authors:

Saeid Naderi

Affiliations Laboratoire d'Ecologie Alpine, CNRS-UMR 5553, Université Joseph Fourier, Grenoble, France, Natural Resources Faculty, University of Guilan, Guilan, Iran

Hamid-Reza Rezaei

Affiliations Laboratoire d'Ecologie Alpine, CNRS-UMR 5553, Université Joseph Fourier, Grenoble, France, Environmental Sciences Department, Gorgan University of Agriculture and Natural Resources, Gorgan, Iran

Pierre Taberlet

Affiliation: Laboratoire d'Ecologie Alpine, CNRS-UMR 5553, Université Joseph Fourier, Grenoble, France

Stéphanie Zundel

Affiliation: Laboratoire d'Ecologie Alpine, CNRS-UMR 5553, Université Joseph Fourier, Grenoble, France

Seyed-Abbas Rafat

Affiliation: Animal Science Department, Faculty of Agriculture, University of Tabriz, Tabriz, Iran

Hamid-Reza Naghash

Affiliation: Laboratoire d'Ecologie Alpine, CNRS-UMR 5553, Université Joseph Fourier, Grenoble, France

Mohamed A. A. El-Barody

Affiliation: Animal Production Department, Faculty of Agriculture, Minia University, Minia, Egypt

Okan Ertugrul

Affiliation: Department of Genetics, Faculty of Veterinary Medicine, Ankara University, Ankara, Turkey

François Pompanon

*To whom correspondence should be addressed. E-mail: [email protected]

Affiliation: Laboratoire d'Ecologie Alpine, CNRS-UMR 5553, Université Joseph Fourier, Grenoble, France

for the Econogene Consortium

Introduction

More than 10,000 years ago, the transition of humans from hunting to the manipulation of the behavior of certain animals lead to the process of domestication [1]. This process contributed to the rise of human civilization by enabling people to settle into a sedentary lifestyle. The goat was one of the first domesticated animals [2]–[4]. It was a source of milk, meat, dung for fuel and materials for clothing and building such as hair, bone and skin [1], [5]. Archaeological studies suggested that the domestic goat Capra hircus was domesticated from the bezoar Capra aegagrus in the Fertile Crescent [e.g. 6]–[8]. This origin was confirmed by genetic studies based on mitochondrial [e.g. 9], [10] and nuclear DNA [11].

From the beginning of the domestication process, the exchange and transportation of domestic animals has been related to human migration and trade. This resulted in genetic (e.g., selection, gene flow) and demographic processes that explain the present worldwide distribution of more than 300 different breeds of Capra hircus and their genetic structure [2]. Thus, the present genetic diversity bears the molecular signature of past events, such as rapid demographic expansions. Therefore, the study of this diversity helps to reconstitute the evolutionary history of the goat [12] and could bring new facts that help to understand the history of domestication.

Mitochondrial DNA is commonly used for the study of domesticated species. The control region has been especially used for describing the genetic polymorphism of goats [13], because it is variable and structured enough across the geographical range of the species, and evolves at a constant rate [12]. Moreover, it allows maternal lineages to be followed and is less sensitive to introgression from wild species than nuclear DNA [13]. However, studies on nuclear genes are needed because they give information on gene flow and selection processes that had a great influence on the evolution of livestock species [12].

Luikart et al. [13] conducted the first study of the overall genetic structure of domestic goats at the worldwide scale. They analyzed 406 individuals representing 88 breeds from the old world. They found three mitochondrial haplogroups (A, B and C) that diverged more than 200,000 years ago and have undergone demographic expansion at different times. This would suggest multiple maternal origins of domestic goats or introgression of other haplotypes after the first domestication event. Moreover, they showed that most of genetic diversity occurred within breeds, and interpreted the very weak geographic structure as the result of the extensive transportation of goats among continents.

The initial global survey by Luikart et al. [13] has been followed by regional studies describing more precisely the genetic diversity of goat breeds. However, these studies were always realized in restricted geographic regions corresponding to different countries such as Pakistan [14], India [15], China [16], South Korea [17], Sicily [18], Spain [19], [20] and Portugal [21]. The existence of three new haplogroups has been suggested [14], [15], [18]. However, this has sometimes been based only on a few individuals, and without comparing the new divergent haplotypes to a sample representative of the worldwide haplotype diversity. In general, the identification of a new haplogroup might be controversial in the absence of standardized criteria. All previous studies describing the mitochondrial polymorphism of domestic animals use the term of “maternal lineage” for characterizing a group of closely related haplotypes. However, this term is ambiguous as it usually corresponds to many haplotypes, and thus to many maternal lineages sensu stricto. As a consequence, we propose to use “mitochondrial haplogroup” instead of “maternal lineage”, a term that is already in common use in genetic studies.

In this context, the goals of the present study are (i) to characterize the domestic goat mtDNA diversity based on a worldwide sampling and make a global synthesis including previous studies, (ii) to establish the relationships between mitochondrial haplogroups and to propose a clear nomenclature, and (iii) to give standard criteria for the definition of mitochondrial haplogroups. For this purpose we used data from previous studies (1484 sequences retrieved from GeneBank), and we analyzed 946 new samples from all over the old world. New samples were especially taken from localities that have not been adequately sampled before, and that may have played an important role in the history of goat domestication (i.e., Middle East and especially the Fertile Crescent).

Results

Sequence polymorphism

The HVI fragment of the control region shows a high polymorphism with 336 variable sites over the 558 bp of the alignment. We observed 285 substitutions (226 transitions and 59 transversions) and 110 insertions/deletions (from 1 to 76 bp). The 2430 individuals correspond to 1540 different haplotypes.

Phylogenetic analysis and genetic structure of domestic goats

The Neighbor-joining tree of the 2430 domestic goats (Figure 1) shows 6 highly divergent groups corresponding to different mitochondrial haplogroups called A, B, C, D, F (according to previous studies) and G. Each group has high haplotype diversity (Table 1), and has been defined by high bootstrap values (except a bootstrap of 53 % for A ; Figure 1A), and by high mean pairwise distance with all other groups (see below). The A haplogroup is the most represented when considering either the number of individuals or the number of haplotypes and is highly dominant all over the old world (Table 1 and Figure 2). Except for two individuals situated at the base of the B group, this clade is composed of two sub-groups, B1 (35 haplotypes) and B2 (9 haplotypes), as previously defined by Chen et al. [16]. The B group is mostly found in whole Asia, with a few individuals from the Sub-Saharan Africa and one European goat from Greece. The B2 individuals are restricted to China and Mongolia. Goats from the C group are from whole Asia and Europe and the D group is present in the whole Asia and Northern Europe. The three goats from the F group are from Sicily. The G group has not been reported until now, and is present in Middle East and Northern Africa, near the Fertile Crescent.

[Figure omitted. See PDF.]

Figure 1. Neighbor-joining trees of domestic goat based on 1540 mtDNA haplotypes (A) and on the 22 reference mtDNA haplotypes (B).

Distances were calculated using the Kimura 2-Parameter model with gamma correction (alpha = 0.28). On the (A) tree, the numbers on the branches represent bootstrap values out of 1000 replications, and the stars point out the position of reference individuals for each haplogroup used to construct the (B) tree (see Table 5).

https://doi.org/10.1371/journal.pone.0001012.g001

[Figure omitted. See PDF.]

Figure 2. Geographic distribution of domestic goat mtDNA haplogroups.

The size of each circle is proportional to the sample size and each specific haplotype is represented by a different colour.

https://doi.org/10.1371/journal.pone.0001012.g002

[Figure omitted. See PDF.]

Table 1. Genetic diversity of goat mtDNA haplogroups

https://doi.org/10.1371/journal.pone.0001012.t001

The haplotype diversity is very high all over the Eurasia and Africa with a value above 0.97 in 39 of the 54 studied countries (Table 2). More than 77% of the mtDNA variation is distributed within breeds while about 11% is found among breeds within geographic regions and 12 % among geographic regions (Table 3). Nevertheless this low but significant geographic structure is coherent with the fact that all breeds are composed of individuals from the A group, with eventually a lower percentage of individuals from other haplogroups (for about 25% of the breeds). This low geographic structure is also confirmed by the distribution of all haplogroups that are present in several regions (except for F). Most of the mtDNA diversity is distributed among groups and within geographic regions, while less than 4% of this variability is found among regions within groups (Table 3).

[Figure omitted. See PDF.]

Table 2. Geographic origin and characteristics of the studied domestic goat

https://doi.org/10.1371/journal.pone.0001012.t002

[Figure omitted. See PDF.]

Table 3. Partition of the genetic variance among haplogroups, breeds and continental regions revealed by hierarchical AMOVAs

https://doi.org/10.1371/journal.pone.0001012.t003

Demography of mitochondrial haplogroups

Because of the low number of goats in the F group, demographic parameters were not estimated for this group. The overall mismatch distribution shows a multi-modal distribution (Figure 3). The first peak with a maximum of 10 pairwise differences corresponds to the differences between haplotypes from the same group. Two other peaks with maxima at 27 and 39 pairwise differences correspond to differences between haplotypes from different groups. The distributions of within-groups and between-groups pairwise differences have a very thin overlap around 20 mismatches. The mismatch distribution analysis reveals a unimodal bell-shaped distribution of pairwise sequence differences for all haplogroups (Figure 3), except for B that is bimodal (data not shown). B1 and B2 are unimodal, and individuals from these sub- groups generally differ by 8 or 9 mismatches (always less than 14 mismatches). This unimodal pattern that is less clear for the D group, perhaps because of the low sample size (n = 13), would be coherent with recent demographic expansions. The time of expansion would differ according to the group, as suggested by the different means of pairwise distribution (Figure 3) and the estimations made under a model of pure demographic expansion [22] (Table 4). However, the validity of the expansion model used for estimating the expansion time is only accepted for the A, C groups (SSD P-Values <0.00001 and ≤0.05 respectively see Table 4). All groups have high growth rates indicating high demographic expansion (Table 4). The estimates differ according to the groups, but the overlapping of confidence intervals, as well as the different sample sizes, preclude further interpretation.

[Figure omitted. See PDF.]

Figure 3. Mismatch distributions for mtDNA haplogroups of domestic goats.

For the overall dataset, the distribution of pairwise differences were realized separately for comparisons between and within haplogroups.

https://doi.org/10.1371/journal.pone.0001012.g003

[Figure omitted. See PDF.]

Table 4. Estimation of demographic parameters from genetic data

https://doi.org/10.1371/journal.pone.0001012.t004

Discussion

High mtDNA diversity in domestic goat

The very high mt DNA diversity may partly result from a high mutation rate of the control region. Higher pedigree divergence rates than phylogenetic divergence rates have been shown for the control region in human [23] and other animals (e.g.,[24], [25]). This could explain that we observe a higher diversity than the one expected with the phylogenetic mutation rate estimated for Bovidae (i.e., 30.1 % of divergence per Myr on the total control region sequence based on the Bos/Bison divergence [26]). Such high variability could also result from the selection of polymorphism but, to our knowledge, this has never been shown for the control region. Another explanation would be the capture of a large part of the diversity of the wild ancestor (i.e., the bezoar) during the domestication, with a large maternal effective population size. Testing this last hypothesis requires comparing the diversity of goats to that of the bezoar [27].

Characteristics and nomenclature of mitochondrial haplogroups

Five reliable mitochondrial haplogroups have previously been described in domestic goats [13]–[15], [18]. However, most of the previous studies were based on local samples and thus only considered a part of the whole haplotype variability. Therefore, it may be difficult to assess the pertinence of defining a new group on the base of few haplotypes. It would also be difficult to make the correspondence between several studies analyzing samples from different geographic origins. Our study can lead to a clear nomenclature of goat mitochondrial haplogroups, because we analyzed 2430 goats representing 1540 different haplotypes from all over Africa, Asia and Europe (946 new sequences mainly from the region of domestication and 1484 sequences from previous studies). We revealed the existence of 6 highly divergent groups. Five of them (A, B, C, D and F) have already been described, and one (G) is a new group. The two sequences that have been previously used to define the E group [15] now fall within the A haplogroup. This is partly due to the finding of new haplotypes, which are intermediate between those from A and E used by Joshi et al. [15]. Therefore, the E group cannot be considered as a mitochondrial haplogroup anymore. The B clade is composed of two groups (B1 and B2) that have previously been described as “sub-lineages” by Chen et al. [16]. We agree that the B1 and B2 are part of the same haplogroup because the genetic divergence between them (pairwise differences always lower than 14 mismatches) is lower than the divergence between all pairs of haplogroups (more than 20 mismatches). They must be considered as two subgroups because even with a low divergence they are supported by robust bootstrap values.

Standard criteria for defining goat mitochondrial haplogroups

More generally, previous considerations point out the problem of defining groups and sub-groups. A new haplogroup is defined when it highly diverges from all other haplotypes. However, the haplogroups may change over time, as more and more haplotypes will be available. We faced this situation for the E haplogroup that is no valid any more. There is therefore a need for standard and easy-to-use criteria in order to assign new goat haplotypes to existing haplogroups or to define new haplogroups. A haplotype can be related to an existing group if it presents a moderate genetic divergence from this group. The difficulty may be to define what is a “moderate” divergence. It can be deduced from the distributions of pairwise sequence differences within and between haplogroups. For goats, almost all haplotypes from the same group differ by less than 20 mismatches (whatever the group) while haplotypes from different groups usually present more than 20 mismatches (Figure 3). This threshold value would give a quick and easy indication for almost all studied haplotypes. However, it may be inadequate for some haplotypes (about 1% in our study) because the two mismatch distributions overlap.

Given the increasing number of sequences available, analyzing new haplotypes together with all previously published sequences will be time consuming and will require huge computational resources. Moreover several programs cannot be used because the algorithm complexity does not allow managing such datasets. Especially when a few haplotypes from restricted localities are studied, their assignation to haplogroups should be quick and easy. For a first approach, an accurate solution would be to place the new different haplotypes in a phylogenetic tree containing sequences of reference representative of the diversity of C. hircus mitochondrial DNA. For this purpose we have selected 22 haplotypes representing the variability of the 6 present goat mitochondrial haplogroups (Table 5 and Figure 1B).

[Figure omitted. See PDF.]

Table 5. The 22 reference individuals of the 6 domestic goat haplogroups

https://doi.org/10.1371/journal.pone.0001012.t005

Four of the 1540 haplotypes present a tandemly repeated sequence of 76 bp. Three individuals are from the A group (from Iran, Morocco and India) and one from the B1 sub-group (Malaysia). Such tandem repeats are common in vertebrate species [28] and have already been found in the Bovidae family [29]. They are attributed to slippage-mispairing events that are more likely to appear in regions where the polymerase activity is interrupted [28]. This phenomenon corresponding to a single duplication event is found in a few individuals from different haplogroups, and has occurred more than once in the history of goats.

Genetic structure of domestic goats

Our results show that most of the genetic variation is found among goat haplogroups, with a weak phylogeographic structure. The strongly dominant A group (91 % of the goats) is distributed worldwide, and even if the other groups have more restricted distributions they still occupy large geographic areas (Figure 2). The F group is the exception, with three haplotypes restricted to a single locality (Sicily) that could have been brought along from recently captured wild goats. However, the sampling effort may still be insufficient to see the whole distribution of haplogroups other than A, because of their low frequency. The differences among geographic regions at the worldwide scale are low (about 12%) but significant. This is concordant with previous results showing a very weak phylogeographic structure of goats [13] and sheep [30], [31] compared to cattle [32], [33]. The genetic differences among continental regions could partly result from the differential geographic distribution of mitochondrial haplogroups. However, there is still a low but significant genetic variation (3.5%) among region within groups, indicating regional differentiations of haplotypes. At the regional scale, the lack of geographic structure has also been reported in several places [16], [19], [21] while a structure has been found in India [15]. The weak phylogeographic structure found today in goats has been explained by a high mobility of this species in relation to human migration and commercial trade [12], [13], [34]. This mobility would have been higher than those of cattle due to their versatility in feeding habits and ability to live under extreme conditions [1]. However, the mixing of goat haplogroups could have existed before the worldwide translocation of goats. The presence of goats in Cyprus 10,000 years ago [35] suggests that goats could have been translocated within the domestication area since the first domestication events. Moreover, we cannot exclude that the mt-haplogroups were already mixed in the wild ancestor before domestication. When considering the local scale, the genetic pattern of domestic goats also seems related to human history. For instance, the geographic structure found in Indian goats would have a common historical basis in the sequential migrations of human populations with different cultural and linguistic characteristics [15].

However, the information given by mitochondrial markers is limited because it does not detect male-mediated gene flow and does not predict the nuclear genomic diversity [12]. In particular, the breeds cannot be distinguished on the base of mtDNA [16], [19], [20] while nuclear markers show a genetic structure [36]–[38]. Our study confirmed that more than 77 % of the mtDNA variation is found within breeds and that nearly 25% of the breeds are composed of at least 2 haplogroups. This is in accordance with the recent fragmentation of local goat populations into discrete breeds about 200 years ago, under strong selection pressures on a few phenotypic traits [39]. This structure can be seen on nuclear markers linked to selected parts of the genome, but not on mitochondrial markers. Then, looking at the evolutionary history of breeds using mtDNA markers could lead to misinterpretation. For example, a breed composed of two mitochondrial haplogroups would have a bimodal mismatch distribution due to within- and between-breeds pairwise differences, and should not be interpreted in term of demographic history of the breed. Thus, fully understanding the evolutionary history of domestic goats would also require the use of nuclear markers.

Demography of mitochondrial haplogroups

The present structure of the genetic diversity retains the signature of past demographic events and helps reconstitute the evolutionary history [40]. The estimation of demographic parameters remains difficult because of the difficulties of verifying the hypothesis of the models used, of estimating accurate initial parameters (e.g., absolute date of domestication) and sometimes because of low sample sizes. However, rough estimations from the present work and previous studies [13], [15], [16] are concordant and agree on the same scenario. All haplogroups had a recent demographic expansion corresponding roughly to the period when domestication took place about 10,000 years ago. It is difficult to give relative dates of expansion because of large confidence intervals, especially for D and G groups, but our results confirm that the expansions of B and C groups were more recent than that of A [13]. Also, our results show that all groups had high growth rates, with a tendency for slower growth in B2 sub-group and C and G. A faster growth of A relative to C is in accordance with archaeozoological data: the genotyping of fossil goats showed that about 7000 years ago A and C were equally represented in Southern France [34] while A is strongly dominant in Southern Europe now.

Limits of genetic data from domestic goats for reconstituting the history of domestication

Divergence time between haplogroups has been estimated on adequate molecular markers (mainly cytochrome b) between 103,000 and 597,800 years [13]–[16]. All these values are far greater than the domestication time, showing that most of goat genetic diversity existed before domestication, and that several haplogroups were domesticated in one or several events. However, the genetic data available for domestic goats does not permit furthering our understanding of the domestication process and identifying potential domestication centre(s). A higher genetic diversity would have been expected near the Fertile Crescent where the goat domestication took place according to archaeological data, and where extensive sampling has been done. But the haplotype diversity is similar all over the world (more than 80% of the countries with a haplotype diversity greater than 0.9), because of the high migration rates in domestic goats due to human migration and commercial trade.

Moreover, the presence of a possible ancestral haplotype in a particular area does not prove that this is a domestication centre, since many events could have occurred to mask the real history (e.g., coalescence or founder effects). For instance the domestication of a B sub-group in China supported by genetic data [16] is doubtful since the wild ancestor of the domestic goat (i.e. the bezoar Capra aegagrus) has credibly never been present in this area [11], [41]. Overall, in order to fully understand the domestication of goats it is necessary to characterize the genetic diversity of wild goat species, and to establish the evolutionary relationships between wild and domesticated haplotypes.

Materials and Methods

Sampling and DNA extraction

Samples were collected from 946 individuals from 42 countries of the old world (See Table 2) from which 569 individuals were studied within the Econogene project (www.econogene.eu). Samples consisted of ear tissue preserved in ethanol 95% until extraction, or of blood collection. DNA was extracted from tissue using the Qiagen DNeasy tissue kit following the manufacturer's instructions, and from blood samples using QIAamp DNA blood kit.

To have a good coverage of the goat breeds, the dataset was completed with 1484 sequences containing the Capra hircus HVI control region (450 to 1200 bp long) retrieved from GenBank (Table 2).

DNA amplification and sequencing

The HVI segment of the control region was sequenced for all blood and tissue DNA extracts. Using the primers CAP-F (5′-CGTGTATGCAAGTACATTAC-3′) and CAP-R (5′-CTGATTAGTCATTAGTCCATC-3′), we amplified a fragment of 598 bp (without primers) that corresponds to the positions 15,653 to 16,250 on the complete goat mitochondrial sequence of reference ([42]; accession number AF533441). PCR amplifications were conducted in a 25 µl volume with 2 mM MgCl2, 200 µM of each dNTP, 1 µM of each primer and 1 unit of AmpliTaq Gold Polymerase (Applied Biosystems). After a 10 min period at 95°C for polymerase activation, 35 cycles were run with the following steps: 95°C: 30 s, 55°C: 30 s, 72°C: 1 min. PCR products were purified using the Qiaquick PCR purification kit (Qiagen). 35 ng of purified DNA from this PCR product was used for sequencing with the CAP-F or CAP-R primer. Sequence reactions were performed for both DNA strands by using the ABI PRISM Dye Terminator Cycle Sequencing Reaction Kit (Applied Biosystems) in a 20 µl volume with 2 µM of each primer. 25 cycles were run with the following steps 96°C: 30 s, 55°C: 30 s, 60°C: 4 min. Excess dye terminators were removed by spin-column purification and the products were electrophorezed on an ΑΒΙ 3700 PRISM DNA sequencer (Applied Biosystems) using the POP 7 polymer.

Sequences were edited for correction with the SeqScape v2.5 software (Applied Biosystems). All sequences were deposited in GenBank (Accession Numbers EF617601- EF618546, Table 2 and Table S1).

Sequences from GenBank and from our dataset were aligned with Mega v3.1 [43], and then adjusted by eye. For further analyses, we only kept the region used by Luikart et al. [13] because this is the part of the sequence available for most of the GenBank records, and also the most informative one. This region is 481 bp long and corresponds to the positions 15,707 to 16,187 on the Capra hircus reference sequence (mtDNA complete sequence of C. hircus, Accession number AF533441 [42]). According to the insertion/deletion events, the analyzed sequences ranged from 481 to 558 bp. For Indian goats a shorter fragment of 453 bp has been sequenced [15] and the 28 missing nucleotides were treated as missing data. The alignment of the 2430 sequences used in this study is provided as supplementary information (Table S1).

Data analysis

The substitution model used for the HVI region was the Kimura 2-parameters model, as previously used on several subsets of the present dataset (e.g., [13], [15]). The heterogeneity in substitution rates among nucleotide sites was modelled by a gamma distribution. The alpha parameter was estimated as the mean of 10 estimations by a maximum-likelihood method under the Kimura 2-parameters model using PAML v 2.0.2 [44]. Each estimation was based on the analysis of 1000 individuals randomly chosen in the dataset of 2430 individuals. The estimated value (alpha = 0.28) was similar to that estimated for the same region on a smaller sample of domestic and wild goats by Luikart et al. [13]. These settings were used for further phylogenetic reconstruction and analysis of genetic diversity. We used 1484 published sequences for checking the validity of the halpogroups previously defined (see Table 2 and Table S1 for references and GenBank accession numbers).

Given the very high number of sequences analyzed, the phylogenetic tree was constructed using the Neighbor-joining method using PAUP* v 4.0 [45], with 1000 bootstraps for measuring branch robustness. The ARLEQUIN v 3.0 software [46] was used for estimating haplotype and nucleotide diversity, for analyzing mismatch distribution within mitochondrial haplogroups, and for estimating the parameters of demographic expansion. Four individuals that showed a 76 bp insertion were discarded for mismatch analyses and the analyses were thus performed on 481 bp long sequences. The expansion time was estimated under a model of pure demographic expansion [22] with parameters set to default values in ARLEQUIN 3.0. The parameter of demographic expansion τ was estimated according to the method of Schneider and Excoffier [47]. The validity of the expansion model was tested using the sum of square deviations (SSD) between the observed and expected mismatches [47]. Growth rates of mitochondrial haplogroups were estimated with Lamarc v2.2 [48] using a bayesian framework allowing migrations across haplogroups (with a maximum of 10000 migration events, default priors used for migration rates estimation). The estimation of growth rates was done with linear prior (upper bound of 1000 and lower bound of −500), 10 initial chains (500 samples, sampling interval of 20 and burn-in period of 1000) and 2 final chains (10000 samples, sampling interval of 20 and burn-in period of 1000).

In order to test the geographic structure of the mtDNA haplotype diversity, the goat distribution has been partitioned in 7 geographic regions (Northern Europe, Southern Europe, Northern Africa, Sub-Saharan Africa, Middle East, Western Asia and Eastern Asia, see Table 2). Two hierarchical AMOVA were performed using ARLEQUIN v3.0 to test the partition of the genetic variance among haplogroups and among continents within haplogroups, as well as among continents and among breeds within continents. This second AMOVA was performed on the 1602 goats for which the breeds were known.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Alignements of the 2430 control region sequences of domestic goat. The code for geographic regions are defined in Table 2. Missing data are coded as ‘?’.

https://doi.org/10.1371/journal.pone.0001012.s001

(1.24 MB XLS)

Acknowledgments

We would like to thank Prof. Oscar Gaggiotti for fruitful discussions and constructive comments, Cécile Albert for help with figures, and Jessica Scriven and Alison Cleary for help with the English. We are grateful to Ebrahim Ghaderi and Prof. Mahamoud Abo-Shehada for their help for collecting samples in Iran and Jordan. We also thank two anonymous reviewers for constructive and stimulating comments.

Institution at which research was done: Laboratoire d'Ecologie Alpine, Université Joseph Fourier, Grenoble, France.

Members of the Econogene Consortium (www.econogene.eu):

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, Aristotle University of Thessaloniki, Fac. of Agriculture, Dept. of Animal Production, Thessaloniki, Greece,

Angiolillo Antonella, Dipartimento di Scienze Animali Vegetali e dell'Ambiente, Università del Molise, Campobasso, Italy,

Baret Philip, Department of Animal Genetics, Université cath. 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, University 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, Ecole Polytechnique Fédérale de Lausanne (EPFL), Laboratoire de Systèmes d'Information Géographique (LASIG), 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 University of Giessen, Giessen, 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, University Joseph Fourier et CNRS, Grenoble, France,

Erhardt Georg, Institut für Tierzucht und Haustiergenetik, Justus-Liebig University of Giessen, Giessen, Germany,

Ertugrul Okan, Veterinary Faculty, Ankara University, Ankara, Turkey,

Eva-Maria Prinzenberg, Institut für Tierzucht und Haustiergenetik, Justus-Liebig University of Giessen, Giessen, Germany,

Eveline Ibeagha-Awemu, Institut für Tierzucht und Haustiergenetik, Justus-Liebig University of Giessen, Giessen, Germany,

Ewa Strzelec, Department of Sheep and Goat Breeding, Warsaw Agricultural University, Warsaw, Poland,

Fadlaoui Aziz, Department of Animal Genetics, Cath. 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, Aristotle University of Thessaloniki, Fac. of Agriculture, Dept. of Animal Production, Thessaloniki, Greece,

Gesine Lühken, Institut für Tierzucht und Haustiergenetik, Justus-Liebig University of Giessen, Giessen, Germany,

Giovenzana Stefano, Istituto di Zootecnia Generale, Università degli Studi di Milano, Milano, Italy,

Gutscher Katja, Institut für Tierzucht und Haustiergenetik, Justus-Liebig University of Giessen, Giessen, 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 University of Giessen, Giessen, Germany,

Istvan Anton, Research Institute for Animal Breedeing and Nutrition, Department of Genetics, Herceghalom, Hungary,

Juma Gabriela, School of Biosciences, Cardiff University, Cardiff, UK,

Joost Stéphane, Ecole Polytechnique Fédérale de Lausanne (EPFL), Laboratoire de Systèmes d'Information Géographique (LASIG), Lausanne, Switzerland,

Jones Sam, The Sheep Trust, Biology Area-8, University of York, York, England,

Karetsou Katerina, Aristotle University of Thessaloniki, Fac. of Agriculture, Dept. of Animal Production, Thessaloniki, Greece,

Kliambas Georgios, Aristotle University of Thessaloniki, Fac. of Agriculture, Dept. of Animal Production, Thessaloniki, Greece,

Koban Evren, Department of Biology, Middle East Technical University, Ankara, Turkey,

Kutita Olga, Aristotle University of Thessaloniki, Fac. of Agriculture, Dept. of Animal Production, Thessaloniki, Greece,

Lazlo Fesus, Research Institute for Animal Breedeing and Nutrition, Department of Genetics, Herceghalom, Hungary,

Lenstra Johannes A., Faculty of Veterinary Medicine, Utrecht University, Utrecht, Netherlands,

Ligda Christina, National Agricultural Research Foundation, Thessaloniki, Greece,

Lipsky Shirin, Institut für Tierzucht und Haustiergenetik, Justus-Liebig University of Giessen, Giessen, Germany,

Luikart Gordon, Laboratoire de Biologie des Populations d'Altitude, University 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, 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 University of Giessen, Giessen, 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 recherche 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, 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 Stéphanie , Laboratoire d'Ecologie Alpine (LECA), Université Joseph Fourier, Grenoble, France,

Author Contributions

Conceived and designed the experiments: PT FP SN HR SR. Performed the experiments: SN HR SZ. Analyzed the data: FP SN HR. Contributed reagents/materials/analysis tools: SN HR SZ SR HN ME OE. Wrote the paper: PT FP SN HR.

Citation: Naderi S, Rezaei H-R, Taberlet P, Zundel S, Rafat S-A, Naghash H-R, et al. (2007) Large-Scale Mitochondrial DNA Analysis of the Domestic Goat Reveals Six Haplogroups with High Diversity. PLoS ONE2(10): e1012. https://doi.org/10.1371/journal.pone.0001012

References

1. Clutton-Brock J (1999) Domesticated animals from early times. The Natural History Museum: Cambridge University Press. J. Clutton-Brock1999Domesticated animals from early times.The Natural History Museum: Cambridge University Press238

2. Porter V (1996) Goats of the world. Ipswich, Uk: Farming Press. V. Porter1996Goats of the world.Ipswich, UkFarming Press172

3. Pringle H (1998) Reading the signs of ancient animal domestication. Science 282: 1448–1450.H. Pringle1998Reading the signs of ancient animal domestication.Science28214481450

4. Zeder M, Emshwiller E, Smith B, Bradley D (2006) Documenting domestication: the intersection of genetics and archaeology. Trends Genet 22: 139–155.M. ZederE. EmshwillerB. SmithD. Bradley2006Documenting domestication: the intersection of genetics and archaeology.Trends Genet22139155

5. MacHugh D, Bradley D (2001) Livestock genetic origins: goats buck the trend. Proc Natl Acad Sci USA 98: 5382–5384.D. MacHughD. Bradley2001Livestock genetic origins: goats buck the trend.Proc Natl Acad Sci USA9853825384

6. Meadow RH (1996) The Origins and Spread of Agriculture and Pastoralism in northwestern South Asia. In: Harris DR. The Origins and Spread of Agriculture and Pastoralism in Eurasia. London: UCL Press. pp. 390–412.RH Meadow1996The Origins and Spread of Agriculture and Pastoralism in northwestern South Asia. In: Harris DR. The Origins and Spread of Agriculture and Pastoralism in Eurasia.LondonUCL Press390412

7. Zeder M, Hesse B (2000) The initial domestication of goats (Capra hircus) in the Zagros Mountains 10,000 years ago. Sciences 287: 2254–2257.M. ZederB. Hesse2000The initial domestication of goats (Capra hircus) in the Zagros Mountains 10,000 years ago.Sciences28722542257

8. Zeder M, Vigne JD, Helmer D (2005) A view from the Zagros: new perspectives on livestock domestication in the Fertile Crescent. In: Vigne JD, Peters J, Helmer D, editors. The first steps of animal domestication: Oxbow Books. pp. 125–146.M. ZederJD VigneD. Helmer2005A view from the Zagros: new perspectives on livestock domestication in the Fertile Crescent.JD VigneJ. PetersD. HelmerThe first steps of animal domestication: Oxbow Books125146

9. Takada T, Kikkawa Y, Yonekawa H, Kawakami S, Amano T (1997) Bezoar (Capra aegagrus) is a matriarchal candidate for ancestor of domestic goat (Capra hircus): evidence from the mitochondrial DNA diversity. Biochem Genet 35(9–10): 315–326.T. TakadaY. KikkawaH. YonekawaS. KawakamiT. Amano1997Bezoar (Capra aegagrus) is a matriarchal candidate for ancestor of domestic goat (Capra hircus): evidence from the mitochondrial DNA diversity.Biochem Genet35(9–10)315326

10. Manceau V, Despres L, Bouvet J, Taberlet P (1999) Systematics of the genus Capra inferred from mitochondrial DNA sequence data. Mol Phylogenet Evol 13(3): 504–510.V. ManceauL. DespresJ. BouvetP. Taberlet1999Systematics of the genus Capra inferred from mitochondrial DNA sequence data.Mol Phylogenet Evol13(3)504510

11. Pidancier N, Jordan S, Luikart G, Taberlet P (2006) Evolutionary history of the genus Capra (Mammalia, Artiodactyla): discordance between mitochondrial DNA and Y-chromosome phylogenies. Mol Phylogenet Evol 40: 739–749.N. PidancierS. JordanG. LuikartP. Taberlet2006Evolutionary history of the genus Capra (Mammalia, Artiodactyla): discordance between mitochondrial DNA and Y-chromosome phylogenies.Mol Phylogenet Evol40739749

12. Bruford M, Bradley D, Luikart G (2003) DNA markers reveal the complexity of livestock domestication. Nat Rev Genet 3: 900–910.M. BrufordD. BradleyG. Luikart2003DNA markers reveal the complexity of livestock domestication.Nat Rev Genet3900910

13. Luikart G, Gielly L, Excoffier L, Vigne JD, Bouvet J, et al. (2001) Multiple maternal origins and weak phylogeographic structure in domestic goats. Proc Natl Acad Sci USA 98: 5927–5932.G. LuikartL. GiellyL. ExcoffierJD VigneJ. Bouvet2001Multiple maternal origins and weak phylogeographic structure in domestic goats.Proc Natl Acad Sci USA9859275932

14. Sultana S, Mannen H, Tsuji S (2003) Mitochondrial DNA diversity of Pakistani goats. Anim Genet 34(6): 417–421.S. SultanaH. MannenS. Tsuji2003Mitochondrial DNA diversity of Pakistani goats.Anim Genet34(6)417421

15. Joshi MB, Rout PK, Mandal AK, Tyler-Smith C, Singh L, et al. (2004) Phylogeography and origin of Indian domestic goats. Mol Biol Evol 21(3): 454–462.MB JoshiPK RoutAK MandalC. Tyler-SmithL. Singh2004Phylogeography and origin of Indian domestic goats.Mol Biol Evol21(3)454462

16. Chen SY, Su YH, Wu SF, Sha T, Zhang YP (2005) Mitochondrial diversity and phylogeographic structure of Chinese domestic goats. Mol Phylogenet Evol 37(3): 804–814.SY ChenYH SuSF WuT. ShaYP Zhang2005Mitochondrial diversity and phylogeographic structure of Chinese domestic goats.Mol Phylogenet Evol37(3)804814

17. Odahara S, Chung H, Choi S, Yu S, Sasazaki S, et al. (2006) Mitochondrial DNA diversity of Korean native goats. Asian Australas J Anim Sci 19: 482–485.S. OdaharaH. ChungS. ChoiS. YuS. Sasazaki2006Mitochondrial DNA diversity of Korean native goats.Asian Australas J Anim Sci19482485

18. Sardina MT, Ballester M, Marmi J, Finocchiaro R, van Kaam JB, et al. (2006) Phylogenetic analysis of Sicilian goats reveals a new mtDNA lineage. Anim Genet 37(4): 376–378.MT SardinaM. BallesterJ. MarmiR. FinocchiaroJB van Kaam2006Phylogenetic analysis of Sicilian goats reveals a new mtDNA lineage.Anim Genet37(4)376378

19. Amills M, Capote J, Tomas A, Kelly L, Obexer-Ruff G, et al. (2004) Strong phylogeographic relationships among three goat breeds from the Canary Islands. J Dairy Res 71(3): 257–262.M. AmillsJ. CapoteA. TomasL. KellyG. Obexer-Ruff2004Strong phylogeographic relationships among three goat breeds from the Canary Islands.J Dairy Res71(3)257262

20. Azor PJ, Monteagudo LV, Luque M, Tejedor MT, Rodero E, et al. (2005) Phylogenetic relationships among Spanish goats breeds. Anim Genet 36(5): 423–425.PJ AzorLV MonteagudoM. LuqueMT TejedorE. Rodero2005Phylogenetic relationships among Spanish goats breeds.Anim Genet36(5)423425

21. Pereira F, Pereira L, Van Asch B, Bradley D, Amorim A (2005) The mtDNA catalogue of all Portuguese autochthonous goat (Capra hircus) breeds: high diversity of female lineages at the western fringe of European distribution. Mol Ecol 14: 2313–2318.F. PereiraL. PereiraB. Van AschD. BradleyA. Amorim2005The mtDNA catalogue of all Portuguese autochthonous goat (Capra hircus) breeds: high diversity of female lineages at the western fringe of European distribution.Mol Ecol1423132318

22. Rogers A, Harpending H (1992) Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol 9: 552–569.A. RogersH. Harpending1992Population growth makes waves in the distribution of pairwise genetic differences.Mol Biol Evol9552569

23. Howell N, Smejkal CB, Mackey DA, Chinnery PF, Turnbull DM, et al. (2003) The Pedigree Rate of Sequence Divergence in the Human Mitochondrial Genome: There Is a Difference Between Phylogenetic and Pedigree Rates. Am J Hum Genet 72: 659–670.N. HowellCB SmejkalDA MackeyPF ChinneryDM Turnbull2003The Pedigree Rate of Sequence Divergence in the Human Mitochondrial Genome: There Is a Difference Between Phylogenetic and Pedigree Rates.Am J Hum Genet72659670

24. Denver DR, Morris K, Lynch M, Vassilieva LL, Thomas WK (2000) High Direct Estimate of the Mutation Rate in the Mitochondrial Genome of Caenorhabditis elegans. Science 289: 2342.DR DenverK. MorrisM. LynchLL VassilievaWK Thomas2000High Direct Estimate of the Mutation Rate in the Mitochondrial Genome of Caenorhabditis elegans.Science2892342

25. Lambert DM, Ritchie PA, Millar CD, Holland B, Drummond AJ, et al. (2002) Rates of Evolution in Ancient DNA from Adélie Penguins. Science 295: 2270.DM LambertPA RitchieCD MillarB. HollandAJ Drummond2002Rates of Evolution in Ancient DNA from Adélie Penguins.Science2952270

26. Bradley D, MacHugh D, Cunningham P, Loftus R (1996) Mitochondrial diversity and the origins of African and European cattle. Proc Natl Acad Sci USA 93: 5131–5135.D. BradleyD. MacHughP. CunninghamR. Loftus1996Mitochondrial diversity and the origins of African and European cattle.Proc Natl Acad Sci USA9351315135

27. Naderi S, Rezaei HR, Pompanon F, Blum MGB, Negrini R, et al. (In Prep.) Genetic evidence for a large scale domestication in goatsS. NaderiHR RezaeiF. PompanonMGB BlumR. NegriniIn Prep.Genetic evidence for a large scale domestication in goats

28. Fumagalli L, Taberlet P, Favre L, Hausser J (1996) Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews. Mol Biol Evol 13: 31–46.L. FumagalliP. TaberletL. FavreJ. Hausser1996Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.Mol Biol Evol133146

29. Hiendleder S, Lewalski H, Wassmuth R, Janke A (1998) The complete mitochondrial DNA sequence of the domestic sheep (Ovis aries) and comparison with the other major ovine haplotype. J Mol Evol 47: 441–448.S. HiendlederH. LewalskiR. WassmuthA. Janke1998The complete mitochondrial DNA sequence of the domestic sheep (Ovis aries) and comparison with the other major ovine haplotype.J Mol Evol47441448

30. Meadows JRS, Li K, Kantanen J, Tapio M, Sipos W, et al. (2005) Mitochondrial Sequence Reveals High Levels of Gene Flow Between Breeds of Domestic Sheep from Asia and Europe. J Hered 96(5): 494–501.JRS MeadowsK. LiJ. KantanenM. TapioW. Sipos2005Mitochondrial Sequence Reveals High Levels of Gene Flow Between Breeds of Domestic Sheep from Asia and Europe.J Hered96(5)494501

31. Meadows JRS, Cemal I, Karaca O, Gootwine E, Kijas JW, et al. (2007) Five Ovine Mitochondrial Lineages Identified From Sheep Breeds of the Near East. Genetics 175: 1371–1379.JRS MeadowsI. CemalO. KaracaE. GootwineJW Kijas2007Five Ovine Mitochondrial Lineages Identified From Sheep Breeds of the Near East.Genetics17513711379

32. Mannen H, Kohno M, Nagata Y, Tsuji S, Bradley DG, et al. (2004) Independent mitochondrial origin and historical genetic differentiation in North Eastern Asian cattle. Mol Phylogenet Evol 32: 539–544.H. MannenM. KohnoY. NagataS. TsujiDG Bradley2004Independent mitochondrial origin and historical genetic differentiation in North Eastern Asian cattle.Mol Phylogenet Evol32539544

33. Freeman AR, Meghan CM, Machugh DE, Loftus RT, Achukwi MD, et al. (2004) Admixture and diversity in West African cattle populations. Mol Ecol 13: 3477–3487.AR FreemanCM MeghanDE MachughRT LoftusMD Achukwi2004Admixture and diversity in West African cattle populations.Mol Ecol1334773487

34. Fernández H, Hugues S, Vigne JD, Helmer D, Hodgins G, et al. (2006) Divergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas. Proc Natl Acad Sci USA 103: 15375–15379.H. FernándezS. HuguesJD VigneD. HelmerG. Hodgins2006Divergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas.Proc Natl Acad Sci USA1031537515379

35. Vigne J.-D, et al. pp. 52–75.J.-D Vigne5275in Archaeozoology of the Near East IV, Proc. 4th int. Symp. Archaeozoology of Southwestern Asia and adjacent areas (ASWA; Paris, June 1998) M. Mashkour, A. M. Choyke, H. Buitenhuis, F. Poplin, Eds. (Archaeological Research and Consultancy, Groningen, 2000), vol. Publicaties 32. in Archaeozoology of the Near East IV, Proc. 4th int. Symp. Archaeozoology of Southwestern Asia and adjacent areas (ASWA; Paris, June 1998) M. Mashkour, A. M. Choyke, H. Buitenhuis, F. Poplin, Eds. (Archaeological Research and Consultancy, Groningen, 2000), vol. Publicaties 32.

36. Tuñon M, Gonzalez JP, Vallejo M (1989) Genetic relationships between 14 native spanish breeds of goat. Anim Genet 20: 205–212.M. TuñonJP GonzalezM. Vallejo1989Genetic relationships between 14 native spanish breeds of goat.Anim Genet20205212

37. Cañon J, García D, García-Atance M, Obexer-Ruff AG, Lenstra J, et al. (2006) Geographical partitioning of goat diversity in Europe and the Middle East. Anim Genet 37: 327–334.J. CañonD. GarcíaM. García-AtanceAG Obexer-RuffJ. Lenstra2006Geographical partitioning of goat diversity in Europe and the Middle East.Anim Genet37327334

38. Pariset L, Cappuccio I, Ajmone Marsan P, Dunner S, Luikart G, et al. (2006) Assessment of population structure by single nucleotide polymorphisms (SNPs) in goat breeds. J Chromatography B. L. ParisetI. CappuccioP. Ajmone MarsanS. DunnerG. Luikart2006Assessment of population structure by single nucleotide polymorphisms (SNPs) in goat breeds.J Chromatography B

39. Taberlet P, Ajmone-Marsan P, Valentini A, Rezaei HR, Naderi S, et al. (In Press) Are cattle, sheep, and goats endangered species? Mol Ecol. P. TaberletP. Ajmone-MarsanA. ValentiniHR RezaeiS. NaderiIn PressAre cattle, sheep, and goats endangered species?Mol Ecol

40. Luikart G, England P, Tallmon D, Jordan S, Taberlet P (2003) The power and promise of population genomics: from genotyping to genome typing. Nat Rev Genet 4: 981–994.G. LuikartP. EnglandD. TallmonS. JordanP. Taberlet2003The power and promise of population genomics: from genotyping to genome typing.Nat Rev Genet4981994

41. Shackleton DM (1997) Wild sheep and goats and their relatives: Status survey and conservation action plan for Caprinae. Gland, Switzerland: IUCN. DM Shackleton1997Wild sheep and goats and their relatives: Status survey and conservation action plan for Caprinae.Gland, SwitzerlandIUCN390

42. Parma P, Feligini M, Greeppi G, Enne G (2003) The complete nucleotide sequence of goat (Capra hircus) mitochondrial genome. Goat mitochondrial genome. DNA Seq 14(3): 199–203.P. ParmaM. FeliginiG. GreeppiG. Enne2003The complete nucleotide sequence of goat (Capra hircus) mitochondrial genome. Goat mitochondrial genome.DNA Seq14(3)199203

43. Kumar S, Tamura K, Nei M (2004) MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignement. Brief Bioinform 5: 150–163.S. KumarK. TamuraM. Nei2004MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignement.Brief Bioinform5150163

44. Yang Z (1997) PAML: a program package for phylogenetic analysis by maximum likelihood. CABIOS 13: 555–556.Z. Yang1997PAML: a program package for phylogenetic analysis by maximum likelihood.CABIOS13555556

45. Swofford DL (2003) PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods), Version 4. Sunderland, Massachusetts: Sinauer Associates. DL Swofford2003PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods), Version 4.Sunderland, MassachusettsSinauer Associates

46. Excoffier L, Laval G, Schneider S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online 1: 47–50.L. ExcoffierG. LavalS. Schneider2005Arlequin ver. 3.0: An integrated software package for population genetics data analysis.Evolutionary Bioinformatics Online14750

47. Schneider S, Excoffier L (1999) Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics 152: 1079–1089.S. SchneiderL. Excoffier1999Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA.Genetics15210791089

48. Kuhner M (2006) LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics 22: 768–770.M. Kuhner2006LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters.Bioinformatics22768770