Complex Inheritance of Melanoma and Pigmentation of Coat and Skin in Grey Horses

Abstract 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 (h2), and we decomposed this parameter into polygenic heritability (h2 POLY), heritability due to the Grey (STX17) mutation (h2 STX17), and heritability due to agouti (ASIP) locus (h2 ASIP). A high heritability was found for greying (h2 = 0.79), vitiligo (h2 = 0.63), and speckling (h2 = 0.66), while a moderate heritability was estimated for melanoma (h2 = 0.37). The additive component of ASIP was significantly different from zero only for melanoma (h2 ASIP = 0.02). STX17 controlled large proportions of phenotypic variance (h2 STX17 = 0.18–0.55) and overall heritability (h2 STX17/h2 = 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.


Introduction
Recent developments in molecular genetics have enabled molecular dissection of quantitative traits in humans [1], model organisms [2] and domestic animals [3]. Genome-wide association (GWA) studies suggest that variability of complex traits is caused by many loci, most exerting tiny effects, whereas loci exerting moderate-tolarge effects or loci that explain more than 5-10% of phenotypic variation are rare [3,4,5]. In human populations, those genes with moderate-to-large effects do appear in low frequency as rare or ''private'' mutations [6]. In contrast, appearance of moderate-to-large effect mutations at intermediate or high frequencies is documented in domestic animal populations, perhaps, as consequence to change in selection pressure caused by domestication [7].
The genetic variation of several complex traits in Grey horses is considerably affected by at least one gene of moderate-to-large effect. Grey horses are born with their base colour (e.g. black, bay, chestnut), which then greys early in life due to loss of melanocytes, a process similar to the greying of hair in humans, which typically occurs much later in life. The mode of inheritance of the greying phenotype is autosomal dominant. A grey horse will be either GG or Gg, non-grey horses carry the gg genotype. Previously we reported that greying with age is caused by a 4.6-kb duplication in intron 6 of STX17, which encodes syntaxin 17 [8]. In addition to grey level, the Grey mutation was found to strongly influence melanoma, vitiligo and speckling grade as well as to capacitate the effects of the loss-of-function mutation at the agouti (ASIP) locus on melanoma grade. Recently, we have shown that the duplicated region contains a melanocyte-specific enhancer that becomes much stronger when duplicated [9]. Examination of melanoma tissue revealed copy number variation of the 4.6 kb intron in STX17, with a difference between blood and tumor DNA. The highest copy numbers occurred in tumors classified as aggressive [10]. Grey horses are the only animals to present progressive coat greying, melanoma, vitiligo-like skin depigmentation and coat speckling, making them an excellent model for studying the inheritance of traits with complex genetic background.
Studies of coat colour in horses led to the definition of discrete phenotypes (e.g. bay, black, brown, grey) controlled by a few genes showing epistatic interaction [11]. This typology, although practical for breeders to some extent, could not account for variation in greying level. Curik et al. [12] provided a quantitative description of progressive greying, while Toth et al. [13] found relationship between that quantitative greying measure and total melanin content of horse hair.
Melanoma occurs frequently in grey horses, in three histopathologically defined clinical patterns. In the first pattern, which describes the majority of cases, the melanoma grows slowly over many years without evidence of regional or distant metastases. In the second pattern, the melanoma results from malignant transformation of a benign melanoma (melanocytoma; [14]). In the third pattern, which is rare, the melanoma is malignant from onset [15]. Although most melanomas in Grey horses present benign features at diagnosis, nearly 66% become malignant later [16]. Melanoma in Grey horses occurs most frequently underneath the tail, in the peri-anal region, and around the lips and eyelids [17,18]. The frequency of melanoma occurrence is around 80% in horses older than 15 years [19].
The hereditary component of melanoma in Grey horses was first studied by Rieder et al. [20]. Using segregation analysis, they were unable to establish whether the mode of inheritance was monogenic, polygenic or mixed, because of the relatively low number of horses examined (n = 71). Nevertheless, models including a polygenic component fitted the data significantly better than did a non-genetic model. In a more recent study involving 296 grey Lipizzan horses, Seltenhammer et al. [21] estimated a heritability of 0.36 for melanoma grade.
Grey horses also show vitiligo-like depigmentation. In fact, Seltenhammer et al. [21] found equine vitiligo in 50% of older grey Lipizzaner horses. In humans, vitiligo is a chronic disorder characterized by sharply delimited, progressive, patchy loss of pigmentation resulting from death or dysfunction of the cells responsible for skin and hair pigmentation [22]. While the pathogenesis of human vitiligo is unknown, it is considered a complex hereditary disease, and several theories have been proposed to explain it [23,24]. Similarly, although linked to progressive greying, the pathogenesis of equine vitiligo remains largely unknown. Like melanoma, vitiligo occurs far more frequently in Grey horses than for those with solid colours, indicating that it has a hereditary component that may be associated with progressive greying. Although rare, studies in humans seem to indicate that melanoma and vitiligo might be genetically linked [25,26].
When this study was conceived, we hypothesized that while the effects of STX17 are large for all four traits and ASIP has a considerable effect of melanoma, a large part of the genetic variation remains unexplained. We further hypothesized that the genetic correlations among those traits are strongly influenced by STX17. Here we have extended our previous studies by increasing the number of genotyped horses and parameters examined, including genotype at the STX17 and ASIP loci. Firstly, for all four traits (melanoma, grey level, vitiligo, speckling), we decomposed estimated repeatability into a component due to permanent environmental effects and a component describing narrow sense heritability. We further decomposed estimated narrow sense heritability into a component describing background polygenic heritability and a components describing heritability caused by STX17 and ASIP. We also estimated genetic correlations among all four traits using various models designed to evaluate the contribution from STX17 and ASIP genes and the contribution from pleiotropy due to additional polygenic background versus that of STX17 and ASIP genes. Finally, we graphically illustrated the genetic relationship among melanoma grade, grey level, vitiligo grade and speckling grade with respect to STX17 genotype.

Results
Estimated fixed effects and variance components for melanoma grade, grey level, vitiligo grade and speckling grade are shown in Table 1, together with means and standard errors. Linear effects of age were significant (P,0.001) for all four traits. Thus, for example, a 6-year-old grey horse will present a grey colour that is 5.88 L* units darker than that of a 7-year-old. Linear regression coefficients relating each trait to STX17 were highly significant (P,0.001) for all traits studied, while the corresponding coefficients involving ASIP were significant only for melanoma grade (P,0.001). In quantitative genetics, estimated linear regression coefficients can be interpreted as gene substitution effects of STX17 (a STX17 ) and ASIP (a ASIP ). Thus, the regression coefficient of 213.78 for the Grey mutation in relation to the trait Greying indicates that heterozygotes Gg on average will be 13.78 L* units darker than homozygous GG horses. ASIP has three genotypes, so the melanoma grade is expected to change by 0.38 units (2a ASIP ) when moving from an AA to an aa genotype.
The frequency of the Grey allele was between 0.85 and 0.89 for the various data sets. Univariate estimates were obtained for heritability (h 2 ) and for repeatability (R), which here corresponds to the upper bound of broad sense heritability, together with estimates of phenotypic variance (V P ) and variance due to permanent environmental effects (c 2 ), here variance introduced by repeated measurement ( Table 2). Proportions of variance in grey level and vitiligo grade due to permanent environmental effects were 0.03 and 0.00, respectively, whereas these proportions were much larger for melanoma grade (0.26) and speckling grade (0.12). Very high h 2 values, that explained 63-79% of the phenotypic variance, were estimated for grey level, vitiligo grade and speckling grade. In contrast, moderate h 2 was estimated for melanoma grade. Across all models standard errors of hertiability

Author Summary
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.
estimates ranged from 0.08 to 0. 12 Multivariate models were used to estimate phenotypic (r P ) and genetic correlations (r POLY ) according to three scenarios (Table 3). In the first scenario, we estimated genetic correlations with a standard polygenic model, ignoring the effects of STX17 and ASIP.
In the second scenario, we used a model that includes the fixed effects of these two mutations. For the third scenario, we applied a polygenic model only to the data for GG horses. In the first scenario, the genetic correlations among the traits followed the same pattern as the phenotypic ones, with the genetic values occasionally somewhat higher (e.g. 0.67 vs. 0.52 in the case of grey level and vitiligo grade). When the genetic effects were included (second scenario) or when only homozygous grey horses were considered (third scenario), correlations dropped substantially. This is a strong indication of the pleiotropic effects of the Grey mutation. Only the genetic relationship between grey level and vitiligo grade remained significant and considerable across all three scenarios (0.6760.06 in the first, 0.4860.09 in the second, and 0.5060.10 in the third). Visualization of the genetic correlations due to polygenic effects is presented by scatter plots of estimated breeding values ( Figure 1). Distances from a homozygous (GG) genotype of STX17 to a heterozygous one (Gg) have a diagonal shift in the contour plot centres. In addition to these pleiotropic effects of the STX17 mutation, the fact that the contour plots for breeding values were elliptical specifically for the GG and Gg genotypes indicates that genetic correlations due to polygenic background additive effects remained even after accounting for the STX17 mutation. This was particularly obvious for the genetic correlation between grey level and vitiligo grade, for which the estimated genetic correlation after accounting for major genes remained significant (r POLY = 0.4860.09).

Discussion
We present a case in which the genetic components of four complex traits could be decomposed into the effects of polygenic additive effects and the monogenic effects of STX17 and ASIP mutations. The data available was still insufficient to allow accurate estimation of polygenic dominance variance or of higher-order, non-additive genetic variances. Although melanoma and vitiligo have been intensively studied in humans, few estimates . **Genetic parameters were calculated as follows: V P = V POLY +V STX17 +V ASIP +V pe +V e ; R = (V POLY +V STX17 +V ASIP +V pe ,)/V P ; c 2 = V pe /V P ; h 2 = (V POLY +V STX17 +V ASIP )/V P ; h 2 POLY = V POLY /V P ; h 2 STX17 = V STX17 /V P ; h 2 ASIP = V ASIP / V P . For all traits except melanoma, we assumed V ASIP = 0. The frequencies ranged from 0.85 to 0.89 for the STX17 G allele and from 0.51 to 0.53 for the ASIP A allele in the various data sets. doi:10.1371/journal.pgen.1003248.t002 of their quantitative inheritance are available. The present study provides some of the first quantitative insights into their inheritance in grey horses. While STX17 explains large proportions of the phenotypic variance, we were surprised to see that, when only the data from homozygous GG horses were analysed, the residual polygene component still explained a large part of the variation in melanoma and vitiligo grade (Table 3). These results argue for the need to search for more genes involved in the expression of these traits.
The STX17 mutation accounts for considerable amount of the estimated genetic correlations among all traits analysed; in other words, it exhibits strong pleiotropic effects. While precise estimates of genetic correlation require extremely large data sets, our findings are supported by scatter plots of estimated breeding values ( Figure 1). Graphical illustrations also indicated that polygenic additive effects were causing negative genetic correlations between melanoma grade and speckling grade as well as between grey level and speckling grade. However, we were not able to confirm those evidences numerically since related estimates of genetic correlations had very wide confidence intervals (Table 3). On the other side, we were able to show that estimated genetic correlation between melanoma and vitiligo was negligible after accounting for the pleiotropic effect of STX17 mutation (Table 3 and Figure 1). However, it is highly speculative to conclude the same pattern, just one gene linking melanoma and vitiligo, is applying to humans.
The results obtained here, further, support 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 at intermediate to high frequencies. For example, a single nucleotide substitution in intron 3 of IGF2 in pigs explain about 30% of the residual phenotypic variance for lean meat in ham in a wild boar/Large White intercross [27,28]. Grisart et al. [29] have also shown that the mutation in DGAT1, appearing in different breeds at varying   [30]. Moderate-to large effects have been also demonstrated for the myostatin mutation on the meat and carcass quality in beef cattle [31]. The genetic architecture of complex traits presented here, in which quantitative variation is explained by a strong contribution from a single gene and/or mutation together with a substantial polygenic additive component, may be frequent in traits related to colouration and pigmentation of coat and skin as successful adaptation requires fast and slow responses. A similar genetic architecture was presented by Hayes et al. [32] in a study where KIT, MITF, and a locus on chromosome 8 together explained 24% of the variation in the proportion of black coat in cattle. Recently, Liu et al. [33] also indicated that a few genes played a major role in human eye colour. Coat and skin colouration, the consequences of melanindependent pigmentation, are evolutionarily important traits because they are involved in several aspects of survival, such as hiding from predators and thermoregulation [34]. Fang et al. [35] concluded that coat colour variation in domestic animals has been shaped by artificial selection. The increase in the frequency of mutated allele G, here up to 89%, is strongly associated with human appreciation of grey colour in horses. The high incidence of melanoma and vitiligo in Grey horses is most likely the result of negative pleiotropic effects of the mutations contributing to the grey phenotype, in particular the Grey mutation itself. The presence of this prime example provides, thus, an additional insight into genetic architecture of complex traits and their biology.

Quantitative description of the melanoma phenotype
Melanoma was detected by adspection and palpation. Melanoma grade was defined according to a modified classification system on a scale from 0.0 to 5.0 that allows for intermediate grades (e.g. 0.0, 0.5, 4.5), described in [17] (Table 4). Adspection was conducted at sites melanomas typically occur, such as the perianal and anal region, the perineal region, udder, and praeputium. The tail was bent upwards in order to detect even the smallest, plaque-like lesions. Lips and eyelids, the parotis, the peri-ocular region and ears were also examined. Finally, the whole integument was checked for potential tumours. Variation of melanoma grade in four grey Lipizzan horses is illustrated in Figure 2.

Quantitative description of the progressive greying phenotype
Progressive greying was quantified with a Minolta Chromameter CR210 using the CIE L*a*b* colour space. In this system, colour is quantified by its reflection along three axes: white-black (L*), red-green (a*) and yellow-blue (b*). We referred only to the L* parameter defined on a scale from 0 (black -total absorption) to 100 (white -total reflection). Various grey levels in seven Lipizzan horses are illustrated in Figure 3. Each horse was measured at four places (neck, shoulder, belly and croup), and the average of the four measurements was used to represent the grey level.

Quantitative description of the vitiligo phenotype
Vitiligo was graded on a scale from 0 (no vitiligo) to 3 (severe vitiligo) based on adspection of typical sites, such as the peri-anal and anal region, the perineal region, udder, praeputium and the face, especially around the nostrils and eyes. Vitiligo grade was evaluated simultaneously with melanoma. The overall vitiligo grade was the average of included vitiligo grades from patches in the perianal and the facial regions. Variation of vitiligo grade across four grey Lipizzan horses is illustrated in Figure 4.

Quantitative description of the speckling phenotype
Adult grey horses may present coloured specks or spots on a grey background. The amount of speckling was graded on a scale from 0 (not speckled) to 3 (heavily speckled). Variation in speckling grade across four grey Lipizzan horses is illustrated in Figure 5.

Genotyping
Genotyping for the STX17 mutation was carried out on 760 horses and for the ASIP gene on 667 horses; the agouti locus controls the inheritance of black and bay base colour [36]. Genotyping procedures were performed as described in Pielberg et al. [8]. The genotypes of additional horses were deduced by combining genotyping and pedigree information. For instance, a single coloured offspring qualifies a grey parent as heterozygous, grey offspring with one coloured parent must be heterozygous. In this way we extended the number of genotyped horses to 966 for the STX17 mutation and 873 for the ASIP. However, deducing genotypes may introduce bias in the estimation of major gene effects. To avoid potential bias, we performed analyses on data sets with and without deduced genotypes. The results were very similar, those presented here are therefore based on data sets containing deduced genotypes.

Estimating (co)variance components and gene effects
All traits analysed (grey level and grades of melanoma, vitiligo and speckling) strongly depend on the age of the horse. Grey level can vary widely in animals during the first 6-10 years; from 10 years onwards, this variation drops considerably and all horses reach their final coat colour. Both melanoma and vitiligo show an age of onset of 6-8 years, though they do appear earlier in rare cases. Speckling has been defined only for adult horses, because the density of spots or speckles is difficult to assess while horses are still dark. Thus, we analysed grey level data only for horses seven years old and younger. Conversely, data on melanoma, vitiligo and speckling grade were used only for horses seven years old and older. The dynamics of melanoma, grey level, vitiligo and speckling with respect to age has been shown in Pielberg et al. [8].
Statistical analyses were performed with univariate and multivariate general mixed linear models using the ASReml package (version 3; [37]). These models are known as (Individual) Animal models; see [38] for more detailed explanation. For all of the variables analysed, the statistical models showed no serious statistical violations; residuals from the models followed a normal distribution and variances were roughly homogeneous. The following modeling strategy was applied. First, we fitted a model that included the effects of the stud by sex by year of measurement (fixed), age at measurement in months (covariate), STX17 mutation effects (indicator variables, with GG = 0 and Gg = 1 treated as covariates), ASIP additive effects (indicator variables, with AA = 0, Aa = 1 and aa = 2 treated as covariates) and ASIP dominance effects (indicator variables, with AA = 0, Aa = 1 and aa = 0 treated as covariates) and all two-factor interactions. Two factor interactions and the ASIP dominance effects were nonsignificant for any trait and were dropped from further analyses. Initial and all subsequent models also included additive genetic animal effects (random) and permanent environmental effects (random), since measurements were taken repeatedly.
The following final model was used to analyse the data: y~XbzZazWpeze where the vector y represents the phenotypic values; b is the vector of fixed effects; a, pe and e are vectors of, respectively, polygenic additive effects (breeding values), permanent environmental (repeated measurement) effects and residual error values; X is the design matrix for fixed effects; and Z and W are design matrices for random effects. For all four traits analysed, the final models accounted for the stud by sex by year of measurement effect, for the effect of age at measurement and for the effect of STX17, while the effect of ASIP was also accounted for only in the case of melanoma grade. In addition, additive polygenic (V POLY ), permanent environmental (V pe ) and residual variances (V e ), respectively, were defined as where A is an additive genetic relationship matrix and I is an identity matrix. In multivariate models we respected the modeling results obtained in univariate models. Based on variance components and effects estimated from statistical analyses, together with allele frequencies obtained by counting, we derived quantitative genetic parameters to explain the complex inheritance of the traits analysed. Additive genetic variances of STX17 mutation (V STX17 ) and ASIP (V ASIP ) were calculated as V STX17 or V ASIP = 2pqa 2 , where estimated linear regression coefficients are equal to gene substitution effects a (a = a for STX17 and a = a+(p2q)d for ASIP), p and q are allele frequencies of STX17 or/and ASIP, a is an additive value and d is a dominance value. Calculation of single gene (mutation) variances enabled us to calculate additive genetic variance as a sum of V POLY , V STX17 and V ASIP , and phenotypic variance as a sum of V POLY , V STX17 , V ASIP , V pe , and V e . In further analysis, we calculated proportions of genetically explained variations from the phenotypic variance as follows: repeatability, defined as correlation among observations within individuals that provides an upper limit to broad sense heritability, R = (V POLY +V STX17 + V ASIP +V pe ,)/V P ; permanent environmental effects, c 2 = V pe /V P ; narrow sense heritability, h 2 = (V POLY +V A-STX17 +V A-ASIP )/V P ; background polygenic heritability, h 2 POLY = V POLY /V P ; STX17 mutation heritability, h 2 STX17 = V STX17 /V P ; and ASIP heritability, h 2 ASIP = V ASIP /V P . For grey level, vitiligo grade and speckling grade, a ASIP was not significant, so we assumed V ASIP was zero. The analyses were performed for three scenarios: the first and third ignored STX17 and ASIP effects, using the full data set and a  data set containing GG animals only; the second applied the full model as described above for melanoma and excluded ASIP (because not significant) for the other traits.
Bivariate models were run for all pairs of traits applying the single trait model specifications from above. Genetic (r POLY ) and environmental (r e ) correlations between traits x 1 and x 2 were calculated as r POLY = Cov POLY_ (x 1 , x 2 )/[(V POLY_ x 1 )(V POLY_ x 2 )] 0.5 and r e = Cov e (x 1 , x 2 )/[(V e x 1 )(V e x 2 )] 0.5 , while phenotypic correlations (r P ) were calculated as r P = r POLY (h 2 x1 ) 0.5 (h 2 x2 ) 0.5 + r e [(12h 2 x1 )] 0.5 [(12h 2 x2 )] 0.5 . Contour plots with 95% confidence intervals, which present breeding values standardized to a mean of zero and a standard deviation of one, were created for two traits obtained from bivariate animal models with respect to STX17 genotype. These plots illustrate the decomposition of genetic correlations into additive polygenic background effects and effects of a single mutation. More detailed derivations of quantitative genetic parameters are provided in [39].