A Dig into the Past Mitochondrial Diversity of Corsican Goats Reveals the Influence of Secular Herding Practices

Abstract

The goat (Capra hircus) is one of the earliest domesticated species ca. 10,500 years ago in the Middle-East where its wild ancestor, the bezoar (Capra aegagrus), still occurs. During the Neolithic dispersal, the domestic goat was then introduced in Europe, including the main Mediterranean islands. Islands are interesting models as they maintain traces of ancient colonization, historical exchanges or of peculiar systems of husbandry. Here, we compare the mitochondrial genetic diversity of both medieval and extant goats in the Island of Corsica that presents an original and ancient model of breeding with free-ranging animals. We amplified a fragment of the Control Region for 21 medieval and 28 current goats. Most of them belonged to the A haplogroup, the most worldwide spread and frequent today, but the C haplogroup is also detected at low frequency in the current population. Present Corsican goats appeared more similar to medieval goats than to other European goat populations. Moreover, 16 out of the 26 haplotypes observed were endemic to Corsica and the inferred demographic history suggests that the population has remained constant since the Middle Ages. Implications of these results on management and conservation of endangered Corsican goats currently decimated by a disease are addressed.

About the Authors:

Sandrine Hughes

* E-mail: [email protected]

Affiliation: Paléogénomique et Evolution Moléculaire, Institut de Génomique Fonctionnelle de Lyon, Université de Lyon, Université Lyon 1, CNRS UMR 5242, INRA, Ecole Normale Supérieure de Lyon, 46 allée d'Italie, 69364 Lyon Cedex 07, France

Helena Fernández

Affiliation: Laboratoire d'Ecologie Alpine, CNRS UMR 5553, Université Joseph Fourier, B.P. 53, 38041 Grenoble Cedex 9, France

Thomas Cucchi

Affiliations Centre National de la Recherche Scientifique, UMR 7209, Muséum National d'Histoire Naturelle, «Archéozoologie, Archéobotanique: Sociétés, Pratiques et Environnements», Département “Ecologie et Gestion de la Biodiversité” CP 56, 75005 Paris, France, Department of Archaeology, University of Aberdeen, Aberdeen, United Kingdom

Marilyne Duffraisse

Affiliation: Paléogénomique et Evolution Moléculaire, Institut de Génomique Fonctionnelle de Lyon, Université de Lyon, Université Lyon 1, CNRS UMR 5242, INRA, Ecole Normale Supérieure de Lyon, 46 allée d'Italie, 69364 Lyon Cedex 07, France

François Casabianca

Affiliation: Institut National de la Recherche Agronomique, UR 045 Laboratoire de Recherches sur le Développement de l'Elevage, Quartier Grossetti, 20250 Corte, France

Daniel Istria

Affiliation: Laboratoire d'Archéologie Médiévale Méditerranéenne, CNRS UMR 6572, 5 rue du château de l'Horloge, BP 647, 13094 Aix-en-Provence, France

François Pompanon

Affiliation: Laboratoire d'Ecologie Alpine, CNRS UMR 5553, Université Joseph Fourier, B.P. 53, 38041 Grenoble Cedex 9, France

Jean-Denis Vigne

Affiliation: Centre National de la Recherche Scientifique, UMR 7209, Muséum National d'Histoire Naturelle, «Archéozoologie, Archéobotanique: Sociétés, Pratiques et Environnements», Département “Ecologie et Gestion de la Biodiversité” CP 56, 75005 Paris, France

Catherine Hänni

Affiliation: Paléogénomique et Evolution Moléculaire, Institut de Génomique Fonctionnelle de Lyon, Université de Lyon, Université Lyon 1, CNRS UMR 5242, INRA, Ecole Normale Supérieure de Lyon, 46 allée d'Italie, 69364 Lyon Cedex 07, France

Pierre Taberlet

Affiliation: Laboratoire d'Ecologie Alpine, CNRS UMR 5553, Université Joseph Fourier, B.P. 53, 38041 Grenoble Cedex 9, France

Introduction

In the Near East, cradle of the domestication process, goats were among the first to be domesticated ca. 10,500 years ago [1]–[4]. Several thousands years later, domestic goats (Capra hircus) were dispersed beyond the natural distribution of its wild ancestor (Capra aegagrus). They spread in Anatolia and Europe (starting from 8,800 calBP) throughout the Neolithic dispersal, along with pigs, cattle and sheep [5], [6]. Today, goats are present all over the world with more than 867 million of individuals [7]. In order to better assess the historical processes of domestication, goat mitochondrial genetic diversity has been largely studied across the old world (Europe, Asia, and more recently Africa). It is structured in six different haplogroups A, B, C, D, F and G [3], [8]–[16], with more than 90% of goats solely from the A haplogroup [14]. Moreover, a very weak phylogeographic structure is observed at the continent scale [8] contrary to other domestic species, such as sheep. Thus, A and C haplogroups have a worldwide distribution although B is mostly present in Asia. Some genetic structure is suspected however at a more restricted geographical scale [16] and some haplogroups, such as G and F, are now restricted to small regions (Middle-East [14] and Sicily [12] respectively).

If many studies have been dedicated to the characterization of the mainland diversity few were pursued on islands [12], [17]. However, large Mediterranean islands are of particular importance in the description of the genetic diversity of domestic species since they are considered as biodiversity hot spots, have a high degree of endemism and present a reservoir of cultural practices that have disappeared from the Mainland. In the case of goats, the study of genetic diversity on large Mediterranean islands is highly relevant for several reasons. First, goats are found on most of these islands from the beginning of the Neolithic diffusion and can serve as testimony of this spread [18]. Second, imported domestic goats, are physically isolated from their Mainland relatives. Third, breeding and husbandry practices on the Mediterranean islands are usually different to those on the continent because islands present large but restricted areas and preserve traditional practices mentioned previously. From these observations, we expect that mitochondrial diversity observed on islands can present evidence of historical events or ancient diffusion that would be lost elsewhere, such as the presence of the F haplogroup in Sicily that is unique outside the wild ancestor's area of distribution [3], [12], [14].

