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.

Early life investment in antlers and body growth reflects adult performance in roe deer population under supplementary feeding conditions.

Recent research has challenged the idea that cervid antlers are such costly traits, supporting the assertion early-life antler investment is an honest signal of adult phenotypic quality. However, inferences were made based on antler measurements while growing (velvet) and thus, studies on fully-formed clean antlers are needed to avoid possible bias caused by the inter-individual variation in antler growth phenology. We studied a semi-captive population of European roe deer inhabiting a sub-Mediterranean area (Valsemana research station) and living under optimal conditions (ad libitum food supply and veterinary care). Based on repeated measurements taken from 146 individuals, we assessed whether allocation to secondary sexual traits during early life affected body mass or antler development during adulthood. Furthermore, we evaluated which body measurements better predicted future male quality. Additionally, using 488 individuals, we described age-class-specific variation in male body measurements and assessed the decline in antler size between adult and senescent stages (i.e. magnitude of senescence). Results agree with the assertion that there is no negative effect of a high investment in antler during early life on adult antler size or body mass, but we suggest the use of clean antlers as a more robust and reliable measure. The variables that better predicted body mass during adulthood were yearling body mass and body size at 6 months. Antler score between 10 and 18 months resulted in the best indicator of adult antler size. Finally, we support the idea that the magnitude of senescence in antler size is smaller in environments with higher resource availability during winter.

Comparative antler proteome of sika deer from different developmental stages.

Antler is a special bone tissue that has the ability to regenerate completely periodically. It is the fastest growing bone in the animal kingdom. Antler provides a valuable research model for bone growth and mineralization. Antler grows longitudinally by endochondral ossification with their growth center located in its tip. Many scholars have carried out detailed studies on morphology and gene expression of antler tip. However, few scholars have analyzed the protein expression patterns of antler tip at different development stages. This study used label-free proteomics approach to analyze the protein expression dynamics of the antler tip in six developmental periods (15, 25, 45, 65, 100 and 130 days after the previous antler cast) and costal cartilage. In result, 2052 proteins were confidently quantified, including 1937 antler proteins and 1044 costal cartilage proteins. Moreover, 913 antler core proteins and 132 antler-special proteins were obtained. Besides, the stages special proteins and differentially expressed proteins (DEPs) in different development stages were analyzed. A total of 875 DEPs were determined by one-way AVOVA. It is found that the growth period (15, 25, 45 and 65 days) showed more up-regulated protein including several chondrogenesis-associated proteins (collagen types II, collagen types XI, HAPLN1, PAPSS1 and PAPSS2). In ossification stages, the up-regulated proteins related with lysosome (CTSD, CTSB, MMP9, CAII) indicated that the antler has higher bone remodeling activity. Given the up-regulated expression of immune-related molecules (S100A7, CATHL7, LTF, AZU1, ELANE and MPO), we speculate that the local immune system may contribute to the ossification of antler tip. In conclusion, proteomics technology was used to deeply analyze the protein expression patterns of antler at different development stages. This provides a strong support for the research on the molecular regulation mechanism of rapid growth and ossification of velvet antler.

Anti-tumour activity of deer growing antlers and its potential applications in the treatment of malignant gliomas.

A recent study showed that antlers have evolved a high rate of growth due to the expression of proto-oncogenes and that they have also evolved to express several tumour suppressor genes to control the risk of cancer. This may explain why deer antler velvet (DAV) extract shows anti-tumour activity. The fast growth of antler innervation through the velvet in close association to blood vessels provides a unique environment to study the fast but non-cancerous proliferation of heterogeneous cell populations. We set out to study the anti-cancer effect of DAV in glioblastoma (GB) cell lines in comparison with temozolomide, a chemotherapeutic drug used to treat high-grade brain tumours. Here we report, for the first time, that DAV extract from the tip, but not from mid-parts of the antler, exhibits an anti-tumour effect in GB cell lines (T98G and A172) while being non-toxic in non-cancerous cell lines (HEK293 and HACAT). In T98G cells, DAV treatment showed reduced proliferation (37.5%) and colony-formation capacity (84%), inhibited migration (39%), induced changes in cell cycle progression, and promoted apoptosis. The anticancer activity of DAV extract as demonstrated by these results may provide a new therapeutic strategy for GB treatment.

