Association of inbreeding and regional equine leucocyte antigen homozygosity with the prevalence of insect bite hypersensitivity in Old Kladruber horse.

Inbreeding depression is the reduction of performance caused by mating of close relatives. In livestock populations, inbreeding depression has been traditionally estimated by regression of phenotypes on pedigree inbreeding coefficients. This estimation can be improved by utilising genomic inbreeding coefficients. Here we estimate inbreeding depression for insect bite hypersensitivity (IBH) prevalence, the most common allergic horse disease worldwide, in Old Kladruber horse. In a deep pedigree with 3214 horses (187 genotyped), we used a generalised linear mixed model with IBH phenotype from 558 horses examined between 1996 and 2009 (1368 records). In addition to the classical pedigree information, we used the single-step approach that enabled joint use of pedigree and genomic information to estimate inbreeding depression overall genome and equine leucocyte antigen (ELA) class II region. Significant inbreeding depression was observed in all models fitting overall inbreeding coefficients (odds ratio between 1.018 and 1.074, P < 0.05) with the exception of Kalinowski's new inbreeding (P = 0.0516). The increase of ELA class II inbreeding was significantly associated with increased prevalence of IBH (odds ratio 1.018; P = 0.027). However, when fitted jointly with the overall inbreeding coefficient, the effect of ELA class II inbreeding was not significant (odds ratio 1.016; P = 0.062). Overall, the higher ELA class II and/or overall inbreeding (pedigree or genomic) was associated with increased prevalence of IBH in Old Kladruber horses. The single-step approach provides an efficient use of all the available pedigree, genomic, and phenotype information for estimation of overall and regional inbreeding effects.

Genome-wide mapping of the dominance effects based on breed ancestry for semen traits in admixed Swiss Fleckvieh bulls.

Heterosis is the beneficial deviation of crossbred progeny from the average of parental lines for a particular trait. Heterosis is due to nonadditive genetic effects with dominance and epistatic components. Recent advances in genotyping technology have encouraged researchers to estimate and scan heterosis components for a range of traits in crossbred populations, applying various definitions of such components. In this study, we defined the intralocus (dominance) component of heterosis using local genetic ancestry and performed genome-wide association analysis for admixed Swiss Fleckvieh bulls and their parental populations, Red Holstein Friesian and Swiss Simmental, for semen traits. A linear mixed model for 41,824 SNP, including SNP additive genetic, breed additive, and breed dominance effects on 1,178 bulls (148 Red Holstein Friesian, 213 Swiss Simmental, and 817 Swiss Fleckvieh) with a total of 43,782 measurements was performed. In total, 19 significant regions for breed dominance were identified for volume (2 regions on Bos taurus autosome 10 and 22) and percentage of live spermatozoa (17 regions on Bos taurus autosome 3, 4, 5, 7, 13, 14, and 17), and genes associated with spermatogenesis, sperm motility, and male fertility traits were located there. No significant region for breed dominance was detected for total number of spermatozoa. The signals for breed dominance were relatively wide, most likely due to limited numbers of recombination events in a small number of generations (10-15 generations) of crossbreeding in the recent Swiss Fleckvieh composite.

Effects of breed proportion and components of heterosis for semen traits in a composite cattle breed.

The aim of this study was to estimate the non-additive genetic effects of the dominance component of heterosis as well as epistatic loss on semen traits in admixed Swiss Fleckvieh, a composite of Simmental (SI) and Red Holstein Friesian (RHF) cattle. Heterosis is the additional gain in productivity or fitness of cross-bred progeny over the mid-purebred parental populations. Intralocus gene interaction usually has a positive effect, while epistatic loss generally reduces productivity or fitness due to lack of evolutionarily established interactions of genes from different breeds. Genotypic data on 38,205 SNP of 818 admixed, as well as 148 RHF and 213 SI bulls as the parental breeds were used to predict breed origin of alleles. The genomewide locus-specific breed ancestries of individuals were used to calculate effects of breed difference as well as the dominance component of heterosis, while proxies for two definitions of epistatic loss were derived from 100,000 random pairs of loci. The average Holstein Friesian ancestry in admixed bulls was estimated 0.82. Results of fitting different linear mixed models showed including the dominance component of heterosis considerably improved the model adequacy for three of the four traits. Inclusion of epistatic loss increased the accuracy of the models only for our new definition of the epistatic effect for two traits, while the other definition was so highly correlated with the dominance component that statistical separation was impossible.