Here we characterize goat mitochondrial DNA (mtDNA) diversity through time by studying modern, but above all historical, goat populations from Corsica. We compared them with current continental or island breeds of the Mediterranean Basin to document the microevolution and the influence of insularity and husbandry practices on their genetic diversity. Corsica is an 183 km long and 83 km wide island located in the Northwestern part of the Mediterranean Basin. It formed a unique block with another island, Sardinia, during most of the Pleistocene until the land masses split approximately around 11,000 years ago. The presence of domestic Caprinae is attested in Corsica from the beginning of the Neolithic, ca. 7,600 calBP [19]. The herding system in Corsica, although close to other free-ranging and seasonal transhumance characteristic of other Mediterranean islands, displays a very interesting peculiarity namely the “wandering” of flocks in the mountains for weeks, between the end of the lactation period and the beginning of the births, under very loose surveillance from the herder [20]. These free-ranging practices in Corsica are not recent as Polybius had already mentioned them in his book XII during the II^nd century Before Christ. Nonetheless, being able to ascertain that transhumance and free ranging were common herding practices in Corsica since the Neolithic is a difficult task for zooarchaeologists [19]. Goats are known to be hardy, tough and able to adapt to very difficult habitats compared to other livestock. More than 200,000 goats were still present in Corsica less than 80 years ago but this number has decreased to 30,000 individuals during the last decades. Nevertheless, goats have retained a particular status in Corsica where pastoralism is still strongly established [21]. Besides, a Corsican breed has been recently recognized by the French CNAG (Commission Nationale d'Amélioration Génétique depending on Ministry of Agriculture) in 2003 and by decree of the French Ministry of Agriculture in 2007. It is a dairy breed, with relatively long hairs, either of uniform colour or multicoloured, characterized by its rustic character and ability to adapt to difficult grounds. This was possible because efforts were made to protect and promote local breeds by avoiding mixings with commercial breeds [22]. However, the traditional Corsican breed is now endangered because the traditional husbandry system seems difficult to maintain. The strains on this traditional system are not only economic but also social. Traditional goat husbandry practices implicate daily mobility from the herder, which most of the new generation does not wish to pursue [23]. More recently, the occurrence of Johne's disease (i.e., paratuberculosis) presents another very serious threat and has decimated flocks.

This study aims to (i) better characterize the current and past (medieval) mitochondrial genetic diversity of Corsican goats using a Control Region (CR) fragment; (ii) bring information about goat dispersal in the Neolithic by testing the congruence between scenario proposed and data observed on this island; (iii) explain the maintenance of endemic variability in Corsica; (iv) discuss implications for conservation of the Corsican breed. For these purposes, we gathered samples of 28 present-day individuals and 29 bones dated from the Middle Ages (XII^th and XIV^th centuries).

Methods

The archeological site of Rostino

The castrum (strong hold) of Rostino [24] is situated in the North East of Corsica (Figure 1). Occupied between the XII^th and XIV^th centuries AD, this late medieval site has yielded the largest assemblage of Caprinae in Corsica in a good state of preservation [25], which is rare in the acidic soil of Corsica [19]. Among the domestic species of the XIV^th century deposit, Caprinae represent more than 70% of the identified mammal remains [25]. The economy of the castrum relies on specialized caprine exploitation where the production of sheep and goat complement each other: sheep for meat production and goat for milk and hair production [25]. During the XIV^th century, the caprine exploitation specialized in goats with Capra hircus representing more than 70% of the total caprine remains [25]. This large amount of late medieval Corsican Capra hircus represents a great opportunity to investigate the consequence of the secular herding practices and the selective choices made by herders to renew their flock given the genetic diversity of goats in large Mediterranean islands.

[Figure omitted. See PDF.]

Figure 1. Ancient and present-day sampling of goats in Corsica.

The medieval samples are all located in the archeological site of Rostino, in the northern part of the island. The present-day samples come from 5 different localities identified by numbers on the map.

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

Archeological samples

We analysed 29 bone fragments from Rostino: 17 were excavated from a XII^th century deposit and 12 from a deposit dated to the XIV^th century. These bones have been identified as Capra hircus using morphoscopic criteria on both dental [26], [27] and appendicular [28]–[30] characters (Table S1). There is no ambiguity about the origin of bones from domestic animals as the wild ancestor (Capra aegagrus) has never been present in Europe [31]. According to the type of bone (mandible, radius, humerus …), the laterality (left, right) and detailed information from the excavation, we were able to clearly identify 25 different individuals from the 29 fragments. The four other fragments very probably came from previously identified individuals (Table S1). Two samples were subsequently identified as sheep by the molecular analyses. This is not surprising as the inter-specific distinction between sheep and goat on fragmented bones cannot be ascertained with a 100% reliability [29]. Molecular analyses have precisely proven to be useful in this case [32] as has, more recently, the analysis of collagen by mass spectrometry [33].

Ancient DNA extracts and PCR amplifications

Retrieval of the DNA preserved in bones was performed in ultra clean rooms dedicated to ancient DNA experiments (French National Platform of Paleogenetics PALGENE, CNRS, ENS Lyon). No more than 4 capra samples were treated in the same session of DNA extraction along with a bone from another species (cervids, ursids) and a blank control [34], [35]. 100 to 500 milligrams of each bone fragments were reduced in powder and suspended in 5 or 10 ml of EDTA buffer as described in [36]. We extracted the DNA using one, or both, of the two following protocols: a classical phenol-chloroform approach [36], [37] or direct purification using Qiaquick column (Qiagen kit) [34]. 25 out of the 29 samples were independently re-extracted in a second laboratory dedicated to ancient DNA in another city, Grenoble, using the Qiaquick protocol.

A 130 bp fragment of the CR (HVI) was amplified with the CapFII and CapRII primers with conditions identical to those described in [38]. At least 2 or 3 independent positive amplifications per sample were obtained, cloned and sequenced following protocols described in [34], [35]. The final sequence of one individual was obtained by making the majority-rule consensus of all consensus of all different clones obtained from each of the independent amplification. More than 80% of the differences observed between clones were G to A or C to T punctual substitutions. This result is consistent with ancient DNA degradation profiles where deamination of cytosines is known as the major factor of artifactual substitutions [39].

Present-day Corsican goat sequences

We sampled 28 goats in 5 different localities in Corsica (7 from Moltifao, 6 from Tralonca, 11 from Corte, 3 from Altiani, 1 from Quenza; Figure 1, Table S1). 14 sequences were already published [14] and 14 were produced for this study. To amplify the CR fragment we used the same primers (CapFII and CapRII) than for medieval sequences with slightly modified conditions: 40 cycles were performed instead of 50–60 for ancient DNA and products were directly sequenced on both DNA strands.

Sequence analyses

All the sequences obtained were aligned (Seaview v4 [40]) and the primers removed leading to a fragment of 130 base pairs. Four data sets were constituted: (i) all medieval sequences, (ii) the XII^th century sequences, (iii) the XIV^th century sequences and (iv) the present-day sequences.

Firstly, the mitochondrial haplogroup of each Corsican sequence found was determined by performing phylogenetic analyses. The different Corsican haplotypes were analyzed together with 20 haplotypes of known haplogroups. These 20 haplotypes corresponded to the 20 reference sequences that were different for the 130 bp fragment under study, in the dataset selected by Naderi et al. [14] to represent the worldwide variability of the whole HVI-control region (558 bp). Identical sequences, or haplotypes, were identified in the Corsican dataset by using Fabox [41]. After estimating the better model of evolution using jModeltest program [42] and the Akaike Information Criterion (AIC), we performed Bayesian analyses (BA) with MrBayes v3.1.2 [43], [44] (independently confirmed by Maximum Likelihood analyses, not shown). The parameters used for BA were the following: GTR+I+G (nst = 6 and rates = invgamma), 5,000,000 generations sampled every 1000^th generation, 4 chains, a burn-in period of 500 trees (i.e. 10% of generations) visually confirmed using Tracer v1.4.1 (developed by Rambaut A and Drummond A; available from http://beast.bio.ed.ac.uk/Tracer), allcompat option. Two independent runs were performed with an average standard deviation of split frequencies at completion of 0.006598. The average values obtained for alpha and proportion of invariable sites parameters were 0.287 and 0.323 respectively.

Secondly, we assessed the relationships between only medieval or all Corsican sequences, by performing median-joining networks, using the Network software ([45], available at fluxus-engineering.com) with default parameters. Network Publisher was used to manipulate the networks. To compare the Corsican mitochondrial genetic diversity with the Mediterranean or worldwide diversity, we defined supplementary datasets, one by haplogroup observed, gathering all the sequences published and covering the 130 bp fragment. These datasets combined with the Corsican sequences were used to draw median-joining networks (references are given in the legends of Figures S1 and S2).

Thirdly, to assess the past demographic history of the Corsican goats we performed a Bayesian Skyline Plot (BSP) using BEAST v.1.5.4 software [46]. All the Corsican sequences were used and average tip dates were given for all medieval sequences: 850 years BP for the XII^th century and 650 years BP for those of the XIV^th century. BEAUti v1.5.3. was used to build the xml file by using the following parameters: HKY+I+G4 (best model for this dataset assessed by jModeltest and AIC criteria [42]); uncorrelated lognormal relaxed clock; 5 groups and 100 millions of iterations with parameters saved every 10 000 iterations; Burn-in: the first 10% were discarded. The results of 4 independent runs were analysed and the Bayesian Skyline Plot reconstructed with Tracer v.1.4.1 (Figure S3).

Finally, to describe the genetic diversity observed for the Corsican goats, we computed different classical parameters using either DnaSP v5 [47] or Arlequin v3.5.1.2 [48] including haplotype diversity and frequencies, sequence diversity, pairwise comparisons, and population comparisons (F_ST, Fu tests). We compared the different Corsican datasets with each other but also with other datasets corresponding to mainland or island populations (see [16] for accession numbers and geographical origin). We focused in particular on Sardinia's island (75 sequences, accession numbers FJ571522 to FJ571596, [17]; Corsica's closest island) and Portugal (the biggest dataset generated for mainland goats, 288 sequences, accession numbers AY961629 to AY961916, [11]). We also considered 4 different datasets (see Supporting Information for details and references) corresponding to the southern or the northern area of the Mediterranean Sea, the Mediterranean islands (Sicily and Sardinia), and finally 8 Neolithic goats of Baume d'Oullens [38]. To reduce possible biases due to large differences in the size of the datasets (low number of Corsican sequences), we randomly sampled 49 sequences of the non-Corsican sequences and repeated this operation at least three times. The genetic parameters estimated by DnaSP v5 [47] or Arlequin v3.5.1.2 [48] were then computed on these resampled datasets of equal size.

Results

Medieval goat samples

We analyzed 29 bones and obtained reproducible and congruent sequences of Capra hircus for 25 of them (Table S1). We are confident that these results are authentic as we obtained the same results in both laboratories. We also took the ancient DNA precautions recommended by the community as we are used to do (e.g. [34], [35], [38], [49], [50]). When different bones were supposed to be from a single individual (Ro-5, Ro-10 and Ro-22), systematically we obtained the same sequence confirming the first assessments. Finally, we determined 21 sequences coming from different individuals: 10 dated to the XII^th century and 11 dated to the first half of the XIV^th century. All medieval sequences were from the A haplogroup as shown by the phylogenetic analysis (Figure 2). The six haplogroups appeared monophyletic and were supported by posterior probabilities (pp) higher than 0.9 for 3 of them B, C, G (F is not concerned as only one sequence is used). The A haplogroup that had the highest number of sequences and that was the more diverse received the lowest support (pp<0.5). Nevertheless, the clustering of the newly determined sequences inside the A haplogroup raises no doubt as confirmed by subsequent network analyses (see also mismatch distributions, Figure S4). Substantial diversity is observed among medieval haplotypes as seen on the network (Figure 3A). Among the 21 sequences, we detected 14 unique haplotypes (Table S1) with haplotype Ha 04 being the most frequent (5 individuals). Six out of the 14 haplotypes have never been described before. On average, the mean number of pairwise differences observed between sequences reaches 4.65±2.37 (Table 1, Figure S4). The diversity appeared not significantly different between the XII^th century with 9 haplotypes for 10 sequences (13 polymorphic sites) and the XIV^th century with 8 haplotypes for 11 sequences (15 polymorphic sites). According to the network performed on all medieval sequences (14 haplotypes, 19 polymorphic sites), only 3 haplotypes were shared between both periods (Figure 3A). Two of them had a central position in the network (Ha 04 and 09; Figure 3A).

[Figure omitted. See PDF.]

Figure 2. Mitochondrial diversity of Corsican goats compared to world-wide goats diversity.

Ancient (circles) and current (stars) Corsican sequences have been used with 20 sequences of reference corresponding to all haplogroups described up to now [14] to reconstruct a Bayesian phylogenetic tree using a GTR+I+G model of evolution (see text). Only a single sequence by haplotype has been used, so haplotypes in common between medieval and today periods are indicated by both a circle and star. Only posterior probabilities higher than 0.6 are indicated.

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

[Figure omitted. See PDF.]

Figure 3. Networks generated with CR sequences of Corsican goats from A) medieval samples only or B) combined with present-day samples.