Comparative transcriptome analysis of the main beam and brow tine of sika deer antler provides insights into the molecular control of rapid antler growth.

BACKGROUND: Deer antlers have become a valuable model for biomedical research due to the capacities of regeneration and rapid growth. However, the molecular mechanism of rapid antler growth remains to be elucidated. The aim of the present study was to compare and explore the molecular control exerted by the main beam and brow tine during rapid antler growth. METHODS: The main beams and brow tines of sika deer antlers were collected from Chinese sika deer (Cervus nippon) at the rapid growth stage. Comparative transcriptome analysis was conducted using RNA-Seq technology. Differential expression was assessed using the DEGseq package. Functional Gene Ontology (GO) enrichment analysis was accomplished using a rigorous algorithm according to the GO Term Finder tool, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis was accomplished with the R function phyper, followed by the hypergeometric test and Bonferroni correction. Quantitative real-time polymerase chain reaction (qRT-PCR) was carried out to verify the RNA levels for differentially expressed mRNAs. RESULTS: The expression levels of 16 differentially expressed genes (DEGs) involved in chondrogenesis and cartilage development were identified as significantly upregulated in the main beams, including transcription factor SOX-9 (Sox9), collagen alpha-1(II) chain (Col2a1), aggrecan core protein (Acan), etc. However, the expression levels of 17 DEGs involved in endochondral ossification and bone formation were identified as significantly upregulated in the brow tines, including collagen alpha-1(X) chain (Col10a1), osteopontin (Spp1) and bone sialoprotein 2 (Ibsp), etc. CONCLUSION: These results suggest that the antler main beam has stronger growth capacity involved in chondrogenesis and cartilage development compared to the brow tine during rapid antler growth, which is mainly achieved through regulation of Sox9 and its target genes, whereas the antler brow tine has stronger capacities of endochondral bone formation and resorption compared to the main beam during rapid antler growth, which is mainly achieved through the genes involved in regulating osteoblast and osteoclast activities. Thus, the current research has deeply expanded our understanding of the intrinsic molecular regulation displayed by the main beam and brow tine during rapid antler growth.

To beat or not to beat: Behavioral plasticity during the antler growth period affects cortisol but not testosterone concentrations in red deer (Cervus elaphus) males.

Out of rut, male red deer (Cervus elaphus) associate themselves in bachelor groups where animals compete for rank position via agonistic interactions. In a previous study on red deer, males were recognized either as "Non-Fighters" (NF, low frequency of attacks) or "Fighters" (F, high frequency of attacks). This study, therefore, aims to verify the consistency of the inter-individual differences in fighting attitude across different social contexts and investigate whether they could be considered an individual characteristic. Behavioral consistency was presumed across three different sampling seasons, assuming that NF would have lower cortisol (C) and testosterone (T) concentrations than the F males. In 2015 the males were kept in one large group and labelled NF and F. In 2016, the herd was divided into two subgroups ("NF" and "F") based on the frequency of attacks. Finally, in 2017, the males were divided into two randomly composed subgroups. Data about agonistic behavior and concentration of C and T were collected during each season. In 2015 the individuals differed only for the fighting attitude. After the division, the frequency of the attacks always increased, being consistently lower in NF than in F. Unexpectedly, a slight increase in the concentration of C was detected in the NF in 2016, compared to the F who experienced no difference neither in 2015 nor 2017. No significant differences were found in T. We concluded that, even though the males had shown behavioral plasticity, their diversified interaction-prone attitude had been maintained despite the modifications of the social environment.

The Antler Cycle and Fecal Testosterone of Male Sambar Deer Rusa unicolor unicolor at the Horton Plains National Park in Sri Lanka.

