Northern Breeds Research Articles

Comparative locomotor costs of domestic dogs reveal energetic economy of wolf-like breeds.

The broad diversity in morphology and geographic distribution of the 35 free-ranging members of the family Canidae is only rivaled by that of the domesticated dog, Canis lupus familiaris. Considered to be among nature's most elite endurance athletes, both domestic and wild canids provide a unique opportunity to examine the variability in mammalian aerobic exercise performance and energy expenditure. To determine the potential effects of domestication and selective breeding on locomotor gait and economy in canids, we measured the kinematics and mass-specific metabolism of three large (>20 kg) dog breed groups (northern breeds, retrievers and hounds) of varying morphological and genomic relatedness to their shared progenitor, the gray wolf. By measuring all individuals moving in preferred steady-state gaits along a level transect and on a treadmill, we found distinct biomechanical, kinematic and energetic patterns for each breed group. While all groups exhibited reduced total cost of transport (COT) at faster speeds, the total COT and net COT during trotting and galloping were significantly lower for northern breed dogs (3.0 and 2.1 J kg-1m-1, respectively) relative to hound (4.2 and 3.4 J kg-1m-1, respectively) and retriever dogs (3.8 and 3.0 J kg-1m-1, respectively) of comparable mass. Similarly, northern breeds expended less energy per stride (3.5 J kg-1stride-1) than hounds or retrievers (5.0 and 4.0 J kg-1stride-1, respectively). These results suggest that, in addition to their close genetic and morphological ties to gray wolves, northern breed dogs have retained highly cursorial kinematic and physiological traits that promote economical movement across the landscape.

Genetic diversity of Italian goat breeds assessed with a medium-density SNP chip.

BackgroundAmong 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.ResultsAnalysis 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).ConclusionsOur 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-015-0140-6) contains supplementary material, which is available to authorized users.

Loess and Bee-Eaters II: The ‘loess’ of North Africa and the nesting behaviour of the Northern Carmine Bee-Eater (Merops nubicus Gmelin 1788)

The Northern Carmine Bee-Eater (Merops nubicus) lives and breeds in a well demarcated region stretching across Africa close to the 15°N line of latitude. The Bee-Eater zone appears to be associated with a band of loess, defined by Scheidig on his 1934 map as second-order loess. Bee-eaters are known to favour loess for nesting tunnels and it appears that the 15°N material is sufficiently loess-like. Obvious sources for particulate materials for the 15°N band are the Fonta-Djalon highlands which supply sedimentary material to the River Niger; the Bodele Depression, the deepest part of Lake Megachad, source of dust for the World; the Ethiopian highlands at the eastern end of 15°N which supply silt to the Nile system and particulates to the 15°N region. In soil moisture terms the region is ustic, which is possibly a necessary condition for bee-eater nests. The clastic material requires an ustic environment. The River Niger can be seen as a loess river; in some senses a mirror-image of a major loess river like the Danube; but where a restricted range of particle inputs leads to a restricted range of loess deposit outputs. Nevertheless loess river considerations can be applied. The Niger delivers second-order loess and an important loessic admixture to the landscape. Enough loess for selective nesters like the Carmine Bee-Eaters to build their nest tunnels in it. It seems likely that climate change will cause a change in bee-eater distribution; it seems unlikely that they will abandon their nesting regions, the living and wintering zones may shift.

A genome-wide scan for signatures of selection in Chinese indigenous and commercial pig breeds

BackgroundModern breeding and artificial selection play critical roles in pig domestication and shape the genetic variation of different breeds. China has many indigenous pig breeds with various characteristics in morphology and production performance that differ from those of foreign commercial pig breeds. However, the signatures of selection on genes implying for economic traits between Chinese indigenous and commercial pigs have been poorly understood.ResultsWe identified footprints of positive selection at the whole genome level, comprising 44,652 SNPs genotyped in six Chinese indigenous pig breeds, one developed breed and two commercial breeds. An empirical genome-wide distribution of Fst (F-statistics) was constructed based on estimations of Fst for each SNP across these nine breeds. We detected selection at the genome level using the High-Fst outlier method and found that 81 candidate genes show high evidence of positive selection. Furthermore, the results of network analyses showed that the genes that displayed evidence of positive selection were mainly involved in the development of tissues and organs, and the immune response. In addition, we calculated the pairwise Fst between Chinese indigenous and commercial breeds (CHN VS EURO) and between Northern and Southern Chinese indigenous breeds (Northern VS Southern). The IGF1R and ESR1 genes showed evidence of positive selection in the CHN VS EURO and Northern VS Southern groups, respectively.ConclusionsIn this study, we first identified the genomic regions that showed evidences of selection between Chinese indigenous and commercial pig breeds using the High-Fst outlier method. These regions were found to be involved in the development of tissues and organs, the immune response, growth and litter size. The results of this study provide new insights into understanding the genetic variation and domestication in pigs.

