Assessing the contribution of breeds to genetic diversity in conservation schemes

The quantitative assessment of genetic diversity within and between populations is important for decision making in genetic conservation plans. In this paper we define the genetic diversity of a set of populations, S, as the maximum genetic variance that can be obtained in a random mating population that is bred from the set of populations S. First we calculated the relative contribution of populations to a core set of populations in which the overlap of genetic diversity was minimised. This implies that the mean kinship in the core set should be minimal. The above definition of diversity differs from Weitzman diversity in that it attempts to conserve the founder population (and thus minimises the loss of alleles), whereas Weitzman diversity favours the conservation of many inbred lines. The former is preferred in species where inbred lines suffer from inbreeding depression. The application of the method is illustrated by an example involving 45 Dutch poultry breeds. The calculations used were easy to implement and not computer intensive. The method gave a ranking of breeds according to their contributions to genetic diversity. Losses in genetic diversity ranged from 2.1% to 4.5% for different subsets relative to the entire set of breeds, while the loss of founder genome equivalents ranged from 22.9% to 39.3%.     

Contributions of Portuguese cattle breeds to genetic diversity using marker-estimated kinships

The quantitative assessment of genetic diversity within and between populations is important for decision-making in genetic conservation plans. In our study, we applied the livestock core set method to define the contribution of 15 cattle breeds, 11 of which are Portuguese indigenous cattle breeds, to genetic diversity. In livestock core set theory genetic diversity is defined as the maximum genetic variance that can be obtained in a random-mating population that is bred from the populations present in that core set. Two methods to estimate marker-estimated kinships to obtain the contributions to the core set were used in this study: the weighted log-linear model (WLM) and the weighted log-linear mixed model (WLMM). The breeds that contributed most to diversity in the core set were Holstein-Friesian followed by the Portuguese Mertolenga and Cachena for both WLM and WLMM methods. The ranking of relative contributions of cattle breeds was maintained when we considered only the Portuguese cattle breeds. Furthermore, we were able to identify the marginal contributions and respective losses of diversity for each of the 11 Portuguese cattle breeds when we considered a subset of populations that are not threatened of being lost (the Safe set composed of the four exotic breeds present in this study). When WLM was used losses in genetic diversity ranged from 2.68 to 0.65% while the loss in founder genome equivalents ranged from 37.37 to 8.43% for Mertolenga and Brava de Lide breeds respectively. When WLMM was used losses in genetic diversity and founder genome equivalents were less extreme than for the WLM method, ranging from 1.27 to 0.69 and 26.8 to 12.99 respectively.     

Estimation of breed contributions to present and future genetic diversity of 44 North Eurasian cattle breeds using core set diversity measures

Extinction of breeds threatens genetic diversity of livestock species. The need to conserve genetic diversity is widely accepted but involves in general two questions: (i) is the expected loss of diversity in a set of breeds within a defined future time horizon large enough to establish a conservation plan, and if so (ii) which breeds should be prioritised for such a conservation plan? The present study uses a marker assisted methodology to address these questions. The methodology combines core set diversity measures with a stochastic method for the estimation of expected future diversity and breed marginal diversities. The latter is defined as the change in the total diversity of all breeds caused by a one unit decrease in extinction probability of a particular breed. The stochastic method was validated by means of simulations. A large field data set consisting of 44 North Eurasian cattle breeds was analysed using simplified determined extinction probabilities. The results show that the expected loss of diversity in this set within the next 20 to 50 years is between 1 and 3% of the actual diversity, provided that the extinction probabilities which were used are approximately valid. If this loss is to be reduced, it is sufficient to include those three to five breeds with the highest marginal diversity in a conservation scheme.     

Conservation priorities for Ethiopian sheep breeds combining threat status,breed merits and contributions to genetic diversity

Prioritizing livestock breeds for conservation needs to incorporate both genetic and non-genetic aspects important for the survival of the breeds. Here, we apply a maximum-utility-strategy to prioritize 14 traditional Ethiopian sheep breeds based on their threat status, contributions to farmer livelihoods (current breed merits) and contributions to genetic diversity. Contributions of the breeds to genetic diversity were quantified using Eding''s marker-estimated kinship approaches. Non-genetic aspects included threats (e.g. low population size, low preferences by farmers) and current merits (economic, ecological and cultural merits). Threat analysis identified eight of the 14 breeds as threatened. Analysis of current merits showed that sub-alpine and arid-lowland breeds contribute most to farmer livelihoods in comparison to other breeds. The highest contribution to the genetic diversity conserved was from the Simien breed. Simien showed high between-breed (low between-breed kinship = 0.04) as well as high within-breed diversity (low within-breed kinship = 0.09 and high H_E = 0.73 and allelic richness = 6.83). We combined the results on threat status, current breed merits and contributions to genetic diversity to produce a ranking of the 14 breeds for conservation purposes. Our results balance the trade-offs between conserving breeds as insurance against future uncertainties and current sustainable utilization. The ranking of breeds provides a basis for conservation strategies for Ethiopian sheep and contributes to a regional or global conservation plan.     