This study is aimed at evaluating the relationship between endogenous testosterone levels and antler development in male sambar deer (Rusa unicolor unicolor) inhabiting the Horton Plains National Park, Sri Lanka. Seven antler growth stages of sambar were documented based on phenotypic observations for the first time in Sri Lanka as (a) cast, (b) growing 1-single spike, (c) growing 2-antler fork into a Y as the second tine appears, (d) growing 3-velvet begins to harden as the third tine appears, (e) growth completed-velvet shedding begins, (f) hard antler, and (g) casting. Fecal samples were collected every month for a period of eighteen months from male sambar deer in different stages of the antler growth cycle, feeding in the wet patana grasslands of the park, and the fecal testosterone level was estimated by radioimmunoassay. Ten animals were randomly selected from each antler stage for the experiment. The results disclose that the highest concentrations of testosterone were recorded in the hard antler stage. Velvet shedding was preceded by an increase in the testosterone level, and it is the sudden drop in the testosterone concentration which triggers the antler casting. The casting stage corresponded with the lowest mean testosterone concentration. Although the study was able to conclude a clear relationship between the fecal testosterone levels of the male sambar deer in the Horton Plains National Park and their antler stages, there is no clear seasonality for the antler cycle.

Social environment modulates investment in sex trait versus lifespan: red deer produce bigger antlers when facing more rivalry.

Theory predicts that the plastic expression of sex-traits should be modulated not only by their production costs but also by the benefits derived from the presence of rivals and mates, yet there is a paucity of evidence for an adaptive response of sex-trait expression to social environment. We studied antler size, a costly and plastic sex trait, and tooth wear, a trait related to food intake and longevity, in over 4,000 male Iberian red deer (Cervus elaphus hispanicus) from 56 wild populations characterized by two contrasting management practices that affect male age structure and adult sex-ratio. As a consequence, these populations exhibit high and low levels of male-male competition for mating opportunities. We hypothesized that males under conditions of low intra-sexual competition would develop smaller antlers, after controlling for body size and age, than males under conditions of high intra-sexual competition, thus reducing energy demands (i.e. reducing intake and food comminution), and as a consequence, leading to less tooth wear and a concomitant longer potential lifespan. Our results supported these predictions. To reject possible uncontrolled factors that may have occurred in the wild populations, we carried out an experimental design on red deer in captivity, placing males in separate plots with females or with rival males during the period of antler growth. Males living with rivals grew larger antlers than males living in a female environment, which corroborates the results found in the wild populations. As far as we know, these results show, for the first time, the modulation of a sexual trait and its costs on longevity conditional upon the level of intra-sexual competition.

Development of novel EST microsatellite markers for genetic diversity analysis and correlation analysis of velvet antler growth characteristics in Sika deer.

BACKGROUND: Sika deer is one of the most popular and valued animals in China. However, few studies have been conducted on the microsatellite of Sika deer, which has hampered the progress of genetic selection breeding. To develop and characterize a set of microsatellites for Sika deer which provide helpful information for protection of Sika deer natural resources and effectively increase the yield and quantity of velvet antler. RESULTS: We conducted a transcriptome survey of Sika deer using next-generation sequencing technology. One hundred eighty-two thousand two hundred ninety-five microsatellite markers were identified in the transcriptome, 170 of 200 loci were successfully amplified across panels of 140 individuals from Shuangyang Sika deer population. And 29 loci were found to be obvious polymorphic. Number of alleles is from 3 to 14. The expected heterozygosity ranged from 0.3087 to 0.7644. The observed heterozygosity ranged from 0 to 0.7698. The polymorphism information content values of those microsatellites varied ranged from 0.2602 to 0.7507. The marker-trait association was tested for 6 important and kernel characteristics of two-branched velvet antler in Shuangyang Sika deer through one-way analysis of variance. The results showed that marker-trait associations were identified with 8 different markers, especially M009 and M027. CONCLUSIONS: This study not only provided a large scale of microsatellites which were valuable for future genetic mapping and trait association in Sika deer, but also offers available information for molecular breeding in Sika deer.

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.

Analysis of genetic information from the antlers of Rangifer tarandus (reindeer) at the rapid growth stage.

