Gene expression supports a single origin of horns and antlers in hoofed mammals.

Horns, antlers, and other bony cranial appendages of even-toed hoofed mammals (ruminant artiodactyls) challenge traditional morphological homology assessments. Cranial appendages all share a permanent bone portion with family-specific integument coverings, but homology determination depends on whether the integument covering is an essential component or a secondary elaboration of each structure. To enhance morphological homology assessments, we tested whether juvenile cattle horn bud transcriptomes share homologous gene expression patterns with deer antlers relative to pig outgroup tissues, treating the integument covering as a secondary elaboration. We uncovered differentially expressed genes that support horn and antler homology, potentially distinguish them from non-cranial-appendage bone and other tissues, and highlight the importance of phylogenetic outgroups in homology assessments. Furthermore, we found differentially expressed genes that could support a shared cranial neural crest origin for horns and antlers and expression patterns that refine our understanding of the timing of horn and antler differentiation.

Evolution of horn length and lifting strength in the Japanese rhinoceros beetle Trypoxylus dichotomus.

What limits the size of nature's most extreme structures? For weapons like beetle horns, one possibility is a tradeoff associated with mechanical levers: as the output arm of the lever system-the beetle horn-gets longer, it also gets weaker. This "paradox of the weakening combatant" could offset reproductive advantages of additional increases in weapon size. However, in contemporary populations of most heavily weaponed species, males with the longest weapons also tend to be the strongest, presumably because selection drove the evolution of compensatory changes to these lever systems that ameliorated the force reductions of increased weapon size. Therefore, we test for biomechanical limits by reconstructing the stages of weapon evolution, exploring whether initial increases in weapon length first led to reductions in weapon force generation that were later ameliorated through the evolution of mechanisms of mechanical compensation. We describe phylogeographic relationships among populations of a rhinoceros beetle and show that the "pitchfork" shaped head horn likely increased in length independently in the northern and southern radiations of beetles. Both increases in horn length were associated with dramatic reductions to horn lifting strength-compelling evidence for the paradox of the weakening combatant-and these initial reductions to horn strength were later ameliorated in some populations through reductions to horn length or through increases in head height (the input arm for the horn lever system). Our results reveal an exciting geographic mosaic of weapon size, weapon force, and mechanical compensation, shedding light on larger questions pertaining to the evolution of extreme structures.

A deletion at the polled PC locus alone is not sufficient to cause a polled phenotype in cattle.

Dehorning is a common practice in the dairy industry, but raises animal welfare concerns. A naturally occurring genetic mutation (PC allele) comprised of a 212 bp duplicated DNA sequence replacing a 10-bp sequence at the polled locus is associated with the hornless phenotype (polled) in cattle. To test the hypothesis that the 10 bp deletion alone is sufficient to result in polled, a CRISPR-Cas9 dual guide RNA approach was optimized to delete a 133 bp region including the 10 bp sequence. Timing of ribonucleoprotein complex injections at various hours post insemination (hpi) (6, 8, and 18 hpi) as well as in vitro transcribed (IVT) vs synthetic gRNAs were compared. Embryos injected 6 hpi had a significantly higher deletion rate (53%) compared to those injected 8 (12%) and 18 hpi (7%), and synthetic gRNAs had a significantly higher deletion rate (84%) compared to IVT gRNAs (53%). Embryo transfers were performed, and bovine fetuses were harvested between 3 and 5 months of gestation. All fetuses had mutations at the target site, with two of the seven having biallelic deletions, and yet they displayed horn bud development indicating that the 10 bp deletion alone is not sufficient to result in the polled phenotype.

A variance component estimation approach to infer associations between Mendelian polledness and quantitative production and female fertility traits in German Simmental cattle.

