Mapping the four‐horned locus and testing the polled locus in three Chinese sheep breeds

Four‐horned sheep are an ideal animal model for illuminating the genetic basis of horn development. The objective of this study was to locate the genetic region responsible for the four‐horned phenotype and to verify a previously reported polled locus in three Chinese breeds. A genome‐wide association study (GWAS) was performed using 34 two‐horned and 32 four‐horned sheep from three Chinese indigenous breeds: Altay, Mongolian and Sishui Fur sheep. The top two significant single nucleotide polymorphisms (SNPs) associated with the four‐horned phenotype were both located in a region spanning positions 132.6 to 132.7 Mb on sheep chromosome 2. Similar locations for the four‐horned trait were previously identified in Jacob, Navajo‐Churro, Damara and Sishui Fur sheep, suggesting a common genetic component underlying the four‐horned phenotype. The two identified SNPs were both downstream of the metaxin 2 (MTX2) gene and the HOXD gene cluster. For the top SNP—OAR2:g.132619300G>A—the strong associations of the AA and AG genotypes with the four‐horned phenotype and the GG genotype with the two‐horned phenotype indicated the dominant inheritance of the four‐horned trait. No significant SNPs for the polled phenotype were identified in the GWAS analysis, and a PCR analysis for the detection of the 1.8‐kb insertion associated with polled sheep in other breeds failed to verify the association with polledness in the three Chinese breeds. This study supports the hypothesis that two different loci are responsible for horn existence and number. This study contributes to the understanding of the molecular regulation of horn development and enriches the knowledge of qualitative traits in domestic animals.     

Genome‐wide association reveals the locus responsible for four‐horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo‐Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome‐wide association mapping was performed using 125 animals from the two breeds that contain two‐ and four‐horned individuals. A case–control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9–132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two‐horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four‐horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.     

Genetic control of alopecia in sheep

A hypothesis is suggested which turns down the opinion that alopecia or pathological wool shedding in sheep is due to some pathological processes or disease. We suppose that this trait is controlled by a gene of natural wool shedding inherited from sheep far ancestors, and we propose to designate the recessive alopecia gene as Alp3 and its dominant allele, a wool alopecia resistance gene, as AlpA. Alopecia will only be revealed in homozygotic Alpa/Alpa state, due to recessivity of alopecia gene, but heterozygotic AlpA/Alpa and homozygotic AlpA/AlpA sheep will have a normal wool coat. A new polymorphic trait is supposed to exist which is conditioned by simultaneous presence of two genetic alleles in a sheep population with descrete phenotypic manifestation.     

A genome‐wide association study reveals a quantitative trait locus for days open on chromosome 2 in Japanese Black cattle

Days open (DO), which is the interval from calving to conception, is an important trait related to reproductive performance in cattle. To identify quantitative trait loci for DO in Japanese Black cattle, we conducted a genome‐wide association study with 33 303 single nucleotide polymorphisms (SNPs) using 459 animals with extreme DO values selected from a larger group of 15 488 animals. We identified a SNP on bovine chromosome 2 (BTA2) that was associated with DO. After imputation using phased haplotype data inferred from 586 812 SNPs of 1041 Japanese Black cattle, six SNPs associated with DO were located in an 8.5‐kb region of high linkage disequilibrium on BTA2. These SNPs were located on the telomeric side at a distance of 177 kb from the parathyroid hormone 2 receptor (PTH2R) gene. The association was replicated in a sample of 1778 animals. In the replicated population, the frequency of the reduced‐DO allele (Q) was 0.63, and it accounted for 1.72% of the total genetic variance. The effect of a Q‐to‐q allele substitution on DO was a decrease of 3.74 days. The results suggest that the Q allele could serve as a marker in Japanese Black cattle to select animals with superior DO performance.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.     

Genome-wide association study to identify genetic loci associated with gastrointestinal nematode resistance in Katahdin sheep

Resistance to gastrointestinal nematodes has previously been shown to be a moderately heritable trait in some breeds of sheep, but the mechanisms of resistance are not well understood. Selection for resistance currently relies upon faecal egg counts (FEC), blood packed cell volumes and FAMACHA visual indicator scores of anaemia. Identifying genomic markers associated with disease resistance would potentially improve the selection process and provide a more reliable means of classifying and understanding the biology behind resistant and susceptible sheep. A GWAS was conducted to identify possible genetic loci associated with resistance to Haemonchus contortus in Katahdin sheep. Forty animals were selected from the top and bottom 10% of estimated breeding values for FEC from a total pool of 641 sires and ram lambs. Samples were genotyped using Applied Biosystems™ Axiom™ Ovine Genotyping Array (50K) consisting of 51 572 SNPs. Following quality control, 46 268 SNPs were included in subsequent analyses. Analyses were conducted using a linear regression model in plink v1.90 and a single-locus mixed model in snp and variation suite . Genome-wide significance was determined by a Bonferroni correction for multiple testing. Using linear regression, loci on chromosomes 2, 3, 16, 23 and 24 were significantly associated at the genome level with FEC estimated breeding values, and we identified a region on chromosome 2 that was significant using both statistical analyses. We suggest a potential role for the gene DIS3L2 for gastrointestinal nematode resistance in Katahdin sheep, although further research is needed to validate these findings.     

