The 80-kb DNA duplication on BTA1 is the only remaining candidate mutation for the polled phenotype of Friesian origin
© Rothammer et al.; licensee BioMed Central Ltd. 2014
Received: 30 December 2013
Accepted: 17 May 2014
Published: 3 July 2014
The absence of horns, called polled phenotype, is the favored trait in modern cattle husbandry. To date, polled cattle are obtained primarily by dehorning calves. Dehorning is a practice that raises animal welfare issues, which can be addressed by selecting for genetically hornless cattle. In the past 20 years, there have been many studies worldwide to identify unique genetic markers in complete association with the polled trait in cattle and recently, two different alleles at the POLLED locus, both resulting in the absence of horns, were reported: (1) the Celtic allele, which is responsible for the polled phenotype in most breeds and for which a single candidate mutation was detected and (2) the Friesian allele, which is responsible for the polled phenotype predominantly in the Holstein-Friesian breed and in a few other breeds, but for which five candidate mutations were identified in a 260-kb haplotype. Further studies based on genome-wide sequencing and high-density SNP (single nucleotide polymorphism) genotyping confirmed the existence of the Celtic and Friesian variants and narrowed down the causal Friesian haplotype to an interval of 145 kb.
Almost 6000 animals were genetically tested for the polled trait and we detected a recombinant animal which enabled us to reduce the Friesian POLLED haplotype to a single causal mutation, namely a 80-kb duplication. Moreover, our results clearly disagree with the recently reported perfect co-segregation of the POLLED mutation and a SNP at position 1 390 292 bp on bovine chromosome 1 in the Holstein-Friesian population.
We conclude that the 80-kb duplication, as the only remaining variant within the shortened Friesian haplotype, represents the most likely causal mutation for the polled phenotype of Friesian origin.
The bovine polled phenotype, i.e. the absence of horns, has huge practical importance for breeders and is of special biological interest to geneticists. Currently, the world cattle population is estimated to be about 1.3 billion heads  of which a large proportion is horned. Although historical records indicate the presence of naturally polled cattle as far back as ancient Egypt [2, 3], until recently, horned cattle remained the desired phenotype because it simplified tethering and attachment to harnesses. However, in modern husbandry systems, such practices are no longer used and the presence of horns increases the risk of injuries to handlers and animals, especially as the housing densities have increased. Thus, the presence of horns induces considerable economic losses in the cattle industry due to dehorning practices and the treatment of subsequent secondary infections, but also due to carcass and leather deterioration [2, 4–7]. Today, dehorning of cattle at an early age is common and although records are not exhaustive, there are probably more than 100 million calves that are dehorned each year. Since all dehorning methods are invasive and thus raise animal welfare issues, the possibility of breeding genetically polled cattle is a promising alternative [4, 8]. In addition, understanding the molecular mechanisms that are involved in the inhibition of the bovine horn bud differentiation will contribute knowledge on the mechanisms that underlie ectopic expression in mammals .
Since 1906, the polled phenotype is known to be inherited as an autosomal dominant trait  and the POLLED locus was mapped to bovine chromosome 1 (BTA1 for Bos taurus) in 1993 . In the past years, the position of the POLLED locus on BTA1 was refined [11, 12] and candidate causal mutations were identified [2, 4, 5]. Recently, Medugorac et al.  demonstrated the existence of at least two different alleles at the POLLED locus in cattle. While a complex 202-bp insertion-deletion (InDel), referred to as P202ID, was identified in various cattle breeds from Scandinavia, Scotland, England, the Channel Islands and France down to the Alpine region (thus also referred to as the Celtic allele PC), a 260-kb haplotype including five candidate mutations (three SNP: PG1654405A, PC1655463T, PC1768587A; two InDel: P5ID, P80kbID) was reported to be perfectly associated with the polled phenotype in the Holstein-Friesian breed (thus also referred to as the Friesian allele PF) . Further studies based on genome-wide sequencing and high-density SNP genotyping confirmed the existence of these Celtic and Friesian variants and narrowed the causal Friesian haplotype down to an interval of 145 kb  that included only two of the five initial candidate mutations i.e. SNP PC1768587A and InDel P80kbID. However, high-throughput sequencing of this 145-kb interval identified two new candidate variants i.e. SNPs PT1764239C and PG1855898A, and genotyping of an enlarged cattle breed panel revealed the presence of the Celtic allele in new breeds from Iceland, Scandinavia, UK, Ireland and France . In 2013, Allais-Bonnet et al.  provided evidence that introgression of the Friesian allele was ongoing in some French and German cattle breeds for which the Celtic allele was predominant and Glatzer et al.  reported that SNP AC000158:g.1390292G > A, hereinafter referred to as rs134968598, co-segregated perfectly with the POLLED locus in a set of 443 Holstein-Friesian animals. This SNP is located within intron 3 of the interferon gamma receptor 2 (IFNGR2) gene and is close (1113 bp) to a SNP for which the most significant selection signature was observed in a group of polled cattle . However, these two SNPs are located 378 kb from the above mentioned Friesian haplotype.
This study aimed at (1) refining the 260-kb haplotype which, based on previous results, is associated with the polled phenotype of Friesian origin  and (2) examining the co-segregation of SNP rs134968598 and the POLLED locus in Holstein-Friesian cattle. These aims are highly interesting from a practical and economic point of view for the most important global cattle breed Holstein-Friesian. Identification of the causal candidate mutation in the Friesian haplotype will also contribute to decipher the molecular mechanisms that underlie this complex developmental phenotype with allelic heterogeneity.
No formal ethical approval was required, since no new tissue samples were collected for this study. Indeed, for most animals, the DNA samples from previous studies [4, 12] were in stock in our laboratory and the remaining DNA samples were obtained from routine testing (paternity and polled phenotype test) by appropriate cattle breeding organizations.
To test the previously reported  co-segregation of the POLLED locus with SNP rs134968598, which is located outside of the 260-kb haplotype, we used 130 animals that were representative of the global polled and horned Holstein-Friesian cattle population. These 130 animals included animals that we also used to refine the Friesian haplotype and bulls that had already been analyzed in our previous study . Of the 130 animals, 69 were horned (pp), 54 were heterozygous polled (Pp) and seven were homozygous polled (PP). The POLLED genotype of each polled bull based on progeny testing was available from the respective breeding organizations. The progeny test considered phenotypes of at least nine progeny per bull mated to horned dams.
Sequencing of the PT1764239C candidate mutation
In our previous study, the region captured around the SNP PT1764239C, one of the two new candidate SNPs reported by , was sequenced by high-throughput sequencing, but the sequence coverage was low . Therefore, this region was re-sequenced by conventional Sanger sequencing  using DNA samples from four horned animals for which, in the previous study, the sequence capture data suggested a discrepancy between the known phenotype and the expected genotype, although the coverage was too low to infer reliable genotypes.
Refining the Friesian haplotype
In order to refine the 260-kb haplotype that perfectly associates with the POLLED locus of Friesian origin, we genotyped 13 Holstein-Friesian individuals (that formed a small pedigree with two recombinant bulls and 11 relatives) for the Celtic allele P202ID and the four Friesian variants P5ID, PG1654405A, PC1768587A and P80kbID as described in Medugorac et al. . Furthermore, of the two additional Friesian SNPs PG1855898A and PT1764239C reported in Allais-Bonnet et al. , we included only the candidate causal mutation PG1855898A because preliminary Sanger sequencing results indicated that PT1764239C segregated in horned animals. Thus, the 13 Holstein-Friesian individuals were also genotyped for PG1855898A using PCR-RFLP. Therefore, PCR products of the region around PG1855898A were amplified using the primers reported in , digested with the restriction enzyme PstI and the digestion products were size-separated by electrophoresis on a 2% ethidium-bromide stained agarose gel.
Genotyping of SNP rs134968598 to test co-segregation with the POLLED locus
The 130 animals, including those used to refine the haplotype, were genotyped for the recently published mutation within the gene IFNGR2 (rs134968598) as described in Glatzer et al. .
Results and discussion
Genotyping of SNP rs134968598 to test co-segregation with the POLLED locus
Seven bulls for which the progeny-based POLLED genotype and the genotype at SNP rs134968598 did not match
Refining the Friesian haplotype
After our genotyping data excluded the rs134968598 SNP on BTA1 at position 1 390 292 bp as a candidate mutation, we re-analysed our sequence capture data and partly re-sequenced a region around the PT1764239C candidate mutation of the 145-kb haplotype of Friesian origin reported in Allais-Bonnet et al. . Several horned animals were found heterozygous C/T at SNP PT1764239C and thus this SNP was also excluded as a candidate mutation. In order to refine the Friesian haplotype, which contains the causal mutation for the polled trait in most Holstein-Friesian breeding lines, a small pedigree including two recombinant bulls and relatives (Figure 1B) was genotyped for the six known polymorphisms in this region P202ID, P5ID, PG1654405A, PC1768587A, PA1855898G and P80kbID (for physical positions (according to UMD3.1) see Figure 1A).
ELBOY RED, a horned male offspring of LAWN BOY P-RED, is a recombinant between PG1654405A and PC1768587A, which excludes P5ID and PG1654405A as well as the fairly distant rs134968598 SNP as candidate mutations for the polled phenotype. Previously, we reported and confirmed this recombination in 20 horned offspring of ELBOY RED . A second horned offspring of LAWN BOY P-RED, the recombinant bull ROSSI, carries a large region of the Friesian haplotype that is associated with the polled trait and that comprises P5ID, PG1654405A, PC1768587A and PG1855898A (Figure 1A and B). The genotyping data of ROSSI located the recombination event in an interval of 53 kb between PG1855898A and P80kbID (beginning at g.1909352) and clearly excluded P5ID, PG1654405A, PC1768587A and PG1855898A as candidates for the polled phenotype of Friesian origin.
Finally, based on the combined sequencing and genotyping data from our previous studies [2, 4] and this study, the only remaining causal candidate within the Friesian POLLED haplotype is the 80-kb duplication P80kbID. This conclusion is inconsistent with that of Glatzer et al.  who reported the existence of a progeny-proven heterozygous polled Holstein-Friesian sire that carried neither the P80kbID mutation nor a Celtic mutation. Unfortunately, since the identity of this bull was not provided, we could not verify these results. However, our personal experience from different previous analyses that we carried out is, that reliable genotyping of P202ID and P80kbID is not trivial and that the genotyping-procedure should be validated using anonymous samples. Thus, there is a possibility that the bull reported in  does carry either P202ID or P80kbID. In addition, even if we had demonstrated the absence of both P202ID and P80kbID in the genome of this bull, further analyses would still have been necessary to show that it does not carry another spontaneous mutation, which could be responsible for a dominant horn defect syndrome. Such mutations have been reported to occur and can be confused with the regular POLLED locus [2, 15, 16]. Therefore, the existence of a new candidate mutation at the POLLED locus of this animal should first be investigated and supported by linkage analyses on its relatives before refuting the highly probable causality of P80kbID for the polled phenotype in Holstein-Friesian cattle.
Based on the genotypes of almost 6000 animals, we identified two bulls with a recombination within the 260-kb Friesian haplotype and narrowed this haplotype down to an interval of 53 kb on BTA1. Moreover, the results from these two recombinant bulls and related individuals disagree with the perfect association between the POLLED locus and a SNP within intron 3 of the IFNGR2 gene previously reported by Glatzer et al. . Combining the exhaustive sequencing data and high-density SNP genotyping results from previous studies [2, 4] and those reported here, we provide evidence that supports the P80kbID mutation as the only remaining and most likely causal mutation for the polled phenotype of Friesian origin. In addition to P80kbID, P202ID (the POLLED mutation of Celtic origin) segregates in some Holstein-Friesian families and sporadically combines in some homozygous polled but heterogeneous (PC/PF) animals (Figure 1B). Neither P80kbID nor P202ID are located in any known coding regions [2, 4] but these mutations, even as single copy, result in the inhibition of horn development at a very early stage in cattle [2, 4]. Therefore, future investigations of the molecular mechanisms that underlie the inhibition of horn development in cattle will enrich our knowledge not only on the mechanisms that are involved in bovine horn bud differentiation but, also, more generally, on ectopic expression in mammals.
We thank breeders, breeding organizations and insemination stations, which provided samples for this study. In particular, we thank KH Göpel (Göpel Genetik GmbH, Germany), Dr J Potthast (RUW, Germany), Dr B Weber (Masterrind GmbH, Germany), B Brentrup (WWS Germany GmbH, Germany), C Grohs (INRA, France) and J Thomin (Evolution, France) for sharing samples and phenotypes. We thank the reviewers and editor for their helpful comments, which have greatly improved language, content and clarity of the paper.
- Brown LR: Plan B 4.0: Mobilizing to Save Civilization. 2009, New York: WW Norton & CompanyGoogle Scholar
- Allais-Bonnet A, Grohs C, Medugorac I, Krebs S, Djari A, Graf A, Fritz S, Seichter D, Baur A, Russ I, Bouet S, Rothammer S, Wahlberg P, Esquerre D, Hoze C, Boussaha M, Weiss B, Thepot D, Fouilloux MN, Rossignol MN, van Marle-Koster E, Hreietharsdottir GE, Barbey S, Dozias D, Cobo E, Reverse P, Catros O, Marchand JL, Soulas P, Roy P: Novel insights into the bovine polled phenotype and horn ontogenesis in Bovidae. PLoS ONE. 2013, 8: e63512-10.1371/journal.pone.0063512.PubMed CentralView ArticlePubMedGoogle Scholar
- Roman A: L'élevage bovine en Egypte antique. Bull Soc Fr Hist Méd Sci Vét. 2004, 3: 35-45.Google Scholar
- Medugorac I, Seichter D, Graf A, Russ I, Blum H, Gopel KH, Rothammer S, Förster M, Krebs S: Bovine polledness–an autosomal dominant trait with allelic heterogeneity. PLoS ONE. 2012, 7: e39477-10.1371/journal.pone.0039477.PubMed CentralView ArticlePubMedGoogle Scholar
- Glatzer S, Merten NJ, Dierks C, Wöhlke A, Philipp U, Distl O: A Single nucleotide polymorphism within the Interferon Gamma Receptor 2 gene perfectly coincides with polledness in Holstein cattle. PLoS ONE. 2013, 8: e67992-10.1371/journal.pone.0067992.PubMed CentralView ArticlePubMedGoogle Scholar
- Misch LJ, Duffield TF, Millman ST, Lissemore KD: An investigation into the practices of dairy producers and veterinarians in dehorning dairy calves in Ontario. Can Vet J. 2007, 48: 1249-1254.PubMed CentralPubMedGoogle Scholar
- Prayaga KC: Genetic options to replace dehorning in beef cattle - a review. Aust J Agr Res. 2007, 58: 1-8. 10.1071/AR06044.View ArticleGoogle Scholar
- Graf B, Senn M: Behavioural and physiological responses of calves to dehorning by heat cauterization with or without local anaesthesia. Appl Anim Behav Sci. 1999, 62: 153-171. 10.1016/S0168-1591(98)00218-4.View ArticleGoogle Scholar
- Spillman WJ: Mendel's law in relation to animal breeding. J Hered. 1906, 1: 171-177.Google Scholar
- Georges M, Drinkwater R, King T, Mishra A, Moore SS, Nielsen D, Sargeant LS, Sorensen A, Steele MR, Zhao X, Womack JE, Hetzel J: Microsatellite mapping of a gene affecting horn development in Bos taurus. Nat Genet. 1993, 4: 206-210. 10.1038/ng0693-206.View ArticlePubMedGoogle Scholar
- Drögemuller C, Wöhlke A, Mömke S, Distl O: Fine mapping of the polled locus to a 1-Mb region on bovine chromosome 1q12. Mamm Genome. 2005, 16: 613-620. 10.1007/s00335-005-0016-0.View ArticlePubMedGoogle Scholar
- Seichter D, Russ I, Rothammer S, Eder J, Förster M, Medugorac I: SNP-based association mapping of the polled gene in divergent cattle breeds. Anim Genet. 2012, 43: 595-598. 10.1111/j.1365-2052.2011.02302.x.View ArticlePubMedGoogle Scholar
- Stella A, Ajmone-Marsan P, Lazzari B, Boettcher P: Identification of selection signatures in cattle breeds selected for dairy production. Genetics. 2010, 185: 1451-1461. 10.1534/genetics.110.116111.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanger F, Nicklen S, Coulson AR: DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci USA. 1977, 74: 5463-5467. 10.1073/pnas.74.12.5463.PubMed CentralView ArticlePubMedGoogle Scholar
- Capitan A, Allais-Bonnet A, Pinton A, Marquant-Le Guienne B, Le Bourhis D, Grohs C, Bouet S, Clement L, Salas-Cortes L, Venot E, Chaffaux S, Weiss B, Delpeuch A, Noe G, Rossignol MN, Barbey S, Dozias D, Cobo E, Barasc H, Auguste A, Pannetier M, Deloche MC, Lhuilier E, Bouchez O, Esquerre D, Salin G, Klopp C, Donnadieu C, Chantry-Darmon C, Hayes H: A 3.7 Mb deletion encompassing ZEB2 causes a novel polled and multisystemic syndrome in the progeny of a somatic mosaic bull. PLoS ONE. 2012, 7: e49084-10.1371/journal.pone.0049084.PubMed CentralView ArticlePubMedGoogle Scholar
- Capitan A, Grohs C, Weiss B, Rossignol MN, Reversé P, Eggen A: A newly described bovine type 2 scurs syndrome segregates with a frame-shift mutation in TWIST1. PLoS ONE. 2011, 6: e22242-10.1371/journal.pone.0022242.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.