BACKGROUND: Managing beneficial Mendelian characteristics in dairy cattle breeding programs implies that the correlated genetic effects are considered to avoid possible adverse effects in selection processes. The Mendelian trait polledness in cattle is traditionally associated with the belief that the polled locus has unfavorable effects on breeding goal traits. This may be due to the inferior breeding values of former polled bulls and cows in cattle breeds, such as German Simmental, or to pleiotropic or linkage effects of the polled locus. METHODS: We focused on a variance component estimation approach that uses a marker-based numerator relationship matrix reflecting gametic relationships at the polled locus to test for direct pleiotropic or linked quantitative trait loci (QTL) effects of the polled locus on relevant traits. We applied the approach to performance, health, and female fertility traits in German Simmental cattle. RESULTS: Our results showed no evidence for any pleiotropic QTL effects of the polled locus on test-day production traits milk yield and fat percentage, on the mastitis indicator 'somatic cell score', and on several female fertility traits, i.e. 56 days non return rate, days open and days to first service. We detected a significant and unfavorable QTL effect accounting for 6.6% of the genetic variance for protein percentage only. CONCLUSIONS: Pleiotropy does not explain the lower breeding values and phenotypic inferiority of polled German Simmental sires and cows relative to the horned population in the breed. Thus, intensified selection in the polled population will contribute to increased selection response in breeding goal traits and genetic merit and will narrow the deficit in breeding values for production traits.

Epithelial folding determines the final shape of beetle horns.

The elaborate ornaments and weapons of sexual selection, such as the vast array of horns observed in scarab beetles, are some of the most striking outcomes of evolution. How these novel traits have arisen, develop, and respond to condition is governed by a complex suite of interactions that require coordination between the environment, whole-animal signals, cell-cell signals, and within-cell signals. Endocrine factors, developmental patterning genes, and sex-specific gene expression have been shown to regulate beetle horn size, shape, and location, yet no overarching mechanism of horn shape has been described. Recent advances in microscopy and computational analyses combined with a functional genetic approach have revealed that patterning genes combined with intricate epithelial folding and movement are responsible for the final shape of a beetle head horn.

A universal power law for modelling the growth and form of teeth, claws, horns, thorns, beaks, and shells.

BACKGROUND: A major goal of evolutionary developmental biology is to discover general models and mechanisms that create the phenotypes of organisms. However, universal models of such fundamental growth and form are rare, presumably due to the limited number of physical laws and biological processes that influence growth. One such model is the logarithmic spiral, which has been purported to explain the growth of biological structures such as teeth, claws, horns, and beaks. However, the logarithmic spiral only describes the path of the structure through space, and cannot generate these shapes. RESULTS: Here we show a new universal model based on a power law between the radius of the structure and its length, which generates a shape called a 'power cone'. We describe the underlying 'power cascade' model that explains the extreme diversity of tooth shapes in vertebrates, including humans, mammoths, sabre-toothed cats, tyrannosaurs and giant megalodon sharks. This model can be used to predict the age of mammals with ever-growing teeth, including elephants and rodents. We view this as the third general model of tooth development, along with the patterning cascade model for cusp number and spacing, and the inhibitory cascade model that predicts relative tooth size. Beyond the dentition, this new model also describes the growth of claws, horns, antlers and beaks of vertebrates, as well as the fangs and shells of invertebrates, and thorns and prickles of plants. CONCLUSIONS: The power cone is generated when the radial power growth rate is unequal to the length power growth rate. The power cascade model operates independently of the logarithmic spiral and is present throughout diverse biological systems. The power cascade provides a mechanistic basis for the generation of these pointed structures across the tree of life.

Computational analyses decipher the primordial folding coding the 3D structure of the beetle horn.

The beetle horn primordium is a complex and compactly folded epithelial sheet located beneath the larval cuticle. Only by unfolding the primordium can the complete 3D shape of the horn appear, suggesting that the morphology of beetle horns is encoded in the primordial folding pattern. To decipher the folding pattern, we developed a method to manipulate the primordial local folding on a computer and clarified the contribution of the folding of each primordium region to transformation. We found that the three major morphological changes (branching of distal tips, proximodistal elongation, and angular change) were caused by the folding of different regions, and that the folding mechanism also differs according to the region. The computational methods we used are applicable to the morphological study of other exoskeletal animals.

Back to the basics: allometric growth by the horns of bovid mammals.

I used the equivalent of nonlinear analysis of covariance (ANCOVA) to re-examine relative growth by the horns on males and females of alpine ibex (Capra ibex) and mouflon sheep (Ovis gmelini). A prior study of allometric growth by the horns on these animals described a pattern of biphasic allometry for both sexes, with two different mathematical equations being required to capture the pattern of variation over the full range in body size. However, the investigation in question used conventional analytical methods based on logarithmic transformations, which alter bivariate distributions and commonly introduce problems with analysis and interpretation. My new analyses of data for both species revealed that untransformed observations for both males and females are monophasic and that they are described quite well by three-parameter power equations with negative intercepts. Equations for males follow a steep upward trajectory whereas those for females follow much shallower paths. The negative intercepts indicate that males and females of both species must attain a minimum body size before horns begin to develop. Conclusions from the earlier investigation were based on inaccurate perceptions of pattern in the data. Future studies should be based on graphical and analytical analysis of observations expressed on the original arithmetic scale.