Genetic Diversity in Exon 2 of the Major Histocompatibility Complex Class II DQB1 Locus in Nigerian Goats

The DQB1 locus is located in the major histocompatibility complex (MHC) class II region and involved in immune response. We identified 20 polymorphic sites in a 228 bp fragment of exon 2, one of the most critical regions of the MHC DQB1 gene, in 60 Nigerian goats. Four sites are located in the peptide binding region, and 10 amino acid substitutions are peculiar to Nigerian goats, compared with published sequences. A significantly higher ratio of nonsynonymous/synonymous substitutions (dN/dS) suggests that allelic sequence evolution is driven by balancing selection (P < 0.01). In silico functional analysis using PANTHER predicted that substitution P56R, with a subPSEC score of -4.00629 (Pdeleterious = 0.73229), is harmful to protein function. The phylogenetic tree from consensus sequences placed the two northern breeds closer to each other than either was to the southern goats. This first report of sequence diversity at the DQB1 locus for any African goat breed may be useful in the search for disease-resistant genotypes.

Identification of Genetic Variation on the Horse Y Chromosome and the Tracing of Male Founder Lineages in Modern Breeds

The paternally inherited Y chromosome displays the population genetic history of males. While modern domestic horses (Equus caballus) exhibit abundant diversity within maternally inherited mitochondrial DNA, no significant Y-chromosomal sequence diversity has been detected. We used high throughput sequencing technology to identify the first polymorphic Y-chromosomal markers useful for tracing paternal lines. The nucleotide variability of the modern horse Y chromosome is extremely low, resulting in six haplotypes (HT), all clearly distinct from the Przewalski horse (E. przewalskii). The most widespread HT1 is ancestral and the other five haplotypes apparently arose on the background of HT1 by mutation or gene conversion after domestication. Two haplotypes (HT2 and HT3) are widely distributed at high frequencies among modern European horse breeds. Using pedigree information, we trace the distribution of Y-haplotype diversity to particular founders. The mutation leading to HT3 occurred in the germline of the famous English Thoroughbred stallion “Eclipse” or his son or grandson and its prevalence demonstrates the influence of this popular paternal line on modern sport horse breeds. The pervasive introgression of Thoroughbred stallions during the last 200 years to refine autochthonous breeds has strongly affected the distribution of Y-chromosomal variation in modern horse breeds and has led to the replacement of autochthonous Y chromosomes. Only a few northern European breeds bear unique variants at high frequencies or fixed within but not shared among breeds. Our Y-chromosomal data complement the well established mtDNA lineages and document the male side of the genetic history of modern horse breeds and breeding practices.

Genetic Variability of PRNP in Chinese Indigenous Goats

Polymorphism of the prion protein gene (PRNP) is usually associated with scrapie susceptibility or resistance. To determine the variability of PRNP in Chinese indigenous goat breeds, we isolated genomic DNA from goat blood and amplified and sequenced the coding region of the gene. We identified 10 polymorphic sites that gave rise to 28 haplotypes. Clear frequency differences were found between northern and southern breeds and confirmed by genetic distance analysis, except for the Tangshan dairy goat. Phylogeographic analysis supported the idea that northern and southern breeds might be considered separate clusters, except for the Tangshan dairy goat. The finding of significant differences in allele distribution in northern and southern goats, especially if involved in modulating resistance/susceptibility, needs to be carefully considered for the feasibility of selection plans for resistance to scrapie.

Regional occurrence, high frequency but low diversity of mitochondrial DNA haplogroup d1 suggests a recent dog-wolf hybridization in Scandinavia