The circles are proportional to the number of sequences obtained and colours indicated the period: yellow corresponds to XII^th century, red to XIV^th century, orange to all medieval and blue to present-day sequences. See Table 1 for details about haplotypes and archeological samples. Haplotypes indicated in italics are those only observed in Corsica up to now. Haplotype Ha 26 is the only one from the haplogroup C and differs from the closer haplotype by 17 positions that is materialized by a square with this number.

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

[Figure omitted. See PDF.]

Table 1. Genetic diversity parameters.

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

Present-day goat samples

28 individuals from 5 different localities were studied (Table S1 for details). Two different mitochondrial haplogroups were observed (A and C; Figure 2) with a higher proportion of A sequences (92.8% i.e. 26 out of 28 sequences). Considering all A and C sequences, the mean number of pairwise differences reaches 7.01±3.39 (29 polymorphic sites; Figure S4). However, when only sequences of the A haplogroup were considered, this value drops to 4.90±2.47 (19 polymorphic sites, Table 1) which is close to the one observed for A medieval goats. The two sequences from the C haplogroup came from the same locality, Tralonca (Figure 1), and shared the same haplotype Ha 26 that has not been described elsewhere (Figure 3B). 15 haplotypes were obtained for the 26 sequences of the A haplogroup among which, 10 were only observed in Corsica. The most frequent haplotypes were, like for the late medieval goats, Ha 04 and Ha 09, with 4 individuals from 3 different localities in both cases (Figure 3B).

Comparison of Corsican goat diversity through time

