Single nucleotide polymorphisms concordant with the horned/polled trait in Holsteins
© Cargill et al; licensee BioMed Central Ltd. 2008
Received: 27 August 2008
Accepted: 08 December 2008
Published: 08 December 2008
Cattle that naturally do not grow horns are referred to as polled, a trait inherited in a dominant Mendelian fashion. Previous studies have localized the polled mutation (which is unknown) to the proximal end of bovine chromosome 1 in a region approximately 3 Mb in size. While a polled genetic test, Tru-Polled™, is commercially available from MetaMorphix Inc., Holsteins are not a validated breed for this test.
Approximately 160 kb were sequenced within the known polled region from 12 polled and 12 horned Holsteins. Analysis of the polymorphisms identified 13 novel single nucleotide polymorphisms (SNPs) that are concordant with the horned/polled trait. Three of the 13 SNPs are located in gene coding or regulatory regions (e.g., the untranslated region, or UTR) where one is located in the 3'UTR of a gene and the other two are located in the 5'UTR and coding region (synonymous SNP) of another gene. The 3'UTR of genes have been shown to be targets of microRNAs regulating gene expression. In silico analysis indicates the 3'UTR SNP may disrupt a microRNA target site.
These 13 novel SNPs concordant with the horned/polled trait in Holsteins represent a test panel for the breed and this is the first report to the authors' knowledge of SNPs within gene coding or regulatory regions concordant with the horned/polled trait in cattle. These SNPs will require further testing for verification and further study to determine if the 3'UTR SNP may have a functional effect on the polled trait in Holsteins.
Cattle that naturally do not grow horns are termed polled, a trait inherited in an autosomal dominant fashion [1, 2]. De-horning is a common practice in the cattle industry as the presence of horns can lead to injuries such as bruised carcasses and hence, economic loss. Polled cattle are desirable; however, the frequency of the trait is minimal due to the management practice of de-horning calves, which prohibits the later selection and breeding of naturally polled individuals. While de-horning is a management solution, the issue ranks as a high concern with producers and packers [3, 4]. In addition, the process of de-horning creates stress for the cattle  and may be viewed as inhumane.
While a polled genetic test is commercially available from MetaMorphix Inc., Holsteins are not listed as a validated breed for the Tru-Polled™ test . The breeds the Tru-Polled™ test is validated for are Charolais, Gelbvieh, Hereford, Limousin, Salers, and Simmental . Creation of a polled genetic test for Holsteins, the major dairy breed, would be valuable to the dairy industry for inclusion of this trait in selection programs utilizing genetic markers.
The polled mutation in Bos taurus, which is unknown, was localized to the proximal end of bovine chromosome 1 (BTA01) with microsatellite markers . More recent efforts to fine-map the polled locus have included additional microsatellite marker and gene mapping [8–11] and the creation of a BAC-based physical map of the polled region . The location of the most proximal gene, ATP5O, and most distal gene, KRTAP8, of the polled region from these cited sources corresponds to approximately 0.6 Mb and 3.9 Mb respectively on the public bovine genome assembly version 4.0 . One study  did fine map the polled region to a 1 Mb segment that corresponds to approximately 0.6 Mb to 1.6 Mb from the proximal end of BTA01.
Genes targeted for polymorphism detection within the polled region on BTA01.1
interferon (alpha, beta and omega) receptor 2
oligodendrocyte transcription factor 1
oligodendrocyte transcription factor 2
hypothetical protein LOC784884
chromosome 21 open reading frame 59
chromosome 21 open reading frame 63
melanocortin 2 receptor accessory protein
superoxide dismutase 1, soluble
keratin associated protein 11-1
keratin associated protein 8-1
Polymorphisms concordant with the polled trait
SNPs concordant with the polled trait in the Holstein panel.1
Genotypes of the 24 bulls in the Holstein panel for SNP bSYNJ1_C3981T.1
PCR primers, forward (F) and reverse (R) listed 5' to 3', PCR product size in base pairs, and base pair position within the PCR product for the 13 SNPs concordant with the polled trait in Holsteins.
Of the 13 SNPs identified as concordant with the polled trait (Table 2), only three are located in a gene's coding or regulatory (i.e., promoter or UTR) region. The other 10 SNPs are located in introns or are inter-genic and not present in a putative regulatory element as analyzed by WWW Promoter Scan . The bSYNJ1_C3981T polymorphism is located in the 3'UTR of the SYNJ1 gene, and SNPs bC2159_C-193T and bC2159_T372C are located in the 5'UTR and coding region of the C21orf59 gene, respectively. The bC2159_T372C polymorphism is a synonymous SNP, with both alleles coding for a serine amino acid.
MicroRNA target detection
As previously mentioned, the polled concordant SNP bSYNJ1_C3981T is located in the 3'UTR of the gene SYNJ1, and 3'UTRs have been shown to be targets of a gene expression regulatory system orchestrated by microRNAs (e.g., microRNA regulation of myostatin gene expression impacts muscularity in sheep) . In order to examine if the 3'UTR of SYNJ1, and specifically the location of SNP bSYNJ1_C3981T, may be a microRNA target, the current collection of bovine microRNA mature sequences (N = 117) was downloaded from the publicly available miRBase . The microRNA sequences were used in target prediction analysis with the on-line resource RNAhybrid  using the parameters described (see Materials and Methods).
MicroRNAs predicted to target the 3'UTR of SYNJ1 at the location of SNP bSYNJ1_C3981T.
Mature Sequence (5'-3')
A total of 13 novel SNPs concordant with the polled trait in Holsteins (Table 2) were identified. Because of the relatively small sample size and partial family structure (Figure 1) no association analysis was performed. These SNPs will require further testing in an unrelated set of Holsteins for verification. It is also not unexpected the Holstein breed is not validated for the commercial Tru-Polled™ test  as the polled trait has various historical origins in different breeds. One hypothesis is that when the causal functional DNA element (i.e., a gene or regulatory region) for polled is identified there may be multiple breed-specific mutations for the trait in that DNA element, if the trait arose separately in more than one breed. Another hypothesis is the trait was selectively introgressed from the first breed exhibiting polled and the causal mutation will be found to be the same in all breeds. The 13 SNPs reported here do represent the first described genetic markers for the polled trait specifically in Holsteins. In addition, these 13 SNPs are novel as they are not found in a search of public databases, such as NCBI's dbSNP  or the Bovine Genome Project's SNP database .
While discovery of polymorphisms concordant with the polled trait creates utility as a genetic test, the ultimate objective is characterization of a polymorphism with a probable functional effect for the polled trait. SNP bSYNJ1_C3981T is located in the 3'UTR of the SYNJ1 gene. The 3'UTR of genes has been shown to potentially be a target for microRNA regulation of gene expression by post-transcriptionally base-pairing to mRNAs . In silico analysis of the SYNJ1 3'UTR allelic variant created by the bSYNJ1_C3981T SNP revealed eight predicted microRNA target interactions from the 117 available bovine microRNAs (Table 5). Over 500 human microRNAs have been characterized to date  indicating many more bovine microRNAs likely exist as well.
Six of the eight microRNAs (Table 5) are members of a microRNA family, bta-let-7. The let-7 microRNA family has been found to function in late development timing in C. elegans, and one of multiple targets of the let-7 microRNA in C. elegans is a nuclear hormone receptor daf-12 in the seam cells of the hypodermis . The hypodermis is the lowermost layer of the integumentary system, and the integumentary system includes stratified squamous and keratinized epithelium. At a fundamental level, horns are epithelial with a bony core, living tissue, cornified, unbranched, permanent, cannot regenerate, and develop through basal growth .
Based on the in silico analysis presented in this report, it is possible to hypothesize the SYNJ1 3'UTR SNP (bSYNJ1_C3981T) alters a microRNA target site. The effect could be predicted to be tissue specific in such a manner that only horn growth is affected (as the mere presence or absence of horns alone has not been correlated with any other defect to the best of the authors' knowledge) by an unknown mechanism of SYNJ1 function. If the SYNJ1 3'UTR SNP does have a functional effect on the polled trait in Holsteins it, and the SYNJ1 gene itself, becomes a candidate for further investigations of potential causal mutations for the polled trait in other breeds.
In summary, the polymorphisms reported here as concordant with the polled trait in Holsteins can readily be used as a genetic test for this breed. This is the first report, to the authors' knowledge, of SNPs within gene coding or regulatory regions (i.e., not introns or inter-genic) predictive of the horned/polled trait in cattle as previous reports localized and fine-mapped the polled region utilizing inter-genic genetic markers such as microsatellites [7–11]. These SNPs will require further testing for verification and further study to determine if the SYNJ1 3'UTR SNP may have a functional effect on the polled trait in Holsteins.
Animals and DNA samples
Twenty-four Holstein bulls were utilized as a polymorphism detection panel (Figure 1). The majority of the 24 Holsteins are directly related as horned sires and polled sons, with the dams expected to be polled (Figure 1). Semen samples from the bulls were purchased on the open market. In addition to de-horning, the scurs phenotype [2, 23] may complicate identification of truly polled animals. Scurs was not investigated in this study as data was not available. The 12 polled Holstein bulls are registered as polled. All DNA samples were extracted from spermatozoa using the Qiagen Biosprint (Qiagen Inc., Valencia, CA) according to the manufacturer's protocol.
All primers for PCR were designed using Primer3 . Optimal primer annealing temperatures were obtained by using gradient PCR thermocycling conditions of 15 minutes at 95°C, 35 cycles of 45 seconds at 94°C, 45 seconds of gradient temperatures starting at 55° to 66° across twelve sample wells, 45 seconds at 72°C, and 10 minutes 72°C. Once an optimal annealing temperature was found, standard thermocycling conditions were used: 15 minutes at 95°C followed by 35 cycles of 45 seconds at 94°C, 45 seconds at optimal annealing temp, and 45 seconds at 72°C, with a final extension step of 10 minutes at 72°C. Concentrations for a 10 μl PCR volume (gradient and standard) were 5 ng/μl of template DNA, 0.5 μM of each primer (forward and reverse), 1× SIGMA JumpStart PCR Mix (Sigma-Aldrich Co., St. Louis, MO), and 1× combinatorial enhancer solution (CES) .
Putative regulatory element prediction
The on-line resource WWW Promoter Scan  was used to scan targeted gene introns and inter-genic sequences for predicted regulatory elements such as promoters and transcription factor binding sites.
Sequencing and analysis
PCR products were purified using the EXO-SAP-IT PCR Product Clean-up (USB corporation, Cleveland, OH) according to the manufacturer's protocol. Direct sequencing of purified PCR products, in both forward and reverse directions, was conducted using ABI BigDye (Applied Biosystems, Foster City, CA) according to the manufacturer's protocol and resolved on an ABI 3730xl Automated Sequencer (Applied Biosystems, Foster City, CA). Sequence trace alignment and polymorphism detection were carried out using recent versions of Phred/Phrap [26, 27] and Consed .
MicroRNA target prediction
The current collection of bovine microRNA mature sequences (N = 117) was downloaded from the publicly available miRBase . The microRNA sequences were used in target prediction analysis with the on-line resource RNAhybrid . Default target prediction parameters were used, specifically setting 1 hit per target and an energy cut-off of -14. A microRNA was considered to target the 3'UTR of a gene if the specified parameters were met and the seed sequence of the microRNA, base positions 2–8 from the 5' end of the microRNA , contained no more than one gap or mismatch with the 3'UTR sequence being tested.
The authors wish to thank Andrew Brown, Lisa Kemp, David McCarter, and Julie Oermann for technical assistance and Fernanda Rodriguez for microRNA analysis consultation.
- White WT, Ibsen HL: Horn inheritance in Galloway-Holstein cattle crosses. J Gen. 1936, 32: 33-49.View ArticleGoogle Scholar
- Long CR, Gregory KE: Inheritance of the horned scurred, and polled condition in cattle. J Hered. 1978, 69: 395-400.Google Scholar
- Roeber DL, Mies PD, Smith CD, Belk KE, Field TG, Tatum JD, Scanga JA, Smith GC: National market-cow and bull beef quality audit 1999: a survey of producer-related defects in market cow and bulls. J Ani Sci. 2001, 79: 658-665.Google Scholar
- McKenna DR, Roeber DL, Bates PK, Schmidt TB, Hale DS, Griffin DB, Savell JW, Brooks JC, Morgan JB, Montgomery TH, Belk KE, Smith GC: National beef quality audit 2000: survey of targeted cattle and carcass characteristics related to quality, quantity, and value of fed steers and heifers. J Ani Sci. 2002, 80: 1212-1222.Google Scholar
- Graf B, Senn M: Behavioural and physiological responses of calves to dehorning by heat cauterization with or without local anaesthesia. App Ani Behav Sci. 1999, 62: 153-171.View ArticleGoogle Scholar
- MetaMorphix Inc. [http://www.metamorphixinc.com/faqtrupolled.html]
- 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 Gen. 1993, 4: 206-210.View ArticleGoogle Scholar
- Schmutz SM, Marquess FLS, Berryere TG, Moker JS: DNA marker assisted selection of the polled condition in Charolais cattle. Mamm Gen. 1995, 6: 710-713.View ArticleGoogle Scholar
- Brenneman RA, Davis SK, Sander JO, Burns BM, Wheeler TC, Turner JW, Taylor JF: The polled locus maps to BTA1 in a Bos indicus × Bos taurus cross. J Hered. 1996, 87: 156-161.View ArticlePubMedGoogle Scholar
- Harlizius B, Tammen I, Eichler K, Eggen A, Hetzel DJ: New markers on bovine chromosome 1 are closely linked to the polled gene in Simmental and Pinzgauer cattle. Mamm Gen. 1997, 8: 255-257.View ArticleGoogle Scholar
- Drögemüller 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 Gen. 2005, 16: 613-620.View ArticleGoogle Scholar
- Wunderlich KR, Abbey CA, Clayton DR, Song Y, Schein JE, Georges M, Coppieters W, Adelson D, Taylor JF, Davis SL, Gill CA: A 2.5-Mb contig constructed from Angus, Longhorn and horned Hereford DNA spanning the polled interval on bovine chromosome 1. Anim Gen. 2006, 37: 592-594.View ArticleGoogle Scholar
- Bovine Genome Project. [http://www.hgsc.bcm.tmc.edu/projects/bovine/]
- WWW Promoter Scan. [http://www-bimas.cit.nih.gov/molbio/proscan/]
- Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibé B, Bouix J, Caiment F, Elsen JM, Eychenne F, Larzul C, Laville E, Meish F, Milenkovic D, Tobin J, Charlier C, Georges M: A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat Gen. 2006, 38: 813-818.View ArticleGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nuc Acids Res. 2006, D140-D144.miRBase, [http://microrna.sanger.ac.uk/index.shtml]34 Database,
- Rehmsmeier M, Steffen P, Höchsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA. 2004, 10: 1507-1517.RNAhybrid, [http://bibiserv.techfak.uni-bielefeld.de/rnahybrid/]PubMed CentralView ArticlePubMedGoogle Scholar
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798.View ArticlePubMedGoogle Scholar
- National Center for Biotechnology Information (NCBI). [http://www.ncbi.nlm.nih.gov/]
- Wienholds E, Plasterk RHA: MicroRNA function in animal development. FEBS Letters. 2005, 579: 5911-5922.View ArticlePubMedGoogle Scholar
- Esquela-Kerscher A, Johnson SM, Bai L, Saito K, Partridge J, Reinert KL, Slack FJ: Post-embryonic expression of C. elegans microRNAs belonging to the lin-4 and let-7 families in the hypodermis and the reproductive system. Dev Dyn. 2005, 234: 868-877.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall BK: Horns and Ossicones. Bones and Cartilage. 2005, Elsevier Academic Press, San Diego, CA, 95-101.View ArticleGoogle Scholar
- Scurs phenotype listing in the Online Mendelian Inheritance in Animals database. [http://www.ncbi.nlm.nih.gov/sites/entrez?db=omia&cmd=search&term=scurs]
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Meth Molec Bio. 2000, 132: 365-86.Google Scholar
- Ralser M, Querfurth R, Warnatz HJ, Lehrach H, Yaspo ML, Krobitsch S: An efficient and economic enhancer mix for PCR. Biochem Bioph Res Comm. 2006, 347: 747-751.View ArticleGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.View ArticlePubMedGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.View ArticlePubMedGoogle Scholar
- Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202.View ArticlePubMedGoogle Scholar
- Gu Z, Eleswarapu S, Jiang H: Identification and characterization of microRNAs from the bovine adipose tissue and mammary gland. FEBS Letters. 2007, 581: 981-988.View ArticlePubMedGoogle Scholar
- Coutinho LL, Matukumalli LK, Sonstegard TS, Van Tassell CP, Gasbarre LC, Capuco AV, Smith TP: Discovery and profiling of bovine microRNAs from immune-related and embryonic tissues. Phys Gen. 2007, 29: 35-43.View ArticleGoogle 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 cited.