The domestic dog mitochondrial DNA (mtDNA)-gene pool consists of a homogenous mix of haplogroups shared among all populations worldwide, indicating that the dog originated at a single time and place. However, one small haplogroup, subclade d1, found among North Scandinavian/Finnish spitz breeds at frequencies above 30%, has a clearly separate origin. We studied the genetic and geographical diversity for this phylogenetic group to investigate where and when it originated and whether through independent domestication of wolf or dog-wolf crossbreeding. We analysed 582 bp of the mtDNA control region for 514 dogs of breeds earlier shown to harbour d1 and possibly related northern spitz breeds. Subclade d1 occurred almost exclusively among Swedish/Finnish Sami reindeer-herding spitzes and some Swedish/Norwegian hunting spitzes, at a frequency of mostly 60–100%. Genetic diversity was low, with only four haplotypes: a central, most frequent, one surrounded by two haplotypes differing by an indel and one differing by a substitution. The substitution was found in a single lineage, as a heteroplasmic mix with the central haplotype. The data indicate that subclade d1 originated in northern Scandinavia, at most 480–3000 years ago and through dog-wolf crossbreeding rather than a separate domestication event. The high frequency of d1 suggests that the dog-wolf hybrid phenotype had a selective advantage.

Matrilineal genetic inter-introgression of Bos taurus and Bos indicus in China

The genetic resources of cattle breeds in China are very abundant, with the existence of 28 indigenous and other cattle breeds. To further understand the population structure of these indigenous breeds from different geographical locations, we carried out a phylogenetic analysis in maternal lineages by employing mtDNA polymorphisms. The entire mtDNA D-loop sequences of 73 individuals from 14 representative native cattle breeds were determined and analyzed in conjunction with the published 144 D-loop sequences for 19 Chinese cattle breeds retrieved from GenBank. These sequences belonged to 100 haplotypes defined by polymorphisms at 122 polymorphic sites. These haplotypes were clustered into two distinct lineages: taurine Bos taurus and zebu Bos indicus. The two main peaks of the poly-peaked mismatch curve revealed the inter-mismatch distribution between the Bos taurus and Bos indicus sequence types within the cattle samples. The zebu mtDNA dominated in the southern breeds (62.5%), while taurine mtDNA dominated in the northern breeds (96.7%). The geographical distribution of these frequencies revealed an inter-introgression pattern of Bos taurus and Bos indicus in maternal lineages in China. The significant genetic divergence between the northern and the southern breeds resulted from the dominant geographical segregation of the Qinling Moutains and temperature weather conditions. The heterogeneous distribution of taurine and zebu types in other regions could be attributed to the historical events and introgression of taurine cattle imported from home and abroad.

Molecular genetic variation in sheep of the central Volga area inhabited by Finno-Ugric peoples

Based on morphology, native northern European sheep breeds belong to the short tailed type, of which the Romanov is the only native example still distributed in northwest Russia. Besides this, there exist local sheep populations kept by Finno-Ugric peoples in the central Volga region, which represent additional genetic resources in the area. Four sheep populations from the central Volga region were genotyped for 20 microsatellites and compared with geographically proximate breeds (Estonian Whitehead and Blackhead, the Finnsheep and an exported and a native population of Russian Romanov) and with local populations in Estonia, Finland and Russian Karelia. Between-breed analyses including admixture analysis using molecular genetic markers and the phenotypic characteristics indicated that the Volgaic populations have not remained pure. The Viena population from Russian Karelia, the Romanov breed and, to some extent, the Komi population, have escaped extensive mixing, making them most attractive for conservation programmes. The study compared imported and native Romanov breed populations and the results suggest that the diversity parameters are markedly similar in these two populations.;

Genetic variability of the PRNP gene in goat breeds from Northern and Southern Italy