No significant difference was observed between the goats of the XII^th and the XIV^th centuries (non significant F_ST-value; Table S2). The difference became significant when medieval and present goats were compared (0.036, p-value 0.027). This was due to the presence of the C haplogroup in the present-day sequences since the test was no more significant when only A sequences were taken into account (0.027, p-value 0.099). Finally, there were no significant differences between the medieval and present-day sequences, with or without considering the C sequences. Medieval and present-day goats shared 4 haplotypes (Ha 04, 06, 07 and 09, Figure 3B) but only one was specific to Corsica (Ha 06). The 3 remaining haplotypes were generally frequent in other populations.

No significant changes were observed in the demographic history of the Corsican population using the Bayesian Skyline Plot (Figure S3). Neither a sign of expansion nor of a crash were observed from the medieval times to date, as rather a constant population size pattern was obtained. This is congruent with the Tajima's D values obtained (Table 1) but not with all the Fu's Fs values computed. However, this latter parameter can reflect other factors than population growth (selection, bottleneck, … [16]).

Comparison of Corsican goat diversity with other geographical places

We compared the diversity observed in Corsica with other Mediterranean island populations and mainland breeds. The 75 sequences of Sardinia were clustered in 46 haplotypes that were all from the A haplogroup [17]. Similarly, the 288 Portuguese sequences were represented by 104 haplotypes, from which only one was from the C haplogroup and all the others from the A haplogroup [11]. By expanding the comparison to larger or other areas, we observed that the C haplogroup is described only in Europe (in present-day Northern Mediterranean area and already at the Neolithic time, Table S3 and not shown). Both the medieval and present-day Corsican groups appeared significantly different from all other groups whatever the resampled datasets taken into account (Table 1, Tables S2 and S3). Similarly, all the non-Corsican groups of goats also appeared significantly different from each other (data not shown). However, four different haplotypes (Ha 05, 08, 09 and 10) were shared between Corsica and Sardinia islands and six (Ha 07, 08, 09, 10, 12 and 21) between Corsican and Portuguese goats.

The median-joining network performed on 39 worldwide sequences of the C haplogroup (Figure S1) revealed a classical expansion structure with European sequences on one side and Asian sequences on the other (with the single exception of one Swiss haplotype; [8]). As expected, the sequences obtained from Neolithic goats of the archeological site of Baume d'Oullens in France [38], among the first goats to have been diffused in Europe, appeared in the central part of the European cluster. The Corsican haplotype showed the highest number of substitutions with this central node (3 substitutions). The same analysis performed on the sequences from the A haplogroup restricted to the Mediterranean area, showed a more complex history with no clear emerging pattern (Figure S2).

Discussion

DNA preservation in medieval samples

From the 29 archeological bones dated back to the XII^th and XIV^th centuries, we obtained 21 sequences from different ancient goat individuals: 10 from the XII^th century deposit and 11 from the XIV^th century deposit, which represents a surprisingly high DNA preservation for remains in Corsica. Such a good preservation of the DNA, which is here correlated to the good preservation of the bones themselves like in most late Medieval Corsican sites [19], [51], is probably due to the recent age of the site.

Characterization through time of mitochondrial genetic diversity in Corsican goats

The comparison of Middle Ages and present mitochondrial diversity was carried out using 28 present-day goats from five different localities. This comparison may be slightly biased by differences existing in the time span and geographic distribution for either ancient or modern samples. Indeed, the sampled geographic area is larger for the modern goats (Figure 1) whereas the time span is longer for the ancient samples.

Our results tend to prove that the mitochondrial diversity of Corsican goats has remained relatively constant since the Middle Ages. Moreover, we detect no significant demographic changes (F_ST, BSP and Table 1) or decrease of genetic diversity. The only difference between both periods is the occurrence of two C haplotypes in the present-day samples, all other goats belonging to the A haplogroup. However, given a dataset of 21 medieval individuals and assuming a constant frequency of the C haplogroup in Corsica (2/28 = 0.0714), we have a 21% of chance of having missed this haplogroup in the Middle Ages sampling.

The occurrence of the C haplogroup in Corsican goats in the context of the Neolithic diffusion

The presence of the A and C haplogroups in Corsica is in agreement with the goat mitochondrial variability observed in European countries. Most of European goats are from the A haplogroup, with C haplotypes found at a rare frequency in Switzerland, Portugal, Spain and Slovenia [3], [8], [11], [14]. This is also consistent with previous paleogenetic studies that already detected both haplogroups in goats from Southern France in the early Neolithic period [38]. The median-joining network performed on all the C haplotypes found worldwide (Figure S1) gave results congruent with the diffusion of goats in the Neolithic. Indeed we observed: i) a clear separation of the European haplotypes from the Asian ones with only one exception, a shared haplotype between China and Switzerland; ii) a star-like pattern for European haplotypes suggesting a population expansion with a central position for the Neolithic haplotypes [38]; iii) a divergent haplotype for Corsica compatible with a subsequent isolation.

Without a doubt, the origin of the C haplogroup in Europe can be traced back to the Neolithic spread where its frequency was probably higher than the one observed now [38]. Today, the C haplotypes are relatively scarce in goats and appear more like reminiscent testimonies of this first diffusion. Interestingly, not one has been detected yet in the Southern Mediterranean area or in other Mediterranean islands than Corsica (see Supporting Information). According to the few data we have, it is difficult to conclude on the ancient origin of the Corsican C haplotype or a more recent origin linked to subsequent exchanges with the Northern Mediterranean mainland. However, its position in the network is compatible with the first explanation. Further ancient DNA studies, for a larger area and older period, would be very interesting to highlight this question.

Variability in Corsica vs other Mediterranean areas

The relative stability through time of goat diversity in Corsica could be explained by regular importation of goats from other continental areas or islands, e.g. for commercial trade, as many contacts by sea have been reported during the last centuries in the Mediterranean area [52]. However, a striking point is that the diversity observed in Corsica is substantial for both periods (Middle Ages and present-day, Figure 3A and Figure 3B respectively) and differs from that of other places (Table 1, data not shown); only about half of the Corsican haplotypes (14/26) are shared with goats from other geographical regions despite our study focused on a short CR fragment (130 bp). This result is not unexpected however. Previous studies - usually targeting the 480 bp fragment first described in [8] and covering our shorter fragment - have already shown that goats were more polymorphic than other livestock (cattle, sheep, pig) on the CR [11], [15], [53]. Moreover, the specific analysis of the A haplogroup (more than 90% of the modern goats, [14]) confirmed high haplotype diversities for 20 populations/countries with values close to one [16]. Pereira et al. [16] showed a strong correspondence between mitochondrial genetic and geographic distances suggesting that after the initial expansion, differentiation among regions has been established and maintained [16]. A similar conclusion was obtained with large-scale nuclear SNP analyses obtained for 16 breeds of goats [54] and from microsatellites analyses including Corsican and other European goat populations [55]. This seems in agreement with what we observed in Corsica when compared with other Mediterranean populations (F_ST comparisons, Table S2, data not shown).