Locus-specific ancestry to detect recent response to selection in admixed Swiss Fleckvieh cattle.

Identification of selection signatures is one of the current endeavors of evolutionary genetics. Admixed populations may be used to infer post-admixture selection. We calculated local ancestry for Swiss Fleckvieh, a composite of Simmental (SI) and Red Holstein Friesian (RHF), to infer such signals. Illumina Bovine SNP50 BeadChip data for 300 admixed, 88 SI and 97 RHF bulls were used. The average RHF ancestry across the whole genome was 0.70. To identify regions with high deviation from average, we considered two significance thresholds, based on a permutation test and extreme deviation from normal distribution. Regions on chromosomes 13 (46.3-47.3 Mb) and 18 (18.7-25.9 Mb) passed both thresholds in the direction of increased SI. Extended haplotype homozygosity within (iHS) and between (Rsb) populations was calculated to explore additional patterns of pre- and post-admixture selection signals. The Rsb score of admixed and SI was significant in a wide region of chromosome 18 (6.6-24.6 Mb) overlapped with one area of strong local ancestry deviation. FTO, with pleiotropic effect on milk and fertility, NOD2 on dairy and NKD1 and SALL1 on fertility traits are located there. Genetic differentiation of RHF and SI (Fst ), an alternative indicator of pre-admixture selection in pure populations, was calculated. No considerable overlap of peaks of local ancestry deviations and Fst was observed. We found two regions with significant signatures of post-admixture selection in this very young composite, applying comparatively stringent significance thresholds. The signals cover relatively large genomic areas and did not allow pinpointing of the gene(s) responsible for the apparent shift in ancestry proportions.

Genetic relationships among American donkey populations: insights into the process of colonization.

This study presents the first insights into the genetic diversity and structure of the American donkey metapopulation. The primary objectives were to detect the main structural features underlying variability among American donkey populations, identify boundaries between differentiated gene pools, and draw the main colonization pathways since the introduction of donkeys into America in the 15th century. A panel of 14 microsatellite markers was applied for genotyping 350 American donkeys from 13 countries. The genetic structure of this metapopulation was analysed using descriptive statistics and Bayesian model-based methods. These populations were then compared to a database containing information on 476 individuals from 11 European breeds to identify the most likely ancestral donor populations. Results showed the presence of two distinct genetic pools, with confluence of the two in Colombia. The southern pool showed a unique genetic signature subsequent to an older founder event, but lacked any significant influence of modern gene flow from Europe. The northern pool, conversely, may have retained more ancestral polymorphisms and/or have experienced modern gene flow from Spanish breeds. The Andalusian and, to a lesser extent, the Catalan breeds have left a more pronounced footprint in some of the American donkey populations analysed.

GRAIN: a computer program to calculate ancestral and partial inbreeding coefficients using a gene dropping approach.

GRain is freely available software intended to enable and promote testing of hypotheses with respect to purging and heterogeneity of inbreeding depression. The program is based on a stochastic approach, the gene dropping method, and calculates various coefficients from large and complex pedigrees. GRain calculates, together with the 'classical' inbreeding coefficient, ancestral inbreeding coefficients proposed by Ballou, (1997) J. Hered., 88, 169 and Kalinowski et al., (2000) Conserv. Biol., 14, 1375 as well as an ancestral history coefficient (AHC ), defined here for the first time. AHC is defined as the number that tells how many times during pedigree segregation (gene dropping) a randomly taken allele has been in IBD status. Furthermore, GRain enables testing of heterogeneity and/or purging of inbreeding depression with respect to different founders/ancestors by calculating partial coefficients for all previously obtained coefficients.

Origin of porcine circovirus type 2 (PCV2) from swine affected by PCV2-associated diseases in Croatia.