The role of common ancestry and gene flow in the evolution of human‐directed play behaviour in dogs

Among‐population variance of phenotypic traits is of high relevance for understanding evolutionary mechanisms that operate in relatively short timescales, but various sources of nonindependence, such as common ancestry and gene flow, can hamper the interpretations. In this comparative analysis of 138 dog breeds, we demonstrate how such confounders can independently shape the evolution of a behavioural trait (human‐directed play behaviour from the Dog Mentality Assessment project). We combined information on genetic relatedness and haplotype sharing to reflect common ancestry and gene flow, respectively, and entered these into a phylogenetic mixed model to partition the among‐breed variance of human‐directed play behaviour while also accounting for within‐breed variance. We found that 75% of the among‐breed variance was explained by overall genetic relatedness among breeds, whereas 15% could be attributed to haplotype sharing that arises from gene flow. Therefore, most of the differences in human‐directed play behaviour among breeds have likely been caused by constraints of common ancestry as a likely consequence of past selection regimes. On the other hand, gene flow caused by crosses among breeds has played a minor, but not negligible role. Our study serves as an example of an analytical approach that can be applied to comparative situations where the effects of shared origin and gene flow require quantification and appropriate statistical control in a within‐species/among‐population framework. Altogether, our results suggest that the evolutionary history of dog breeds has left remarkable signatures on the among‐breed variation of a behavioural phenotype.     

Conservation priorities of genetic diversity in domesticated metapopulations: a study in taurine cattle breeds

We estimated neutral diversity of 21 European cattle breeds with 105 microsatellites. Nine of them resembled unselected Balkan Buša strains with diffuse breeding barriers and the 12 others were strongly differentiated, isolated breeds. Because of the impact of neutral genetic diversity on long-term population adaptive capacity, we discuss the long-term outcome of different conservation priorities in a subdivided metapopulation of the investigated cattle breeds. The optimal contribution to a pool of total genetic diversity allocated more than 95% of long-term relevant neutral diversity to virtually unselected strains of the Balkan Buša, while the maximization of total variance preferred inbred breeds. Current artificial selection methods, such as genomic selection sped up and a recovery of underestimated traits becomes quickly impossible. We emphasize that currently neutral and even deleterious alleles might be required for future genotypes in sustainable and efficient livestock breeding and production systems of a 21st century. We provide cumulative evidences that long-term survival relies on genetic complexity and complexity relies on allelic diversity. Our results suggest that virtually unselected, nonuniform strains harbor a crucial proportion of neutral diversity and should be conserved with high global priority. As one example, we suggest a cooperative maintenance of the nondifferentiated, highly fragmented, and fast vanishing metapopulation of Balkan Buša.     

Marker-assisted conservation of European cattle breeds: An evaluation

Two methods have been developed for the assessment of conservation priorities on the basis of molecular markers. According to the Weitzman approach, contributions to genetic diversity are derived from genetic distances between populations. Alternatively, diversity within and across populations is optimized by minimizing marker-estimated kinships. We have applied, for the first time, both methods to a comprehensive data set of 69 European cattle breeds, including all cosmopolitan breeds and several local breeds, for which genotypes of 30 microsatellite markers in 25–50 animals per breed have been obtained. Both methods were used to calculate the gain in diversity if a breed was added to a set of nine non-endangered breeds. Weitzman-derived diversities were confounded by genetic drift in isolated populations, which dominates the genetic distances but does not necessarily increase the conservation value of a breed. Marker-estimated kinships across populations were less disturbed by genetic drift than the Weitzman diversities and assigned high conservation values to Mediterranean breeds, which indeed have genetic histories that differ from the non-endangered breeds. Prospects and limitations of marker-assisted decisions on conservation priorities are discussed.     

Relative breed contributions to neutral genetic diversity of a comprehensive representation of Iberian native cattle