The effect of early burn injury on sensitivity to future painful stimuli in dairy heifers.

Animals that experience painful procedures as neonates are more sensitive to pain later in life. We evaluated whether disbudding with a heated iron at 3 (n = 12), 35 (n = 9), or 56 (n = 20) d of age affected heifers' pain responses to vaccine injections at 11 mo of age. Heifers responded to the injection procedure with struggling and changes in eye temperature and heart rate variability compared to a sham procedure the day before, and still had a heightened response 6 d later, regardless of disbudding age. However, some heart rate variability indices suggested increased sympathetic dominance in heifers disbudded at 35 d, compared to the other 2 age groups, independent of the injection procedure. We also found that heifers disbudded at 3 or 35 d had a higher mean heart rate after the injection procedure compared to those disbudded at 56 d. We conclude that: (1) heifers find injections aversive; and (2) there is some evidence that disbudding age influences autonomic nervous system activity later in life.

Formation, structure, and function of extra-skeletal bones in mammals.

This review describes the formation, structure, and function of bony compartments in antlers, horns, ossicones, osteoderm and the os penis/os clitoris (collectively referred to herein as AHOOO structures) in extant mammals. AHOOOs are extra-skeletal bones that originate from subcutaneous (dermal) tissues in a wide variety of mammals, and this review elaborates on the co-development of the bone and skin in these structures. During foetal stages, primordial cells for the bony compartments arise in subcutaneous tissues. The epithelial-mesenchymal transition is assumed to play a key role in the differentiation of bone, cartilage, skin and other tissues in AHOOO structures. AHOOO ossification takes place after skeletal bone formation, and may depend on sexual maturity. Skin keratinization occurs in tandem with ossification and may be under the control of androgens. Both endochondral and intramembranous ossification participate in bony compartment formation. There is variation in gradients of density in different AHOOO structures. These gradients, which vary according to function and species, primarily reduce mechanical stress. Anchorage of AHOOOs to their surrounding tissues fortifies these structures and is accomplished by bone-bone fusion and Sharpey fibres. The presence of the integument is essential for the protection and function of the bony compartments. Three major functions can be attributed to AHOOOs: mechanical, visual, and thermoregulatory. This review provides the first extensive comparative description of the skeletal and integumentary systems of AHOOOs in a variety of mammals.

Are scurs in heterozygous polled (Pp) cattle a complex quantitative trait?

BACKGROUND: Breeding genetically hornless, i.e. polled, cattle provides an animal welfare-friendly and non-invasive alternative to the dehorning of calves. However, the molecular regulation of the development of horns in cattle is still poorly understood. Studying genetic characters such as polledness and scurs, can provide valuable insights into this process. Scurs are hornlike formations that occur occasionally in a wide variety of sizes and forms as an unexpected phenotype when breeding polled cattle. METHODS: We present a unique dataset of 885 Holstein-Friesian cattle with polled parentage. The horn phenotype was carefully examined, and the phenotypic heterogeneity of the trait is described. Using a direct gene test for polledness, the polled genotype of the animals was determined. Subsequently, the existence of a putative scurs locus was investigated using high-density genotype data of a selected subset of 232 animals and two mapping approaches: mixed linear model-based association analyses and combined linkage disequilibrium and linkage analysis. RESULTS: The results of an exploratory data analysis indicated that the expression of scurs depends on age at phenotyping, sex and polled genotype. Scurs were more prevalent in males than in females. Moreover, homozygous polled animals did not express any pronounced scurs and we found that the Friesian polled allele suppresses the development of scurs more efficiently than the Celtic polled allele. Combined linkage and linkage disequilibrium mapping revealed four genome-wide significant loci that affect the development of scurs, one on BTA5 and three on BTA12. Moreover, suggestive associations were detected on BTA16, 18 and 23. The mixed linear model-based association analysis supports the results of the combined linkage and linkage disequilibrium analysis. None of the mapping approaches provided convincing evidence for a monogenic inheritance of scurs. CONCLUSIONS: Our results contradict the initial and still broadly accepted model for the inheritance of horns and scurs. We hypothesise an oligogenetic model to explain the development of scurs and polledness.