To determine the variability of the prion protein gene (PRNP) in goats from Northern and Southern Italy. Genomic DNA isolated from goat blood was polymerase chain reaction (PCR)-amplified for the coding region of the PRNP gene and then sequenced. In total, 13 polymorphic sites were identified: G37V, T110P, G127S, M137I, I142M, I142T, H143R, R154H, P168Q, T194P, R211Q, Q222K and S240P (substitutions I142T and T194P are novel) giving rise to 14 haplotypes. Clear frequency differences between Northern and Southern breeds were found and confirmed by genetic distance analysis. Differences in allele distribution were found between Northern and Southern goats, in particular regarding the M142 and K222 alleles, possibly associated to scrapie resistance; philogeographical analysis supported the idea that Northern and Southern breeds may be considered as separate clusters. In Italy only limited studies have been carried out on caprine PRNP genotype distribution; this study is important to fill this lack of information. Moreover the finding of significant differences among allele distributions in Northern and Southern goats, especially if involved in modulating resistance/susceptibility, need to be carefully considered for the feasibility of selection plans for resistance to scrapie.

Barking up the wrong tree: Modern northern European dogs fail to explain their origin

BackgroundGeographic distribution of the genetic diversity in domestic animals, particularly mitochondrial DNA, has often been used to infer centers of domestication. The underlying presumption is that phylogeographic patterns among domesticates were established during, or shortly after the domestication. Human activities are assumed not to have altered the haplogroup frequencies to any great extent. We studied this hypothesis by analyzing 24 mtDNA sequences in ancient Scandinavian dogs. Breeds originating in northern Europe are characterized by having a high frequency of mtDNA sequences belonging to a haplogroup rare in other populations (HgD). This has been suggested to indicate a possible origin of the haplogroup (perhaps even a separate domestication) in central or northern Europe.ResultsThe sequences observed in the ancient samples do not include the haplogroup indicative for northern European breeds (HgD). Instead, several of them correspond to haplogroups that are uncommon in the region today and that are supposed to have Asian origin.ConclusionWe find no evidence for local domestication. We conclude that interpretation of the processes responsible for current domestic haplogroup frequencies should be carried out with caution if based only on contemporary data. They do not only tell their own story, but also that of humans.

Genetic structure of European sheep breeds

Large-scale evaluations of genetic diversity in domestic livestock populations are necessary so that region-specific conservation measures can be implemented. We performed the first such survey in European sheep by analysing 820 individuals from 29 geographically and phenotypically diverse breeds and a closely related wild species at 23 microsatellite loci. In contrast to most other domestic species, we found evidence of widespread heterozygote deficit within breeds, even after removing loci with potentially high frequency of null alleles. This is most likely due to subdivision among flocks (Wahlund effect) and use of a small number of rams for breeding. Levels of heterozygosity were slightly higher in southern than in northern breeds, consistent with declining diversity with distance from the Near Eastern centre of domestication. Our results highlight the importance of isolation in terms of both geography and management in augmenting genetic differentiation through genetic drift, with isolated northern European breeds showing the greatest divergence and hence being obvious targets for conservation. Finally, using a Bayesian cluster analysis, we uncovered evidence of admixture between breeds, which has important implications for breed management.

Selection on MHC-linked microsatellite loci in sheep populations

Microsatellites within the major histocompatibility complex (MHC) region have received increasing attention as proxy measures of the level of polymorphism at the Mhc genes themselves. We assessed the diversity of microsatellite loci within or in close proximity of the Mhc genes in several breeds of domestic sheep (Ovis aries) and the wild Mouflon (Ovis orientalis musimon). This was compared to variation at other microsatellite loci scattered throughout the sheep genome. Significantly higher number of alleles were observed at the MHC microsatellites. The sheep breeds studied fell into high- and low-diversity group. This grouping is not related to the agricultural use of the breeds, whether for milk, meat or wool. It is, however, correlated with the geographic origins of the breeds. Southern breeds are genetically more diverse than northern breeds. The observed heterozygosity was in most cases lower than Hardy-Weinberg expectations. The potential impact of selective breeding by man on this is discussed. Neutrality tests indicated that for most of the breeds, the distribution of alleles at the MHC-linked microsatellites are more even than would be expected if the genes were neutral and sampled from populations under drift-mutation equilibrium. Hitchhiking due to tight linkage with alleles at the MHC loci that are under balancing selection is proposed as a possible explanation for this pattern.

MtDNA Diversity and genetic lineages of eighteen cattle breeds from Bos taurus and Bos indicus in China