This study is aimed at establishing priorities for the optimal conservation of genetic diversity among a comprehensive group of 40 cattle breeds from the Iberian Peninsula. Different sets of breed contributions to diversity were obtained with several methods that differ in the relative weight attributed to the within- and between-breed components of the genetic variation. The contributions to the Weitzman diversity and the expected heterozygosity (He) account for between- and within-breed variation only, respectively. Contributions to the core set obtained for several kinship matrices, incorporate both sources of variation, as well as the combined contributions of Ollivier and Foulley and those of Caballero and Toro. In general, breeds that ranked high in the different core set applications also ranked high in the contribution to the global He, for example, Sayaguesa, Retinta, Monchina, Berrenda en Colorado or Marismeña. As expected, the Weitzman method prioritised breeds with low contributions to the He, like Mallorquina, Menorquina, Berrenda en Negro, Mostrenca, Vaca Palmera or Mirandesa, all showing highly negative contributions to He – that is, their removal would significantly increase the average He. Weighing the within- and between-breed components with the F_ST produced a balanced set of contributions in which all the breeds ranking high in both approaches show up. Unlike the other methods, the contributions to the diversity proposed by Caballero and Toro prioritised a good number of Portuguese breeds (Arouquesa, Barrosã, Mertolenga and Preta ranking highest), but this might be caused by a sample size effect. Only Sayaguesa ranked high in all the methods tested. Considerations with regard to the conservation scheme should be made before adopting any of these approaches: in situv. cryoconservation, selection and adaptation within the breeds v. crossbreeding or the creation of synthetic breeds. There is no general consensus with regard to balancing within- and between-breed diversity and the decision of which source to favour will depend on the particular scenario. In addition to the genetic information, other factors, such as geographical, historical, economic, cultural, etc., also need to be considered in the formulation of a conservation plan. All these aspects will ultimately influence the distribution of resources by the decision-makers.     

Exploration of conservation and development strategies with a limited stakeholder approach for local cattle breeds

Many local breeds have become endangered due to their substitution by high-yielding breeds. To conserve local breeds, effective development strategies need to be investigated. The aim of this study was to explore conservation and development strategies based on quantified strengths, weaknesses, opportunities and threats (SWOT) for two local cattle breeds from Northern Germany, namely the German Angler (GA) and Red Dual-Purpose cattle (RDP). The data comprised 158 questionnaires regarding both breeds’ SWOT, which were answered by 78 farmers of GA and 80 farmers of RDP. First, data were analysed using the SWOT-Analytic Hierarchy Process (AHP) method, which combines the qualitative strategic decision tool of SWOT analysis and the quantitative tool of AHP. Second, prioritised SWOT factors were discussed with stakeholders in order to form final conservation and development strategies at breed level. For GA prioritised strengths were daily gain, meat quality, milk production and the usage of new biotechnologies, weaknesses were genetic gain in milk production and inbreeding, opportunities were organic farming and breed-specific characteristics and threats were milk prices and dependency regarding the dairy business. Consequently, three conservation and development strategies were formed: (1) changing relative weights and the relevant breeding goal to drift from milk to meat, (2) increasing genetic gain and control the rate of inbreeding by the implementation of specific selection programs and (3) selection of unique and breed characteristic components on product level, that is, milk-fat and fine muscle fibers. For RDP defined strengths were robustness, high adaptability for different housing systems and a balanced dual-purpose of milk and meat, weaknesses were inbreeding, breed extinction, genomic selection with young bulls and milk yield, opportunities were organic farming and dual-purpose aspects and threats were milk and decreasing beef cattle prices. Thus, three conservation and development strategies were identified: (1) adjust relative weights and the relevant breeding goal to balance milk and meat yield, (2) increasing genetic gain and avoid extinction by implementing targeted selection programs and (3) selection of unique and breed characteristic traits on breed level, that is, environmental robustness. Quantified SWOT establish a basis for the exploration of conservation and development strategies at breed level. Explored strategies are promising even if the stakeholder approach was limited for small populations regarding a small number of stakeholder groups. The used approach reflects farmers’ individual convenience better than existing quantitative strategy decision tools on their own.     

猪a1-岩藻糖转移酶基因(FUT1)在26个中外猪种中的遗传变异研究