Porcine circovirus type 2 (PCV2) causes some of the most significant economic losses in pig production. Several multisystemic syndromes have been attributed to PCV2 infection, which are known as PCV2-associated diseases (PCVDs). This study investigated the origin and evolution of PCV2 sequences in domestic pigs and wild boars affected by PCVDs in Croatia. Viral sequences were recovered from three wild boars diagnosed with PCV2-systemic disease (PCV2-SD), 63 fetuses positive for PCV2 DNA as determined by PCR, 14 domestic pigs affected with PCV2-SD (displaying severe interstitial nephritis) and five domestic pigs with proliferative and necrotising pneumonia. Seventeen complete PCV2 genomes were recovered. Phylogenetic and evolutionary analyses based on median-joining phylogenetic networks, amino acid alignments and principal coordinate analysis were performed using complete genomes, as well as complete and partial ORF sequences for ORF1 and ORF2. Two of the 17 PCV2 sequences belonged to PCV2a, 14 to PCV2b and one was unclustered. PCV2b was the predominant genotype in Croatia and has been linked to international trade as a route of introduction. Correlation between particular viral strains with PCVDs is lacking.

Lack of mitochondrial DNA structure in Balkan donkey is consistent with a quick spread of the species after domestication.

A total of 132 mtDNA sequences from 10 Balkan donkey populations were analysed to ascertain their regional genetic structure and to contribute to the knowledge of the spreading of the species after domestication. The Balkan donkey sequences were compared with those from 40 Burkina Faso donkeys as an African outgroup to account for possible local Balkan scenarios. The 172 sequences gave 62 different haplotypes (55 in Balkan donkey). Virtually all the analysed populations had haplotypes assigned to either Clade 1 or Clade 2 even though the relative proportion of Clade 1 or 2 haplotypes differed across populations. Geographical maps constructed using factors computed via principal component analysis showed that the Balkan donkey populations are not spatially structured. AMOVA confirmed a lack of genetic structure in Balkan donkey mtDNA. Balkan populations were poorly differentiated (ΦST  = 0.071). Differentiation between the Balkan donkey and the African outgroup also was low. The lack of correspondence between geographical areas and maternal genetic structure is consistent with the hypothesis suggesting a very quick spread of the species after domestication. The current research illustrates the difficulties to trace routes of expansion in donkey, as the species has no geographical structure.

Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations.

Using genome-wide SNP data, we calculated genomic inbreeding coefficients (FROH  > 1  Mb , FROH  > 2 Mb , FROH  > 8 Mb and FROH  > 16 Mb ) derived from runs of homozygosity (ROH) of different lengths (>1, >2, >8 and > 16 Mb) as well as from levels of homozygosity (FHOM ). We compared these values of inbreeding coefficients with those calculated from pedigrees (FPED ) of 1422 bulls comprising Brown Swiss (304), Fleckvieh (502), Norwegian Red (499) and Tyrol Grey (117) cattle breeds. For all four breeds, population inbreeding levels estimated by the genomic inbreeding coefficients FROH  > 8 Mb and FROH  > 16 Mb were similar to the levels estimated from pedigrees. The lowest values were obtained for Fleckvieh (FPED  = 0.014, FROH  > 8 Mb  = 0.019 and FROH  > 16 Mb  = 0.008); the highest, for Brown Swiss (FPED  = 0.048, FROH  > 8 Mb  = 0.074 and FROH  > 16 Mb  = 0.037). In contrast, inbreeding estimates based on the genomic coefficients FROH  > 1 Mb and FROH  > 2 Mb were considerably higher than pedigree-derived estimates. Standard deviations of genomic inbreeding coefficients were, on average, 1.3-1.7-fold higher than those obtained from pedigrees. Pearson correlations between genomic and pedigree inbreeding coefficients ranged from 0.50 to 0.62 in Norwegian Red (lowest correlations) and from 0.64 to 0.72 in Tyrol Grey (highest correlations). We conclude that the proportion of the genome present in ROH provides a good indication of inbreeding levels and that analysis based on ROH length can indicate the relative amounts of autozygosity due to recent and remote ancestors.