Reindeer is the only deer species in which both males and females regularly grow antlers, providing an excellent model for studying the rapid growth and annual regeneration of antlers. The study of genetic information from reindeer is the basis for revealing the unique mechanism of antler growth. In the present study, we obtained 18.86 GB of clean reads, which were assembled to obtain 94,575 unigenes (average length: 704.69). Among these reads, 30,980 sequences were identified by searching a database of known proteins and then annotated with Gene Ontology (GO) terms, Clusters of Orthologous Groups (COG) classifications and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. All 7,480 simple sequence repeats (SSRs) were detected. A total of 84,435 and 82,226 high-quality single-nucleotide polymorphisms (SNPs) were identified in male and female reindeer, respectively. We identified 31 genes that were highly expressed in reindeer antlers. These genes regulate cell activities that are closely associated with the process of rapid tissue growth. Our results provide a basis for studying reindeer antlers and for further studying the molecular genetics, population genetics, and functional genomics of reindeer.

Novel insights into Dhh signaling in antler chondrocyte proliferation and differentiation: Involvement of Foxa.

The desert hedgehog (Dhh) is crucial for spermatogenesis and Leydig cell differentiation, but little is known regarding its physiological function in cartilage. In this study, Dhh mRNA was abundant in antler chondrocytes, where it advanced cell proliferation concomitant with accelerated transition from the G1 to the S phase and induced elevation of the hypertrophic chondrocyte markers, Col X and Runx2. Silencing of Ptch1 resulted in appreciable Smo accumulation and enhanced rDhh stimulation of Smo, whose impediment by cyclopamine obscured the proliferative function of Dhh and alleviated its guidance of chondrocyte differentiation. Further analysis evidenced the noteworthy positive action of Smo in the bridging between Dhh and Gli transcription factors. Obstruction of Gli1 by GANT58 caused the failed stimulation of Col X and Runx2 by rDhh. Analogously, siRNA against Gli1-3 hindered chondrocyte differentiation in the context of rDhh. Simultaneously, Gli transcription factors mediated the regulation of Dhh on Foxa1, Foxa2, and Foxa3, whose knockdown impaired chondrocyte differentiation. Attenuation of Foxa antagonized the augmentation of Col X and Runx2 generated by rDhh. Collectively, Dhh signaling through its target Foxa appears to induce antler chondrocyte proliferation and differentiation.

Biological adaptations in the Arctic cervid, the reindeer (Rangifer tarandus).

The reindeer is an Arctic species that exhibits distinctive biological characteristics, for which the underlying genetic basis remains largely unknown. We compared the genomes of reindeer against those of other ruminants and nonruminant mammals to reveal the genetic basis of light arrhythmicity, high vitamin D metabolic efficiency, the antler growth trait of females, and docility. We validate that two reindeer vitamin D metabolic genes (CYP27B1 and POR) show signs of positive selection and exhibit higher catalytic activity than those of other ruminants. A mutation upstream of the reindeer CCND1 gene endows an extra functional binding motif of the androgen receptor and thereby may result in female antlers. Furthermore, a mutation (proline-1172→threonine) in reindeer PER2 results in loss of binding ability with CRY1, which may explain circadian arrhythmicity in reindeer.

Transcriptomic analysis of different tissue layers in antler growth Center in Sika Deer (Cervus nippon).

BACKGROUND: With the unprecedented rapid growth rate (up to 2.75 cm/day), velvet antler is an invaluable model for the identification of potent growth factors and signaling networks for extremely fast growing tissues, mainly cartilage. Antler growth center (AGC) locates in its tip and consists of five tissue layers: reserve mesenchyme (RM), precartilage (PC), transition zone (TZ), cartilage (CA) and mineralized cartilage (MC). The aim of this study was to investigate the transcription dynamics in the AGC using RNA-seq technology. RESULTS: Five tissue layers in the AGC were collected from three 3-year-old male sika deer using our previously reported sampling method (morphologically distinguishable). After sequencing (15 samples; triplicates/tissue layer), we assembled a reference transcriptome de novo and used RNA-seq to measure gene expression profiles across these five layers. Nine differentially expressed genes (DEGs) were selected from our data and subsequently verified using qRT-PCR. The results showed a high consistency with the RNA-seq results (R2 = 0.80). Nine modules were constructed based on co-expression network analysis, and these modules contained 370 hub genes. These genes were found to be mainly involved in mesenchymal progenitor cell proliferation, chondrogenesis, osteogenesis and angiogenesis. Combination of our own results with the previously published reports, we found that Wnt signaling likely plays a key role not only in stimulating the antler stem cells or their immediate progeny, but also in promoting chondrogenesis and osteogenesis during antler development. CONCLUSION: We have successfully assembled a reference transcriptome, generated gene expression profiling across the five tissue layers in the AGC, and identified nine co-expressed modules that contain 370 hub genes and genes predorminantly expressed in and highly relevant to each tissue layer. We believe our findings have laid the foundation for the identification of novel genes for rapid proliferation and chondrogenic differentiation of antler cells.

