The Southwestern fringe of Europe as an important reservoir of caprine biodiversity
© Martínez et al. 2015
Received: 11 May 2015
Accepted: 28 October 2015
Published: 5 November 2015
Portugal and Spain, with six and 22 officially recognized caprine breeds, encompass 25 % of the European Union goat census. Many of these populations have suffered strong demographic declines because of competition with exotic breeds and the phasing-out of low income rural activities. In this study, we have investigated the consequences of these and other demographic processes on the genetic diversity, population structure and inbreeding levels of Iberian and Atlantic goats.
A sample of 975 individuals representing 25 officially recognized breeds from Portugal and Spain, two small populations not officially recognized (Formentera and Ajuí goats) and two ecotypes of the Tinerfeña and Blanca Celtibérica breeds were genotyped with a panel of 20 microsatellite markers. A wide array of population genetics methods was applied to make inferences about the genetic relationships and demography of these caprine populations.
Genetic differentiation among Portuguese and Spanish breeds was weak but significant (FST = 0.07; P < 0.001), which is probably the consequence of their short splitting times and extensive gene flow due to transhumance. In contrast, Canarian goats were strongly differentiated because of prolonged geographic isolation. Most populations displayed considerable levels of diversity (mean He = 0.65).
High diversity levels and weak population structures are distinctive features of Portuguese and Spanish breeds. In general, these local breeds have a reduced census, but are still important reservoirs of genetic diversity. These findings reinforce the need for the implementation of management and breeding programs based on genetic data in order to minimize inbreeding, maintain overall genetic and allelic diversities and breed identities, while at the same time taking into account the within-breed genetic structure.
Goat production is a major economic activity in the Mediterranean basin, with key cultural and environmental implications, especially for smallholders and producers in marginal regions. Many local breeds that are well adapted to harsh conditions have been developed throughout the centuries in order to produce milk, meat, leather and fiber products. However, intensification of production systems since the second part of the 20th century has resulted in the demographic regression of many local populations and the concurrent expansion of a few high-producing exotic breeds . Intensive selection programs and artificial insemination have also contributed to erode the genetic reservoir represented by local breeds. Moreover, enhanced transport and communication systems have generated more uniform production environments in which a few transboundary breeds are clearly predominant. These developments have led to growing concerns about the uncontrolled loss of local goat breeds and the erosion of genetic resources that unavoidably result from census decline and increased inbreeding . In Europe, about 7 % of caprine breeds have already disappeared and many more are at the verge of extinction . In the end, the disappearance of these breeds may result in the loss of traits that are essential for adaptation to extensive farming, such as resistance to various diseases and ability to graze on poor pastures .
The importance of goat production in the Iberian Peninsula is well illustrated by the fact that it holds nearly 25 % of the caprine census of the European Union . There are six and 22 caprine breeds officially recognized in Portugal and Spain, respectively, which are raised mostly on marginal and forest lands under extensive conditions. Although a few genetic diversity studies have been carried out for some of these breeds [4–7], no comprehensive analysis has been published regarding their overall diversity and genetic structure. In the current study, we aimed at investigating the amount of diversity, population structure, level of inbreeding and genetic relationships across a broad array of local breeds from the Iberian Peninsula and Atlantic and Balearic archipelagos.
Genetic diversity parameters estimated with 20 microsatellite loci in 29 Portuguese and Spanish goat populations
He ± SD
Ho ± SD
6.60 ± 2.50
0.695 ± 0.043
0.654 ± 0.027
7.10 ± 3.09
0.687 ± 0.048
0.626 ± 0.020
6.70 ± 2.81
0.659 ± 0.039
0.634 ± 0.017
Blanca de Rasqueraa
6.25 ± 2.81
0.634 ± 0.050
0.586 ± 0.017
4.50 ± 2.21
0.611 ± 0.062
0.551 ± 0.038
5.61 ± 2.40
0.688 ± 0.042
0.677 ± 0.029
6.50 ± 2.54
0.654 ± 0.045
0.537 ± 0.022
6.65 ± 2.62
0.665 ± 0.041
0.631 ± 0.017
7.15 ± 2.66
0.663 ± 0.042
0.621 ± 0.017
6.55 ± 2.26
0.652 ± 0.044
0.574 ± 0.021
6.80 ± 2.88
0.683 ± 0.041
0.627 ± 0.017
6.6 0 ± 2.37
0.655 ± 0.049
0.615 ± 0.017
7.25 ± 3.01
0.697 ± 0.036
0.666 ± 0.017
6.40 ± 3.14
0.672 ± 0.041
0.674 ± 0.017
6.35 ± 2.62
0.657 ± 0.037
0.605 ± 0.017
4.20 ± 1.61
0.598 ± 0.051
0.558 ± 0.034
6.35 ± 2.39
0.652 ± 0.045
0.573 ± 0.017
6.80 ± 2.61
0.649 ± 0.040
0.600 ± 0.017
5.85 ± 2.28
0.651 ± 0.029
0.618 ± 0.017
6.60 ± 2.89
0.637 ± 0.038
0.612 ± 0.017
4.15 ± 1.53
0.497 ± 0.038
0.507 ± 0.018
5.30 ± 2.47
0.603 ± 0.037
0.588 ± 0.018
6.00 ± 2.58
0.595 ± 0.037
0.590 ± 0.017
5.95 ± 2.54
0.632 ± 0.046
0.621 ± 0.017
6.80 ± 3.21
0.671 ± 0.045
0.618 ± 0.020
6.35 ± 2.62
0.681 ± 0.036
0.647 ± 0.020
6.40 ± 2.33
0.685 ± 0.034
0.655 ± 0.020
6.90 ± 2.83
0.674 ± 0.042
0.601 ± 0.020
Preta de Montesinho
6.75 ± 2.79
0.669 ± 0.043
0.563 ± 0.019
Genomic DNA was obtained from hair or blood samples using the Chelex 100 chelating resin (Bio Rad Laboratories, Hercules, CA, USA) according to the methodology described by Walsh et al. . Twenty microsatellite loci were chosen according to the recommendations of the FAO/ISAG . More specifically, the following loci were analyzed: BM1329, BM6506, BM6526, BM8125, CRSM60, CSRD247, ETH010, ETH225, HAUT27, ILSTS011, INRA063, MAF065, MAF209, McM527, MM12, OarFCB048, OarFCB304, SPS115, SRCRSP08, and TGLA122 (see Additional file 1: Table S1). Each multiplex PCR was carried out in 25 µL reaction tubes containing 30–60 ng of genomic DNA, 5 µL of 2× PCR mastermix and 2 pmol of each primer (forward primer labeled at the 5′ end with HEX, FAM or NED fluorescent dyes). The PCR mix was subjected to an initial denaturation step at 95 °C for 5 min, followed by 35 cycles of 45 s at 95 °C, 35 s at 55 °C and 35 s at 72 °C with a final extension step of 10 min at 72 °C. PCR products were separated by electrophoresis using an ABI377 equipment (Applied Biosystems, Madrid, Spain) according to the manufacturer’s recommendations. Allele sizes were determined with the internal size standard GeneScan-400HD ROX (Applied Biosystems, Madrid, Spain). Reference samples were used in each assay to ensure the consistency of allele assignments.
Population genetics analyses
Allele frequencies for each locus, total number of alleles per locus (NA), observed (Ho) and unbiased expected (He) heterozygosities, and mean number of alleles (MNA) per population were calculated using the MICROSATELLITE TOOLKIT software . Weir and Cockerham  F-statistics and allelic richness across all loci per breed were estimated with FSTAT . A graphic representation of the matrix containing pairwise F ST distances derived from the 29 goat populations studied was generated with R-lequin . The mean F ST and their 95 % confidence interval across loci after 1000 bootstraps were calculated using the GENETIX 4.05 software . The same program was used to calculate F ST P values after 1000 permutations of individuals across the entire population. A Mantel test  was performed using the statistical software R version 2.15.0  to test for isolation-by-distance, by investigating the correlation that exists between genetic and geographical distances between pairs of breeds. Fisher’s exact tests for ascertaining Hardy–Weinberg (HW) equilibrium across loci and populations and estimation of the frequencies of null alleles were performed with GENEPOP 1.2 . The relationship between the percentage of missing data for each locus and F IS was examined  with the Spearman’s rank correlation using the SAS statistical package v. 9.1.3 . A positive relationship between F IS and missing data for a given locus indicates that amplification failure is due to individuals that carry a null allele in the homozygous state.
The model-based clustering program STRUCTURE  was used to investigate population structure and to estimate the proportions of individual genotypes derived from the inferred clusters. We considered an admixture model and uncorrelated allele frequencies. The putative number of clusters (K) ranged from 2 to 35, and five independent runs were performed with a number of Monte-Carlo Markov chain iterations that ranged from 250,000 (50,000 burn-in length) for K between 2 and 10, to 500,000 (150,000 burn-in length) for K between 11 and 35. To identify the most probable number of ancestral populations, we investigated the mean and variance of the likelihood plots of the data for different K-values (i.e., plot of ln Pr(X/K) vs. K). Bar charts representing the proportions of the genotype membership coefficient of each individual (q) in a given ancestral population, as obtained with STRUCTURE, were visualized using DISTRUCT . The degree of admixture or ancestry diversity of each breed was calculated as 1 − Σ(qk)2, where q k is the average fraction of the genetic ancestry of a given breed that belongs to the kth ancestral population, estimated by STRUCTURE analyses . The correlation between the proportion of mixed ancestry and the expected heterozygosity was calculated to evaluate the relationship between admixture and within-breed diversity. Average genotype membership coefficients in each cluster (Q) were converted into genetic distances among breeds, following the methodology described by Cañón et al. . In this case, a FORTRAN program that implements the computation of Weitzman’s diversity was used, and with the ultra-metric distance matrix obtained, a hierarchical tree resulting from Weitzman’s algorithm  was constructed with MEGA 5 . Information on genetic markers was combined with spatial data to draw synthetic contour maps of the Iberian Peninsula, representing the geographical patterns of genetic variability [27, 28]. These maps were drawn based on the interpolation of genetic contributions to each breed, as computed in the analysis with STRUCTURE for K = 3. The kriging interpolation method was used  and the graphical library of statistical software R  was used to display the maps. Each breed was represented on the map by the coordinates of its center of geographical dispersion.
Genetic structure was further investigated by a factorial analysis of correspondence (FAC) using the function “AFC 3D by populations” implemented in GENETIX 4.05 , which is analogous to a principal component analysis. In this approach, allele frequencies are used to infer the relative position of each breed based on Chi square distances.
The 20 microsatellite markers analyzed were highly variable, with a total of 236 alleles across all populations. The number of alleles per locus ranged from 3 (MAF209) to 22 (OarFCB304), with an average of 11.8 (see Additional file 2: Table S2). Expected (He) and observed (Ho) heterozygosities per locus across all breeds ranged from 0.199 (ETH225) to 0.86 (MM12), and from 0.20 (ETH225) to 0.80 (MM12), respectively. The overall means across loci and breeds for He and Ho were equal to 0.65 and 0.61, respectively (see Additional file 2: Table S2). Among the breeds studied, the smallest MNA per locus (<5.0) was found for the Palmera, Formentera and Guadarrama goat breeds, and the largest MNA (>7.0) for Celtibérica, Florida and Moncaína (Table 1). Palmera had the lowest mean allelic richness (3.2), while the highest estimate was found for Pirenaica (5.5) and Moncaína (5.3). The amount of data was not sufficient to calculate allelic richness for the Guadarrama and Retinta breeds because some individuals were not successfully amplified and genotyped at all microsatellite markers and thus, there were many blanks in the dataset. The He and Ho by breed, averaged across loci, and the measure of within-breed discrepancy among loci (F IS) are in Table 1. He per breed ranged from 0.50 (Palmera) to 0.70 (Florida and Pirenaica). Of the 29 populations studied, 12 exhibited a highly significant deficiency in heterozygotes (P < 0.001) and, in this regard, the Verata (0.18) and Preta de Montesinho (0.16) breeds showed the highest within-breed F IS coefficients (Table 1). Only two breeds (Verata and Negra Serrana) had more than two loci that deviated from HW equilibrium (Table 1), and only two microsatellite markers (OarFCB304 and SPS115) did not adjust to the HW equilibrium model for more than one breed (P < 0.001), as shown in Additional file 2: Table S2.
Frequencies of null alleles ranged from 1.2 % for Tenerife Sur to 7.7 % for Verata breeds (Table 1). In a locus-by-locus analysis, the markers with the highest frequencies of null alleles were SPS115 (9.7 %) and INRA063 (6.6 %) (see Additional file 2: Table S2). Moreover, the relationship between the percentage of missing data for each locus and F IS was not significant (P > 0.01).
Measurement of Wright’s F ST coefficient
The Mantel test showed a strong correlation (R2 = 0.72) between F ST genetic distances and geographic distances among breed pairs (P < 0.0001 after 10,000 permutations), although this value dropped to 0.029 when Canarian breeds were removed from the analysis (see Additional file 3: Figure S1). Thus, for breeds located in the Iberian Peninsula there was no clear association between genetic and geographical distances.
Structure analysis and tree building
Microsatellite data for the goat populations studied were analyzed with the Bayesian model-based clustering implemented in STRUCTURE  and the plot of average likelihoods of the data for different values of K (ln Pr(X/K) vs. K) was used to infer the most likely number of genetic clusters, which was determined to be 12 (see Additional file 4: Figure S2). The proportional contributions of the inferred ancestral populations per breed are in Additional file 5: Table S3 for K = 12.
To gain additional insights into population structure, the fractional contributions of ancestral populations to each breed were converted into genetic distances among breeds and represented in UPGMA trees for K = 3 (see Additional file 7: Figure S3) and K = 12 (Fig. 3), as described by Cañón et al.  and García et al. . These analyses showed three main groups for K = 3: (1) Canarian breeds; (2) Blanca de Rasquera, Payoya, Pirenaica, Formentera, Moncaína, and Blanca Celtibérica; and (3) a main group including the remaining breeds (see Additional file 7: Figure S3).
The tree resulting from the ultra-metric distance matrix obtained with the Weitzman’s algorithm for K = 12 (Fig. 3) indicates that only four groups can be clearly identified, i.e.: (1) Canary Islands breeds; (2) a large group including the Spanish breeds Verata, Florida, Malagueña, Retinta, Blanca Andaluza, Celtibérica, and Guadarrama, and the Portuguese breeds Serrana, Preta de Montesinho, Serpentina, and Charnequeira; (3) the Pirenaica/Moncaína group; and (4) the Payoya/Formentera group. The populations that are in each group (iii) and (iv) shared the same genetic background, with average membership proportions of about 0.5 and 0.7, respectively (see Additional file 5: Table S3). It should be noted that the remaining breeds (Blanca de Rasquera, Murciano-Granadina, Negra Serrana, Pitiusa, Mallorquina and Algarvia) had an estimated membership fraction greater than 0.4 in their main respective inferred clusters (see Additional file 5: Table S3). The clustering of the Azpi Gorri breed was cryptic and shared its main genetic background inconsistently with either Blanca Celtibérica or Bravia goats depending on which STRUCTURE run was considered.
A highly significant positive correlation (r = 0.71, P < 0.0001) was estimated between breed He and levels of admixture based on global STRUCTURE results for K = 12 (see Additional file 5: Table S3), which suggested that the levels of admixture of breeds could be responsible for a large proportion of their estimated genetic variability.
Factorial analysis of correspondence
A weak population structure in goat populations from the Iberian Peninsula and Balearic archipelago
We studied 29 autochthonous goat populations that were sampled in Portugal and Spain by using 20 microsatellite markers in a survey that included endangered populations (e.g., Pitiusa and Blanca Andaluza), commercial goat breeds with larger census (e.g., Murciano-Granadina and Malagueña) and a feral population from the Canary Islands (Ajuí). Estimates of genetic diversity and analyses of population structure with microsatellites were previously reported for certain specific Iberian goat populations [6, 7, 30] but, to our knowledge, this study is the most comprehensive analysis of the genetic diversity of native goats from these two countries. This broad representativeness offers an unprecedented perspective about the variation and population structure of Portuguese and Spanish goats. One of the major features revealed by our study is that the population structure of Iberian goats is weak, and that populations from Portugal and Spain (with the exception of Canarian goats) are, in general, poorly differentiated. This observation can be explained by two factors. First, the splitting times of Portuguese and Spanish breeds are probably very short (in the range of a few hundred years) precluding a strong genetic differentiation. This is well-known for sheep breeds such as the Merino, which was maintained as a joint population exclusively in the Iberian Peninsula until the 18th century . Certain goat populations may have been subject to similar handling procedures. Second, there are no strong geographical barriers between Portugal and Spain that would hinder a bidirectional gene flow. Indeed, transhumance trails crossed the border that separates both countries , particularly in the southern part of the Iberian Peninsula, and small ruminants would be periodically taken from one side to the other, providing the opportunity for genetic exchange to occur. Within Portugal and Spain, seasonal pastoralism and transhumance, which were typical of most goat production systems in the Iberian Peninsula until the mid-20th century, contributed to weaken the population structure . It has also been claimed that long-distance cyclic migrations and the great mobility of goats are the main causal factors that explain the poor phylogeographic structure detected with mitochondrial markers in the Iberian Peninsula  and at a worldwide scale . In our study, it is interesting to highlight that the genetic affinity between the breeds in Cluster II of Fig. 4 (which reveals the existence of two clouds of affinity in the center and west of the Iberian Peninsula, respectively) follows a clear north to south dispersion path. Since it is known that nearly all the historical transhumance trails also followed the same direction , it is expected that the observed geographic relationship among breeds probably reflects this gene flow and admixture with migrants along the routes of transhumance.
In several instances, the between-breed genetic affinities reported in our study were quite unexpected. For example, the Payoya breed (from the southern Mediterranean coast of the Peninsula) clustered with the Formentera population (from the Balearic Islands) in the absence of any historical record that would justify such a relationship. However, an in-depth analysis of this cluster allowed for a clear separation of the two breeds, thus suggesting that migration to the Balearic Islands may have occurred in the past. A close relationship was expected to exist between the Moncaína and Guadarrama breeds, given their very close geographic distribution and common origin, and the fact that these two populations are officially considered as varieties of the Pirenaica breed. However, the analysis with STRUCTURE revealed a close relationship between Pirenaica and Moncaína goats, while the Guadarrama breed was clearly separated from both populations.
We were unable to detect a clear differentiation between several of the breeds analyzed (Figs. 3, 5), an outcome that may have demographic and technical causes. The availability of high-density panels of single nucleotide polymorphisms (SNPs), which until now have been very little used for population genetic studies in goats, will allow a more refined analysis of goat population structure, and contribute to a better understanding of the recent evolutionary history of domestic goats [35–37].
Canarian goats are strongly differentiated from their Iberian and Balearic counterparts
The Canarian breeds showed a strong genetic differentiation from their Iberian and Balearic counterparts, probably because they descend from the North African goat populations that were introduced by the first Imazighen settlers of the Canary Islands [38, 39]. After the conquest by the Spaniards in the 15th century, the Canary archipelago became an important maritime trading platform between Europe and the Americas , providing the opportunity for Canarian goats to hybridize with foreign breeds from Spain and other European and African countries. Notably, the intensity of this hybridization process may have varied from one island to another .
Existence of genetic substructure at the within-breed level
The existence of within-breed genetic heterogeneity should be understood as an essential part of their history, instead of considering it as a negative feature . Several examples of within-breed genetic heterogeneity were identified in our study. The Murciano-Granadina breed is one of the most important dairy goat populations in Spain and results from the administrative unification, in the 1970’s, of two well differentiated genetic groups . This event may explain the substructure revealed by the significant F IS estimate obtained for this breed (Table 1). A somewhat similar situation was evidenced in the Serrana breed, which has three distinct ecotypes that are raised in different geographical regions i.e. Transmontana, Serra and Ribatejana. This geographical isolation may have facilitated the emergence of a population substructure resulting in the significant deficiency in heterozygotes observed in our study. Another case is represented by the Formentera goats, which are classified as part of the Pitiüsa breed in the Spanish Official Breed Catalogue . However, our findings indicate that Formentera and Pitiüsa goats are clearly differentiated (Figs. 3, 5), probably because of genetic drift and inbreeding. In other breeds such as Moncaína, Verata, Serpentina and Preta de Montesinho (Table 1), factors such as inbreeding, population substructure and genetic drift resulting from small population sizes could also have contributed to the observed deficiency in heterozygotes, although the possibility of population substructure cannot be excluded. Finally, the Celtibérica and Blanca Celtibérica populations are considered as two geographic varieties of the same breed, but Celtibérica is reared mainly in the Southern Central part of the Iberian Peninsula while Blanca Celtibérica is raised in Eastern Spain (Castellón province), along the Mediterranean coast. Some decades ago, Blanca Celtibérica was considered extinct  but a few breeders still maintain this breed in the Castellón region. All our analyses indicated a clear separation between these two populations, and Blanca Celtibérica goats were more differentiated from the remaining Iberian breeds than Celtibérica (Figs. 3, 5), which showed some level of genetic affinity with Andalusian goats (Blanca Andaluza, Malagueña and Florida).
Portuguese and Spanish goats display high levels of diversity
In the last decades, many caprine local breeds have suffered a strong demographic decline because of the phasing-out of low income farming activities, replacement of local breeds by cosmopolitan high-producing breeds, intensification of agricultural practices, and the broad use of artificial insemination (AI) as well as other factors . AI is only used for a few Spanish breeds such as the Murciano-Granadina, Malagueña, Florida breeds and the endangered Payoya and Pitiusa breeds and it is always applied within the framework of the selection or conservation programs of these breeds (see Additional file 8: Table S5). In Portugal, AI is used on a small scale only in the Serrana breed, but all goat breeds have germplasm conservation programs, including semen cryopreservation (see Additional file 8: Table S5).
Overall, we found that genetic diversity was high within most goat breeds studied here, with an overall mean of ~12 alleles/locus and an average He of 0.65 for all microsatellite loci analyzed. At the breed level, the mean number of alleles per locus was about 6, with an allelic richness corrected for sample size of about 5 and a He of about 0.65. These values are greater than those found for Asian goats , and similar to those reported for European, Chinese, Indian, and Brazilian breeds [24, 43–46]. However, these comparisons are not straightforward, because they depend on the panel of genetic markers used, the breeds analyzed and sample size (which has a strong effect on the number of alleles detected).
Among the Portuguese and Spanish goat breeds, 12 out of 29 showed significantly positive F IS values, which indicates a deficiency in heterozygotes that could be due to inbreeding, the Wahlund effect (population substructure), and other causes . The Verata, Preta de Montesinho, Pitiüsa and Serpentina breeds exhibited the highest F IS values, but at present it is difficult to infer if this observation is due to population structure, inbreeding or a combination of both factors. As previously said, in the Pitiüsa population, the census is rather small, so inbreeding may have contributed significantly to the high F IS that we detected (Table 1). In contrast, in the Verata breed (F IS = 0.183), there are over 8000 breeding females and the rather high F IS observed may reflect the existence of population substructure.
Overall, we found that the levels of diversity were high, although this finding should not give rise to an excessive optimism because many of the breeds that we analyzed are undergoing a steady and sustained demographic decline that may lead, in the worst case scenario, to their disappearance. One of the most compelling cases is the Blanca de Rasquera population, which in the first half of the 20th century had a census of 30,000 individuals , whereas today it has decreased to 5000 goats. In Portugal and Spain, as in many other Western countries, the proportion of the economic active population dedicated to farming activities is decreasing at a rapid pace, and sheep and goats have been progressively displaced to marginal areas, while local genetic resources are being replaced by more productive industrial and transboundary breeds. Moreover, intensive selection schemes, population fragmentation into discrete subpopulations and artificial insemination contribute to reduce the variability of local as well as transboundary breeds. The joint effects of all these threats could significantly shrink the current gene pool of goat breeds and cause a loss of biodiversity that may have deleterious effects on adaptive traits such as resistance to various diseases, adaptation to harsh environments and the ability to cope with climate changes.
Our results indicate that the levels of genetic diversity are high in the native goat breeds from Portugal and Spain. Indeed, local breeds, often with a small census, are important reservoirs of genetic diversity. With regard to population structure, Portuguese and Spanish breeds were weakly differentiated. In contrast, Canarian breeds showed a strong genetic differentiation from their Iberian and Balearic counterparts, probably due to influences from their North African counterparts combined with continued geographic isolation. Among the Portuguese and Spanish populations, 12 out of 29 showed significantly positive F IS values, which could be due to a variety of reasons, such as inbreeding and the Wahlund effect. Although, overall, the levels of diversity that we observed are high, these results should be taken with caution, because many of the breeds analyzed are undergoing a steady and sustained demographic decline that could lead to their disappearance. Management and breeding programs based on genetic data should be undertaken in order to minimize inbreeding, maintain overall genetic and allelic diversities and breed identities, and also account for within-breed genetic structure.
AMM, LTG, JVD and NS conceived and designed the experiments; AMM, VL and CBS performed the experiments; NS, JC, MA and AM analyzed the data; AMM, LTG and NS wrote the manuscript; JVD, JC, MA, CBS, CG, AM, PZ and VL reviewed and edited the manuscript. Members of the BIOGOAT Consortium provided biological samples, populations’ information and logistic support. All authors read and approved the final manuscript.
The authors wish to express thanks to breeders associations and research groups who kindly provided biological samples used in this study. Members of the CYTED XII-H and CONBIAND networks are thanked for valuable cooperation over the years. This study has been partially supported by the Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA) RZ2004-00037-00-00, RZ2006- 00005-C02-00, RZ2007-00005-C02-01 and RZ2007-00005-C02-02 projects. CG was supported by a grant from the Fundação para a Ciência e a Tecnologia (Ref. SFRH/BCC/52687/2014). CBS was supported by a grant from the Fundação para a Ciência e a Tecnologia (Ref. SFRH/BM/1770/2004).
List of the members of the BioGoat Consortium
IFAPA Alameda del Obispo, Córdoba, Spain
Instituto Canario de Investigaciones Agrarias, La Laguna, Tenerife, Spain
Departamento de Producción Animal, Universidad Complutense de Madrid, Spain
Departamento de Producción Animal, Universidad Complutense de Madrid, Spain
Departament de Ciència Animal i dels Aliments. Universitat Autònoma de Barcelona, Bellaterra, Spain
Servicio de Ganadería, Diputación Foral de Bizkaia, Spain
Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, Bellaterra, Spain
Inmaculada Martin-Burriel, Laboratorio de Genética Bioquímica, Facultad de Veterinaria, Universidad de Zaragoza, Spain
Asociación Nacional de Criadores de Caprino de Raza Murciano-Granadina (CAPRIGRAN), Granada, Spain
Laboratorio de Genética Bioquímica, Facultad de Veterinaria, Universidad de Zaragoza, Spain
Daniel Martín Santana
Departamento de Genética, Universidad de Córdoba, Córdoba, Spain
Águeda Pons Barro
Servei de Millora Agraria i Pesquera (SEMILLA EPE), Mallorca, Islas Baleares, Spain
INIAV, 2005-048 Vale de Santarém, Portugal
Departament de Biologia, Universitat de Girona, Girona, Spain
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- FAO Commission on genetic resources for food and agriculture. The state of the world’s animal genetic resources for food and agriculture—in brief. Pilling D, Rischkowsky B, editors. Rome; 2007. ftp://fao.org/docrep/fao/010/a1260e/a1260e00.pdf. Accessed 16 Oct 2014.
- Taberlet P, Valentini A, Rezaei HR, Naderi S, Pompanon F, Negrini R, et al. Are cattle, sheep, and goats endangered species? Mol Ecol. 2008;17:275–84.View ArticlePubMedGoogle Scholar
- European Union—Directorate-general for agriculture and rural development. Agriculture in the European Union: statistical and economic information. 2013. http://ec.europa.eu/agriculture/statistics/agricultural/2013/pdf/full-report_en.pdf. Accessed 25 June 2014.
- Amills M, Capote J, Tomas A, Kelly L, Obexer-Ruff G, Angiolillo A, et al. Strong phylogeographic relationships among three goat breeds from the Canary Islands. J Dairy Res. 2004;71:257–62.View ArticlePubMedGoogle Scholar
- Azor PJ, Monteagudo LV, Luque M, Tejedor MT, Rodero E, Sierra I, et al. Phylogenetic relationships among Spanish goats breeds: phylogenetic relationships among goats breeds. Anim Genet. 2005;36:423–5.View ArticlePubMedGoogle Scholar
- Martínez AM, Acosta J, Vega-Pla JL, Delgado JV. Analysis of the genetic structure of the canary goat populations using microsatellites. Livest Sci. 2006;102:140–5.View ArticleGoogle Scholar
- Bruno-de-Sousa C, Martinez AM, Ginja C, Santos-Silva F, Carolino MI, Delgado JV, et al. Genetic diversity and population structure in Portuguese goat breeds. Livest Sci. 2011;135:131–9.View ArticleGoogle Scholar
- Martinez AM, Vega-Pla JL, Leon JM, Camacho ME, Delgado JV, Ribeiro MN. Is the Murciano-Granadina a single goat breed? A molecular genetics approach. Arq Bras Med Vet Zootec. 2010;62:1191–8.View ArticleGoogle Scholar
- Walsh PS, Metzger DA, Higushi R. Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques. 1991;10:506–13.PubMedGoogle Scholar
- FAO Commission on genetic resources for food and agriculture. Molecular genetic characterization of animal genetic resources. FAO Animal Production and Health Guidelines. Secondary guidelines: measurement of domestic animal diversity (MoDAD), Rome. 1993. http://www.fao.org/docrep/014/i2413e/i2413e00.pdf. Accessed 14 Sep 2015.
- Park SDE. Trypanotolerance in West African cattle and the population genetic effects of selection. PhD thesis, University of Dublin. 2002.Google Scholar
- Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.View ArticleGoogle Scholar
- Goudet J. FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Lausanne University. 2001. http://www2.unil.ch/popgen/softwares/fstat.htm. Accessed 11 Oct 2012.
- Lisher H. Population genetics data acquisition and visualization—R-lequin. University of Bern. 2008. http://heidi.chnebu.ch/doku.php?id=r-lequin. Accessed 15 April 2014.
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Université de Montpellier II. 2004. http://kimura.univ-montp2.fr/genetix/. Accessed 04 April 2012.
- Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27:209–20.PubMedGoogle Scholar
- R Development Core Team. R. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2012.Google Scholar
- Raymond M, Rousset F. GENEPOP (Version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.Google Scholar
- Beaumont M, Barratt EM, Gottelli D, Kitchener AC, Daniels MJ, Pritchard JK, et al. Genetic diversity and introgression in the Scottish wildcat. Mol Ecol. 2001;10:319–36.View ArticlePubMedGoogle Scholar
- SAS. Statistical analysis with SAS/STAT® Software V9.1. SAS Institute Inc. 2009.Google Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMed CentralPubMedGoogle Scholar
- Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4:137–8.View ArticleGoogle Scholar
- Tapio M, Ozerov M, Tapio I, Toro MA, Marzanov N, Cinkulov M, et al. Microsatellite-based genetic diversity and population structure of domestic sheep in northern Eurasia. BMC Genet. 2010;11:76.PubMed CentralView ArticlePubMedGoogle Scholar
- Cañon J, Garcia D, Garcia-Atance MA, Obexer-Ruff G, Lenstra JA, Ajmone-Marsan P, et al. Geographical partitioning of goat diversity in Europe and the Middle East. Anim Genet. 2006;37:327–34.View ArticlePubMedGoogle Scholar
- García D, Corral N, Cañón J. Combining inter- and intrapopulation information with the Weitzman approach to diversity conservation. J Hered. 2005;96:704–12.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9.View ArticlePubMedGoogle Scholar
- Manel S, Schwartz MK, Luikart G, Taberlet P. Landscape genetics: combining landscape ecology and population genetics. Trends Ecol Evol. 2003;18:189–97.View ArticleGoogle Scholar
- Storfer A, Murphy MA, Evans JS, Goldberg CS, Robinson S, Spear SF, et al. Putting the “landscape” in landscape genetics. Heredity. 2007;98:128–42.View ArticlePubMedGoogle Scholar
- Hengl T. A practical guide to geostatistical mapping of environmental variables. Luxembourg: Office for official publications of the European Communities; 2007.Google Scholar
- Martínez A, Ferrando A, Manunza A, Gómez M, Landi V, Jordana J, et al. Inferring the demographic history of a highly endangered goat breed through the analysis of nuclear and mitochondrial genetic signatures. Small Ruminant Res. 2012;104:78–84.View ArticleGoogle Scholar
- Diez-Tascón C, Littlejohn RP, Almeida PA, Crawford AM. Genetic variation within the Merino sheep breed: analysis of closely related populations using microsatellites. Anim Genet. 2000;31:243–51.View ArticlePubMedGoogle Scholar
- Manzano P, Casas R. Past, present and future of trashumancia in Spain: nomadism in a developed country. Pastoralism. 2010;1:72–90.Google Scholar
- Rigueiro-Rodríguez A, McAdam J, Mosquera-Losada MR. Agroforestry in Europe—current status and future prospects. Santiago de Compostela: Springer; 2009.Google Scholar
- Luikart G, Gielly L, Excoffier L, Vigne JD, Bouvet J, Taberlet P. Multiple maternal origins and weak phylogeographic structure in domestic goats. Proc Natl Acad Sci USA. 2001;98:5927–32.PubMed CentralView ArticlePubMedGoogle Scholar
- Pariset L, Cuteri A, Ligda C, Ajmone-Marsan P, Valentini A, ECONOGENE Consortium. Geographical patterning of sixteen goat breeds from Italy, Albania and Greece assessed by single nucleotide polymorphisms. BMC Ecol. 2009;9:20.PubMed CentralView ArticlePubMedGoogle Scholar
- Brito LF, Jafarikia M, Grossi DA, Kijas JW, Porto-Neto LR, Ventura RV, Salgorzaei M, Schenkel F, et al. Characterization of linkage disequilibrium, consistency of gametic phase and admixture in Australian and Canadian goats. BMC Genet. 2015;16:67.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicoloso L, Bomba L, Colli L, Negrini R, Milanesi M, Mazza R, et al. Genetic diversity of Italian goat breeds assessed with a medium-density SNP chip. Genet Sel Evol. 2015;47:62.PubMed CentralView ArticlePubMedGoogle Scholar
- Catálogo oficial de razas de Ganado de España. Accessed 2 Dec 2014.Google Scholar
- Rodero A, Delgado JV, Rodero E. Primitive Andalusian livestock and their implications in the discovery of America. Arch Zootec. 1992;41:383–400.Google Scholar
- Capote J, Delgado JV, Fresno M, Camacho ME, Molina A. Morphological variability in the Canary goat population. Small Ruminant Res. 1998;27:167–72.View ArticleGoogle Scholar
- European Cattle Genetic Diversity Consortium. Marker-assisted conservation of European cattle breeds: an evaluation. Anim Genet. 2006;37:475–81.View ArticleGoogle Scholar
- Barker JSF, Tan SG, Moore SS, Mukherjee TK, Matheson JL, Selvaraj OS. Genetic variation within and relationships among populations of Asian goats (Capra hircus). J Anim Breed Genet. 2001;118:213–33.View ArticleGoogle Scholar
- Fatima S, Bhong CD, Rank DN, Joshi CG. Genetic variability and bottleneck studies in Zalawadi, Gohilwadi and Surti goat breeds of Gujarat (India) using microsatellites. Small Ruminant Res. 2008;77:58–64.View ArticleGoogle Scholar
- Glowatzki-Mullis ML, Muntwyler J, Bäumle E, Gaillard C. Genetic diversity measures of Swiss goat breeds as decision-making support for conservation policy. Small Ruminant Res. 2008;74:202–11.View ArticleGoogle Scholar
- Ribeiro MN, Bruno-de-Sousa C, Martinez-Martinez A, Ginja C, Menezes MPC, Pimenta-Filho EC, et al. Drift across the Atlantic: genetic differentiation and population structure in Brazilian and Portuguese native goat breeds: Brazilian and Portuguese goat breeds. J Anim Breed Genet. 2012;129:79–87.View ArticlePubMedGoogle Scholar
- Wei C, Lu J, Xu L, Liu G, Wang Z, Zhao F, et al. Genetic structure of Chinese indigenous goats and the special geographical structure in the Southwest China as a geographic barrier driving the fragmentation of a large population. PLoS One. 2014;9:e94435.PubMed CentralView ArticlePubMedGoogle Scholar
- Hedrick PW. High inbreeding in sheep or erroneous estimation? J Hered. 2013;104:298–9.View ArticlePubMedGoogle Scholar
- Carné I Fructuoso S, Roig N, Jordana I Vidal J. The Blanca de Rasquera goat breed: morphological and phaneroptical characterization. Arch Zootec. 2007;56:319–30.Google Scholar