Background: Runs of homozygosity are long, uninterrupted stretches of homozygous genotypes that enable reliable estimation of levels of inbreeding (i.e., autozygosity) based on high-throughput, chip-based single nucleotide polymorphism (SNP) genotypes. While the theoretical definition of runs of homozygosity is straightforward, their empirical identification depends on the type of SNP chip used to obtain the data and on a number of factors, including the number of heterozygous calls allowed to account for genotyping errors. We analyzed how SNP chip density and genotyping errors affect estimates of autozygosity based on runs of homozygosity in three cattle populations, using genotype data from an SNP chip with 777 972 SNPs and a 50 k chip. Results: Data from the 50 k chip led to overestimation of the number of runs of homozygosity that are shorter than 4 Mb, since the analysis could not identify heterozygous SNPs that were present on the denser chip. Conversely, data from the denser chip led to underestimation of the number of runs of homozygosity that were longer than 8 Mb, unless the presence of a small number of heterozygous SNP genotypes was allowed within a run of homozygosity. Conclusions: We have shown that SNP chip density and genotyping errors introduce patterns of bias in the estimation of autozygosity based on runs of homozygosity. SNP chips with 50 000 to 60 000 markers are frequently available for livestock species and their information leads to a conservative prediction of autozygosity from runs of homozygosity longer than 4 Mb. Not allowing heterozygous SNP genotypes to be present in a homozygosity run, as has been advocated for human populations, is not adequate for livestock populations because they have much higher levels of autozygosity and therefore longer runs of homozygosity. When allowing a small number of heterozygous calls, current software does not differentiate between situations where these calls are adjacent and therefore indicative of an actual break of the run versus those where they are scattered across the length of the homozygous segment. Simple graphical tests that are used in this paper are a current, yet tedious solution.

The dominant phenotype of greying with age in horses, caused by a 4.6-kb duplication in intron 6 of STX17, is associated with a high incidence of melanoma and vitiligo-like skin depigmentation. However, the progressive greying and the incidence of melanoma, vitiligo-like depigmentation, and amount of speckling in these horses do not follow a simple inheritance pattern. To understand their inheritance, we analysed the melanoma grade, grey level, vitiligo grade, and speckling grade of 1,119 Grey horses (7,146 measurements) measured in six countries over a 9-year period. We estimated narrow sense heritability (h²), and we decomposed this parameter into polygenic heritability (h²_POLY), heritability due to the Grey (STX17) mutation (h²_STX17), and heritability due to agouti (ASIP) locus (h²_ASIP). A high heritability was found for greying (h² = 0.79), vitiligo (h² = 0.63), and speckling (h² = 0.66), while a moderate heritability was estimated for melanoma (h² = 0.37). The additive component of ASIP was significantly different from zero only for melanoma (h²_ASIP = 0.02). STX17 controlled large proportions of phenotypic variance (h²_STX17 = 0.18–0.55) and overall heritability (h²_STX17/h² = 0.28–0.83) for all traits. Genetic correlations among traits were estimated as moderate to high, primarily due to the effects of the STX17 locus. Nevertheless, the correlation between progressive greying and vitiligo-like depigmentation remained large even after taking into account the effects of STX17. We presented a model where four traits with complex inheritance patterns are strongly influenced by a single mutation. This is in line with evidence of recent studies in domestic animals indicating that some complex traits are, in addition to the large number of genes with small additive effects, influenced by genes of moderate-to-large effect. Furthermore, we demonstrated that the STX17 mutation explains to a large extent the moderate to high genetic correlations among traits, providing an example of strong pleiotropic effects caused by a single gene. Clarifying the genetic architecture of complex traits is a problem with profound implications for agriculture, biology, and medicine. Using data from Lipizzan horses with the grey coat phenotype, we present an example of a single mutation (intronic duplication in STX17) that explains 18%–55% of phenotypic variation in four complex traits, while polygenic background additive effects also explain 11%–57% of phenotypic variation. This study provides a prime example of complex traits being influenced by genes of moderate-to-large effect and supports further the evidence of recent studies in domestic animals that some complex traits are, in addition to the large number of genes with small additive effects, influenced by genes of moderate-to-large effect. We further show that the STX17 mutation accounts for a large proportion of the estimated genetic correlations between the traits. This case of strong pleiotropic effects of a single mutation on complex traits makes this work of significant general interest for biology and medicine.