Genome‐wide analysis of polymorphisms identified domestication‐associated long low‐diversity region carrying important rice grain size/weight quantitative trait loci

Rice grain size and weight are major determinants of grain quality and yield and so have been under rigorous selection since domestication. However, the genetic basis for contrasting grain size/weight trait among Indian germplasms and their association with domestication‐driven evolution is not well understood. In this study, two long (LGG) and two short grain (SGG) genotypes were resequenced. LGG (LGR and PB 1121) differentiated from SGG (Sonasal and Bindli) by 504 439 single nucleotide polymorphisms (SNPs) and 78 166 insertion‐and‐deletion polymorphisms. The LRK gene cluster was different and a truncation mutation in the LRK8 kinase domain was associated with LGG. Phylogeny with 3000 diverse rice accessions revealed that the four sequenced genotypes belonged to the japonica group and were at the edge of the clades indicating them to be the potential source of genetic diversity available in Indian rice germplasm. Six SNPs were significantly associated with grain size/weight and the top four of these could be validated in mapping a population, suggesting this study as a valuable resource for high‐throughput genotyping. A contiguous long low‐diversity region (LDR) of approximately 6 Mb carrying a major grain weight quantitative trait loci (harbouring OsTOR gene) was identified on Chromosome 5. This LDR was identified as an evolutionary important site with significant positive selection and multiple selection sweeps, and showed association with many domestication‐related traits, including grain size/weight. The aus population retained more allelic variations in the LDR than the japonica and indica populations, suggesting it to be one of the divergence loci. All the data and analyses can be accessed from the RiceSzWtBase database.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin,USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine‐scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine‐scale adaptation.     

Evaluation of polymorphisms in microRNA‐binding sites and pancreatic cancer risk in Chinese population

As promising biomarkers and therapy targets, microRNAs (miRNAs) are involved in various physiological and tumorigenic processes. Genetic variants in miRNA‐binding sites can lead to dysfunction of miRNAs and contribute to disease. However, systematic investigation of the miRNA‐related single nucleotide polymorphisms (SNPs) for pancreatic cancer (PC) risk remains elusive. We performed integrative bioinformatics analyses to select 31 SNPs located in miRNA‐target binding sites using the miRNASNP v2.0, a solid database providing miRNA‐related SNPs for genetic research, and investigated their associations with risk of PC in two large case‐control studies totally including 1847 cases and 5713 controls. We observed that the SNP rs3802266 is significantly associated with increased risk of PC (odds ratio (OR) = 1.21, 95% confidence intervals (CI) = 1.11‐1.31, P = 1.29E‐05). Following luciferase reporter gene assays show that rs3802266‐G creates a stronger binding site for miR‐181a‐2‐3p in 3′ untranslated region (3′UTR) of the gene ZHX2. Expression quantitative trait loci (eQTL) analysis suggests that ZHX2 expression is lower in individuals carrying rs3802266‐G with increased PC risk. In conclusion, our findings highlight the involvement of miRNA‐binding SNPs in PC susceptibility and provide new clues for PC carcinogenesis.     

Evidence for the effect of serotonin receptor 1A gene (HTR1A) polymorphism on tractability in Thoroughbred horses

Tractability, or how easily animals can be trained and controlled, is an important behavioural trait for the management and training of domestic animals, but its genetic basis remains unclear. Polymorphisms in the serotonin receptor 1A gene (HTR1A) have been associated with individual variability in anxiety‐related traits in several species. In this study, we examined the association between HTR1A polymorphisms and tractability in Thoroughbred horses. We assessed the tractability of 167 one‐year‐old horses reared at a training centre for racehorses using a questionnaire consisting of 17 items. A principal components analysis of answers contracted the data to five principal component (PC) scores. We genotyped two non‐synonymous single nucleotide polymorphisms (SNPs) in the horse HTR1A coding region. We found that one of the two SNPs, c.709G>A, which causes an amino acid change at the intracellular region of the receptor, was significantly associated with scores of four of five PCs in fillies (all Ps < 0.05) and one PC in colts (P < 0.01). Horses carrying an A allele at c.709G>A showed lower tractability. This result provides the first evidence that a polymorphism in a serotonin‐related gene may affect tractability in horses with the effect partially different depending on sex.     