肠毒素大肠杆菌F18(ECF18)是引起仔猪断奶后水肿和腹泻病的主要病原菌,a1—岩藻糖转移酶基因(FUT1)是ECF18侵染猪小肠的受体蛋白候选基因。通过采用PCR—RFLP方法检测了5个西方商业猪种以及21个中国地方猪种(群)1458个个体在FUT1基因开放阅读框架的307核苷酸位点的G-A点突变(M307^G-A)遗传变异。结果表明：5个外来猪种以及中国地方猪种中的临高猪在该FUT1基因位点存在多态性,其他中国地方猪种均表现为极端的单态分布,只有易感的GG基因型,没有多态性。由此提示：1)如果猪FUT1 M307^G-A点突变是决定猪小肠上ECF18受体表达与否的关键因素,则绝大部分中国地方猪种均不具备抵抗ECF18的遗传基础,这除了表明ECF18抗性基因有可能起源于西方猪种外,同时也表明对中国地方猪种中在这个位点惟一存在多态性的海南临高猪的品种资源保存具有非常重要的意义。2)一般而言,在中国的养猪生产实践中,中国地方猪种的仔猪抗水肿与腹泻病能力普遍强于外来猪种,研究的结果提示有必要对中国地方猪种所具备的上述遗传抗性做更深入的研究,寻找、定位其相应的QTL或／和抗性基因。     

Using molecular markers and multivariate methods to study the genetic diversity of local European and Asian chicken breeds

French and Asian subsets of chicken breeds were first analysed using 22 microsatellites and then compared to the AVIANDIV European set using 14 loci. Positive correlations were observed between F _IT or F _ST and typological values or variance of markers using the multivariate analysis mcoa . The first axis of the multivariate representation separated Asian from European breeds, revealing breeds with Asian ancestor. Using all or 14 loci, correct assignation rate was always higher than 93%. The Weitzman index and the aggregate diversity D were calculated using 22 loci within French and Asian breeds. The French breed Coucou de Rennes and the Hua-Tung breed seemed to contribute the most to the global diversity of each subset. This approach on French-only breeds and then on French with AVIANDIV domestic breeds (14 loci) showed that the Marans breed contributed the most. The AVIANDIV framework could be useful to evaluate the genetic diversity of local breeds and to help in connecting national and regional conservation policies.     

Comparison of molecular breeding values based on within- and across-breed training in beef cattle

Background

Although the efficacy of genomic predictors based on within-breed training looks promising, it is necessary to develop and evaluate across-breed predictors for the technology to be fully applied in the beef industry. The efficacies of genomic predictors trained in one breed and utilized to predict genetic merit in differing breeds based on simulation studies have been reported, as have the efficacies of predictors trained using data from multiple breeds to predict the genetic merit of purebreds. However, comparable studies using beef cattle field data have not been reported.

Methods

Molecular breeding values for weaning and yearling weight were derived and evaluated using a database containing BovineSNP50 genotypes for 7294 animals from 13 breeds in the training set and 2277 animals from seven breeds (Angus, Red Angus, Hereford, Charolais, Gelbvieh, Limousin, and Simmental) in the evaluation set. Six single-breed and four across-breed genomic predictors were trained using pooled data from purebred animals. Molecular breeding values were evaluated using field data, including genotypes for 2227 animals and phenotypic records of animals born in 2008 or later. Accuracies of molecular breeding values were estimated based on the genetic correlation between the molecular breeding value and trait phenotype.

Results

With one exception, the estimated genetic correlations of within-breed molecular breeding values with trait phenotype were greater than 0.28 when evaluated in the breed used for training. Most estimated genetic correlations for the across-breed trained molecular breeding values were moderate (> 0.30). When molecular breeding values were evaluated in breeds that were not in the training set, estimated genetic correlations clustered around zero.

Conclusions

Even for closely related breeds, within- or across-breed trained molecular breeding values have limited prediction accuracy for breeds that were not in the training set. For breeds in the training set, across- and within-breed trained molecular breeding values had similar accuracies. The benefit of adding data from other breeds to a within-breed training population is the ability to produce molecular breeding values that are more robust across breeds and these can be utilized until enough training data has been accumulated to allow for a within-breed training set.     

Single-nucleotide-polymorphism-based association mapping of dog stereotypes