In order to clarify the origin and genetic diversity of indigenous cattle breeds in China, we carried out phylogenetic analysis of representatives of those breeds by employing mitochondrial gene polymorphism. Complete cyt b gene sequences, 1140 bp in length, were determined for a total of 136 individuals from 18 different breeds and these sequences were clustered into two distinct genetic lineages: taurine (Bos taurus) and zebu (Bos indicus). In analysis of the cyt b gene diversity, Chinese cattle showed higher nucleotide (0.00923) and haplotype diversity (0.848) than the reports from other studies, and the animals from the taurine lineage indicated higher nucleotide diversity (0.00330) and haplotype diversity (0.746) than the ones from the zebu lineage (0.00136; 0.661). The zebu mtDNA dominated in the southern breeds (63.3-100%), while the taurine dominated in the northern breeds (81.8-100%). Six cattle breeds from the central area of China exhibited intermediate frequencies of zebu mtDNA (25-71.4%). This polymorphism revealed a declining south-to-north gradient of female zebu introgression and a geographical hybrid zone of Bos taurus and Bos indicus in China.

Polymorphism of goat PrP gene in Northern and Southern Italian breeds

Scrapie in sheep is an infectious disease with genetic susceptibility playing an important role. The European Union applied genetics to control scrapie in the sheep population. In goats no clear association between PrP genotypes and resistance to scrapie has been demonstrated yet. Beside investigating this further, it is important to know the PrP genotype distribution in goat breeds. Actually, to estimate the frequency of interesting alleles in a population is a necessary preliminary step to understand how feasible it could be to select them. Aim of this work is to analyse the variability of the PrP gene in goat breeds from a Northern and a Southern Italian regions.Blood samples were collected from 300 goats of the four most common goat breeds reared in Piedmont region, belonging to 81 herds, and from 178 goats of 7 breeds from Basilicata, belonging to 19 herds. PrP gene polymorphisms were detected by direct DNA sequencing on both strands of the PCR products.In total, thirteen point mutations were identified: G37V, T110P, G127S, M137I, I142M, I142T, H143R, R154H, P168Q, T194P, R211Q, Q222K and S240P (the substitutions I142T and T194P are novel). This variability gave rise to 14 alleles. The allele with S127 mutation and the allele with Q211 mutation were found only in Northern goats while V37, I137, T142, R143 and P194 alleles were detected only in Southern goats.Regarding the alleles that were present in both groups, the M142 allele was statistically associated with Northern goats (χ2=39.84, p-value=<0.001) while K222 and H154 alleles were significantly associated with breeds from the South (χ2=54.75, p-value=<0.001; χ2=5.68, p-value=<0.05, respectively). The substitutions T110P, M137I, I142T, P168Q, and T194P were found at a low rate in both Southern and Northern goats and no statistically significant differences were observed between the groups.The present study revealed a high variability in the caprine PrP gene and clear frequency differences between Northern and Southern breeds, with mutations exclusively or significantly associated with one group or the other. M142 and K222 are known to be related to protection against scrapie in goats. More case-control studies are necessary to confirm the protective role of these variants; in fact if this role is confirmed it will allow the introduction of breeding programme to select for resistance to TSEs also in goats.

Prioritization for Conservation of Northern European Cattle Breeds Based on Analysis of Microsatellite Data

Northern European indigenous cattle breeds are currently endangered and at a risk of becoming extinct. We analyzed variation at 20 microsatellite loci in 23 indigenous, 3 old imported, and 9 modern commercial cattle breeds that are presently distributed in northern Europe. We measured the breeds' allelic richness and heterozygosity, and studied their genetic relationships with a neighbor-joining tree based on the Chord genetic distance matrix. We used the Weitzman approach and the core set diversity measure of Eding et al. (2002) to quantify the contribution of each breed to the maximum amount of genetic diversity and to identify breeds important for the conservation of genetic diversity. We defined 11 breeds as a "safe set" of breeds (not endangered) and estimated a reduction in genetic diversity if all nonsafe (endangered) breeds were lost. We then calculated the increase in genetic diversity by adding one by one each of the nonsafe breeds to the safe set (the safe-set-plus-one approach). The neighbor-joining tree grouped the northern European cattle breeds into Black-and-White type, Baltic Red, and Nordic cattle groups. Väne cattle, Bohus Poll, and Danish Jersey had the highest relative contribution to the maximum amount of genetic diversity when the diversity was quantified by the Weitzman diversity measure. These breeds not only showed phylogenetic distinctiveness but also low within-population variation. When the Eding et al. method was applied, Eastern Finncattle and Lithuanian White Backed cattle contributed most of the genetic variation. If the loss of the nonsafe set of breeds happens, the reduction in genetic diversity would be substantial (72%) based on the Weitzman approach, but relatively small (1.81%) based on the Eding et al. method. The safe set contained only 66% of the observed microsatellite alleles. The safe-set-plus-one approach indicated that Bohus Poll and Väne cattle contributed most to the Weitzman diversity, whereas the Eastern Finncattle contribution was the highest according to the Eding et al. method. Our results indicate that both methods of Weitzman and Eding et al. recognize the importance of local populations as a valuable resource of genetic variation.