Implication of the traditional Corsican husbandry in the maintenance of the variability

Except for some haplotypes that are common in many countries (Figure S2) and probably constitute traces of the initial diffusion, 46% of the Corsican haplotypes found have not been previously described elsewhere. The preservation through time of this endemic genetic diversity and its constant level since medieval time could suggest that relatively large effective population sizes have been maintained in domestic goats through exchanges of animals. But ethnographic insights into the herding practices carried out in the Niolu [20] provide another possible or complementary explanation for the preserved genetic diversity of the Corsican goat breed.

Typically in the past for the Mediterranean area, goats were usually moved according to the seasonal changes (transhumance) to gain access to more reliable food. Because fodder resources fluctuate in the wild according to different factors (e.g. annual climatic conditions), the system developed in Corsica has been extensive with goats left free-ranging most of the times although under the careful control of the herders. This particular system led to the characteristics observed in the Corsican goats [21], [56]: i) the high diversity of coat colours is encouraged as individuals can be more easily identified by sight; ii) only the strongest and toughest goats can generally survive in this relatively hostile environment, explaining why the introduction of goats from industrial breeds usually failed; iii) large herds are usually managed to maximize the herd's productivity overall instead of individual productivity. Indeed, selection is not performed to optimize for instance, either milk productivity or the fecundity, but instead to obtain a constant productivity of the herd by year whatever the weather or difficulties encountered [57]. All these points lead to a strategy where the phenotypic diversity, and the underlying genetic diversity, is maximized in order to obtain herds that are highly adaptable, rustic and robust.

In order to do that, herders exert strong selection pressure while forming their herds to gain in productivity but also to favour behavioural traits. Indeed, along with the search for constant productivity, one of the main objectives of herders is to design flocks that will maintain “families” from the same maternal lineage [20]. A young descendant female is usually chosen according to its mother's and grandmother's family, taking into account its desirable productive traits as well as its abilities to endure the tough conditions. Animals coming from the same family/pool will live together and move together more easily and naturally. Such “familial” behaviour is extremely beneficial to the cohesion of herds under free-ranging exploitation and especially during the movements on the pasturelands (transhumance). On the contrary, introduction of new animals from other herds can induce significant disturbance in the herd movements by breaking the cohesion of the group. This practice of herders of course has the effect of reducing the genetic diversity. Hence, to insure the persistence of the genetic diversity and to “change the blood”, “cambià u sangui” in the language of traditional Corsican herders, selected males of these different maternal lineages or “families” will be exchanged between flocks to avoid inbreeding and, without control, renew the pool of the mating males within a seasonal time span. As mitochondria are inherited maternally, this system will naturally lead to maintaining high mitochondrial genetic diversity between herds or “families”. So, the diversity we observed would be not promoted by large herds with many exchanges of females between them but on the contrary by traditional practices. These ones rely on very few introductions of females coming from other herds to keep the kinship within the herd as tight as possible to reinforce cohesive behaviour during the ranging but also on a mixing of the genetic diversity through the exchanges of males only.

Such entangled practices in the Corsican husbandry system can explain why we do not observe change of the mitochondrial diversity in Corsican goats since the late Medieval period, and probably earlier if we had been able to investigate the genetic diversity of goats from earlier periods. It's thus highly probable that sustainable husbandry practices in today's Corsica that are so well adapted to their environment are, at least partly, the result of practices over millenaries contributing to the maintenance of a relatively high level of genetic diversity. Finally, analyses of additional samples and genetic markers coupled with population simulations using varying population genetic parameters should help to test between both hypotheses, regular importation of goats or millenaries husbandry practices.

Towards a protection of the Corsican goat breed

In this study, we observed that the mitochondrial diversity of goats in Corsica Island has been maintained since the Middle Ages to date. In a time where rustic breeds are endangered and industrial breeds tend to reduce the genetic resources [22], Corsican goats, by the model of husbandry used and the high diversity observed constitute an interesting pool to preserve for the future management of domestic genetic resources. All the actors concerned by the Corsican goat (regional association in charge of its management, public authorities, researchers and extension services) should so pay attention to make their preservation successful all the more because paratuberculosis has started to decimate flocks.

Ethic statement

The medieval bones of goats were excavated from the Rostino archeological site. Daniel Istria, in charge of the excavation, authorized their analyses. Tissue samples were collected in Corsica in the framework of the ECONOGEN project (European Union contract QLK5-CT-2001-02461), following the European ethical rules implemented in all European projects.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Network generated with CR sequences of the C haplogroup (130 bp). Only domestic goats have been considered (35 sequences; Luikart et al. 2001 [8], Sultana et al. 2003, Joshi et al. 2004, Pereira et al. 2005 [11], Chen et al. 2005, Liu et al. unpublished, Naderi et al. 2007 [14]). Neolithic sequences coming from the archeological site of Baume d'Oullen (Fernández et al. 2006 [38]) were also taking into account (4 sequences).

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

(DOC)

Figure S2.

Network generated with CR sequences of the A haplogroup (130 bp) for goats coming from the Mediterranean Sea around Corsica. 584 sequences were used for the median-joining network analysis coming from: North Mediterranean area (Italy, France, Spain), Mediterranean Islands (Malta, Sicily, Sardinia), South Mediterranean area (Morocco, Algeria, Tunisia). See Pereira et al. 2009 [16] for the accession numbers of the sequences used and their geographical origin. Positions were weighted inversely to the number of mutations observed by position on a first run-test.

https://doi.org/10.1371/journal.pone.0030272.s002

(DOC)

Figure S3.

Bayesian Skyline Plot. The analyses were performed using the 49 Corsican sequences with date-tips for the medieval sequences (see text for details). X-axis: Time in years; Y-axis: Population size (Neτ) in log-scale. Mean is plotted with the 95% HPD. 3 runs of 100 M of iterations were performed (ESS >200 for all parameters).

https://doi.org/10.1371/journal.pone.0030272.s003

(DOC)

Figure S4.

Mismatch distributions for medieval and extant Corsican goat populations. The numbers of pairwise differences are given on the x-axis and their frequency on the y-axis.

https://doi.org/10.1371/journal.pone.0030272.s004

(PDF)

Table S1.

Corsican goats sampling and mitochondrial genotyping results. A star indicates archeological samples for which the molecular identification was Ovis aries and a dash when no amplification was obtained.

https://doi.org/10.1371/journal.pone.0030272.s005

(DOC)

Table S2.

Population pairwise F_ST. Pairwise difference was used as distance method. F_ST values are given at upper right corner and corresponding p-values at bottom left. Significant p-values are highlighted in color (orange for intra-Corsica comparisons, pink for inter comparisons). NS: Non-significant, * p-value between 0.05 and 0.01, ** p-value between 0.01 and 0.001, *** p-value < to 0.001.

https://doi.org/10.1371/journal.pone.0030272.s006

(DOC)

Table S3.

Genetic diversity parameters. Comparison for A and C haplogroups between Corsican and other Mediterranean or Portuguese datasets (see Figures S1, S3 and text).