While the negative effects of inbreeding and reduced heterozygosity on fecundity and survival are well established, only a few investigations have been carried out concerning their influence on morphological traits. This topic is of particular interest for a small and closed population such as the Lipizzan horse. Thus, 27 morphological traits were measured in 360 Lipizzan mares and were regressed on the individual inbreeding coefficients, as well as on the individual heterozygosity and mean squared distances (mean d²) between microsatellite alleles within an individual. Both individual heterozygosity and mean d² were based on 17 microsatellite loci dispersed over 14 chromosomes. The results obtained by multivariate analysis reveal significant effects of stud (P <.0001), age at measurement (P <.0001), and mean d² (P =.0143). In univariate analyses, significant associations were obtained between length of pastern-hindlimbs and inbreeding coefficient (P <.01), length of cannons-hindlimb and mean d² (P <.01), and length of neck and mean d² (P <.001). After adjustment of single-test P values for multiple tests (Hochberg's step-up Bonferroni method), only the association of the length of neck and mean d² remained significant (P =.0213). Thus, no overall large effects of inbreeding, microsatellite heterozygosity, and mean d² on morphological traits were observed in the Lipizzan horse.

Background: Greying with age in horses is an autosomal dominant trait, associated with loss of hair pigmentation, melanoma and vitiligo-like depigmentation. We recently identified a 4.6 kb duplication in STX17 to be associated with the phenotype. The aims of this study were to investigate if the duplication in Grey horses shows copy number variation and to exclude that any other polymorphism is uniquely associated with the Grey mutation. Results: We found little evidence for copy number expansion of the duplicated sequence in blood DNA from Grey horses. In contrast, clear evidence for copy number expansions was indicated in five out of eight tested melanoma tissues or melanoma cell lines. A tendency of a higher copy number in aggressive tumours was also found. Massively parallel resequencing of the ~350 kb Grey haplotype did not reveal any additional mutations perfectly associated with the phenotype, confirming the duplication as the true causative mutation. We identified three SNP alleles that were present in a subset of Grey haplotypes within the 350 kb region that shows complete linkage disequilibrium with the causative mutation. Thus, these three nucleotide substitutions must have occurred subsequent to the duplication, consistent with our interpretation that the Grey mutation arose more than 2,000 years before present. Conclusions: These results suggest that the mutation acts as a melanoma-driving regulatory element. The elucidation of the mechanistic features of the duplication will be of considerable interest for the characterization of these horse melanomas as well as for the field of human melanoma research.

Runs of homozygosity (ROH) islands are stretches of homozygous sequence in the genome of a large proportion of individuals in a population. Algorithms for the detection of ROH depend on the similarity of haplotypes. Coverage gaps and copy number variants (CNV) may result in incorrect identification of such similarity, leading to the detection of ROH islands where none exists. Misidentified hemizygous regions will also appear as homozygous based on sequence variation alone. Our aim was to identify ROH islands influenced by marker coverage gaps or CNV, using Illumina BovineHD BeadChip (777 K) single nucleotide polymorphism (SNP) data for Austrian Brown Swiss, Tyrol Grey and Pinzgauer cattle. ROH were detected using clustering, and ROH islands were determined from population inbreeding levels for each marker. CNV were detected using a multivariate copy number analysis method and a hidden Markov model. SNP coverage gaps were defined as genomic regions with intermarker distances on average longer than 9.24 kb. ROH islands that overlapped CNV regions (CNVR) or SNP coverage gaps were considered as potential artefacts. Permutation tests were used to determine if overlaps between CNVR with copy losses and ROH islands were due to chance. Diversity of the haplotypes in the ROH islands was assessed by haplotype analyses. In Brown Swiss, Tyrol Grey and Pinzgauer, we identified 13, 22, and 24 ROH islands covering 26.6, 389.0 and 35.8 Mb, respectively, and we detected 30, 50 and 71 CNVR derived from CNV by using both algorithms, respectively. Overlaps between ROH islands, CNVR or coverage gaps occurred for 7, 14 and 16 ROH islands, respectively. About 37, 44 and 52% of the ROH islands coverage in Brown Swiss, Tyrol Grey and Pinzgauer, respectively, were affected by copy loss. Intersections between ROH islands and CNVR were small, but significantly larger compared to ROH islands at random locations across the genome, implying an association between ROH islands and CNVR. Haplotype diversity for reliable ROH islands was lower than for ROH islands that intersected with copy loss CNVR. Our findings show that a significant proportion of the ROH islands in the bovine genome are artefacts due to CNV or SNP coverage gaps.