Geographical partitioning of goat diversity in Europe and the Middle East

Thirty microsatellite markers were analysed in 1426 goats from 45 traditional or rare breeds in 15 European and Middle Eastern countries. In all populations inbreeding was indicated by heterozygosity deficiency (mean FIS = 0.10). Genetic differentiation between breeds was moderate with a mean FST value of 0.07, but for most (c. 71%) northern and central European breeds, individuals could be assigned to their breeds with a success rate of more than 80%. Bayesian-based clustering analysis of allele frequencies and multivariate analysis revealed at least four discrete clusters: eastern Mediterranean (Middle East), central Mediterranean, western Mediterranean and central/northern Europe. About 41% of the genetic variability among the breeds could be explained by their geographical origin. A decrease in genetic diversity from the south-east to the north-west was accompanied by an increase in the level of differentiation at the breed level. These observations support the hypothesis that domestic livestock migrated from the Middle East towards western and northern Europe and indicate that breed formation was more systematic in north-central Europe than in the Middle East. We propose that breed differentiation and molecular diversity are independent criteria for conservation.

Native breeds demonstrate high contributions to the molecular variation in northern European sheep

Population contribution to genetic diversity can be estimated using neutral variation. However, population expansion or hybridization of diverged ancestries may weaken correlation between neutral and non-neutral variation. Microsatellite variation was studied at 25 loci in 20 native and 12 modern or imported northern European sheep breeds. Breed contributions to total gene diversity, allelic richness and mean allele-sharing distance between individuals were measured. Indications of changes in population size and admixtures of divergent ancestries were investigated and the extent of inbreeding was estimated. The northern European sheep demonstrated signs of reduction in effective population size. Many old, small populations made a substantial positive contribution to total molecular variation, but populations with several divergent major ancestries did not contribute substantially to molecular variation, with the exception of the Norwegian Rygja sheep. However, several diverged major ancestries may cause it to contribute less to non-neutral variation than expected from the microsatellite data. Breed uniqueness and within-breed variability generally had opposite effects on breed contributions to molecular diversity. The degree of inbreeding did not reflect the breed contribution to total gene diversity or allelic richness, but inbred populations increased the mean allele-sharing distance between individuals. Our study indicates breed conservation to be especially important in maintaining allelic variation in northern European sheep and supports the evolutionary importance of peripheral populations.

Casein Haplotype Structure in Five Italian Goat Breeds

The aim of this work was to investigate the genetic structure of the casein gene cluster in 5 Italian goat breeds and to evaluate the haplotype variability within and among populations. A total of 430 goats from Vallesana, Roccaverano, Jonica, Garganica, and Maltese breeds were genotyped at αs1-casein (CSN1S1), αs2-casein, (CSN1S2), β-casein (CSN2), and κ-casein (CSN3) loci using several genomic techniques and milk protein analysis. Casein haplotype frequencies were estimated for each breed. Principal component analysis was carried out to highlight the relationship among breeds. Allele and haplotype distributions indicated considerable differences among breeds. The haplotype CSN1S1*F- CSN1S2*F-CSN3*D occurred in all breeds with frequencies >0.100 and was the most common haplotype in the Southern breeds. A high frequency of CSN1S1*0-CSN1S2*C-CSN3*A haplotype was found in Vallesana population (0.162). Principal component analysis clearly separated the Northern and Southern breeds by the first component. The variability of the caprine casein loci and variety of resulting haplotypes should be exploited in the future using specific breeding programs aiming to preserve biodiversity and to select goat genetic lines for specific protein production.