The contribution of dominance and inbreeding depression in estimating variance components for litter size in Pannon White rabbits.

In a synthetic closed population of Pannon White rabbits, additive (VA ), dominance (VD ) and permanent environmental (VPe ) variance components as well as doe (bF d ) and litter (bF l ) inbreeding depression were estimated for the number of kits born alive (NBA), number of kits born dead (NBD) and total number of kits born (TNB). The data set consisted of 18,398 kindling records of 3883 does collected from 1992 to 2009. Six models were used to estimate dominance and inbreeding effects. The most complete model estimated VA and VD to contribute 5.5 ± 1.1% and 4.8 ± 2.4%, respectively, to total phenotypic variance (VP ) for NBA; the corresponding values for NBD were 1.9 ± 0.6% and 5.3 ± 2.4%, for TNB, 6.2 ± 1.0% and 8.1 ± 3.2% respectively. These results indicate the presence of considerable VD . Including dominance in the model generally reduced VA and VPe estimates, and had only a very small effect on inbreeding depression estimates. Including inbreeding covariates did not affect estimates of any variance component. A 10% increase in doe inbreeding significantly increased NBD (bF d  = 0.18 ± 0.07), while a 10% increase in litter inbreeding significantly reduced NBA (bF l  = -0.41 ± 0.11) and TNB (bF l  = -0.34 ± 0.10). These findings argue for including dominance effects in models of litter size traits in populations that exhibit significant dominance relationships.

Mitochondrial DNA and Y-chromosome diversity in East Adriatic sheep.

Variation in mitochondrial DNA (mtDNA) and Y-chromosome haplotypes was analysed in nine domestic sheep breeds (159 rams) and 21 mouflon (Ovis musimon) sampled in the East Adriatic. Mitochondrial DNA analyses revealed a high frequency of type B haplotypes, predominantly in European breeds, and a very low frequency of type A haplotypes, which are more frequent in some Asian breeds. Mitochondrial haplotype Hmt-3 was the most frequent (26.4%), and 37.1%, 20.8% and 7.6% of rams had haplotypes one, two and three mutations remote from Hmt-3 respectively. In contrast, Y-chromosome analyses revealed extraordinary paternal allelic richness: HY-6, 89.3%; HY-8, 5.0%; HY-18, 3.1%; HY-7, 1.3%; and HY-5, 1.3%. In fact, the number of haplotypes observed is comparable to the number found in Turkish breeds and greater than the number found in European breeds so far. Haplotype HY-18 (A-oY1/135-SRYM18), identified here for the first time, provides a link between the haplotype HY-12 (A-oY1/139-SRYM18) found in a few rams in Turkey and haplotype HY-9 (A-oY1/131-SRYM18) found in one ram in Ethiopia. All mouflons had type B mtDNA haplotypes, including the private haplotype (Hmt-55), and all were paternally monomorphic for haplotype HY-6. Our data support a quite homogeneous maternal origin of East Adriatic sheep, which is a characteristic of European breeds. At the same time, the high number of haplotypes found was surprising and intriguing, and it begs for further analysis. Simultaneous analysis of mtDNA and Y-chromosome information allowed us to detect a large discrepancy between maternal and paternal lineages in some populations. This is most likely the result of breeder efforts to 'upgrade' local populations using rams with different paternal origins.

Prediction of breed composition in an admixed cattle population.