https://doi.org/10.1371/journal.pone.0030272.s007

(DOC)

Acknowledgments

We thank Hamid Reza Rezaei and Saeid Naderi for having given us details about samples of Corsican goats they have published. We also thank members of the paleogenetics lab and the Palgene Platform (CNRS, ENS Lyon) for help during the experiments. We are also grateful to Javier Oliver and Benjamin Gillet for critical reading of the manuscript. Finally, we express gratitude to Joanne Burden who kindly proofread the English.

Author Contributions

Conceived and designed the experiments: PT HF TC SH CH. Performed the experiments: MD HF. Analyzed the data: SH TC. Contributed reagents/materials/analysis tools: DI J-DV PT CH. Wrote the paper: SH FP TC. Provided expertise for Corsican husbandry, ancient DNA and archeozoological data: FC CH J-DV.

Citation: Hughes S, Fernández H, Cucchi T, Duffraisse M, Casabianca F, Istria D, et al. (2012) A Dig into the Past Mitochondrial Diversity of Corsican Goats Reveals the Influence of Secular Herding Practices. PLoS ONE 7(1): e30272. https://doi.org/10.1371/journal.pone.0030272

References

1. Zeder MA, Hesse B (2000) The initial domestication of goats (Capra hircus) in the Zagros mountains 10,000 years ago. Science 287: 2254–2257.MA ZederB. Hesse2000The initial domestication of goats (Capra hircus) in the Zagros mountains 10,000 years ago.Science28722542257

2. Peters J, von den Driesch A, Helmer D (2005) The First Steps of Animal Domestication. New Archaeological Approaches. Vigne JDPJ, Helmer D, editors. Oxford, UK: Oxbow Books. J. PetersA. von den DrieschD. Helmer2005JDPJ VigneD. HelmerThe First Steps of Animal Domestication. New Archaeological ApproachesOxford, UKOxbow Books

3. Naderi S, Rezaei HR, Pompanon F, Blum MG, Negrini R, et al. (2008) The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals. Proc Natl Acad Sci U S A 105: 17659–17664.S. NaderiHR RezaeiF. PompanonMG BlumR. Negrini2008The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals.Proc Natl Acad Sci U S A1051765917664

4. Vigne JD (2011) The origins of animal domestication and husbandry: a major change in the history of humanity and the biosphere. C R Biol 334: 171–181.JD Vigne2011The origins of animal domestication and husbandry: a major change in the history of humanity and the biosphere.C R Biol334171181

5. Guilaine J (2003) De la vague à la tombe. J. Guilaine2003De la vague à la tombe.La conquête néolithique de la Méditerranée. Paris. La conquête néolithique de la Méditerranée. Paris.

6. Tresset A, Bollongino R, Edwards CJ, Hughes S, Vigne JD (2009) Early diffusion of domestic bovids in Europe: an indicator for human contacts, exchanges and migrations? In: Hombert J-MadE, F. , editors. The Origin of Man, Language and Languages. Cambridge: McDonald Institute Monograph Series. pp. 69–90.A. TressetR. BollonginoCJ EdwardsS. HughesJD Vigne2009Early diffusion of domestic bovids in Europe: an indicator for human contacts, exchanges and migrations?J-MadE HombertF.The Origin of Man, Language and LanguagesCambridgeMcDonald Institute Monograph Series6990

7. FAOSTAT (2009) Preliminary data for production of Live Animals in 2009 in the World. FAOSTAT website. FAOSTAT2009Preliminary data for production of Live Animals in 2009 in the World. FAOSTAT website.Available: http://faostat.fao.org/ Accessed 2011 May 10. Available: http://faostat.fao.org/ Accessed 2011 May 10.

8. 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 U S A 98: 5927–5932.G. LuikartL. GiellyL. ExcoffierJD VigneJ. Bouvet2001Multiple maternal origins and weak phylogeographic structure in domestic goats.Proc Natl Acad Sci U S A9859275932

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

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

11. Pereira F, Pereira L, Van Asch B, Bradley DG, 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 AschDG 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

12. 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: 376–378.MT SardinaM. BallesterJ. MarmiR. FinocchiaroJB van Kaam2006Phylogenetic analysis of Sicilian goats reveals a new mtDNA lineage.Anim Genet37376378

13. Fan B, Chen S-L, Kijas JH, Liu B, Yu M, et al. (2007) Phylogenetic relationships among Chinese indigenous goat breds inferred from mitochondrial control region sequence. Small Ruminant Research 73: 262–266.B. FanS-L ChenJH KijasB. LiuM. Yu2007Phylogenetic relationships among Chinese indigenous goat breds inferred from mitochondrial control region sequence.Small Ruminant Research73262266

14. Naderi S, Rezaei HR, Taberlet P, Zundel S, Rafat SA, et al. (2007) Large-scale mitochondrial DNA analysis of the domestic goat reveals six haplogroups with high diversity. PLoS One 2: e1012.S. NaderiHR RezaeiP. TaberletS. ZundelSA Rafat2007Large-scale mitochondrial DNA analysis of the domestic goat reveals six haplogroups with high diversity.PLoS One2e1012

15. Wang J, Chen YL, Wang XL, Yang ZX (2008) The genetic diversity of seven indigenous Chinese goat breds. Small Ruminant Research 74: 231–237.J. WangYL ChenXL WangZX Yang2008The genetic diversity of seven indigenous Chinese goat breds.Small Ruminant Research74231237

16. Pereira F, Queiros S, Gusmao L, Nijman IJ, Cuppen E, et al. (2009) Tracing the history of goat pastoralism: new clues from mitochondrial and Y chromosome DNA in North Africa. Mol Biol Evol 26: 2765–2773.F. PereiraS. QueirosL. GusmaoIJ NijmanE. Cuppen2009Tracing the history of goat pastoralism: new clues from mitochondrial and Y chromosome DNA in North Africa.Mol Biol Evol2627652773

17. Vacca GM, Daga C, Pazzola M, Carcangiu V, Dettori ML, et al. (2010) D-loop sequence mitochondrial DNA variability of Sarda goat and other goat breeds and populations reared in the Mediterranean area. Journal of animal breeding and genetics = Zeitschrift fur Tierzuchtung und Zuchtungsbiologie 127: 352–360.GM VaccaC. DagaM. PazzolaV. CarcangiuML Dettori2010D-loop sequence mitochondrial DNA variability of Sarda goat and other goat breeds and populations reared in the Mediterranean area.Journal of animal breeding and genetics = Zeitschrift fur Tierzuchtung und Zuchtungsbiologie127352360

18. Vigne JD (1999) The large “true” Mediterranean islands as a model for the Holocene human impact on the European vertebrate fauna? Recent data and new reflections. In: Benecke N, editor. The Holocene history of the European vertebrate fauna Modern aspects of research Berlin Deutsches Archäologisches Institut, Eurasien-Abteilung (Archäologie in Eurasien, 6). pp. 295–322.JD Vigne1999The large “true” Mediterranean islands as a model for the Holocene human impact on the European vertebrate fauna? Recent data and new reflections.N. BeneckeThe Holocene history of the European vertebrate fauna Modern aspects of research Berlin Deutsches Archäologisches Institut, Eurasien-Abteilung (Archäologie in Eurasien, 6)295322

