- Research Article
- Open Access
Genetic diversity of Italian goat breeds assessed with a medium-density SNP chip
- Letizia Nicoloso†1Email author,
- Lorenzo Bomba†2,
- Licia Colli†2,
- Riccardo Negrini2, 3,
- Marco Milanesi2,
- Raffaele Mazza4,
- Tiziana Sechi5,
- Stefano Frattini1,
- Andrea Talenti1,
- Beatrice Coizet1,
- Stefania Chessa6,
- Donata Marletta7,
- Mariasilvia D’Andrea8,
- Salvatore Bordonaro7,
- Grazyna Ptak9,
- Antonello Carta5,
- Giulio Pagnacco1,
- Alessio Valentini10,
- Fabio Pilla8,
- Paolo Ajmone-Marsan2,
- Paola Crepaldi1 and
- the Italian Goat Consortium
© Nicoloso et al. 2015
- Received: 23 January 2015
- Accepted: 9 July 2015
- Published: 4 August 2015
Among the European countries, Italy counts the largest number of local goat breeds. Thanks to the recent availability of a medium-density SNP (single nucleotide polymorphism) chip for goat, the genetic diversity of Italian goat populations was characterized by genotyping samples from 14 Italian goat breeds that originate from different geographical areas with more than 50 000 SNPs evenly distributed on the genome.
Analysis of the genotyping data revealed high levels of genetic polymorphism and an underlying North–south geographic pattern of genetic diversity that was highlighted by both the first dimension of the multi-dimensional scaling plot and the Neighbour network reconstruction. We observed a moderate and weak population structure in Northern and Central-Southern breeds, respectively, with pairwise FST values between breeds ranging from 0.013 to 0.164 and 7.49 % of the total variance assigned to the between-breed level. Only 2.11 % of the variance explained the clustering of breeds into geographical groups (Northern, Central and Southern Italy and Islands).
Our results indicate that the present-day genetic diversity of Italian goat populations was shaped by the combined effects of drift, presence or lack of gene flow and, to some extent, by the consequences of traditional management systems and recent demographic history. Our findings may constitute the starting point for the development of marker-assisted approaches, to better address future breeding and management policies in a species that is particularly relevant for the medium- and long-term sustainability of marginal regions.
- Goat Breed
- Admixture Analysis
- Discovery Panel
- Weak Population Structure
- Italian Breed
According to archeozoological and genetic data, goats were domesticated some 10 000 years ago in the geographical region that spans from Eastern Anatolia to the Zagros Mountains in Northern Iran [1, 2]. After domestication, goats quickly spread all over the world following human migrations and commercial trade . They rapidly adapted to a very wide range of environmental conditions and started to play economic, cultural and religious roles in many human cultures. Today, goats represent an important source of milk, meat and fiber (e.g., cashmere wool) especially in marginal rural areas, dry lands and mountains, particularly in developing countries. As a consequence of an increase in farmland abandonment in marginal areas, the genetic diversity of many goat populations is being rapidly eroded or lost, particularly in Europe that counts 200 recorded goat breeds (according to the Food and Agriculture Organization (FAO) .
In modern European agriculture, the economic role of goats is mainly linked to the products of either high-yielding dairy breeds that developed in the central Alps (e.g., the Saanen and Toggenburg breeds of Swiss origin) or of local stocks that were often improved by crossing with more productive dairy or meat breeds (e.g., cosmopolitan Boer). In Europe, 96 % of the 2.8 millions of tons of goat products are dairy products and only 4 % are meat products (see FAOSTAT at faostat.fao.org).
Even if the main hotspots of the world goat diversity are probably in Africa and Asia, among European countries, Italy can be considered as a reservoir of genetic resources for the caprine species with 36 breeds recorded by the National Goat and Sheep Breeder Association (www.assonapa.it). Diversity in orography and climate together with historical factors and traditions, led to the development of a large variety of livestock populations, which were later standardized in modern breeds. Although, since the early 1950s the genetic diversity of farm species has suffered from steady erosion, an opposite trend is observed for goats. A number of local goat breeds still populate rural environments where harsh climate and pastures challenge the diffusion of more productive species and cosmopolite breeds. According to the FAO classification, 58.3 % of Italian goat breeds are considered not to be at risk, each with more than 1200 heads, but there is a declining trend. These breeds include three cosmopolite or exotic breeds (Saanen, Camosciata delle Alpi and Maltese), which represent 24 % of recorded Italian goat heads. Among the remaining 41.7 % of breeds classified as at risk, 11 are endangered (number of heads less than 1200 with a declining trend) including Orobica, Valdostana and Ciociara Grigia breeds, and four are in a critical status (number of heads less than 100) as for example the Di Teramo breed.
Italian goats are mainly reared in the Alps and in the Mediterranean environments that are typical of Southern Italy and of the islands. In 2013, the regions that counted the highest proportion of goats were Calabria (43 % of animals), the two major islands Sardinia (18 %) and Sicily (10 %) in Southern Italy, and Lombardy (10 %) in the north (www.assonapa.com). These data from the national goat and sheep breeders association refer only to registered animals that, according to FAOSTAT, represent only 20 % of the goats reared in Italy.
Goat farming systems in Italy vary widely depending on region and on the breeds raised: in the north, two main different farming systems are present: (1) traditional farming, with an indoor system in the winter, natural grazing in the spring and autumn, and vertical transhumance in the summer; (2) intensive and semi-intensive indoor farming, with animals kept in flocks of medium to large size and reared under controlled feeding that includes hay and concentrate. The farming system for breeds in Central Italy is characterized by small sedentary herds that practice transhumance from spring to fall while that for breeds in Southern Italy is mainly characterized by small- to large-sized farms with either semi-sedentary farming based on natural pasture or free-ranging farming. In Sardinia, goat farming systems differ from the traditional ones and range from small infrastructures with low management costs to a semi-intensive system as for Maltese goats .
Several molecular studies, which often included limited numbers of breeds and loci, have been carried out on local Italian goat breeds [5–7] in attempts to monitor genetic erosion and identify conservation strategies. Large-scale surveys based on nuclear markers, such as microsatellites , AFLP (amplified fragment length polymorphisms)  or small panels of SNPs (single nucleotide polymorphisms) [9–11], detected a remarkable level of genetic diversity in European goats. Recently, a meta-analysis of worldwide goat microsatellite datasets highlighted a decreasing gradient of diversity from the domestication centre towards Europe and Asia, and a clear phylogeographic structure at both the continental and regional levels . In particular, breed formation seems to have been less systematic in the Middle East than in North-Central Europe, where several breeds are recognised as separate gene pools, partly as a consequence of inbreeding and partly of a strong genetic identity . On the contrary, results obtained from mitochondrial DNA analyses  revealed a weak phylogeographic structure, which suggested a long-lasting intercontinental gene flow as a consequence of the frequent translocation of goats along colonization, migration and commercial routes .
Recently, the availability of SNP panels [15–17] allows the investigation of livestock genomic diversity at a level of resolution that is impossible to reach with other types of markers. In this study, we exploited the medium-density (>50 000 SNPs) BeadChip available for goat  to assess the genome-wide diversity of 14 Italian goat breeds, as a contribution to biodiversity conservation and prioritisation actions.
Name of the breeds, breed acronyms, sample size before (n-PreQC) and after (n-PostQC) genotyping quality control procedures, expected heterozygosity (HE), observed heterozygosity (HO), Wright’s inbreeding coefficient (FIS), and proportion of polymorphic SNPs (PN)
Camosciata delle Alpi
Valpassiria or Passeirer Gebirgziege
Maltese from Sicily
Maltese from Sardinia
Blood samples were collected according to the recommendations of the European Council (1986) concerning animal care. Whole blood was collected in Vacutainer tubes with K-EDTA as anticoagulant and stored at −20 °C until genomic DNA was extracted using a commercial kit (NucleoSpin Blood, Macherey-Nagel) according to the manufacturer’s instructions.
Genotyping and SNP quality control
DNA samples were genotyped using the GoatSNP50 BeadChip (Illumina Inc., San Diego, CA) developed by the International Goat Genome Consortium (IGGC) . SNP typing was outsourced at the Associazione Italiana Allevatori - Laboratorio di Genetica e Servizi (http://www.lgscr.it) and Porto Conte Ricerche s.r.l. (Alghero, Sassari, Italy) facilities. Raw signal intensities of the 53 347 SNPs were converted into genotype calls with GenomeStudio software v2011.1 and by using the SNP genomic locations and cluster files made available by the IGGC. GenABEL ver. 1.7-6  was used for quality control (QC) procedures with standard thresholds i.e., a SNP call rate greater than 0.95, a MAF (minor allele frequency) greater than 0.01 and an individual genotype call rate greater than 0.95.
Data parsing, recoding and formatting were performed using PLINK ver. 1.07 suite . Expected (HE) and observed heterozygosities (HO) were calculated with an in-house script. Arlequin software ver. 126.96.36.199  (http://cmpg.unibe.ch/software/arlequin3) was used to (i) calculate population specific inbreeding coefficients (FIS); (ii) compute the FST [23, 24] distance matrix between breeds; and (iii) perform an Analysis of MOlecular VAriance (AMOVA, ) at different hierarchical levels to test the differentiation between breeds and between groups of breeds from distinct geographical areas (i.e., Northern Italy, Central Italy, Southern Italy and islands).
A Neighbor-network graph based on between-breeds FST distances was obtained with the software SplitsTree ver. 4.10 . The R package GenABEL ver. 1.7-6  was used to build a multi-dimensional scaling (MDS) plot, based on a matrix of (1-pairwise genomic kinship) distances between individuals. ADMIXTURE ver. 1.22  software was used for population structure analysis with a number of hypothetical pseudo-populations, K, that ranged from 2 to 25. To evaluate optimal partitioning, cross-validation (CV) error values were computed for each K using a 5-fold cross-validation procedure.
After quality control: (i) 551 markers with a MAF less than 0.01 and (ii) 1660 SNPs and four animals with call rates less than 0.95 were excluded. The final working dataset included 350 animals and 51 136 SNPs of which only 0.46 % carried rare alleles (0.01 ≤ MAF ≤ 0.05).
The percentage of within-breed polymorphic SNPs ranged from 95.13 % to 99.70 % (Table 1) with the highest values found for the SAA and CAM breeds (99.66 % and 99.70 %, respectively), which were both included in the SNP chip discovery panel and the lowest value found for the TER breed (95.13 %). Expected and observed heterozygosities ranged from 0.37 (VAL) to 0.41 (ARG) and from 0.36 (VAL) to 0.41 (SAA), respectively (Table 1). The only significant deviation from Hardy-Weinberg equilibrium was detected in the NIC breed, for which an inbreeding coefficient (FIS) of 0.074 (P < 0.05) revealed a slight excess of homozygotes. The pairwise FST values between breeds ranged from 0.013 to 0.164 [See Additional file 1: Table S1]. As a general trend, the ORO breed from Northern Italy and the TER breed from the centre of the peninsula showed the highest FST values, while ARG, NIC, CGI and ASP breeds from Southern and Central Italy showed the lowest FST values.
Although the 50 K SNP panel was developed from sequence data for goat breeds such as Saanen, Alpine, Creole, Boer, Kacang, and Savanna (http://www.goatgenome.org/), a large number of polymorphic SNPs were detected for the Italian breeds (although not all included in the discovery panel), which suggests that the impact of ascertainment bias is small in our dataset.
The analysis of genotyping data for Italian goats revealed: (i) high levels of genetic polymorphism, (ii) a limited amount of inbreeding, (iii) a geographical pattern that underlies the distribution of genomic diversity, and (iv) a moderate and weak population structure in the Northern and Central-Southern breeds, respectively. The latter two results can also be influenced by traditional management practices, recent demographical events and adaptation to different climatic conditions. However, a more extensive sampling of Italian goat breeds that would cover more evenly the geographical range is necessary to assess the impact of these factors.
Geographical distribution of Italian goat molecular diversity
The presence of a clear North–south geographical distribution of genetic diversity (Fig. 3) along the longitudinal axis of the Italian Peninsula was highlighted by both the first dimension of the MDS plot and the Neighbour-network reconstruction. A similar geographical pattern has been described in previous studies based on dominant or multiallelic markers in European goats and cattle [11, 28, 29] and on the IlluminaOvineChip50 in Italian sheep breeds .
Probably, isolation by geographical distance played a major role in shaping the differentiation of Italian goats, although it may have acted in synergy with other factors. Introgression from gene pools of animals that are native to foreign countries (e.g., Southern Italian breeds may have been crossed with other Southern European or Northern African breeds at the time of the domination of the Arabs) may have also played a role, together with adaptation to local environmental conditions, as already suggested by a genome-wide analysis based on AFLP markers . The effect of geographical isolation at the local scale can be seen in Northern Italian goat breeds that are traditionally reared in the Alps, where geographical barriers such as peaks and steep valleys can strongly limit or prevent gene flow. This is suggested by the barely visible or absent overlap between these breeds on the MDS plot. However, to better clarify the relative roles of all these factors, a comparison with a larger set of European and Mediterranean breeds is necessary.
Although the sizes of the populations analyzed were small, HE values were in line with those calculated for other breeds that are included in the SNP discovery panel [18, 32]. In addition, the difference between HO and HE was not significant for any of the investigated breeds, with the sole exception of NIC for which a slightly positive FIS was found. It is likely that pedigree records and occasional DNA-based controls helped to design mating plans that enabled farmers to control inbreeding. The institution of National Registers and Herd Books (e.g., for CAM in 1973, ORO in 1993, BIO in 1997) further contributed to preserve goat breeds from indiscriminate crossbreeding. This is particularly the case for breeds with distinctive phenotypic traits, such as horn shape in VAL and GIR, and horns and coat colours in the ORO breed.
Population structure of Italian goats
The analysis of population structure highlighted a moderate tendency for clustering for most of the Northern Italian breeds. This was also confirmed by the MDS plot and can be explained by the demographic history of these breeds that have been reared for a long time in geographically and culturally separate valleys. As a consequence, it is likely that these breeds experienced reproductive isolation and reduced gene flow and thus acquired a strong genetic identity . The effect of cultural and geographical separation due to alpine barriers has already been observed in the genetic makeup of some human populations in the Alps, which still today reflect the topographic features of these mountains . The VPS and BIO breeds are the main exception to this trend with a shared strong genetic component. This is probably due to the geographical proximity of the breeding areas of these two breeds, which facilitates the exchange of bucks and does. As already reported in the literature [10, 28, 35], our results from admixture analysis confirmed the strong genetic identity of the ORO breed: indeed at K = 2, this breed is clearly assigned to one of the two gene pools (average Q scores > 0.97; Fig. 4a), and at K = 4 [See Additional file 3: Figure S2], all ORO individuals are assigned to a separate gene pool. There are several possible explanations for this situation. The large number of monomorphic SNPs and rare alleles detected for this breed suggests that a strong drift effect has taken place. In fact, the ORO breed, is among the first populations of goats in Northern Italy to have experienced a strong reproductive isolation because of several specific phenotypic traits i.e., four well characterized coat colour patterns that differ from those of breeds in the nearby regions, long hair and a particular shape of horn) and also because of a dramatic demographic decrease of about 90 % in less than 15 years, in spite of the agricultural policy of the European Union (EU) to support native breeds at risk of extinction. Hopefully, new insights on the ORO origin as well as a renewed concern towards its conservation will come from the results of the complete genome sequencing of several unrelated individuals that is currently in progress (unpublished data).
Among the central Italian breeds investigated, TER had the smallest population size, with only 58 animals recorded by the National Breeders Association (www.assonapa.it). Nevertheless, an excess of TER heterozygotes was observed, which, together with the large number of rare SNPs identified for this breed and the wide area covered by TER individuals on the MDS plot, suggest that crossbreeding with bucks from other breeds occurred. However, since admixture analysis assigned TER individuals to a well-defined cluster at a K value as low as 3, which suggests a low level of admixture, if a recent episode of introgression occurred, it probably involved a donor gene pool from a breed that is not included in our dataset.
CGI from Central Italy is also at risk of extinction; in 2008, when the national Register was created, this breed counted only 181 animals that were distributed across a few farms in the Lazio region. Thanks to local (Regolamento (European Council EC) 1698/2005 - Programma di Sviluppo Rurale 2007–2013) and EU economic support, after five years, the population size had increased by 3.7 folds and reached 674 animals in 2013 (www.assonapa.it). The distribution of individual points on the MDS plot and an HE of 0.402 indicate that there is a certain level of genetic variability within this breed. According to admixture analysis, this breed shares a common genomic background with three Southern Italian populations (NIC, ASP and ARG). Although these three breeds are bred in non-contiguous geographical areas - i.e., Southern Lazio for CGI, Calabria for NIC and ASP and Sicily for ARG –, their pairwise FST distances were the lowest, which can be partly explained by the impact of transhumance in Central and Southern Italy  and the trade or occasional exchange of bucks and does between these breeds, which share a similar grey coat colour.
NIC is the only breed that showed a significant FIS value. In 2013, it counted 4975 heads as a result of the adoption in Southern Italy of the Council Regulation (EC) No 1698/2005 to support rural development that includes the presence of autochthonous breeds in danger of extinction; this is funded by the European Agricultural Fund for Rural Development (EAFRD). The absence of male rotational schemes between farms has probably caused inbreeding at the farm level (www.goatit.eu), which shows that conservation strategies that are based exclusively on remuneration of farmers are not efficient when proper breeding schemes and extension services are lacking .
Besides demographic factors, historical factors are likely to have an impact on the genetic makeup of Southern Italian and insular breeds. In fact, since the Neolithic age, the area that covers North Africa, Malta, Sicily, Sardinia, and Southern Italy has been deeply interconnected by routes of migration, trade and conquest. The tightly interconnected history of these areas was favoured by their role of crossroads of the trading and warring routes through the land-locked Mediterranean Sea .
The differentiation between Maltese goats that were sampled in Sardinia and Sicily as revealed by admixture and Neighbour-network analyses, is probably the consequence of the combined effects of genetic drift, small population size, founder effects and reproductive isolation (since the Southern Tyrrhenian Sea acted as a strong genetic barrier). However, the two populations were not separated on the MDS plot based on a genomic kinship distance matrix. In fact, this measure may be showing that these two populations share a common origin since it can account for older evolutionary relationships and is less affected by recent population dynamics . Anyway, the ancestral origin of the Maltese breed as a whole remains uncertain and previous studies suggested that the distant roots of this breed may be in the Middle-Eastern side of the Mediterranean basin and that it probably derives from crosses between North African and Italian goats [39, 40]. To test this hypothesis, it is necessary to compare the MAL breed with a larger set of goat breeds from the Mediterranean area, Northern Africa and Spain.
According to our results, the axis of the main source of genetic variation for Italian goat populations stretches along the longitudinal axis of the Italian peninsula. Among the major factors which could have acted on the goat genome in Italy, reproductive isolation due to geographical distance, adaptation to local conditions and breeders’ management may all have a role, although the identification of their relative impacts and contribution is not straightforward. To assess the importance of any of these factors, a larger and more comprehensive set of breeds is necessary, in particular to have a more uniform distribution of sampling locations. Nevertheless, some breeds displayed a clear genetic identity that confirms previous findings [9–11, 28, 35], although its source is not always completely understood e.g., the ORO breed. The adoption of European or National conservation policies has boosted the population size of breeds at risk of extinction and probably affected the recent evolutionary history of goat populations in Italy. However, the lack of extensive records on demographic trajectories makes it difficult to either confirm or challenge this hypothesis. Based on these observations, it is clear that understanding the extent, distribution and origin of present-day genetic diversity is a complex task that requires other sources of information than molecular data only.
A thorough genomic characterization of breeds represents a key point to develop efficient conservation strategies, which, to become effective, should also take into account population viability. This variable, in turns, depends on budgetary limitations, management practices and on the existence of services to support shepherds.
Based on the use of a standardized genotyping array, such as the GoatSNP50 BeadChip, it will be possible to combine various datasets and to provide a global picture of goat genetic diversity both at a local and global scale. This will help to understand the origins of genetic diversity and to manage biodiversity of these animal genetic resources that are particularly relevant for poor and marginal rural areas of the world.
This work was funded by the Italian Ministry of Agriculture (grants INNOVAGEN).
- Clutton-Brock J. A natural history of domesticated mammals. Cambridge: Cambridge University Press; 1999.Google Scholar
- Naderi S, Rezaei HR, Pompanon F, Blum MGB, Negrini R, Naghash HR, et al. The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals. Proc Natl Acad Sci U S A. 2008;105:17659–64.PubMed CentralPubMedView ArticleGoogle Scholar
- FAO. Status and trends of animal genetic resources-2012. Commission on genetic resources for food and agriculture. 15–19 April 2013: Rome. http://www.fao.org/docrep/meeting/027/mg046e.pdf. 15 Jan 2015 Accessed
- Usai MG, Casu S, Molle G, Decandia M, Ligios S, Carta A. Using cluster analysis to characterize the goat farming system in Sardinia. Livest Sci. 2006;104:63–76.View ArticleGoogle Scholar
- Iamartino D, Bruzzone A, Lanza A, Blasi M, Pilla F. Genetic diversity of Southern Italian goat population assessed by microsatellite markers. Small Ruminant Res. 2005;57:249–55.View ArticleGoogle Scholar
- Sechi T, Usai MG, Miari S, Mura L, Casu S, Carta A. Identifying native animals in crossbred populations: the case of the Sardinian goat population. Anim Genet. 2007;38:614–20.PubMedView ArticleGoogle Scholar
- Ajmone-Marsan P, Negrini R, Crepaldi P, Milanesi E, Gorni C, Valentini A, et al. Assessing genetic diversity in Italian goat populations using AFLP® markers. Anim Genet. 2001;32:281–8.PubMedView ArticleGoogle Scholar
- Cañón J, García D, García-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.PubMedView ArticleGoogle Scholar
- Pariset L, Cappuccio I, Ajmone-Marsan P, Dunner S, Luikart G, England PR, et al. Assessment of population structure by single nucleotide polymorphisms (SNPs) in goat breeds. J Chromatogr B. 2006;833:117–20.View ArticleGoogle Scholar
- Pariset L, Joost S, Ajmone-Marsan P, Valentini A, Econogene Consortium. Landscape genomics and biased FST approaches reveal single nucleotide polymorphisms under selection in goat breeds of North-East Mediterranean. BMC Genet. 2009;10:7.PubMed CentralPubMedView ArticleGoogle 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 CentralPubMedView ArticleGoogle Scholar
- Ajmone-Marsan P, Colli L, Han JL, Achilli A, Lancioni H, Joost S, et al. The characterization of goat genetic diversity: Towards a genomic approach. Small Ruminant Res. 2014;121:58–72.View ArticleGoogle Scholar
- Amills M, Ramírez O, Tomàs A, Badaoui B, Marmi J, Acosta J, et al. Mitochondrial DNA diversity and origins of South and Central American goats. Anim Genet. 2009;40:315–22.PubMedView ArticleGoogle Scholar
- Luikart G, Gielly L, Excoffier L, Vigne JD, Bouvet J, Taberlet P. Multiple maternal origins and weak phylogeographic structure in domestic goat. Proc Natl Acad Sci U S A. 2001;98:5927–32.PubMed CentralPubMedView ArticleGoogle Scholar
- ovine HapMap Consortium, Gibbs RA, Taylor JF, Van Tassell CP, Barendse W, Eversole KA, et al. Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science. 2009;324:528–32.View ArticleGoogle Scholar
- Kijas JW, Lenstra JA, Hayes B, Boitard S, Porto Neto LR, San Cristobal M, et al. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2012;10, e1001258.PubMed CentralPubMedView ArticleGoogle Scholar
- Muir WM, Wong GK, Zhang Y, Wang J, Groenen MA, Crooijmans RP, et al. Genome-wide assessment of worldwide chicken SNP genetic diversity indicates significant absence of rare alleles in commercial breeds. Proc Natl Acad Sci U S A. 2008;105:17312–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Tosser-Klopp G, Bardou P, Bouchez O, Cabau C, Crooijmans R, Dong Y, et al. Design and characterization of a 52 K SNP chip for goats. PLoS One. 2014;9, e86227.PubMed CentralPubMedView ArticleGoogle Scholar
- FAO. World watch list for domestic animal diversity. Rome: FAO; 2000.Google Scholar
- Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007;23:1294–6.PubMedView ArticleGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.PubMed CentralPubMedView ArticleGoogle Scholar
- Excoffier L, Lischer HEL. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.PubMedView ArticleGoogle Scholar
- Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.View ArticleGoogle Scholar
- Wright S. The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution. 1965;19:395–420.View ArticleGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131:479–91.PubMed CentralPubMedGoogle Scholar
- Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:254–67.PubMedView ArticleGoogle Scholar
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.PubMed CentralPubMedView ArticleGoogle Scholar
- Negrini R, D’Andrea M, Crepaldi P, Colli L, Nicoloso L, Guastella AM, et al. Effect of microsatellite outliers on the genetic structure of eight Italian goat breeds. Small Ruminant Res. 2012;103:99–107.View ArticleGoogle Scholar
- Negrini R, Nijman IJ, Milanesi E, Moazami-Goudarzi K, Williams JL, Erhardt G, et al. Differentiation of European cattle by AFLP fingerprinting. Anim Genet. 2007;38:60–6.PubMedView ArticleGoogle Scholar
- Ciani E, Crepaldi P, Nicoloso L, Lasagna E, Sarti FM, Moioli B, et al. Genome-wide analysis of Italian sheep diversity reveals a strong geographic pattern and cryptic relationships between breeds. Anim Genet. 2014;45:256–66.PubMedView ArticleGoogle Scholar
- Colli L, Joost S, Negrini R, Nicoloso L, Crepaldi P, Ajmone-Marsan P, et al. Assessing the spatial dependence of adaptive loci in 43 European and Western Asian goat breeds using AFLP markers. PLoS One. 2014;9, e86668.PubMed CentralPubMedView ArticleGoogle Scholar
- Kijas JW, Ortiz JS, McCulloch R, James A, Brice B, Swain B, et al. Genetic diversity and investigation of polledness in divergent goat populations using 52 088 SNPs. Anim Genet. 2013;44:325–35.PubMedView ArticleGoogle Scholar
- Bigi D, Zanon A. Atlante delle razze autoctone: bovini, equini, ovicaprini, suini allevati in Italia. Milano: Edagricole; 2008.Google Scholar
- Thomas MG, Barnes I, Weale ME, Jones AL, Forster P, Bradman N, et al. New genetic evidence supports isolation and drift in the Ladin communities of the South Tyrolean Alps but not an ancient origin in the Middle East. Eur J Hum Genet. 2008;16:124–34.PubMedView ArticleGoogle Scholar
- Crepaldi P, Negrini R, Milanesi E, Gorni C, Cicogna M, Ajmone-Marsan P. Diversity in five goat populations of the Lombardy Alps: comparison of estimates obtained from morphometric traits and molecular markers. J Anim Breed Genet. 2001;118:173–80.View ArticleGoogle Scholar
- FAO. The state of the world’s animal genetic resources for food and agriculture. Edited by B Rischkowsky and D Pilling. Rome. 2007. http://www.fao.org/docrep/010/a1250e/a1250e00.htm. 15 Jan 2015 Accessed
- Horden P, Purcell N. The Corrupting Sea: a study of Mediterranean history. London: Blackwell Publishers; 2000.Google Scholar
- Meuwissen TH, Odegard J, Andersen-Ranberg I, Grindflek E. On the distance of genetic relationships and the accuracy of genomic prediction in pig breeding. Genet Sel Evol. 2014;46:49.PubMed CentralPubMedView ArticleGoogle Scholar
- Porter V. Goats of the World. Ipswich: Farming Press; 1996.Google Scholar
- Gall C. Goat breeds of the world. Weikersheim: Margraf Verlag; 1996.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.