Phenotypic stereotypes are traits, often polygenic, that have been stringently selected to conform to specific criteria. In dogs, Canis familiaris, stereotypes result from breed standards set for conformation, performance (behaviors), etc. As a consequence, phenotypic values measured on a few individuals are representative of the breed stereotype. We used DNA samples isolated from 148 dog breeds to associate SNP markers with breed stereotypes. Using size as a trait to test the method, we identified six significant quantitative trait loci (QTL) on five chromosomes that include candidate genes appropriate to regulation of size (e.g., IGF1, IGF2BP2 SMAD2, etc.). Analysis of other morphological stereotypes, also under extreme selection, identified many additional significant loci. Less well-documented data for behavioral stereotypes tentatively identified loci for herding, pointing, boldness, and trainability. Four significant loci were identified for longevity, a breed characteristic not under direct selection, but inversely correlated with breed size. The strengths and limitations of the approach are discussed as well as its potential to identify loci regulating the within-breed incidence of specific polygenic diseases.     

Estimating the genetic architecture of quantitative traits

Understanding and estimating the structure and parameters associated with the genetic architecture of quantitative traits is a major research focus in quantitative genetics. With the availability of a well-saturated genetic map of molecular markers, it is possible to identify a major part of the structure of the genetic architecture of quantitative traits and to estimate the associated parameters. Multiple interval mapping, which was recently proposed for simultaneously mapping multiple quantitative trait loci (QTL), is well suited to the identification and estimation of the genetic architecture parameters, including the number, genomic positions, effects and interactions of significant QTL and their contribution to the genetic variance. With multiple traits and multiple environments involved in a QTL mapping experiment, pleiotropic effects and QTL by environment interactions can also be estimated. We review the method and discuss issues associated with multiple interval mapping, such as likelihood analysis, model selection, stopping rules and parameter estimation. The potential power and advantages of the method for mapping multiple QTL and estimating the genetic architecture are discussed. We also point out potential problems and difficulties in resolving the details of the genetic architecture as well as other areas that require further investigation. One application of the analysis is to improve genome-wide marker-assisted selection, particularly when the information about epistasis is used for selection with mating.     

Marker-Assisted Introgression in Backcross Breeding Programs

The efficiency of marker-assisted introgression in backcross populations derived from inbred lines was investigated by simulation. Background genotypes were simulated assuming that a genetic model of many genes of small effects in coupling phase explains the observed breed difference and variance in backcross populations. Markers were efficient in introgression backcross programs for simultaneously introgressing an allele and selecting for the desired genomic background. Using a marker spacing of 10-20 cM gave an advantage of one to two backcross generations selection relative to random or phenotypic selection. When the position of the gene to be introgressed is uncertain, for example because its position was estimated from a trait gene mapping experiment, a chromosome segment should be introgressed that is likely to include the allele of interest. Even for relatively precisely mapped quantitative trait loci, flanking markers or marker haplotypes should cover ~10-20 cM around the estimated position of the gene, to ensure that the allele frequency does not decline in later backcross generations.     

Genetic diversity and relationships of Portuguese and other horse breeds based on protein and microsatellite loci variation

There are three native Portuguese horse breeds: Lusitano, Sorraia and Garrano. This study compares diversity patterns of 17 protein and 12 microsatellite markers in these three as well as 30 other breeds to infer relationships among the breeds and to compare levels of polymorphism of these breeds for use in conservation efforts. The Garrano and the Lusitano showed a high level of genetic diversity, similar to that observed for most of the other analysed breeds, while the Sorraia and Friesian breeds showed low levels of variation for both genetic marker types. The combined protein and microsatellite data produced a tree that fit historical records well and with greater confidence levels than those for either data set alone. The combined genetic diversity and relationship information provides important baseline data for future breed conservation efforts, especially for a critically endangered breed such as the Sorraia.     

Exploring alternative models for sex-linked quantitative trait loci in outbred populations: application to an iberian x landrace pig intercross

We present a very flexible method that allows us to analyze X-linked quantitative trait loci (QTL) in crosses between outbred lines. The dosage compensation phenomenon is modeled explicitly in an identity-by-descent approach. A variety of models can be fitted, ranging from considering alternative fixed alleles within the founder breeds to a model where the only genetic variation is within breeds, as well as mixed models. Different genetic variances within each founder breed can be estimated. We illustrate the method with data from an F(2) cross between Iberian x Landrace pigs for intramuscular fat content and meat color component a*. The Iberian allele exhibited a strong overdominant effect for intramuscular fat in females. There was also limited evidence of one or more regions affecting color component a*. The analysis suggested that the QTL alleles were fixed in the Iberian founders, whereas there was some evidence of segregation in Landrace for the QTL affecting a* color component.     

Accuracies of genomically estimated breeding values from pure-breed and across-breed predictions in Australian beef cattle

Background