Identification of the miRNA-mRNA regulatory network of antler growth centers.

Antler growth is a unique event compared to other growth and development processes in mammals. Antlers grow extremely fast during the rapid growth stage when growth rate peaks at 2 cm per day. Antler growth is driven by a specific endochondral ossification process in the growth center that is in the distal region of the antler tip. In this study, we used state-of-art RNA-seq technology to analyze the expression profiles of mRNAs and miRNAs during antler growth. Our results indicated that the expression levels of multiple genes involved in chondrogenesis and endochondral ossification, including Fn1, Sox9, Col2a1, Acan, Col9a1, Col11a1, Hapln1, Wwp2, Fgfr3, Comp, Sp7 and Ihh, were significantly increased at the rapid growth stage. Our results also indicated that there were multiple differentially expressed miRNAs interacting with differentially expressed genes with opposite expression patterns. Furthermore, some of the miRNAs, including miR-3072-5p, miR-1600, miR-34-5p, miR-6889-5p and miR-6729-5p, simultaneously interacted with and controlled multiple genes involved in the process of chondrogenesis and endochondral ossification. Therefore, we established a miRNA-mRNA regulatory network by identifying miRNAs and their target genes that were differentially expressed in the antler growth centers by comparing the rapid growth stage and the initial growth stage.

Comparison of linear model and artificial neural network using antler beam diameter and length of white-tailed deer (Odocoileus virginianus) dataset.

Evaluation of harvest data remains one of the most important sources of information in the development of strategies to manage regional populations of white-tailed deer. While descriptive statistics and simple linear models are utilized extensively, the use of artificial neural networks for this type of data analyses is unexplored. Linear model was compared to Artificial Neural Networks (ANN) models with Levenberg-Marquardt (L-M), Bayesian Regularization (BR) and Scaled Conjugate Gradient (SCG) learning algorithms, to evaluate the relative accuracy in predicting antler beam diameter and length using age and dressed body weight in white-tailed deer. Data utilized for this study were obtained from male animals harvested by hunters between 1977-2009 at the Berry College Wildlife Management Area. Metrics for evaluating model performance indicated that linear and ANN models resulted in close match and good agreement between predicted and observed values and thus good performance for all models. However, metrics values of Mean Absolute Error and Root Mean Squared Error for linear model and the ANN-BR model indicated smaller error and lower deviation relative to the mean values of antler beam diameter and length in comparison to other ANN models, demonstrating better agreement of the predicted and observed values of antler beam diameter and length. ANN-SCG model resulted in the highest error within the models. Overall, metrics for evaluating model performance from the ANN model with BR learning algorithm and linear model indicated better agreement of the predicted and observed values of antler beam diameter and length. Results of this study suggest the use of ANN generated results that are comparable to Linear Models of harvest data to aid in the development of strategies to manage white-tailed deer.

Compositional analysis of the glycosaminoglycan family in velvet antlers of Sika deer (Cervus nippon) at different growing stages.

Glycosaminoglycans (GAG) from the velvet antlers of Sika deer (Cervus nippon) at the different growing stages (Fukurozuno, Anshi, and Santajo) of bred and wild deer were isolated and their concentrations and sulfation patterns were analyzed. GAG were digested with chondroitinase ABC, ACI, heparinase-I and -III, and keratanase-II into the corresponding repeating disaccharides of chondroitin sulfate (CS), dermatan sulfate (DS), hyaluronan, heparan sulfate (HS), and keratan sulfate. Cartilaginous tissues contained CS-DS at high concentrations with an almost equal ratio of 4- and 6-sulfates, while 4-sulfate-type CS-DS predominantly occupied ossified tissues, but at low concentrations. High O- and N-sulfation degrees of HS correspond to high ossification. Dynamic quantitative changes in CS-DS and compositional changes in CS-DS and HS were closely associated with the mineralization of deer antlers.