A Genome-Wide Association Study Reveals Dominance Effects on Number of Teats in Pigs

Dominance has been suggested as one of the genetic mechanisms explaining heterosis. However, using traditional quantitative genetic methods it is difficult to obtain accurate estimates of dominance effects. With the availability of dense SNP (Single Nucleotide Polymorphism) panels, we now have new opportunities for the detection and use of dominance at individual loci. Thus, the aim of this study was to detect additive and dominance effects on number of teats (NT), specifically to investigate the importance of dominance in a Landrace-based population of pigs. In total, 1,550 animals, genotyped for 32,911 SNPs, were used in single SNP analysis. SNPs with a significant genetic effect were tested for their mode of gene action being additive, dominant or a combination. In total, 21 SNPs were associated with NT, located in three regions with additive (SSC6, 7 and 12) and one region with dominant effects (SSC4). Estimates of additive effects ranged from 0.24 to 0.29 teats. The dominance effect of the QTL located on SSC4 was negative (−0.26 teats). The additive variance of the four QTLs together explained 7.37% of the total phenotypic variance. The dominance variance of the four QTLs together explained 1.82% of the total phenotypic variance, which corresponds to one-fourth of the variance explained by additive effects. The results suggest that dominance effects play a relevant role in the genetic architecture of NT. The QTL region on SSC7 contains the most promising candidate gene: VRTN. This gene has been suggested to be related to the number of vertebrae, a trait correlated with NT.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

A multi-marker test based on family data in genome-wide association study

Background

Requirements for successful implementation of multivariate animal threshold models including phenotypic and genotypic information are not known yet. Here simulated horse data were used to investigate the properties of multivariate estimators of genetic parameters for categorical, continuous and molecular genetic data in the context of important radiological health traits using mixed linear-threshold animal models via Gibbs sampling. The simulated pedigree comprised 7 generations and 40000 animals per generation. Additive genetic values, residuals and fixed effects for one continuous trait and liabilities of four binary traits were simulated, resembling situations encountered in the Warmblood horse. Quantitative trait locus (QTL) effects and genetic marker information were simulated for one of the liabilities. Different scenarios with respect to recombination rate between genetic markers and QTL and polymorphism information content of genetic markers were studied. For each scenario ten replicates were sampled from the simulated population, and within each replicate six different datasets differing in number and distribution of animals with trait records and availability of genetic marker information were generated. (Co)Variance components were estimated using a Bayesian mixed linear-threshold animal model via Gibbs sampling. Residual variances were fixed to zero and a proper prior was used for the genetic covariance matrix.

Results

Effective sample sizes (ESS) and biases of genetic parameters differed significantly between datasets. Bias of heritability estimates was -6% to +6% for the continuous trait, -6% to +10% for the binary traits of moderate heritability, and -21% to +25% for the binary traits of low heritability. Additive genetic correlations were mostly underestimated between the continuous trait and binary traits of low heritability, under- or overestimated between the continuous trait and binary traits of moderate heritability, and overestimated between two binary traits. Use of trait information on two subsequent generations of animals increased ESS and reduced bias of parameter estimates more than mere increase of the number of informative animals from one generation. Consideration of genotype information as a fixed effect in the model resulted in overestimation of polygenic heritability of the QTL trait, but increased accuracy of estimated additive genetic correlations of the QTL trait.

Conclusion

Combined use of phenotype and genotype information on parents and offspring will help to identify agonistic and antagonistic genetic correlations between traits of interests, facilitating design of effective multiple trait selection schemes.     

A genome‐wide association study reveals a quantitative trait locus for age at first calving in delta/notch‐like EGF repeat containing on chromosome 2 in Japanese Black cattle

Age at first calving (AFC) is an important trait for achieving earlier reproductive performance in cattle. To identify quantitative trait loci for AFC in Japanese Black cattle, we conducted a genome‐wide association study using 866 animals with extreme AFC values selected from a larger group of 52 009 animals. We identified single nucleotide polymorphisms (SNPs) on bovine chromosome 2 that were associated with AFC. These SNPs were located within 112.8‐kbp intronic region of delta/notch‐like EGF repeat containing (DNER) and proved to be in a state of high linkage disequilibrium. The association was replicated in an independent sample of 2963 animals. In the replicated population, the frequency of the reduced AFC allele (Q) was 0.463, and the allele accounts for 8% of the total genetic variance. The effect of allele substitution on AFC was a decrease of 11.54 days. The results suggest that the Q allele could serve as a useful marker in Japanese Black cattle to select animals with superior AFC performance.     

