Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high‐density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (F_ROH), which were compared with inbreeding coefficients estimated from pedigree information (F_PED) and using the genomic relationship matrix (F_GRM). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The F_PED were moderately correlated with F_ROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with F_GRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of F_ROH with F_PED and F_GRM. ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium‐density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.     

Runs of homozygosity reveal genome‐wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium‐density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH = 0.099), whereas Comisana showed the lowest (F_ROH = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.     

Association between population structure and allele frequencies of the glycogen synthase 1 mutation in the Austrian Noriker draft horse

The aim of this study was to determine the allele frequency of the glycogen synthase 1 (GYS1) mutation associated with polysaccharide storage myopathy type 1 in the Austrian Noriker horse. Furthermore, we examined the influence of population substructures on the allele distribution. The study was based upon a comprehensive population sample (208 breeding stallions and 309 mares) and a complete cohort of unselected offspring from the year 2014 (1553 foals). The mean proportion of GYS1 carrier animals in the foal cohort was 33%, ranging from 15% to 50% according to population substructures based on coat colours. In 517 mature breeding horses the mutation carrier frequency reached 34%, ranging on a wider scale from 4% to 62% within genetic substructures. We could show that the occurrence of the mutated GYS1 allele is influenced by coat colour; genetic bottlenecks; and assortative, rotating and random mating strategies. Highest GYS1 carrier frequencies were observed in the chestnut sample comprising 50% in foals, 54% in mares and 62% in breeding stallions. The mean inbreeding of homozygous carrier animals reached 4.10%, whereas non‐carrier horses were characterized by an inbreeding coefficient of 3.48%. Lowest GYS1 carrier frequencies were observed in the leopard spotted Noriker subpopulation. Here the mean carrier frequency reached 15% in foals, 17% in mares and 4% in stallions and inbreeding decreased from 3.28% in homozygous non‐carrier horses to 2.70% in heterozygous horses and 0.94% in homozygous carriers. This study illustrates that lineage breeding and specified mating strategies result in genetic substructures, which affect the frequencies of the GYS1 gene mutation.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

Whole-genome association study for the roan coat color in an intercrossed pig population between Landrace and Korean native pig

The roan coat color is characterized by white hairs intermingled with colored hairs. Candidate genes based on comparative phenotypes in horses and cattle involve the KIT and KIT ligand (MGF) genes. Here, we report the result of the whole genome scanning to detect genomic regions responsible for the roan coat color, using a three-generation pedigree of 62 pigs in an intercross between Landrace and Korean native pig. These pigs were genotyped using the PorcineSNP 60 BeadChip (Illumina, USA). The whole genome scan indicated that three genomic regions, 35～36 Mb, 38～39 Mb, and 58～59 Mb on SSC8, were commonly and highly associated/linked with the roan phenotype in the case/control, sib-pair, and linkage test, respectively. The porcine KIT was selected as a candidate gene, because it is located in one of the three significant regions and its function is related to coat color formation. SNPs and Indels within coding sequence (CDS), promoter, and 3′-UTR of KIT were surveyed. Twenty-two SNPs in the CDS reported previously, as well as nine variations in promoter (2 SNPs) and 3′-UTR (5 SNPs and 2 Indels) were detected. Although no causative mutations were identified, these results will help to elucidate the genetic mechanisms involved in the expression of the roan phenotype and will aid in identifying key mutations responsible for the roan phenotype in further studies.     

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F_ROH) in three local dairy cattle breeds. F_ROH values were compared with F estimated from the genomic relationship matrix (F_GRM), based on the difference between observed v. expected number of homozygous genotypes (F_HOM) and the genomic homozygosity of individual i (F_{MOL i}). The molecular coancestry coefficient (f_{MOL ij}) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F_ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F_HOM and F_ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F_ROH and other estimated F coefficients were generally lower ranged from 0.45 (F_{MOL i}−F_ROH) in Cinisara to 0.17 (F_GRM−F_ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome‐wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.     

The inheritance of the leopard complex of spotting patterns in horses

The leopard complex of white spotting patterns in horses consists of the leopard, few-spot leopard, blanket, blanket with spots, varnish roan (or marble), snowflake, frosted, speckled, and mottled patterns. Horses with any of these patterns can produce the other patterns when mated to nonpatterned horses. Twenty-two horses of the Welsh Pony, Noriker, Appaloosa, and Pony of the Americas breeds produced 270 foals in a distribution consistent with a single dominant allele being responsible for the patterns. The symbol for this dominant allele, Lp, is retained from previous work on the leopard pattern. Heterozygotes are less extensively marked than are homozygotes, but the two classes overlap. The differences in the patterns are related to varying degrees of expression of this allele.     

SNP genotypes reveal breed substructure,selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population-specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.     

A missense mutation in the gene for melanocyte-stimulating hormone receptor (MCIR) is associated with the chestnut coat color in horses

The melanocyte-stimulating hormone receptor gene (MC1R) is the major candidate gene for the chestnut coat color in horses since it is assumed to be controlled by an allele at the extension locus. MC1R sequences were PCR amplified from chestnut (e/e) and non-chestnut (E/−) horses. A single-strand conformation polymorphism was found that showed a complete association to the chestnut coat color among 144 horses representing 12 breeds. Sequence analysis revealed a single missense mutation (83Ser → Phe) in the MC1R allele associated with the chestnut color. The substitution occurs in the second transmembrane region, which apparently plays a key role in the molecule since substitutions associated with coat color variants in mice and cattle as well as red hair and fair skin in humans are found in this part of the molecule. We propose that the now reported mutation is likely to be the causative mutation for the chestnut coat color. The polymorphism can be detected with a simple PCR-RFLP test, since the mutation creates a TaqI restriction site in the chestnut allele. Received: 20 May 1996 / Accepted: 31 July 1996     

Genomic regions influencing coat color saturation and facial markings in Fleckvieh cattle

Genomic regions associated with coat color and pigmented areas of the head were identified for Fleckvieh (dual‐purpose Simmental), a red‐spotted and white‐headed cattle breed. Coat color was measured with a chromameter, implementing the CIELAB color space and resulting in numerical representation of lightness, color intensity, red/green and blue/yellow color components, rather than subjective classification. Single marker regression analyses with fixed effects of the sex and barn were applied, and significant regions were determined with the local false discovery rate methodology. The PMEL and ERBB3 genes on chromosome 5 were in the most significant region for the color measurements. In addition to the blue/yellow color component and color intensity, the AP3B2 gene on chromosome 21 was identified. Its function was confirmed for similar traits in a range of model species. The KIT gene on chromosome 6 was found to be strongly associated with the inhibition of circum‐ocular pigmentation and pigmented spots on the cheek.     

Molecular genetics of coat colour variations in White Galloway and White Park cattle

White Galloway cattle exhibit three different white coat colour phenotypes, that is, well marked, strongly marked and mismarked. However, mating of individuals with the preferred well or strongly marked phenotype also results in offspring with the undesired mismarked and/or even fully black coat colour. To elucidate the genetic background of the coat colour variations in White Galloway cattle, we analysed four coat colour relevant genes: mast/stem cell growth factor receptor (KIT), KIT ligand (KITLG), melanocortin 1 receptor (MC1R) and tyrosinase (TYR). Here, we show that the coat colour variations in White Galloway cattle and White Park cattle are caused by a KIT gene (chromosome 6) duplication and aberrant insertion on chromosome 29 (Cs₂₉) as recently described for colour‐sided Belgian Blue. Homozygous (Cs₂₉/Cs₂₉) White Galloway cattle and White Park cattle exhibit the mismarked phenotype, whereas heterozygous (Cs₂₉/wt₂₉) individuals are either well or strongly marked. In contrast, fully black individuals are characterised by the wild‐type chromosome 29. As known for other cattle breeds, mutations in the MC1R gene determine the red colouring. Our data suggest that the white coat colour variations in White Galloway cattle and White Park cattle are caused by a dose‐dependent effect based on the ploidy of aberrant insertions and inheritance of the KIT gene on chromosome 29.     

The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.     

Low genome‐wide homozygosity in 11 Spanish ovine breeds

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (F_ROH) are generally small (F_ROH < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.     

Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (F_ROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (F_HOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (F_PED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. F_HOM was highly correlated (0.671 to 0.985) with F_ROH measures in all breeds. Apulo-Calabrese and Casertana had the highest F_ROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (F_ROH1) to 16 Mbp (F_ROH16)), whereas the lowest F_ROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). F_ROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all F_ROH measures and F_PED (from 0.514 to 0.523) and between F_HOM and F_PED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between F_ROH and F_PED. This might be due to the imperfect measure of F_PED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.     

Novel insights into Sabino1 and splashed white coat color patterns in horses

Within the framework of genome‐wide analyses using the novel Axiom^® genotyping array, we investigated the distribution of two previously described coat color patterns, namely sabino1 (SBI), associated with the KIT gene (KI16+1037A), and splashed white, associated with the PAX3 gene (ECA6:g.11429753C>T; PAX3^C70Y), including a total of 899 horses originating from eight different breeds (Achal Theke, Purebred Arabian, Partbred Arabian, Anglo‐Arabian, Shagya Arabian, Haflinger, Lipizzan and Noriker). Based on the data we collected we were able to demonstrate that, besides Quarter horses, the PAX3^C70Y allele is also present in Noriker (seven out of 189) and Lipizzan (three out of 329) horses. The SB1 allele was present in three breeds (Haflinger, 14 out of 98; Noriker, four out of 189; Lipizzan one out of 329). Furthermore, we examined the phenotypes of SB1‐ and PAX3^C70Y‐carrier horses for their characteristic white spotting patterns. None of the SB1/sb1‐carrier horses met the criteria defining the Sabino1 pattern according to current applied protocols. From 10 heterozygous PAX3^C70Y‐carrier horses, two had nearly a splashed white phenotype. The results of this large‐scale experiment on the genetic association of white spotting patterns in horses underline the influence of gene interactions and population differences on complex traits such as Sabino1 and splashed white.     

Pedigree estimation of the (sub) population contribution to the total gene diversity: the horse coat colour case

A method to quantify the contribution of subpopulations to genetic diversity in the whole population was assessed using pedigree information. The standardization of between- and within-subpopulation mean coancestries was developed to account for the different coat colour subpopulation sizes in the Spanish Purebred (SPB) horse population. The data included 166264 horses registered in the SPB Studbook. Animals born in the past 11 years (1996 to 2006) were selected as the 'reference population' and were grouped according to coat colour into eight subpopulations: grey (64 836 animals), bay (33 633), black (9414), chestnut (1243), buckskin (433), roan (107), isabella (57) and white (37). Contributions to the total genetic diversity were first assessed in the existing subpopulations and later compared with two scenarios with equal subpopulation size, one with the mean population size (13 710) and another with a low population size (100). Ancestor analysis revealed a very similar origin for the different groups, except for six ancestors that were only present in one of the groups likely to be responsible for the corresponding colour. The coancestry matrix showed a close genetic relationship between the bay and chestnut subpopulations. Before adjustment, Nei's minimum distance showed a lack of differentiation among subpopulations (particularly among the black, chestnut and bay subpopulations) except for isabella and white individuals, whereas after adjustment, white, roan and grey individuals appeared less differentiated. Standardization showed that balancing coat colours would contribute preserving the genetic diversity of the breed. The global genetic diversity increased by 12.5% when the subpopulations were size standardized, showing that a progressive increase in minority coats would be profitable for the genetic diversity of this breed. The methodology developed could be useful for the study of the genetic structure of subpopulations with unbalanced sizes and to predict their genetic importance in terms of their contribution to genetic variability.     

Who's behind that mask and cape? The Asian leopard cat's Agouti (ASIP) allele likely affects coat colour phenotype in the Bengal cat breed

Coat colours and patterns are highly variable in cats and are determined mainly by several genes with Mendelian inheritance. A 2‐bp deletion in agouti signalling protein (ASIP) is associated with melanism in domestic cats. Bengal cats are hybrids between domestic cats and Asian leopard cats (Prionailurus bengalensis), and the charcoal coat colouration/pattern in Bengals presents as a possible incomplete melanism. The complete coding region of ASIP was directly sequenced in Asian leopard, domestic and Bengal cats. Twenty‐seven variants were identified between domestic and leopard cats and were investigated in Bengals and Savannahs, a hybrid with servals (Leptailurus serval). The leopard cat ASIP haplotype was distinguished from domestic cat by four synonymous and four non‐synonymous exonic SNPs, as well as 19 intronic variants, including a 42‐bp deletion in intron 4. Fifty‐six of 64 reported charcoal cats were compound heterozygotes at ASIP, with leopard cat agouti (A^Pbe) and domestic cat non‐agouti (a) haplotypes. Twenty‐four Bengals had an additional unique haplotype (A2) for exon 2 that was not identified in leopard cats, servals or jungle cats (Felis chaus). The compound heterozygote state suggests the leopard cat allele, in combination with the recessive non‐agouti allele, influences Bengal markings, producing a darker, yet not completely melanistic coat. This is the first validation of a leopard cat allele segregating in the Bengal breed and likely affecting their overall pelage phenotype. Genetic testing services need to be aware of the possible segregation of wild felid alleles in all assays performed on hybrid cats.