Understanding the effects of the bovine POLLED variants.

Horns are paired appendages on the head of bovine species, comprising an inner bony core and outer keratin sheath. The horn bud forms during early fetal development but ossification of the developing horn does not occur until approximately 1 month after birth. Little is known about the genetic pathways that lead to horn growth. Hornless, or polled, animals are found in all domestic bovids. Histological studies of bovine fetuses have shown that the horn bud does not form in polled individuals. There are currently four known genetic variants for polledness in cattle on BTA1. All of the variants are intergenic, but probably affect regulation of nearby genes or long non-coding RNAs. Transcriptomic studies suggest that the expression of two nearby long non-coding RNAs are affected by the Celtic POLLED variant, but further studies are required to confirm these data. Candidate genes located elsewhere in the genome are involved in regulating bone formation and epithelial-to-mesenchymal transition. Expression of one of these candidate genes, RXFP2, appears to be reduced in the fetal horn bud of polled animals carrying the Celtic variant compared with horned individuals. Investigating horn ontogenesis and the genetic pathway by which the POLLED variants prevent horn development has implications for cattle breeding. If the genetic basis of horn bud formation and polledness is better understood, then new targets may be identified for precision genome editing to create polled individuals.

Beetle horns evolved from wing serial homologs.

Understanding how novel complex traits originate is a foundational challenge in evolutionary biology. We investigated the origin of prothoracic horns in scarabaeine beetles, one of the most pronounced examples of secondary sexual traits in the animal kingdom. We show that prothoracic horns derive from bilateral source tissues; that diverse wing genes are functionally required for instructing this process; and that, in the absence of Hox input, prothoracic horn primordia transform to contribute to ectopic wings. Once induced, however, the transcriptional profile of prothoracic horns diverges markedly from that of wings and other wing serial homologs. Our results substantiate the serial homology between prothoracic horns and insects wings and suggest that other insect innovations may derive similarly from wing serial homologs and the concomitant establishment of structure-specific transcriptional landscapes.

A study of lead uptake and distribution in horns from lead-dosed goats using synchrotron radiation-induced micro X-ray fluorescence elemental imaging.

OBJECTIVE: The principal goal of this study was to investigate the uptake and distribution of lead (Pb) in the horns of Pb-dosed goats, and to explore possible links to their historical Pb dosing records. Horn is a keratinized material that grows in discrete increments with the potential to preserve the historical record of past environmental exposures. While previous studies have leveraged this potential to examine environmental and biological phenomena in horns, Pb uptake has never been explored. METHODS: Horns were collected post-mortem from three goats that had been previously used to produce blood lead reference materials for the New York State proficiency testing program. The animals were periodically dosed with lead acetate, administered orally in a capsule, over a 5 to 8-year period. Horn cross sections were taken from each animal and analyzed using synchrotron radiation-induced micro X-ray fluorescence spectrometry (SR-µXRF) at the Cornell High Energy Synchrotron Source (CHESS). RESULTS: Elemental distribution maps were obtained by SR-µXRF for Pb, Ca, S, Se, and three other elements (Br, Zn and Cu), with values reported quantitatively as a mass fraction (µg/g for trace elements and mg/g for Ca and S). Accumulations of Pb were clearly visible as a series of narrow "rings" in each of the horn samples analyzed. The elements Ca, S, Br, Zn, and Cu were also detected as discrete rings within each cross-section, with Br strongly correlated with S in the samples examined. A marginal increase in Se may coincide with Pb accumulation in horn cross-sections. Annual mineralization estimates based on the relative distribution of Ca and S were used to establish a tentative timeline for horn growth, with each timeline linked to the pattern of Pb accumulation in the corresponding horn cross-section sample. CONCLUSIONS: Following ingestion, absorbed Pb is eventually deposited into caprine horns, resulting in discrete accumulations or "rings." Elemental mapping by SR-µXRF clearly show Ca-rich layers that vary with annual periodicity, consistent with previous reports of horn mineralization. Localized enrichment of Cu, Zn, Br and S appear to coincide with the keratinized regions related to the annual growth ring pattern in horns. Spatial analysis of horns for Pb accumulation may be useful as a qualitative marker of time-resolved exposures that may reflect specific periods of acute Pb absorption.