The major obstacles for the implementation of genomic selection in Australian beef cattle are the variety of breeds and in general, small numbers of genotyped and phenotyped individuals per breed. The Australian Beef Cooperative Research Center (Beef CRC) investigated these issues by deriving genomic prediction equations (PE) from a training set of animals that covers a range of breeds and crosses including Angus, Murray Grey, Shorthorn, Hereford, Brahman, Belmont Red, Santa Gertrudis and Tropical Composite. This paper presents accuracies of genomically estimated breeding values (GEBV) that were calculated from these PE in the commercial pure-breed beef cattle seed stock sector.

Methods

PE derived by the Beef CRC from multi-breed and pure-breed training populations were applied to genotyped Angus, Limousin and Brahman sires and young animals, but with no pure-breed Limousin in the training population. The accuracy of the resulting GEBV was assessed by their genetic correlation to their phenotypic target trait in a bi-variate REML approach that models GEBV as trait observations.

Results

Accuracies of most GEBV for Angus and Brahman were between 0.1 and 0.4, with accuracies for abattoir carcass traits generally greater than for live animal body composition traits and reproduction traits. Estimated accuracies greater than 0.5 were only observed for Brahman abattoir carcass traits and for Angus carcass rib fat. Averaged across traits within breeds, accuracies of GEBV were highest when PE from the pooled across-breed training population were used. However, for the Angus and Brahman breeds the difference in accuracy from using pure-breed PE was small. For the Limousin breed no reasonable results could be achieved for any trait.

Conclusion

Although accuracies were generally low compared to published accuracies estimated within breeds, they are in line with those derived in other multi-breed populations. Thus PE developed by the Beef CRC can contribute to the implementation of genomic selection in Australian beef cattle breeding.     

Genetic Variance for Body Size in a Natural Population of Drosophila Buzzatii

Previous work has shown thorax length to be under directional selection in the Drosophila buzzatii population of Carboneras. In order to predict the genetic consequences of natural selection, genetic variation for this trait was investigated in two ways. First, narrow sense heritability was estimated in the laboratory F2 generation of a sample of wild flies by means of the offspring-parent regression. A relatively high value, 0.59, was obtained. Because the phenotypic variance of wild flies was 7-9 times that of the flies raised in the laboratory, "natural" heritability may be estimated as one-seventh to one-ninth that value. Second, the contribution of the second and fourth chromosomes, which are polymorphic for paracentric inversions, to the genetic variance of thorax length was estimated in the field and in the laboratory. This was done with the assistance of a simple genetic model which shows that the variance among chromosome arrangements and the variance among karyotypes provide minimum estimates of the chromosome's contribution to the additive and genetic variances of the trait, respectively. In males raised under optimal conditions in the laboratory, the variance among second-chromosome karyotypes accounted for 11.43% of the total phenotypic variance and most of this variance was additive; by contrast, the contribution of the fourth chromosome was nonsignificant. The variance among second-chromosome karyotypes accounted for 1.56-1.78% of the total phenotypic variance in wild males and was nonsignificant in wild females. The variance among fourth chromosome karyotypes accounted for 0.14-3.48% of the total phenotypic variance in wild flies. At both chromosomes, the proportion of additive variance was higher in mating flies than in nonmating flies.     

Genetic diversity measures of local European beef cattle breeds for conservation purposes

This study was undertaken to determine the genetic structure, evolutionary relationships, and the genetic diversity among 18 local cattle breeds from Spain, Portugal, and France using 16 microsatellites. Heterozygosities, estimates of Fst, genetic distances, multivariate and diversity analyses, and assignment tests were performed. Heterozygosities ranged from 0.54 in the Pirenaica breed to 0.72 in the Barrosã breed. Seven percent of the total genetic variability can be attributed to differences among breeds (mean F_st = 0.07; P < 0.01). Five different genetic distances were computed and compared with no correlation found to be significantly different from 0 between distances based on the effective size of the population and those which use the size of the alleles. The Weitzman recursive approach and a multivariate analysis were used to measure the contribution of the breeds diversity. The Weitzman approach suggests that the most important breeds to be preserved are those grouped into two clusters: the cluster formed by the Mirandesa and Alistana breeds and that of the Sayaguesa and Tudanca breeds. The hypothetical extinction of one of those clusters represents a 17% loss of diversity. A correspondence analysis not only distinguished four breed groups but also confirmed results of previous studies classifying the important breeds contributing to diversity. In addition, the variation between breeds was sufficiently high so as to allow individuals to be assigned to their breed of origin with a probability of 99% for simulated samples.