Swiss Fleckvieh was established in 1970 as a composite of Simmental (SI) and Red Holstein Friesian (RHF) cattle. Breed composition is currently reported based on pedigree information. Information on a large number of molecular markers potentially provides more accurate information. For the analysis, we used Illumina BovineSNP50 Genotyping Beadchip data for 90 pure SI, 100 pure RHF and 305 admixed bulls. The scope of the study was to compare the performance of hidden Markov models, as implemented in structure software, with methods conventionally used in genomic selection [BayesB, partial least squares regression (PLSR), least absolute shrinkage and selection operator (LASSO) variable selection)] for predicting breed composition. We checked the performance of algorithms for a set of 40 492 single nucleotide polymorphisms (SNPs), subsets of evenly distributed SNPs and subsets with different allele frequencies in the pure populations, using F(ST) as an indicator. Key results are correlations of admixture levels estimated with the various algorithms with admixture based on pedigree information. For the full set, PLSR, BayesB and structure performed in a very similar manner (correlations of 0.97), whereas the correlation of LASSO and pedigree admixture was lower (0.93). With decreasing number of SNPs, correlations decreased substantially only for 5% or 1% of all SNPs. With SNPs chosen according to F(ST) , results were similar to results obtained with the full set. Only when using 96 and 48 SNPs with the highest F(ST) , correlations dropped to 0.92 and 0.90 respectively. Reducing the number of pure animals in training sets to 50, 20 and 10 each did not cause a drop in the correlation with pedigree admixture.

Hierarchical structure of the Sicilian goats revealed by Bayesian analyses of microsatellite information.

Genetic structure and relationship amongst the main goat populations in Sicily (Girgentana, Derivata di Siria, Maltese and Messinese) were analysed using information from 19 microsatellite markers genotyped on 173 individuals. A posterior Bayesian approach implemented in the program STRUCTURE revealed a hierarchical structure with two clusters at the first level (Girgentana vs. Messinese, Derivata di Siria and Maltese), explaining 4.8% of variation (amovaФ(ST) estimate). Seven clusters nested within these first two clusters (further differentiations of Girgentana, Derivata di Siria and Maltese), explaining 8.5% of variation (amovaФ(SC) estimate). The analyses and methods applied in this study indicate their power to detect subtle population structure.

Modeling variance structure of body shape traits of Lipizzan horses.

Heterogeneity of variance of growth traits over age is a common issue in estimating genetic parameters and is addressed in this study by selecting appropriate variance structure models for additive genetic and environmental variances. Modeling and partitioning those variances connected with analyzing small data sets were demonstrated on Lipizzan horses. The following traits were analyzed: withers height, chest girth, and cannon bone circumference. The measurements were taken at birth, and at approximately 6, 12, 24, and 36 mo of age of 660 Lipizzan horses born in Croatia between 1948 and 2000. The corresponding pedigree file consisted of 1,458 horses. Sex, age of dam, and stud-year-season interaction were considered fixed effects; additive genetic and permanent environment effects were defined as random. Linear adjustments of age at measuring were done within measuring groups. Maternal effects were included only for measurements taken at birth and at 6 mo. Additive genetic variance structures were modeled by using uniform structures or structures based on polynomial random regression. Environmental variance structures were modeled by using one of the following models: unstructured, exponential, Gaussian, or combinations of identity or diagonal with structures based on polynomial random regression. The parameters were estimated by using REML. Comparison and fits of the models were assessed by using Akaike and Bayesian information criteria, and by checking graphically the adequacy of the shape of the overall (phenotypic) and component (additive genetic and environmental) variance functions. The best overall fit was obtained from models with unstructured error variance. Compared with the model with uniform additive genetic variance, models with structures based on random regression only slightly improved overall fit. Exponential and Gaussian models were generally not suitable because they do not accommodate adequately heterogeneity of variance. Using the unstructured error variance model, the heritability estimates ranged from 0.17 to 0.33 for withers height, 0.07 to 0.27 for chest girth, and 0.14 to 0.30 for cannon bone circumference. This study demonstrated the necessity of accounting for heterogeneity of variances and covariances for body shape traits in Lipizzan horses, and possible difficulties in estimating variance and covariance components when applying more complicated structure models on a small data set. The choice of models depends not only on overall fit but also on the fit of genetic and environmental components.

Multiple paternal origins of domestic cattle revealed by Y-specific interspersed multilocus microsatellites.