Precise staging of beetle horn formation in Trypoxylus dichotomus reveals the pleiotropic roles of doublesex depending on the spatiotemporal developmental contexts.

Many scarab beetles have sexually dimorphic exaggerated horns that are an evolutionary novelty. Since the shape, number, size, and location of horns are highly diverged within Scarabaeidae, beetle horns are an attractive model for studying the evolution of sexually dimorphic and novel traits. In beetles including the Japanese rhinoceros beetle Trypoxylus dichotomus, the sex differentiation gene doublesex (dsx) plays a crucial role in sexually dimorphic horn formation during larval-pupal development. However, knowledge of when and how dsx drives the gene regulatory network (GRN) for horn formation to form sexually dimorphic horns during development remains elusive. To address this issue, we identified a Trypoxylus-ortholog of the sex determination gene, transformer (tra), that regulates sex-specific splicing of the dsx pre-mRNA, and whose loss of function results in sex transformation. By knocking down tra function at multiple developmental timepoints during larval-pupal development, we estimated the onset when the sex-specific GRN for horn formation is driven. In addition, we also revealed that dsx regulates different aspects of morphogenetic activities during the prepupal and pupal developmental stages to form appropriate morphologies of pupal head and thoracic horn primordia as well as those of adult horns. Based on these findings, we discuss the evolutionary developmental background of sexually dimorphic trait growth in horned beetles.

The effectiveness of clove oil and two different cautery disbudding methods on preventing horn growth in dairy goat kids.

The effectiveness of clove oil and cautery disbudding on horn growth was evaluated in goat kids. The study used 243 Saanen doe kids (4±1 days old; mean±SD) on two goat farms that were disbudded with either (i) clove oil injection (CLOVE), (ii) a cautery iron and bud removed (BUDOFF), or (iii) a cautery iron with bud left intact (BUDON). Each kid received a different treatment per bud, which were balanced between buds (left/right) and randomly allocated. A trained observer monitored bud growth following treatment for 3 months recording either: N: no growth, H: normal horn, S: abnormal horn (scur), or SC: soft, fibrous lump (scorn). After the final observation, buds were assessed for the probability of detecting (i) success (no growth), (ii) scurs, (iii) horns or (iv) scorns [with 95% CI]. The probability of success for BUDOFF (0.77 [0.63, 0.87]) was higher than for BUDON (0.20 [0.11, 0.34]) and CLOVE (0.09 [0.04, 0.18]; P ≤ 0.05). Furthermore, the probability of success for BUDON was higher than for CLOVE (P ≤ 0.05). The probability of scurs was higher for CLOVE (0.72 [0.63, 0.80]) than BUDOFF (0.25 [0.17, 0.34]) and BUDON (0.30 [0.21, 0.39]; P ≤ 0.05). There was no difference in the probability of scurs for BUDOFF and BUDON (P > 0.05). The probability of horns was higher for CLOVE (0.21 [0.15, 0.29]) than BUDON (0.02 [0.01, 0.06]; P ≤ 0.05); horns were not observed for BUDOFF. The probability of scorns for BUDON, the only treatment that led to scorns, was 0.41 (0.25, 0.60). These results suggest that BUDOFF was more effective at preventing growth than CLOVE and BUDON and appears the most effective method, of the methods tested, for disbudding kids. Future research should explore other alternatives to cautery disbudding that may be both efficacious and cause less pain.

Genome-wide association study provides insights into genes related with horn development in Nelore beef cattle.