Estimating quantitative genetic parameters in wild populations: a comparison of pedigree and genomic approaches

The estimation of quantitative genetic parameters in wild populations is generally limited by the accuracy and completeness of the available pedigree information. Using relatedness at genomewide markers can potentially remove this limitation and lead to less biased and more precise estimates. We estimated heritability, maternal genetic effects and genetic correlations for body size traits in an unmanaged long‐term study population of Soay sheep on St Kilda using three increasingly complete and accurate estimates of relatedness: (i) Pedigree 1, using observation‐derived maternal links and microsatellite‐derived paternal links; (ii) Pedigree 2, using SNP‐derived assignment of both maternity and paternity; and (iii) whole‐genome relatedness at 37 037 autosomal SNPs. In initial analyses, heritability estimates were strikingly similar for all three methods, while standard errors were systematically lower in analyses based on Pedigree 2 and genomic relatedness. Genetic correlations were generally strong, differed little between the three estimates of relatedness and the standard errors declined only very slightly with improved relatedness information. When partitioning maternal effects into separate genetic and environmental components, maternal genetic effects found in juvenile traits increased substantially across the three relatedness estimates. Heritability declined compared to parallel models where only a maternal environment effect was fitted, suggesting that maternal genetic effects are confounded with direct genetic effects and that more accurate estimates of relatedness were better able to separate maternal genetic effects from direct genetic effects. We found that the heritability captured by SNP markers asymptoted at about half the SNPs available, suggesting that denser marker panels are not necessarily required for precise and unbiased heritability estimates. Finally, we present guidelines for the use of genomic relatedness in future quantitative genetics studies in natural populations.     

A genome‐wide association study to detect genetic variation for postpartum dysgalactia syndrome in five commercial pig breeding lines

Postpartum dysgalactia syndrome (PDS) in sows is an important disease after parturition with a relevant economic impact, affecting the health and welfare of both sows and piglets. The genetic background of this disease has been discussed and its heritability estimated, but further genetic analyses are lacking in detail. The aim of the current study was to detect loci affecting the susceptibility to PDS through a genome‐wide association approach. The study was designed as a family‐based association study with matched sampling of affected sows and healthy half‐ or full‐sib control sows on six farms. For the study, 597 sows (322 affected vs. 275 healthy control sows) were genotyped on 62 163 single nucleotide polymorphisms (SNPs) using the Illumina PorcineSNP60 BeadChip. After quality control, 585 sows (314 affected vs. 271 healthy control sows) and 49 740 SNPs remained for further analysis. Statistics were performed mainly with the r package genabel and included a principal component analysis. A statistically significant genome‐wide associated SNP was identified on porcine chromosome (SSC) 17. Further promising results with moderate significance were detected on SSC 13 and on an unplaced scaffold with an older annotation on SSC 15. The PRICKLE2 and NRP2 genes were identified as candidate genes near associated SNPs. Several quantitative trait loci (QTL) have been previously described in these genomic regions, including QTL for mammary gland condition, as teat number and non‐functional nipples QTL, as well as QTL for body temperature and gestation length.     

QTL and association analysis for skin and fibre pigmentation in sheep provides evidence of a major causative mutation and epistatic effects

The pursuits of white features and white fleeces free of pigmented fibre have been important selection objectives for many sheep breeds. The cause and inheritance of non‐white colour patterns in sheep has been studied since the early 19th century. Discovery of genetic causes, especially those which predispose pigmentation in white sheep, may lead to more accurate selection tools for improved apparel wool. This article describes an extended QTL study for 13 skin and fibre pigmentation traits in sheep. A total of 19 highly significant, 10 significant and seven suggestive QTL were identified in a QTL mapping experiment using an Awassi × Merino × Merino backcross sheep population. All QTL on chromosome 2 exceeded a LOD score of greater than 4 (range 4.4–30.1), giving very strong support for a major gene for pigmentation on this chromosome. Evidence of epistatic interactions was found for QTL for four traits on chromosomes 2 and 19. The ovine TYRP1 gene on OAR 2 was sequenced as a strong positional candidate gene. A highly significant association (P < 0.01) of grandparental haplotypes across nine segregating SNP/microsatellite markers including one non‐synonymous SNP with pigmentation traits could be shown. Up to 47% of the observed variation in pigmentation was accounted for by models using TYRP1 haplotypes and 83% for models with interactions between two QTL probabilities, offering scope for marker‐assisted selection for these traits.