In this study, we show how Y-specific interspersed multilocus microsatellites, which are loci that yield several amplified bands differing in size from the same male individual and PCR reaction, are a powerful source of information for tracing the history of cattle. Our results confirm the existence of three main groups of sires, which are separated by evolutionary time and clearly predate domestication. These three groups are consistent with the haplogroups previously identified by Götherström et al. (2005) using five Y-specific segregating sites: Y1 and Y2 in taurine (Bos taurus) cattle and Y3 in zebu (Bos indicus) cattle. The zebu cattle cluster clearly originates from a domestication process that was geographically and temporally separated from that of taurine clusters. Our analyses further suggest that: (i) introgression of wild sire genetic material into domesticated herds may have a significant role in the formation of modern cattle, including the formation of the Y1 haplogroup; (ii) a putative domestication event in Africa probably included local Y2-like wild sires; (iii) the West African zebu cattle Y-chromosome may have partially originated from an ancient introgression of humped cattle into Africa; and (iv) the high genetic similarity among Asian zebu sires is consistent with a single domestication process.

Y-specific microsatellites reveal an African subfamily in taurine (Bos taurus) cattle.

Five cattle Y-specific microsatellites, totalling six loci, were selected from a set of 44 markers and genotyped on 608 Bos taurus males belonging to 45 cattle populations from Europe and Africa. A total of 38 haplotypes were identified. Haplogroups (Y1 and Y2) previously defined using single nucleotide polymorphisms did not share haplotypes. Nine of the 27 Y2-haplotypes were only present in African cattle. Network and correspondence analyses showed that this African-specific subfamily clustered separately from the main Y2-subfamily and the Y1 haplotypes. Within-breed genetic variability was generally low, with most breeds (78%) showing haplotypes belonging to a single haplogroup. AMOVA analysis showed that partitioning of genetic variation among breeds can be mainly explained by their geographical and haplogroup assignment. Between-breed genetic variability summarized via Principal Component Analysis allowed the identification of three principal components explaining 94.2% of the available information. Projection of principal components on geographical maps illustrated that cattle populations located in mainland Europe, the three European Peninsulas and Mediterranean Africa presented similar genetic variation, whereas those breeds from Atlantic Europe and British Islands (mainly carrying Y1 haplotypes) and those from Sub-Saharan Africa (belonging to Y2-haplogroup) showed genetic variation of a different origin. Our study confirmed the existence of two large Y-chromosome lineages (Y1 and Y2) in taurine cattle. However, Y-specific microsatellites increased analytical resolution and allowed at least two different Y2-haplotypic subfamilies to be distinguished, one of them restricted to the African continent.

Genetic diversity and population structure of the synthetic Pannon White rabbit revealed by pedigree analyses.

Demographic history, current status, and efficiency of the mating strategy were analyzed using the pedigree of Pannon White (PW) rabbits born between 1992 and 2007. Potential accumulation of detrimental effects and loss of genetic diversity were also considered. Calculations and estimates were done most often for rabbits born in 2007, whereas other reference populations (REFPOPXXXX) were specified explicitly. The pedigree contained 4,749 individuals and 580 founders, and its completeness was 82.1% up to 10 and 94.5% up to 5 generations, respectively. Generation intervals through different pathways averaged 1.2 yr. When adjusted to the pedigree completeness, the amount of inbreeding (F(i)) of rabbits was comparable (5.54%) with that of other livestock populations, whereas the 10 (30) founders contributing the most to inbreeding explained a large part of the population inbreeding [i.e., 42.24% (73.18%)]. The ancestral inbreeding coefficient of REFPOP2004 (10.67%) was one-half that of REFPOP2007 (20.66%), showing its strong dependence on pedigree length. Family variance, inbreeding, and realized effective population size were 84.18 (REFPOP2006; this variable could not be calculated for the last year examined), 37.19, and 91.08, respectively. The effective numbers of ancestors, founders, and founder genomes were 48, 26, and 7.33, respectively. Although the circular mating scheme applied was generally effective, the large accumulated reduction in genetic variability indicates the need to revise and improve the current breeding strategy.

Individual-based assessment of population structure and admixture in Austrian, Croatian and German draught horses.