Proteomic analysis of the effects of antler extract on chondrocyte proliferation, differentiation and apoptosis.

Deer antlers are unique cranial appendages capable of regeneration and rapid growth. In addition, deer antlers have been widely used in traditional Chinese medicine to promote the function of the kidneys, reproductive system, bones and nervous system. It has been shown that water-soluble substances are the major bioactive components within the deer antlers. In this study, we prepared aqueous extracts from deer antlers during a rapid growth stage. We investigated the effects of antler extracts on primary chondrocytes by analyzing their protein expression patterns using isobaric tags for relative and absolute quantitation technology. We demonstrated that antler extracts promote chondrocyte proliferation and prevent chondrocyte differentiation and apoptosis by controlling multiple cellular processes involved in genomic stability, epigenetic alterations, ribosome biogenesis, protein synthesis and cytoskeletal reorganization. Antler extracts significantly increased the expression levels of proliferation markers Mki67 and Stmn1 and differentiation inhibitor Acp5 as well as cellular apoptosis inhibitors Ndufa4l2 and Rcn1. Thus, this study has greatly expanded our current knowledge of the molecular effects of antler extracts on chondrocytes. It has also shed new light on possible strategies to prevent damage to and to treat cartilage and its related diseases by using aqueous extracts from growing Sika deer antlers.

Full-length transcriptome and microRNA sequencing reveal the specific gene-regulation network of velvet antler in sika deer with extremely different velvet antler weight.

Velvet antler displays the fastest and most robust tissue proliferation in the animal world, it is a model for a complete organ development/regeneration, and alternative medicine, tonic made from velvet antler, was beneficial for human. The weight of velvet antler had high biomedical and economic value, but the related regulation mechanisms controlling velvet antler weight remain unclear. In this study, extremely heavy and light velvet antler groups were selected from a sika deer population of 100 individuals with extreme velvet antler weight. A combination of full-length transcriptome sequencing and microRNA sequencing to the proliferation zone in the tip of velvet antler was applied. A total of 55306 transcripts and 1082 microRNAs were identified. Some highly expressed genes (COL1A1, COL1A2, COL3A1, FN1, and ATP6) and microRNAs (miR-21, let-7i, and miR-27b) were highly correlated with the physiological and growth characteristics of velvet antlers. Among the 334 differentially expressed genes, we found that most of the genes were located in the developmental process, especially animal organ development process. It is exciting to see that more blood vessels were found in the growing tip of heavy velvet antler through histological observation, and GO term of blood vessel development was also significant different between two groups. The combination analysis with mRNA and microRNA data in velvet antler showed a specific regulation network involved in the development of bone, mesenchyme, cartilage, and blood vessel, and helped us clearly find out the candidate 14 genes and 6 microRNAs, which could be used for selecting significant DNA markers of velvet antler weight.

Transcriptomic characterization elucidates a signaling network that controls antler growth.

Deer antlers are amazing appendages with the fastest growth rate among mammalian organs. Antler growth is driven by the growth center through a modified endochondral ossification process. Thus, identification of signaling pathways functioning in antler growth center would help us to uncover the underlying molecular mechanism of rapid antler growth. Furthermore, exploring and dissecting the molecular mechanism that regulates antler growth is extremely important and helpful for identifying methods to enhance long bone growth and treat cartilage- and bone-related diseases. In this study, we build a comprehensive intercellular signaling network in antler growth centers from both the slow growth stage and rapid growth stage using a state-of-art RNA-Seq approach. This network includes differentially expressed genes that regulate the activation of multiple signaling pathways, including the regulation of actin cytoskeleton, calcium signaling, and adherens junction. These signaling pathways coordinately control multiple biological processes, including chondrocyte proliferation and differentiation, matrix homeostasis, mechanobiology, and aging processes, during antler growth in a comprehensive and efficient manner. Therefore, our study provides novel insights into the molecular mechanisms regulating antler growth and provides valuable and powerful insight for medical research on therapeutic strategies targeting skeletal disorders and related cartilage and bone diseases.