Transcription profiling provides insights into gene pathways involved in horn and scurs development in cattle
© Mariasegaram et al; licensee BioMed Central Ltd. 2010
Received: 4 March 2010
Accepted: 11 June 2010
Published: 11 June 2010
Two types of horns are evident in cattle - fixed horns attached to the skull and a variation called scurs, which refers to small loosely attached horns. Cattle lacking horns are referred to as polled. Although both the Poll and Scurs loci have been mapped to BTA1 and 19 respectively, the underlying genetic basis of these phenotypes is unknown, and so far, no candidate genes regulating these developmental processes have been described. This study is the first reported attempt at transcript profiling to identify genes and pathways contributing to horn and scurs development in Brahman cattle, relative to polled counterparts.
Expression patterns in polled, horned and scurs tissues were obtained using the Agilent 44 k bovine array. The most notable feature when comparing transcriptional profiles of developing horn tissues against polled was the down regulation of genes coding for elements of the cadherin junction as well as those involved in epidermal development. We hypothesize this as a key event involved in keratinocyte migration and subsequent horn development. In the polled-scurs comparison, the most prevalent differentially expressed transcripts code for genes involved in extracellular matrix remodelling, which were up regulated in scurs tissues relative to polled.
For this first time we describe networks of genes involved in horn and scurs development. Interestingly, we did not observe differential expression in any of the genes present on the fine mapped region of BTA1 known to contain the Poll locus.
Three genes are believed to control the development of horn phenotypes, namely Poll, African horn and Scurs . The Poll gene, which results in the absence of horns, was mapped to the centromeric region of bovine chromosome 1 (BTA1) by Georges et al.  and has been confirmed by several groups [3–5]. Using 20 microsatellites, the Poll locus has been restricted to a 1 Mb interval between the microsatellite markers RP42-218J17_MS1 and BM6438 in the centromeric region of BTA1 . A report of the location of the Scurs locus responsible for producing a loose or "wobbly horn" to BTA19  has not been confirmed . The postulated African horn locus has not been mapped to a chromosome .
The mapping of the Poll locus has not led to the identification of the causative genetic variant for the trait. Mouse genetic models have in the past yielded useful positional candidate genes for important production traits in livestock species, for example MSTN (muscle hypertrophy) , and DGAT1 (milk yield and composition) . As horn development is unique to the Bovidae, no small animal models for the phenotype exist, making the identification of candidate genes difficult. Additionally, ambiguities in phenotypic characterisation of horns as distinct from scurs has also precluded the precise mapping and identification of candidate genes for scurs [5, 8]. A less constrained approach involves the use of DNA microarrays to uncover transcripts associated with different morphologies, as successfully demonstrated in the identification of calmodulin as a mediator of beak size in Darwin's finches 
It is thought that horn development is primarily controlled by the skin . During the initial stages of horn bud development, the epidermis stops making hair and begins to synthesise horn. Once initiated, the primordium of the bony core of the horn forms as a separate centre of ossification in the dermal connective tissue underneath the horn-forming region, which is later fused to the skull . Consequently, gene expression was evaluated in horn bud tissues obtained from the epidermal and dermal tissues of 1-2 week old Brahman calves.
In this study we report results from transcript profiling to identify key genes and molecular mechanisms underpinning the development of horns and scurs in Brahman cattle. We were interested in describing the network of differential gene expression that occurs during the development of the horn and scurs in cattle by contrasting to polled animals. Additionally we wanted to test the hypothesis of a highly differentially expressed gene on BTA1 and 19, for horn and scurs respectively, corresponding to the reported loci.
RNA extraction, RIN evaluation and hybridization
For RNA preparation, the tissues biopsies were blotted dry to remove excess RNAlater®, trimmed of the topmost keratin layer, wrapped in heavy duty aluminium foil and snap-frozen in liquid nitrogen. The frozen samples were physically disrupted by impact and then transferred to 2 ml screw cap tubes containing 1 ml Tri-Reagent (Applied Biosystems) and a mixture of 2- and 4- mm acid washed glass beads. The samples were homogenised in a Thermo-Savant FastPrep® (Thermo Fisher Scientific Inc., MA, USA) for 25 s and the lysate treated as previously described . Total RNA was isolated using the RiboPure™ kit (Applied Biosystems), according to the manufacturer's instructions. All samples that showed OD260/280 ratios > 1.8 were examined on the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) using an RNA 6000 Nano LabChip (Agilent Technologies) and only those yielding RNA with RNA Integrity Number (RIN) > 6.5 were hybridized to a microarray. From the original samples, one of the horn samples was discarded on account of low RIN resulting in a total of 11 samples (6 male and 5 female) corresponding to the three phenotypic categories of horned (3), polled (4) and scurred (4) from the first week time point. These were hybridized against the Agilent Bovine - Four Plex G2519F 44 k array in a dye swap experiment. Amplification of mRNA, hybridization and data scanning services were performed at the SRC Microarray Facility, University of Queensland http://microarray.imb.uq.edu.au/.
Agilent array annotation
The bovine Agilent array contains 21,475 unique 60-mer probes, designed to represent approximately 19,500 distinct bovine genes. Probe annotation was accomplished through a series of sequence searching steps using the Baylor College of Medicine Bovine Genome (version 4) as the reference sequence. The probes were aligned to Btau4.0. Both human and bovine reference gene sets were also aligned to Btau4.0. A comparison was then made of the locations of the probes with respect to the locations of the genes. It was noted where a probe overlapped with an intron or an exon of a gene. If there was no overlap, the nearest gene was noted. A direct BLAST  of all of the Agilent probe sequences against the most current bovine reference gene set was also performed
The gene prediction for each probe was performed using the following algorithm.
Step 1: If the probe sequence overlapped with a bovine exon and that gene name did not begin with 'LOC or 'MGC' we used that gene name. If there were more than one overlapping bovine gene, the name corresponding to the lowest accession number was used.
Step 2: If the overlapping bovine gene began with either 'LOC' or 'MGC' or when there was no overlapping bovine gene, we used the overlapping human gene name. In the event there was more than one overlapping human gene, the name corresponding to the lowest accession number was used.
Step 3: If there was no overlapping gene, we used the name of the nearest gene based on the alignment to Btau4.0.
Step 4: If there was direct blast hit to the bovine reference set and the name of that gene differed from the one we selected using either of the first two steps, the entry was highlighted to signify uncertainty with the annotation. These annotations were then manually checked to choose the correct gene name or to discard the probe when the issue was difficult to resolve.
Microarray design and Statistical analysis
Gene expression intensity signals were subjected to a series of data acquisition criteria based on signal to noise ratio and mean to median correlation as detailed Tan et al. . In brief, we employed the following two editing criteria for data acquisition: First, we required that the signal to noise ratio s > 1, where s is obtained by dividing the background corrected intensity by the standard deviation of the background pixels. Second, we required that the correlation between the mean and the median signal intensities r > 0.85 where r is obtained by dividing the smaller of the mean or median by the larger. Tran et al.  suggested that a correlation of r ≥ 0.85 not only retains more data than other methods, but the data retained are more accurate than those obtained using traditional thresholds or common spot flagging algorithms. Spots that failed to pass these filtering criteria were assigned a reading of zero. This resulted in a total of 1,375,680 gene expression intensity readings (687,840 of each channel, red and green) that were corrected for background intensity and then log2 transformed. The arithmetic mean and standard deviation for the background corrected log transformed red and green intensities were 6.97, SD = 4.21 and 7.90, SD = 3.25, respectively.
where Y ijkhsmn represents the n-th background-adjusted, normalized base-2 log-intensity from the g-th gene (or probe) at the h-th phenotype variety (horn status) and s-th sex, from the i-th chip, j-th array (i.e., there are four microarrays per chip) and k-th dye channel; μ is the overall mean; C represents a comparison fixed group effects with 32 levels and defined as those intensity measurements from the same chip, array and dye channel; G represent the random gene (or probe) effects with 21,475 levels; AG, DG, HG, SG are the random interaction effects of array by gene, dye by gene, horn-phenotype by gene, and sex by gene, respectively; and finally, ε is the random error term.
In what follows, it is understood that the h-th horn phenotype incorporates the three possible phenotypes. That is: h = 1, 2, and 3 for Polled, Horned, and Scurred, respectively.
where iid denotes independently and identically distributed and N denotes the normal distribution. Variance components between genes (σ2 g ), between genes within array (σ2 ag ), between genes within dye (σ2 dg ), between genes within horn phenotype (σ2 hg ), between genes within sex (σ2 sg ), and within genes (σ2 e ) were estimated using restricted (to zero error contrasts) maximum likelihood based on the analytical gradients option of VCE6 software ftp://ftp.tzv.fal.de/pub/vce6/.
Finally, the DE measurement contrasts in (2, 3, 4) were processed by fitting a two-component normal mixture model and posterior probabilities of belonging to the non-null component were used to identify DE genes with an estimated experiment-wise false discovery rate of < 1% as described by McLachlan et al. .
Gene ontology distribution and hierarchial clustering analysis
For all of the differentially expressed probes, the corresponding human GeneIDs were analysed in the DAVID database  and in Ingenuity Pathway Analysis (IPA, Ingenuity® Systems, http://www.ingenuity.com. This allowed for the identification of enriched biological processes and cellular components specific to the three contrasts of interest, namely, PvH (polled versus horn) and PvS (polled versus scurs). The differentially expressed genes assigned to functional clusters with the DAVID software were also analysed using Permutmatrix . This program performed hierarchial clustering and seriation based on normalised gene expression values to produce a dendrogram. To explore molecular interaction networks among the gene expression profiles, equivalent HUGO gene IDs of the probes were uploaded into IPA Ingenuity Pathway Analysis) software together with the fold changes for each contrast of interest. The program generated networks based on its knowledge base of interactions for the gene list.
RT-PCR confirmation of DE genes
Validation of the differential expression of genes was carried out for 3 selected genes, namely desmocollin 1 (DSC1), desmoglein 1 (DSG1) and dehydrogenase 7C (DHRS7C), using qPCR. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as the internal reference gene based on the observation of stable expression in animals from all 3 phenotypic classes in the microarray analysis. Taqman gene expression assays for each of these genes were purchased from Applied Biosystems (Assay IDs: Bt03212566_m1, Bt03212583_m1, Bt03241459_m1, Bt03210915_g1 corresponding to DSC1, DSG1, DHRS7C and GAPDH respectively).
RNA was isolated from tissue biopsies collected in the second week after birth as described earlier. There were 12 samples in each phenotype category of polled, horned and scurred in the validation experiment with 4 biological replicates and 3 repeats in each sample. All of the animals sampled for the horned phenotype and one of the scurred individual were different from those used for the initial microarray study. In each case, a total of 1 μg of total RNA was treated with Turbo DNA-free kit™ (Applied Biosystems) and 200 ng of the treated RNA was used in cDNA synthesis with random hexamers and the SuperScript™ III First-Strand Synthesis System (Applied Biosystems). The cDNA was diluted one in ten, from which 5 μl was used in subsequent pre-amplification reactions using the Taqman PreAmp Master Mix kit (Applied Biosystems P/N 4384556) according to the manufacturer's instruction. For the Real-time PCR step, amplifications were performed in triplicate on the 7900 Sequence Detection System (Applied Biosystems) using standard cycling parameters.
The relative expression level was determined by subtracting the Ct value (cycle number at the threshold level of log-based fluorescence) of the reference gene from the corresponding target gene across all phenotypic classes of interest. A total of three phenotypic contrasts - polled vs horned (PvH), polled vs scurs (PvS) and horned vs scurs (HvS) - were analysed for statistical significance at each gene. The effects of phenotype, gene and their interaction were included in the generalised linear model with Ct as a dependent variable. The model tested the null hypotheses that Ct differences between target and reference genes were the same between a pair of phenotypes represented as combinational effect (CE): CE1 - CE3 = CE2 - CE4 or CE1 - CE2 - CE3 + CE4 = 0, which yields an estimate of ΔΔCt . For example, in a PvH contrast with DSC1 as target gene and GAPDH as reference gene, CE1 represents polled-DSC1, CE3 represents polled-GAPDH, CE2 represents horned-DSC1, and CE4 represents horned-GAPDH. The Proc GLM procedure in SAS (SAS Institute Inc., Cary, NC) was used for the estimation of ΔΔCt and the evaluation of significance.
(i) Array annotation results
Results from annotation of probes on the Agilent bovine array
(ii) Microarray results
Overrepresented GO terms from functional annotation clustering in DAVID corresponding to those categories of terms with enrichment fold greater than 3.0.
Description of change
Enriched gene category
Genes down regulated in horned versus polled animals
COL17A1, KLK5, SPINK5, GJB5, KLK7, IVL, GRHL3, CALML5, KRT5, SPRR3, TCHH, LAMC2, DSP, KRT17, S100A7, LOR, BARX2, KRT10,
striated muscle contraction
MYOM2, SMPX, MYBPC1, MYH7, PGAM2, MYH2, KBTBD10, TTN
DSC1, KRT40, KRT33A, COL17A1, DSC3, KRT25, GJB5, GJB4, KRT23, KRT5, DSG1, LAMC2, KRT27, KRT17, GJA3, KRT10,
KRT40, KRT33A, KRT25, KRTAP1-1, KRT23, KRTAP1-5, KRT5, KRT26, DSP, KRT27, KRT17, KRT10, PKP1
ACTN2, KRT40, DSC1, MYL1, PKP3, DSC3, LMOD2, MYH7, NEB, KRTAP1-1, KRT25, CSTA, MYH2, MYOZ2, KRT23, KRTAP1-5, SPRR3, KRT26, MYL3, KRT10, PKP1, KRT33A, ACTA1, KBTBD10, CGN, IVL, ELMOD1, DST, PDLIM3, MC1R, MYOM2, ELMO3, KRT5, MYBPC1, MYOT, SPTBN2, TCHH, KRT27, DSP, KRT17, LOR
DSC1, COL17A1, PKP3, DSC3, GJB5, CLDN8, CGN, GJB4, DST, CLDN23, DSG1, DSP, GJA3, PKP1
Genes upregulated in horned versus polled animals
GDF6, IGF1, PROC, ANG, DLK1, PTN, SCUBE1, FMOD, PLA2G7, CHI3L1, SPP2, CCL8, CP, HP, MASP1, SFRP4, SPON1, CFI, COL8A1, APOL3, CRLF1, PRSS2, PI15, FGG, FGF11, CCL19, TF, EPDR1
Genes down regulated in scurred versus polled animals
intermediate filament cytoskeleton
KRT40, KRT33A, KRT15, KRTAP1-5, KRTAP2-1, KRT85, NEFH, KRT26, KRTAP1-1, KRT25, KRT27
Genes upregulated in scurred versus polled animals
COL1A2, COL5A2, LUM, COL1A1, COL5A1
COMP, PHEX, IBSP, DMP1, OSTN, ALX1, MEPE, COL1A2, SPP1, POSTN, SPP2, COL1A1, MMP9
TNC, LAMA1, TNN, COL1A2, COL5A2, SPP1, COL6A6, COL1A1, ITGA11, IBSP, COL5A1
(iii) Validation of the DE genes by quantitative Real-Time PCR
Changes in expression of DSC1, DSG1 and DHRS7C levels between polled, horned and scurs evaluated by q-PCR.
Contrasts of interest
q-PCR Fold change#
Microarray fold change#
In the microarray analyses, DSC1 mRNA showed nearly 9-fold higher expression in polled tissues compared to horned, whereas qPCR found a 11-fold increased expression in polled animals compared to horned individuals. In PvS, the microarray and qPCR showed agreement in level of DSC1 expression, where both results were insignificant. DSG1 mRNA expression in the microarray study was 7.5 fold compared to 7.7 fold from qPCR in the polled versus horned animals. In PvS comparison, neither microarray nor qPCR revealed changes in DSG1 expression levels. DHRS7C mRNA expression was 12.6 fold higher in polled versus horned animals compared to 3.8 fold difference between polled and scurred animals from the microarray experiments. The qPCR results show a similar trend of up-regulation of DHRS7C gene in polled versus horned and scurs animals, but with much greater fold change.
This study represents the first attempt at transcript profiling to uncover gene expression changes associated with the development of horns and scurs in cattle. Towards this end we sampled epidermal and dermal tissues from the skull region of 1-2 week old calves and analysed the extracted mRNA against the Agilent 44 K bovine array.
The fine mapping of the polled locus to a 1 Mb region on BTA1  and presence of a well characterised 2.5 Mb contig spanning this region  have enabled comparative mapping to the human genome, to identify the comprehensive list of genes present in this critical region. Interestingly, none of the genes in the 1 Mb interval encompassing the polled locus were found to be differentially expressed, even though all of them were represented on the microarray. While there were no obvious candidates relating to horn development in this region, it remains uncertain whether the inability to identify differential expression in genes from this region is a consequence of the sampling period. Clearly, if the inductors of horn development preceded the 1st week biopsy or even occurred at the foetal stage, then there would be little chance of detecting their presence in our design. However, a novel regulatory mechanism located on BTA1 igniting these processes involved in cell communication and differentiation cannot be ruled out. In the absence of a clear differentially expressed candidate gene for polled, we analysed the gene list for enriched GO term categories (Table 2) and network relationships between them (Figure 6, Figure 7).
One of the most consistent observations in the differential gene expression results from the PvH comparison was the enrichment for genes relating to skin development as well as striated muscle contraction. Further scrutiny of these genes, indicated that entire modules of genes encoding structural components of cell junctions and attachment show lower levels of expression in horned tissues compared to their age matched, polled counterparts.
Cell junctions are fundamental to cellular attachments and tissue structure, consisting of tight junction made up of cadherins, desmosomes and gap junction proteins [26–28]. The down-regulation of some of these genes or their complexes in horned animals has a striking resemblance to the process of epithelial to mesenchymal transition, also known as EMT. This is both an important developmental and pathological process, whereby the dissociation of the intercellular junction complexes enables epithelial cells to acquire a mesenchymal phenotype [29, 30]. In turn, epithelial cells are able to migrate to other environments where they may be involved in tissue morphogenesis . A hallmark of the EMT process involves the loss of intercellular adhesion, in which E-cadherin is vital for its role in promoting homotypic interactions between neighbouring cells . We observed a 4 fold decrease in the expression of E-cadherin (CDH1) in horned animals compared to polled [Additional file 1].
Another notable feature of the PvH comparison is the lower expression of genes encoding cytokeratin intermediate filaments, principally the keratins and keratin associated proteins in the samples from developing horns, when compared to the samples from prospective polled calves (Table 2). Interestingly many of these keratin genes are present in a cluster on BTA5. Since horns are composed of keratin, at first glance this result does not appear to make sense. Viewed in the context of the EMT, however, the gene expression differences may reflect the remodelling of the epidermis prior to horn formation. Cytokeratins provide the necessary attachment for junction protein to connect to neighbouring cells as well as the extracellular matrix and therefore contributing to the integrity of tissues. In a recent study Liovic et al.  showed that mutations in keratin genes caused noticeable down-regulation in many of the keratin-interacting components. Consistent with this study, our own observations showed significantly lowered expression of desomoscollins (DSC1, DSC3), desmogleins (DSG1), desmoplakins (DSP), plakophilins (PKP1, PKP3), claudins (CLDN 8/11/23) and gap junction protein genes (GJA3, GJB4 and GJB5) in horned compared to polled animals, as can be clearly seen in Figure 5 and Figure 6A. The fact that we were able to confirm the differential expression of DSC1 and DSG1 by q-PCR in 2 week old samples as opposed to week 1 biopsies used in the microarray study suggests that the lowered expression of these adhesion molecules maybe be prolonged for the duration of horn development.
Showing a similar trend of lowered expression in horned compared to polled are two further clusters of genes relating to epidermal differentiation that are present on BTA3 and 18, respectively. The cluster on BTA3 enclodes 4 members of the S100 family of small acidic calcium binding proteins including, loricrin (LOR), small proline rich protein 2A (SPRR2A) and involcurin (IVL). All of these protein products are important constituents of the cornified epithelium. The equivalent human gene cluster, termed "the epidermal differentiation complex" is located on HSA 1q21 . Another cluster on BTA18 codes for the kallikreins, which are serine proteases. Four of the subfamily members (KLK5, KLK6, KLK7, KLK11) are down-regulated in PvH. Kallikreins have been implicated in tissue desquamation as well in EMT [35, 36]. One could thus envision a "jack pot" effect where an overlapping regulator modulates expression through their interaction with these distinct but functionally related genes.
Beyond the contribution of cytokeratins to the intermediate filament net, the other important structural components of living cells are actin filaments which contribute to cell motility, structure and integrity. Cell migration involves dynamic interactions between the extracellular matrix and the actin cytoskeleton mediated through integrin adhesions. Several intergrin-associated proteins are known to connect the integrins with the actin cytoskeleton such as alpha-actinin and talins among others [37–39]. The expression of alpha-actinin is 18-fold lower in horned animals versus polled. An important component of the actin cytoskeleton is myosin II which acts as the bridge that organises actin into higher order structures . Phosphorylation of the myosin II light chain facilitates interaction of the heavy chain with F-actin . Several isoforms of the heavy (e.g. MYH2) and regulatory light chains (e.g MYL1) of myosin II also show lower expression in horned compared to polled tissues (Figure 6B). Other components of the actomyosin complex that show a similar trend in PvH include Cingulin and Nebulin. Cingulin has been shown to interact with myosin to link the tight junction proteins to the actomyosin cytoskeleton . Nebulin has been correlated with actin filament length . These changes underscore the flux in cytoskeletal reassembly underpinning horn morphogenesis.
In the PvS comparison, the lower expression of cytokeratin components in scurred animals relative to polled represents the main commonality with the development of horns. However, the main distinguishable feature of scurs was the elevated expression of many members of fibrillar collagen genes which code for components of the proteinaceous extracellular matrix (ECM). For instance both COL1A1 and COL1A2 show over two-fold increased expression in scurs tissues over polled. Other ECM components showing up-regulation are additional members of the collagen family such as COL5A1, COL5A2 and laminin alpha 1 (LAMA1). Some of the most highly expressed genes in the IPA pathways are the matrix metalloproteinases (MMP9, MMP12 and MMP13), MMP13 expression is nearly 10-fold higher in scurred individuals compared to polled. The difference in the gene expression of the MMP between horned and polled tissues is not statistically significant.
In humans, MMP are represented by 24 genes all of which are collectively capable of cleaving virtually all components of the ECM and the basement membrane . The collagenase activity of MMP  supports the induction of the collagens during scurs development. MMP have been shown to have important role during tissue remodelling, wound healing and repair through regulated degradation of the ECM. The other feature of MMP is their involvement in the regulation of a variety of non-matrix substrates such as chemokines, cytokines and growth factors  (Figure 7B). When degradation is uncoordinated, MMP have been implicated in arthritis, cancer and cardiovascular disease . A study on laminitis in cattle has also noted the up-regulation of several MMP in ulcerated bovine claw in conjunction with the observation of the disruption of the basement membrane .
It remains uncertain whether the lack of greater overlap of differentially expressed genes between PvH and PvS contrasts simply reflects a more advanced state of development in the horn tissue, consistent with the observations that scurs phenotypes in calves takes longer to manifest than horns. The q-PCR validation of greater expression of DHRS7C in scurs versus horned tissues is significant in light of its localization on BTA19 and would thus warrant further investigation.
Although there has been no recent study on the physiology of the horn development in cattle, many of the features of horn morphogenesis maybe viewed as equivalent processes to cattle hoof development. The latter has been studied extensively because of the importance of lameness in both cattle and horses. Hoof development involves the intracellular keratinization of epithelial cells and their progressive migration through the suprabasal layers in a program of terminal differentiation, to finally replace the desquamated cornified cells shed from the outermost layers to produce the hardened horn .
The down-regulation of cell adhesion molecules that was evident in horn development is consistent with findings reported by Donetti et al. . They found for instance that expression of Dsc1 was undetectable in the keratinizing oral mucosa to coincide with keratinocyte terminal differentiation. A similar investigation comparing plakophilin 1 (PKP1) deficient keratinocytes to those expressing recombinant plakophilin 1 found that it had a dramatic effect on desmosomal content . It would thus appear that the breakdown of adhesion and dissolution of the cytoskeleton evident in horn tissues might be setting the stage for keratinocyte transition through the epithelium. Additionally, the enrichment of gene networks relating to EMT may also play a part in the skeletal development of horns. Interestingly, we did not identify changes to cell adhesion molecules in scurs tissues. However, the cytokeratins were down-regulated as in horned tissues. Another feature of scurs was the up-regulation of ECM components which we did not observe in horned tissues.
In conclusion, this study provides new insights into the dynamics of horns and scurs development. Most of the research into EMT processes has been established in vitro, where for instance the actin cytoskeleton and adhesion related molecules have been studied as disparate entities. For the first time we demonstrate all these inter-related molecules working together in a coordinated fashion to produce the horn and scurs phenotypes. Importantly, this work generates new hypothesis on potential regulators of horn development, which linkage studies have mapped to chromosome 1. Although the microarray results indicate a lack of differentially expressed transcript in this region, one cannot rule out a novel regulatory element that was not included in the Agilent array. Recently, Wright et al.  showed that a massive duplication of sequences in the intron of SOX5 caused the pea-comb phenotype in chicken. Like the horn forming layer in cattle, the pea-comb of chickens is composed of layers of epidermis, dermis and connective tissues. This begs the question whether the polled phenotype is caused by a similar CNV or some structural change on BTA1 that affects entire cascades of genes involved in EMT like processes?
This project was funded by Meat and Livestock Australia (MLA) through grant no. AHW.144 to KP. The authors would like to thank Paul Williams, Warren Sim, Rob Young, Brian Burns and Geoffry Fordyce for their help in sample collection and data recording. The contributions of polled cattle breeders of 'Fairysprings' and 'Kenilworth Brahmans' is highly appreciated. The authors are grateful to QFAB (Queensland Facility for Advanced Bioinformatics) for access to the IPA program. We also wish to acknowledge Neil Bagnall and Trudi Flatscher-Bader for sharing their expertise on RNA extraction, Mythily Mariasegaram for helpful discussion on EMT and Bill Barendse for useful suggestions on an earlier draft of this manuscript.
- Long CR, Gregory KE: Inheritance of the horned, scurred and polled condition in cattle. Journal of Heredity. 1978, 69 (6): 395-400.Google Scholar
- Georges M, Drinkwater R, King T, Mishra A, Moore SS, Nielsen D, Sargeant LS, Sorensen A, Steele MR, Zhao XY: Microsatellite mapping of a gene affecting horn development in Bos-taurus. Nature Genet. 1993, 4 (2): 206-210. 10.1038/ng0693-206.PubMedView ArticleGoogle Scholar
- Brockmann GA, Martin J, Teuscher F, Schwerin M: Marker controlled inheritance of the polled locus in Simmental cattle. Archiv Fur Tierzucht-Archives of Animal Breeding. 2000, 43 (3): 207-212.Google Scholar
- Harlizius B, Tammen I, Eichler K, Eggen A, Hetzel DJS: New markers on bovine chromosome 1 are closely linked to the polled gene in Simmental and Pinzgauer cattle. Mamm Genome. 1997, 8 (4): 255-257. 10.1007/s003359900404.PubMedView ArticleGoogle Scholar
- Brenneman RA, Davis SK, Sanders JO, Burns BM, Wheeler TC, Turner JW, Taylor JF: The polled locus maps to BTA1 in a Bos indicus × Bos taurus cross. Journal of Heredity. 1996, 87 (2): 156-161.PubMedView ArticleGoogle Scholar
- Drogemuller C, Wohlke A, Momke S, Distl O: Fine mapping of the polled locus to a 1-Mb region on bovine chromosome 1q12. Mamm Genome. 2005, 16 (8): 613-620. 10.1007/s00335-005-0016-0.PubMedView ArticleGoogle Scholar
- Asai M, Berryere TG, Schmutz SM: The scurs locus in cattle maps to bovine chromosome 19. Animal Genetics. 2004, 35 (1): 34-39. 10.1111/j.1365-2052.2003.01079.x.PubMedView ArticleGoogle Scholar
- Capitan A, Grohs C, Gautier M, Eggen A: The scurs inheritance: new insights from the French Charolais breed. BMC Genet. 2009, 10 (33):Google Scholar
- Prayaga KC: Genetic options to replace dehorning in beef cattle - a review. Australian Journal of Agricultural Research. 2007, 58 (1): 1-8. 10.1071/AR06044.View ArticleGoogle Scholar
- Grobet L, Martin LJR, Poncelet D, Pirottin D, Brouwers B, Riquet J, Schoeberlein A, Dunner S, Menissier F, Massabanda J: A deletion in the bovine myostatin gene causes the double-muscled phenotype in cattle. Nature Genet. 1997, 17 (1): 71-74. 10.1038/ng0997-71.PubMedView ArticleGoogle Scholar
- Grisart B, Farnir F, Karim L, Cambisano N, Kim JJ, Kvasz A, Mni M, Simon P, Frere JM, Coppieters W: Genetic and functional confirmation of the causality of the DGAT1 K232A quantitative trait nucleotide in affecting milk yield and composition. Proc Natl Acad Sci USA. 2004, 101 (8): 2398-2403. 10.1073/pnas.0308518100.PubMed CentralPubMedView ArticleGoogle Scholar
- Abzhanov A, Kuo WP, Hartmann C, Grant BR, Grant PR, Tabin CJ: The calmodulin pathway and evolution of elongated beak morphology in Darwin's finches. Nature. 2006, 442 (7102): 563-567. 10.1038/nature04843.PubMedView ArticleGoogle Scholar
- Hall BK: Bones and cartilage: developmental and evolutionary skeletal biology. 2005, San Diego, CA: Elsevier Academic PressGoogle Scholar
- Dove WF: The physiology of horn growth: A study of the morphogenesis, the interaction of tissues, and the evolutionary processes of a mendelian recessive character by means of transplantation of tissues. J Exp Zool. 1935, 69 (3): 347-405. 10.1002/jez.1400690302.View ArticleGoogle Scholar
- Piper EK, Jackson LA, Bagnall NH, Kongsuwan KK, Lew AE, Jonsson NN: Gene expression in the skin of Bos taurus and Bos indicus cattle infested with the cattle tick, Rhipicephalus (Boophilus) microplus. Vet Immunol Immunopathol. 2008, 126 (1-2): 110-119. 10.1016/j.vetimm.2008.06.011.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. Journal of Molecular Biology. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Tan S, Reverter A, Wang Y, Byrne KA, McWilliam SM, Lehnert SA: Gene expression profiling of bovine in vitro adipogenesis using a cDNA microarray. Functional & Integrative Genomics. 2006, 6 (3): 235-249. 10.1007/s10142-005-0016-x.View ArticleGoogle Scholar
- Tran PH, Peiffer DA, Shin Y, Meek LM, Brody JP, Cho KWY: Microarray optimizations: increasing spot accuracy and automated identification of true microarray signals. Nucleic Acids Res. 2002, 30 (12): 9-10.1093/nar/gnf053.View ArticleGoogle Scholar
- Reverter A, Barris W, McWilliam S, Byrne KA, Wang YH, Tan SH, Hudson N, Dalrymple BP: Validation of alternative methods of data normalization in gene co-expression studies. Bioinformatics. 2005, 21 (7): 1112-1120. 10.1093/bioinformatics/bti124.PubMedView ArticleGoogle Scholar
- Reverter A, Wang YH, Byrne KA, Tan SH, Harper GS, Lehnert SA: Joint analysis of multiple cDNA microarray studies via multivariate mixed models applied to genetic improvement of beef cattle. J Anim Sci. 2004, 82 (12): 3430-3439.PubMedGoogle Scholar
- McLachlan GJ, Bean RW, Jones LBT: A simple implementation of a normal mixture approach to differential gene expression in multiclass microarrays. Bioinformatics. 2006, 22 (13): 1608-1615. 10.1093/bioinformatics/btl148.PubMedView ArticleGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for annotation, visualization, and integrated discovery. Genome Biology. 2003, 4 (9): 11-10.1186/gb-2003-4-9-r60.View ArticleGoogle Scholar
- Caraux G, Pinloche S: PermutMatrix: a graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics. 2005, 21 (7): 1280-1281. 10.1093/bioinformatics/bti141.PubMedView ArticleGoogle Scholar
- Yuan JS, Reed A, Chen F, Stewart CN: Statistical analysis of real-time PCR data. Bmc Bioinformatics. 2006, 7: 10.1186/1471-2105-7-85.Google Scholar
- Wunderlich KR, Abbey CA, Clayton DR, Song Y, Schein JE, Georges M, Coppieters W, Adelson DL, Taylor JF, Davis SL: A 2.5-Mb contig constructed from Angus, Longhorn and horned Hereford DNA spanning the polled interval on bovine chromosome 1. Animal Genetics. 2006, 37 (6): 592-594. 10.1111/j.1365-2052.2006.01538.x.PubMedView ArticleGoogle Scholar
- Stokes DL: Desmosomes from a structural perspective. Curr Opin Cell Biol. 2007, 19 (5): 565-571. 10.1016/j.ceb.2007.09.003.PubMed CentralPubMedView ArticleGoogle Scholar
- Kljuic A, Bazzi H, Sundberg JP, Martinez-Mir A, O'Shaughnessy R, Mahoney MG, Levy M, Montagutelli X, Ahmad W, Alta VM: Desmoglein 4 in hair follicle differentiation and epidermal adhesion: Evidence from inherited hypotrichosis and acquired pemphigus vulgaris. Cell. 2003, 113 (2): 249-260. 10.1016/S0092-8674(03)00273-3.PubMedView ArticleGoogle Scholar
- Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P: Cell junctions, cell adhesion, and the extracellular matrix. Molecular biology of the cell. 2008, New York: Garland Science, 1131-1204.Google Scholar
- Shook D, Keller R: Mechanisms, mechanics and function of epithelial-mesenchymal transitions in early development. Mech Dev. 2003, 120 (11): 1351-1383. 10.1016/j.mod.2003.06.005.PubMedView ArticleGoogle Scholar
- Boyer B, Valles AM, Edme N: Induction and regulation of epithelial-mesenchymal transitions. Biochem Pharmacol. 2000, 60 (8): 1091-1099. 10.1016/S0006-2952(00)00427-5.PubMedView ArticleGoogle Scholar
- Savagner P: Leaving the neighborhood: molecular mechanisms involved during epithelial-mesenchymal transition. Bioessays. 2001, 23 (10): 912-923. 10.1002/bies.1132.PubMedView ArticleGoogle Scholar
- Moreno-Bueno G, Portillo F, Cano A: Transcriptional regulation of cell polarity in EMT and cancer. Oncogene. 2008, 27 (55): 6958-6969. 10.1038/onc.2008.346.PubMedView ArticleGoogle Scholar
- Liovic M, D'Alessandro M, Tomic-Canic M, Bolshakov VN, Coats SE, Lane EB: Severe keratin 5 and 14 mutations induce down-regulation of junction proteins in keratinocytes. Exp Cell Res. 2009, doi:10.1016/j.yexcr.2009.07.013Google Scholar
- Zeeuwen P: Epidermal differentiation: The role of proteases and their inhibitors. Eur J Cell Biol. 2004, 83 (11): 761-773. 10.1078/0171-9335-00388.PubMedView ArticleGoogle Scholar
- Emami N, Diamandis EP: Human tissue kallikreins: A road under construction. International Conference on Enzymes - Old Molecules with New Cllinical Applications: Oct 25 2005; Padova, ITALY. 2005, Elsevier Science Bv, 78-84.Google Scholar
- Whitbread AK, Veveris-Lowe TL, Lawrence MG, Nicol DL, Clements JA: The role of kallikrein-related peptidases in prostate cancer: potential involvement in an epithelial to mesenchymal transition. 1st International Symposium on Kallikreins: Sep 01-03 2005; Lausanne, SWITZERLAND. 2005, Walter De Gruyter & Co, 707-714.Google Scholar
- Cortesio CL, Chan KT, Huttenlocher A: Adhesion dyanimcs in motile cells. Cell Junctions - Adhesion, Development, and Disease. Edited by: LaFlamme SE, Kowalczyk A. 2008, Weinheim: Wiley-VCH Verlag GmbH & Co, 71-88.View ArticleGoogle Scholar
- Hemmings L, Rees DJG, Ohanian V, Bolton SJ, Gilmore AP, Patel B, Priddle H, Trevithick JE, Hynes RO, Critchley DR: Talin contains three actin-binding sites each of which is adjacent to a vinculin-binding site. J Cell Sci. 1996, 109: 2715-2726.PubMedGoogle Scholar
- Otey CA, Pavalko FM, Burridge K: An interaction between alpha-actinin and beta-1 integrin subunit in vitro. J Cell Biol. 1990, 111 (2): 721-729. 10.1083/jcb.111.2.721.PubMedView ArticleGoogle Scholar
- Manser E: Rho GTPases in the organisation of the actin cytoskeleton. Rho GTPases. Edited by: Symons M. 2004, NY, USA: Eurekah.com and Kluwer Academic/Plenum Publishers, 107-114.Google Scholar
- Olson MF, Sahai E: The actin cytoskeleton in cancer cell motility. Clin Exp Metastasis. 2009, 26 (4): 273-287. 10.1007/s10585-008-9174-2.PubMedView ArticleGoogle Scholar
- D'Atri F, Citi S: Cingulin interacts with F-actin in vitro. FEBS Lett. 2001, 507 (1): 21-24. 10.1016/S0014-5793(01)02936-2.PubMedView ArticleGoogle Scholar
- Horowits R: Nebulin regulation of actin filament lengths: new angles. Trends in Cell Biology. 2006, 16 (3): 121-124. 10.1016/j.tcb.2006.01.003.PubMedView ArticleGoogle Scholar
- Pardo A, Selman M: Matrix metalloproteases in aberrant fibrotic tissue remodeling. Proc Am Thorac Soc. 2006, 3 (4): 383-388. 10.1513/pats.200601-012TK.PubMedView ArticleGoogle Scholar
- Gross J, Lapiere CM: Collagenolytic activity in amphibian tissues - A tissue culture assay. Proc Natl Acad Sci USA. 1962, 48 (6): 1014-10.1073/pnas.48.6.1014.PubMed CentralPubMedView ArticleGoogle Scholar
- Rath T, Roderfeld M, Graf J, Roeb E: Matrix Metalloproteinases in Inflammatory Bowel Disease - From Basic Research to Clinical Significance. Z Gastroent. 2009, 47 (8): 758-769. 10.1055/s-0028-1109520.View ArticleGoogle Scholar
- Murphy G, Nagase H: Progress in matrix metalloproteinase research. Molecular Aspects of Medicine. 2008, 29 (5): 290-308. 10.1016/j.mam.2008.05.002.PubMed CentralPubMedView ArticleGoogle Scholar
- Hendry KAK, Knight CH, Galbraith H, Wilde CJ: Basement membrane integrity and keratinization in healthy and ulcerated bovine hoof tissue. Journal of Dairy Research. 2003, 70 (1): 19-27. 10.1017/S0022029902005885.PubMedView ArticleGoogle Scholar
- Tomlinson DJ, Mulling CH, Fakler TM: Invited review: Formation of keratins in the bovine claw: Roles of hormones, minerals, and vitamins in functional claw integrity. J Dairy Sci. 2004, 87 (4): 797-809. 10.3168/jds.S0022-0302(04)73223-3.PubMedView ArticleGoogle Scholar
- Donetti E, Bedoni M, Boschini E, Dellavia C, Barajon I, Gagliano N: Desmocollin 1 and desmoglein 1 expression in human epidermis and keratinizing oral mucosa: a comparative immunohistochemical and molecular study. Arch Dermatol Res. 2005, 297 (1): 31-38. 10.1007/s00403-005-0573-9.PubMedView ArticleGoogle Scholar
- South AP, Wan H, Stone MG, Dopping-Hepenstal PJC, Purkis PE, Marshall JF, Leigh IM, Eady RAJ, Hart IR, McGrath JA: Lack of plakophilin 1 increases keratinocyte migration and reduces desmosome stability. J Cell Sci. 2003, 116 (16): 3303-3314. 10.1242/jcs.00636.PubMedView ArticleGoogle Scholar
- Wright D, Boije H, Meadows JRS, Bed'hom B, Gourichon D, Vieaud A, Tixier-Boichard M, Rubin CJ, Imsland F, Hallbook F: Copy Number Variation in Intron 1 of SOX5 Causes the Pea-comb Phenotype in Chickens. Plos Genetics. 2009, 5 (6): 10-10.1371/journal.pgen.1000512.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.