All over Europe, the number of draught horses has decreased drastically during the last 50 years. As a prerequisite for efficient management decisions, we analysed the conservation status in Austrian (Noriker Carinthia - NC, Noriker Salzburg - NS), Croatian (Croatian Coldblood - C, Posavina horse - P) and German (Altmaerkisch Coldblood - A, Black Forest horse - BF, Mecklenburg Coldblood - M, Rhenish German Draught horse - R, Saxon Thuringa Coldblood - ST, Schleswig Draught horse - Sch, South German Coldblood - SG) draught horses (434) using multilocus genotypic information from 30 (effectively 27) microsatellite loci. Populations located in areas with less intensive agricultural production (C, NC, NS, P and SG) had greater diversity within the population and estimated effective population size than A, BF, Sch, M, R and ST populations. The PCA plots revealed that populations form five separate groups. The 'Noriker' group (NC, NS and SG) and the 'Rhenish' group (A, M, R and ST) were the most distinctive (pairwise F(ST) values ranged from 0.078 to 0.094). The 'Croatian' group (C and P) was in the centre, while the BF and Sch populations formed two out-groups. A posterior Bayesian analysis detected further differentiation, mainly caused by political and geographical factors. Thus, it was possible to separate the South German Coldblood from the Austrian Noriker population where no subpopulation structure was detected. The admixture analysis revealed imprecise classification between C and P populations. A small but notable separation of R from A, M and ST populations was detected, while Sch and BF populations remained as out-groups. The information obtained should aid in making efficient conservation decisions.

Quantitative genetic aspects of coat color in horses.

The aim of this study was to estimate genetic parameters for coat color in horses. Besides defining coat color classes (gray, chestnut, bay, and black), the phenotypes were also measured quantitatively according to standardized international procedures (Commission Internationale de l'Eclairage L*, a*, b*), where L* describes lightness, a* describes color saturation from red to green, and b* describes color saturation from yellow to blue. The total color saturation was derived from a* and b* and referred to as Chroma. A total of 294 horses from the breeds Lipizzan, Nonius, Arabian Pure Bred, Shagya Arabian, and Gidran were measured at neck, shoulder, and belly. Heritabilities (within and between breeds or color classes) and repeatabilities were estimated using REML from univariate animal models defined separately for gray and nongray horses. For gray horses, the estimated within-breed heritabilities for L* ranged from 0.45 to 0.49 and for a*, b*, and Chroma from 0.09 to 0.52, indicating moderate polygenic effect. For nongray horses, between-color class heritabilities were high (0.70 to 0.85) and within-color class heritabilities were negligible (except for L* measured on neck and belly, 0.21 and 0.34, respectively). Additionally, the importance of L* was described by the relation with the total melanin content of horse coat hair; for gray and nongray horses, a strong negative linear relationship was detected (P < 0.01). The spectrometric measures and the results of this study demonstrate a possible approach to the estimation of the polygenic component involved in coat color inheritance.

Microsatellite diversity, population subdivision and gene flow in the Lipizzan horse.

Blood samples of 561 Lipizzan horses from subpopulations (studs) of seven European countries representing a large fraction of the breed's population were used to examine the genetic diversity, population subdivision and gene flow in the breed. DNA analysis based on 18 microsatellite loci revealed that genetic diversity (observed heterozygosity = 0.663, gene diversity = 0.675 and the mean number of alleles = 7.056) in the Lipizzan horse is similar to other horse breeds as well as to other domestic animal species. The genetic differentiation between Lipizzan horses from different studs, although moderate, was apparent (pairwise F(ST) coefficients ranged from 0.021 to 0.080). Complementary findings explaining the genetic relationship among studs were revealed by genetic distance and principal component analysis. One genetic cluster consisted of the subpopulations of Austria, Italy and Slovenia, which represent the classical pool of Lipizzan horse breeding. A second cluster was formed by the Croatian, Hungarian and Slovakian subpopulations. The Romanian subpopulation formed a separate unit. The largest genetic differentiation was found between the Romanian and Italian subpopulation. Genetic results are consistent with the known breeding history of the Lipizzan horse. Correct stud assignment was obtained for 80.9% and 92.1% of Lipizzan horses depending on the inclusion or exclusion of migrant horses, respectively. The results of the present study will be useful for the development of breeding strategies, which consider classical horse breeding as well as recent achievements of population and conservation genetics.