The causal mutation for polledness in Nelore (Bos taurus indicus) breed seems to have appeared first in Brazil in 1957. The expression of the polled trait is known to be ruled by a few groups of alleles in taurine breeds; however, the genetic basis of this trait in indicine cattle is still unclear. The aim of this study was to identify genomic regions associated with the hornless trait in a commercial Nelore population. A total of 107,294 animals had phenotypes recorded and 2,238 were genotyped/imputed for 777k SNP. The weighted single-step approach for genome-wide association study (WssGWAS) was used to estimate the SNP effects and variances accounted for by 1 Mb sliding SNP windows. A centromeric region of chromosome 1 with 3.11 Mb size (BTA1: 878,631-3,987,104 bp) was found to be associated with hornless in the studied population. A total of 28 protein-coding genes are mapped in this region, including the taurine Polled locus and the IFNAR1, IFNAR2, IFNGR2, KRTAP11-1, MIS18A, OLIG1, OLIG2, and SOD1 genes, which expression can be related to the horn formation as described in literature. The functional enrichment analysis by DAVID tool revealed cytokine-cytokine receptor interaction, JAK-STAT signaling, natural killer cell mediated cytotoxicity, and osteoclast differentiation pathways as significant (P < 0.05). In addition, a runs of homozygosity (ROH) analysis identified a ROH island in polled animals with 2.47 Mb inside the region identified by WssGWAS. Polledness in Nelore cattle is associated with one region in the genome with 3.1 Mb size in chromosome 1. Several genes are harbored in this region, and they may act together in the determination of the polled/horned phenotype. Fine mapping the locus responsible for polled trait in Nelore breed and the identification of the molecular mechanisms regulating the horn growth deserve further investigation.

Horn growth variation and hunting selection of the Alpine ibex.

Selective hunting can affect demographic characteristics and phenotypic traits of the targeted species. Hunting systems often involve harvesting quotas based on sex, age and/or size categories to avoid selective pressure. However, it is difficult to assess whether such regulations deter hunters from targeting larger "trophy" animals with longer horns that may have evolutionary consequences. Here, we compile 44,088 annually resolved and absolutely dated measurements of Alpine ibex (Capra ibex) horn growth increments from 8,355 males, harvested between 1978 and 2013, in the eastern Swiss Canton of Grisons. We aim to determine whether male ibex with longer horns were preferentially targeted, causing animals with early rapid horn growth to have shorter lives, and whether such hunting selection translated into long-term trends in horn size over the past four decades. Results show that medium- to longer-horned adult males had a higher probability of being harvested than shorter-horned individuals of the same age and that regulations do affect the hunters' behaviour. Nevertheless, phenotypic traits such as horn length, as well as body size and weight, remained stable over the study period. Although selective trophy hunting still occurs, it did not cause a measurable evolutionary response in Grisons' Alpine ibex populations; managed and surveyed since 1978. Nevertheless, further research is needed to understand whether phenotypic trait development is coinfluenced by other, potentially compensatory factors that may possibly mask the effects of selective, long-term hunting pressure.

The 1.78-kb insertion in the 3'-untranslated region of RXFP2 does not segregate with horn status in sheep breeds with variable horn status.

BACKGROUND: The mode of inheritance of horn status in sheep is far more complex than a superficial analysis might suggest. Observations, which were mostly based on crossbreeding experiments, indicated that the allele that results in horns is dominant in males and recessive in females, and some authors even speculated about the involvement of more than two alleles. However, all recent genome-wide association analyses point towards a very strong effect of a single autosomal locus on ovine chromosome 10, which was narrowed down to a putatively causal insertion polymorphism in the 3'-untranslated region of the relaxin/insulin-like family peptide receptor 2 gene (RXFP2). The main objective of this study was to test this insertion polymorphism as the causal mutation in diverse sheep breeds, including breeds with a variable and/or sex-dependent horn status. RESULTS: After re-sequencing a region of about 246 kb that covered the RFXP2 gene and its flanking regions for 24 sheep from six completely horned and six completely polled breeds, we identified the same insertion polymorphism that was previously published as segregating with horn status in these breeds. Multiplex PCR genotyping of 489 sheep from 34 breeds and some crosses between sheep breeds showed a nearly perfect segregation of the insertion polymorphism with horn status in sheep breeds of Central and Western European origin. In these breeds and their crossings, heterozygous males were horned and heterozygous females were polled. However, this segregation pattern was not, or at least not completely, reproducible in breeds with sex-dependent and/or variable horn status, especially in sheep that originated from even more southern European regions and from Africa. In such breeds, we observed almost all possible combinations of genotype, sex and horn status phenotype. CONCLUSIONS: The 1.78-kb insertion polymorphism in the 3'-untranslated region of RXFP2 and SNPs in the 3'-UTR, exon 14 and intron 11 of this gene that we analyzed in this study cannot be considered as the only cause of polledness in sheep and are not useful as a universal marker to define the genetic horn status in sheep.