19. Vigne JD (1988) JD Vigne1988Les Mammifères post-glaciaires de Corse, étude Archéozoologique (XXVIe suppl. à Gallia Préhistoire); CNRS, editor. Paris. Les Mammifères post-glaciaires de Corse, étude Archéozoologique (XXVIe suppl. à Gallia Préhistoire); CNRS, editor. Paris.

20. Ravis-Giordani G (1983) Bergers corses. G. Ravis-Giordani1983Bergers corses.Communautés villageoises du Niolu. Aix-en-Provence: Edisud. Communautés villageoises du Niolu. Aix-en-Provence: Edisud.

21. Bouche R, Aragni CHJ, Bordeaux C (2009) Caprin extensif en Corse: Savoirs durables en quête de développement. Options Méditerranéennes, A 91: 209–213.R. BoucheCHJ AragniC. Bordeaux2009Caprin extensif en Corse: Savoirs durables en quête de développement.Options Méditerranéennes, A91209213

22. Taberlet P, Valentini A, Rezaei HR, Naderi S, Pompanon F, et al. (2008) Are cattle, sheep, and goats endangered species? Mol Ecol 17: 275–284.P. TaberletA. ValentiniHR RezaeiS. NaderiF. Pompanon2008Are cattle, sheep, and goats endangered species?Mol Ecol17275284

23. Pernet F, Lenclud G (1977) Berger en Corse. Essai sur la question pastorale. Grenoble: Presses universitaires de Grenoble. F. PernetG. Lenclud1977Berger en Corse. Essai sur la question pastoraleGrenoblePresses universitaires de Grenoble

24. Istria D (2005) Pouvoirs et fortifications dans le Nord de la Corse, XIe-XIVe siècle. 517 p.D. Istria2005Pouvoirs et fortifications dans le Nord de la Corse, XIe-XIVe siècle.Ajaccio Editions Alain Piazzola. 517 Ajaccio Editions Alain Piazzola.

25. Cucchi T (2003) Production et distribution de biens d'origine animale en contexte seigneurial de Haute-Corse: la faune mammalienne de la salle 1 du castrum de Rostino (XIIIe-XIVe siècles). Bulletin d'Archéologie et d'Histoire de la Corse 1: 70–93.T. Cucchi2003Production et distribution de biens d'origine animale en contexte seigneurial de Haute-Corse: la faune mammalienne de la salle 1 du castrum de Rostino (XIIIe-XIVe siècles).Bulletin d'Archéologie et d'Histoire de la Corse17093

26. Payne S (1985) Morphological distinction between the mandibular teeth of young sheep, Ovis and goats, Capra. Journal of Archaeological Science 12: 139–147.S. Payne1985Morphological distinction between the mandibular teeth of young sheep, Ovis and goats, Capra.Journal of Archaeological Science12139147

27. Helmer D (2000) Discrimination des genres Ovis et Capra à l'aide des prémolaires inférieures 3 et 4 et interprétation des âges d'abattage: l'exemple de Dikili Tash (Grèce). Anthropozoologica 31: D. Helmer2000Discrimination des genres Ovis et Capra à l'aide des prémolaires inférieures 3 et 4 et interprétation des âges d'abattage: l'exemple de Dikili Tash (Grèce).Anthropozoologica31

28. Boessneck J (1969) Oesteological differences beetween sheep (Ovis aries Linne) and goat (Capra hircus Linne). In: D HEB, editor. Science in Archaeology. Londres: Thames and Huston. pp. 331–358.J. Boessneck1969Oesteological differences beetween sheep (Ovis aries Linne) and goat (Capra hircus Linne).HEB DScience in ArchaeologyLondresThames and Huston331358

29. Clutton-Brock J, Dennis-Brian K, Armitage P-A, Jewel P (1990) Osteology of the Soay sheep. Bull Br Mus Nat Hist 56: 1–56.J. Clutton-BrockK. Dennis-BrianP-A ArmitageP. Jewel1990Osteology of the Soay sheep.Bull Br Mus Nat Hist56156

30. Helmer D, Rocheteau M (1994) D. HelmerM. Rocheteau1994Atlas du squelette appendiculaire des principaux genres holocènes de petits ruminants du nord de la méditerranée et du Proche-Orient; APDCA, editor. Juan-les-Pins. Atlas du squelette appendiculaire des principaux genres holocènes de petits ruminants du nord de la méditerranée et du Proche-Orient; APDCA, editor. Juan-les-Pins.

31. Poplin F (1979) Origine du Mouflon de Corse dans une nouvelle perspective paléontologique: par marronnage. Annales de Génétique et Sélection Animale 11: 133–143.F. Poplin1979Origine du Mouflon de Corse dans une nouvelle perspective paléontologique: par marronnage.Annales de Génétique et Sélection Animale11133143

32. Loreille O, Vigne JD, Hardy C, Callou C, Treinen-Claustre F, et al. (1997) First distinction of sheep and goat archaeological bones by the means of their fossil mtDNA. Journal of Archaeological Science 24: 33–37.O. LoreilleJD VigneC. HardyC. CallouF. Treinen-Claustre1997First distinction of sheep and goat archaeological bones by the means of their fossil mtDNA.Journal of Archaeological Science243337

33. Collins M, Buckley M, Grundy HH, Thomas-Oates J, Wilson J, et al. (2010) ZooMS: the collagen barcode and fingerprints. Spectroscopy Europe 22: 11–13.M. CollinsM. BuckleyHH GrundyJ. Thomas-OatesJ. Wilson2010ZooMS: the collagen barcode and fingerprints.Spectroscopy Europe221113

34. Hughes S, Hayden TJ, Douady CJ, Tougard C, Germonpre M, et al. (2006) Molecular phylogeny of the extinct giant deer, Megaloceros giganteus. Mol Phylogenet Evol 40: 285–291.S. HughesTJ HaydenCJ DouadyC. TougardM. Germonpre2006Molecular phylogeny of the extinct giant deer, Megaloceros giganteus.Mol Phylogenet Evol40285291

35. Calvignac S, Hughes S, Tougard C, Michaux J, Thevenot M, et al. (2008) Ancient DNA evidence for the loss of a highly divergent brown bear clade during historical times. Mol Ecol 17: 1962–1970.S. CalvignacS. HughesC. TougardJ. MichauxM. Thevenot2008Ancient DNA evidence for the loss of a highly divergent brown bear clade during historical times.Mol Ecol1719621970

36. Loreille O, Orlando L, Patou-Mathis M, Philippe M, Taberlet P, et al. (2001) Ancient DNA analysis reveals divergence of the cave bear, Ursus spelaeus, and brown bear, Ursus arctos, lineages. Curr Biol 11: 200–203.O. LoreilleL. OrlandoM. Patou-MathisM. PhilippeP. Taberlet2001Ancient DNA analysis reveals divergence of the cave bear, Ursus spelaeus, and brown bear, Ursus arctos, lineages.Curr Biol11200203

37. Orlando L, Bonjean D, Bocherens H, Thenot A, Argant A, et al. (2002) Ancient DNA and the population genetics of cave bears (Ursus spelaeus) through space and time. Mol Biol Evol 19: 1920–1933.L. OrlandoD. BonjeanH. BocherensA. ThenotA. Argant2002Ancient DNA and the population genetics of cave bears (Ursus spelaeus) through space and time.Mol Biol Evol1919201933

38. Fernández H, Hughes 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 U S A 103: 15375–15379.H. FernándezS. HughesJD VigneD. HelmerG. Hodgins2006Divergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas.Proc Natl Acad Sci U S A1031537515379

39. Hofreiter M, Jaenicke V, Serre D, Haeseler Av A, Paabo S (2001) DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA. Nucleic Acids Res 29: 4793–4799.M. HofreiterV. JaenickeD. SerreA. Haeseler AvS. Paabo2001DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA.Nucleic Acids Res2947934799

40. Gouy M, Guindon S, Gascuel O (2010) SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol 27: 221–224.M. GouyS. GuindonO. Gascuel2010SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building.Mol Biol Evol27221224

41. Villesen P (2007) FaBox: an online toolbox for fasta sequences. Mol Ecol Notes 7: 965–968.P. Villesen2007FaBox: an online toolbox for fasta sequences.Mol Ecol Notes7965968

42. Posada D (2008) jModelTest: phylogenetic model averaging. Mol Biol Evol 25: 1253–1256.D. Posada2008jModelTest: phylogenetic model averaging.Mol Biol Evol2512531256

43. Huelsenbeck JP, Ronquist F (2001) MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754–755.JP HuelsenbeckF. Ronquist2001MRBAYES: Bayesian inference of phylogenetic trees.Bioinformatics17754755

44. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.F. RonquistJP Huelsenbeck2003MrBayes 3: Bayesian phylogenetic inference under mixed models.Bioinformatics1915721574

45. Bandelt HJ, Forster P, Rohl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16: 37–48.HJ BandeltP. ForsterA. Rohl1999Median-joining networks for inferring intraspecific phylogenies.Mol Biol Evol163748

46. Drummond AJ, Rambaut A, Shapiro B, Pybus OG (2005) Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 22: 1185–1192.AJ DrummondA. RambautB. ShapiroOG Pybus2005Bayesian coalescent inference of past population dynamics from molecular sequences.Mol Biol Evol2211851192

47. Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452.P. LibradoJ. Rozas2009DnaSP v5: a software for comprehensive analysis of DNA polymorphism data.Bioinformatics2514511452

48. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Ressources 10: 564–567.L. ExcoffierHEL Lischer2010Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows.Mol Ecol Ressources10564567

49. Orlando L, Calvignac S, Schnebelen C, Douady CJ, Godfrey LR, et al. (2008) DNA from extinct giant lemurs links archaeolemurids to extant indriids. BMC Evol Biol 8: 121.L. OrlandoS. CalvignacC. SchnebelenCJ DouadyLR Godfrey2008DNA from extinct giant lemurs links archaeolemurids to extant indriids.BMC Evol Biol8121

50. Calvignac S, Hughes S, Hänni C (2009) Genetic diversity of endangered brown bear (Ursus arctos) populations at the crossroads of Europe, Asia and Africa. Diversity and Distributions 15: 742–750.S. CalvignacS. HughesC. Hänni2009Genetic diversity of endangered brown bear (Ursus arctos) populations at the crossroads of Europe, Asia and Africa.Diversity and Distributions15742750

51. Ruas M-P, Vigne JD (1995) Perspectives archéobotaniques et archéozoologiques pour les périodes historiques en Corse. pp. 113–126.M-P RuasJD Vigne1995Perspectives archéobotaniques et archéozoologiques pour les périodes historiques en Corse.113126Patrimoine d'une île, 1: Recherches récentes d'archéologie médiévale en Corse: Ajaccio: Soc. Archéol. Corse-du-Sud. Patrimoine d'une île, 1: Recherches récentes d'archéologie médiévale en Corse: Ajaccio: Soc. Archéol. Corse-du-Sud.

52. Charlet P, Le Jaouen JC (1977) Les populations caprines du bassin méditerranéen. Aptitudes et évolution. Options Méditerranéennes, A 35: 45–53.P. CharletJC Le Jaouen1977Les populations caprines du bassin méditerranéen. Aptitudes et évolution.Options Méditerranéennes, A354553

53. Bruford MW, Bradley DG, Luikart G (2003) DNA markers reveal the complexity of livestock domestication. Nat Rev Genet 4: 900–910.MW BrufordDG BradleyG. Luikart2003DNA markers reveal the complexity of livestock domestication.Nat Rev Genet4900910

54. Pariset L, Cuteri A, Ligda C, Ajmone-Marsan P, Valentini A (2009) Geographical patterning of sixteen goat breeds from Italy, Albania and Greece assessed by Single Nucleotide Polymorphisms. BMC Ecol 9: 20.L. ParisetA. CuteriC. LigdaP. Ajmone-MarsanA. Valentini2009Geographical patterning of sixteen goat breeds from Italy, Albania and Greece assessed by Single Nucleotide Polymorphisms.BMC Ecol920

55. Canon J, Garcia D, Garcia-Atance MA, Obexer-Ruff G, Lenstra JA, et al. (2006) Geographical partitioning of goat diversity in Europe and the Middle East. Anim Genet 37: 327–334.J. CanonD. GarciaMA Garcia-AtanceG. Obexer-RuffJA Lenstra2006Geographical partitioning of goat diversity in Europe and the Middle East.Anim Genet37327334

56. Bouche R, Gambotti JY, Maestrini O (2004) L'avenir c'est de durer… Quand le changement en élevage extensif nécessite le maintien ou la réinvention de systèmes collectifs. Cas des caprins corses. L'évolution des systèmes de production ovine et caprine: Avenir des systèmes extensifs face au changement de la société, Alghero (ITA). Montpellier, France: CIHEAM. pp. 319–327.R. BoucheJY GambottiO. Maestrini2004L'avenir c'est de durer… Quand le changement en élevage extensif nécessite le maintien ou la réinvention de systèmes collectifs. Cas des caprins corses. L'évolution des systèmes de production ovine et caprine: Avenir des systèmes extensifs face au changement de la société, Alghero (ITA)Montpellier, FranceCIHEAM319327

57. Hugot S, Bouche R (1999) Regards sur les pratiques de l'élevage extensif de la chèvre corse: Préalable à la mise en place d'un schéma de sélection. Options Méditerranéennes, A 38: 137–143.S. HugotR. Bouche1999Regards sur les pratiques de l'élevage extensif de la chèvre corse: Préalable à la mise en place d'un schéma de sélection.Options Méditerranéennes, A38137143