A consensus genome of sika deer (Cervus nippon) and transcriptome analysis provided novel insights on the regulation mechanism of transcript factor in antler development.

BACKGROUND: Sika deer (Cervus nippon) holds significance among cervids, with three genomes recently published. However, these genomes still contain hundreds of gaps and display significant discrepancies in continuity and accuracy. This poses challenges to functional genomics research and the selection of an appropriate reference genome. Thus, obtaining a high-quality reference genome is imperative to delve into functional genomics effectively. FINDINGS: Here we report a high-quality consensus genome of male sika deer. All 34 chromosomes are assembled into single-contig pseudomolecules without any gaps, which is the most complete assembly. The genome size is 2.7G with 23,284 protein-coding genes. Comparative genomics analysis found that the genomes of sika deer and red deer are highly conserved, an approximately 2.4G collinear regions with up to 99% sequence similarity. Meanwhile, we observed the fusion of red deer's Chr23 and Chr4 during evolution, forming sika deer's Chr1. Additionally, we identified 607 transcription factors (TFs) that are involved in the regulation of antler development, including RUNX2, SOX6, SOX8, SOX9, PAX8, SIX2, SIX4, SIX6, SPI1, NFAC1, KLHL8, ZN710, JDP2, and TWST2, based on this consensus reference genome. CONCLUSIONS: Our results indicated that we acquired a high-quality consensus reference genome. That provided valuable resources for understanding functional genomics. In addition, discovered the genetic basis of sika-red hybrid fertility and identified 607 significant TFs that impact antler development.

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.

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.

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.

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.

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.

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.

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.

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.

Screening weight related genes of velvet antlers by whole genome re-sequencing.

The velvet antler is a special organ that has important biological significance for deer, and its growth is a complicated biological metabolism process. Growing evidence suggests that genetics factors play essential roles in the weight of velvet antlers. In this study, we investigated five sika deer (Cervus nippon) populations under the same feeding condition, and screened genetic variations in the 100 samples (including 50 heavy and 50 light velvet antler weight samples) by whole genome re-sequencing. The results showed that 94 genetic variations were related to the velvet antler weight, among which two single nucleotide polymorphism (SNP) sites were located on the exon regions of OAS2 and ALYREF/THOC4, respectively. Furthermore, ALYREF/THOC4 is highly expressed in the velvet antler. The biological functions of these genetic variations were highly related to the growth and development of deer velvet antlers. Collectively, we screened genes related to the velvet antler weight in sika deer populations by whole genome re-sequencing and identified 94 sites as candidate genetic variations related to the velvet antler weight. We hope that it will contribute to further mechanistic studies of velvet antler development and weight variations.

MHC, parasites and antler development in red deer: no support for the Hamilton & Zuk hypothesis.

The Hamilton-Zuk hypothesis proposes that the genetic benefits of preferences for elaborated secondary sexual traits have their origins in the arms race between hosts and parasites, which maintains genetic variance in parasite resistance. Infection, in turn, can be reflected in the expression of costly sexual ornaments. However, the link between immune genes, infection and the expression of secondary sexual traits has rarely been investigated. Here, we explored whether the presence and identity of functional variants (supertypes) of the highly polymorphic major histocompatibility complex (MHC), which is responsible for the recognition of parasites, predict the load of lung and gut parasites and antler development in the red deer (Cervus elaphus). While we found MHC supertypes to be associated with infection by a number of parasite species, including debilitating lung nematodes, we did not find support for the Hamilton-Zuk hypothesis. On the contrary, we found that lung nematode load was positively associated with antler development. We also found that the supertypes that were associated with resistance to certain parasites at the same time cause susceptibility to others. Such trade-offs may undermine the potential genetic benefits of mate choice for resistant partners.

iTRAQ-Based Quantitative Proteomic Analysis of the Potentiated and Dormant Antler Stem Cells.

As the only known organ that can completely regenerate in mammals, deer antler is of real significance in the field of regenerative medicine. Recent studies have shown that the regenerative capacity of the antlers comes from the pedicle periosteum and the cells resident in the periosteum possess the attributes of stem cells. Currently, the molecular mechanism of antler regeneration remains unclear. In the present study, we compared the potentiated and dormant antler stem cells using isobaric tags for the relative and absolute quantification (iTRAQ) labeling of the peptides, coupled with two-dimensional liquid chromatography-tandem mass spectrometry (LC-MS/MS) to compare the proteome profiles. Proteins were identified by searching against the NCBI nr database and our own Cervine transcriptome database, and bioinformatics analysis was conducted to identify the differentially expressed proteins. Based on this searching strategy, we identified 169 differentially expressed proteins in total, consisting of 70 up- and 99 down-regulated in the potentiated vs. dormant antler stem cells. Reliability of the iTRAQ was confirmed via quantitative real-time polymerase chain reaction (qRT-PCR) to measure the expression of selected genes. We identified transduction pathways through the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, such as HIF-1 and PI3K-AKT signaling pathways that play important roles in regulating the regeneration of antlers. In summary, the initiation stage of antler regeneration, a process from dormant to potentiated states in antler stem cells, is regulated by multiple proteins and complicated signal networks.

Direct localisation of molecules in tissue sections of growing antler tips using MALDI imaging.

The astonishing growth rate of deer antlers offers a valuable model for the discovery of novel factors and regulatory systems controlling rapid tissue growth. Numerous molecules have been identified in growing antlers using a variety of techniques. However, little is known about the spatial distribution of these molecules in situ. A technique that has the potential to help in this regard is direct proteomic analysis of tissue sections by matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI-IMS). The present study applied this technique to spatially map molecules in antler tissue sections. Two protonated molecular ions were selected: m/z 6679 and m/z 6200 corresponding to VEGF and thymosin beta-10, respectively. Superimposition of the respective ion images on to histologically stained samples showed distinct spatial distribution across the antler tissue sections which were consistent with the previous reports using in situ hybridization. Two other molecular ions specifically m/z 8100 and m/z 11,800 were also selected, corresponding to reported masses of urocortin precursor and thioredoxin, respectively. As the spatial distribution of these proteins is not specifically known, MALDI-IMS was used as a potential technique to obtain information on their distribution on antler tips. The presence of all these molecules in deer antlers were further confirmed using LC-MS/MS data. The present study also demonstrated that MALDI-IMS could be further used to image antler sections with an extended ion mass range of up to m/z 45,000, thus potentially increasing the ability to discover the distribution of a larger set of molecules that may play an important role in antler growth. We have thus demonstrated that MALDI-IMS is a promising technique for generating molecular maps with high spatial resolution which can aid in evaluating the function of novel molecules during antler growth.

CONTROLLING ANTLER GROWTH IN A CASTRATED INDOCHINESE SIKA DEER CERVUS NIPPON PSEUDAXIS USING A COMMERCIALLY AVAILABLE TRENBOLONE ACETATE AND ESTRADIOL IMPLANT.

A captive Indochinese sika deer (Cervus nippon pseudaxis) was castrated at the age of 5 yr. The resultant abnormal antler growth over the next few years became difficult to manage from both the veterinary and husbandry standpoints. Using a commercially available trenbolone acetate and estradiol implant marketed for domestic cattle heifers, normal mineralization of the abnormal antlers was achieved along with the expected normal casting. The deer was then maintained for 6 yr using an annual implant regimen.

MicroRNA let-7a and let-7f as novel regulatory factors of the sika deer (Cervus nippon) IGF-1R gene.

MicroRNAs and their roles in rapid antler growth and regeneration have attracted much attention. In the present study, we examined the effects of microRNAs let-7a and let-7f on antler cell proliferation. We used a luciferase reporter screen to demonstrate that insulin-like growth factor 1 receptor (IGF-1R) can be regulated by let-7a and let-7f. MTT assay confirmed that chondrocyte proliferation was inhibited by let-7a and let-7f mimics. In contrast, transfection of let-7a and let-7f inhibitors increased chondrocyte proliferation, indicating that inhibitors can competitively bind to endogenous miRNA, reducing the inhibitory effect of miRNA. Moreover, western blotting analysis further identified that let-7a and let-7f mimics suppressed IGF-1R expression, and that let-7a and let-7f inhibitors increased the expression level of IGF-1R. Taken together, our study demonstrates the important roles of let-7a and let-7f in antler proliferation and its potential application in antler development. let-7a and let-7f may represent novel regulatory factors of IGF-1R expression in deer antler.

The allometry between secondary sexual traits and body size is nonlinear among cervids.

Allometric relationships between sexually selected traits and body size have been extensively studied in recent decades. While sexually selected traits generally display positive allometry, a few recent reports have suggested that allometric relationships are not always linear. In male cervids, having both long antlers and large size provides benefits in terms of increased mating success. However, such attributes are costly to grow and maintain, and these costs might constrain antler length from increasing at the same rate as body mass in larger species if the quantity of energy that males can extract from their environment is limiting. We tested for possible nonlinearity in the relationship between antler size and body mass (on a log-log scale) among 31 cervids and found clear deviation from linearity in the allometry of antler length. Antler length increased linearly until a male body mass threshold at approximately 110 kg. Beyond this threshold, antler length did not change with increasing mass. We discuss this evidence of nonlinear allometry in the light of life-history theory and stress the importance of testing for nonlinearity when studying allometric relationships.

Association between melatonin receptor 1A (MTNR1A) gene single-nucleotide polymorphisms and the velvet antler yield of Sika deer.

Melatonin, a secretion from pineal gland is ambiguously considered as the key hormone involved in regulation of the antler cycle in Sika deer. To find out more about the roles of melatonin and its receptor gene, we carried out current study to investigate the association between polymorphisms in melatonin receptor 1A (MTNR1A) gene and the antler yield from Sika deer. A total of 251 Sika deer were analyzed in this study, of which consisted of Wusan Sika deer (n = 163) and Dongfeng Sika deer (n = 88). MTNRA gene was amplified by PCR and genotyped by Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Three polymorphism loci (C518T, C629G and C635T) were detected in exon2 of MTNR1A gene. The restriction site Ecol881 was used for C518T while a C629G polymorphism locus was digested with Mval restriction endonucleases. In Wusan Sika deer the allele frequencies of C and T were 0.637 and 0.363 for C518T, Also C and G alleles in C629G locus were 0.206 and 0.794. Genotypic frequencies of allele CC, CT and TT were 33.7, 59.9 and 6.4 % respectively, It showed 1.8, 37.4 and 60.7 % for frequencies of genotypes CC, CG and GG. In Dongfeng Sika deer the allele frequencies of C and T were 0.518 and 0.482 for C518T, C and G alleles were 0.375 and 0.625 for C629G. Genotypic frequencies were 10.6, 82.4 and 7.1 % for genotypes CC, CT and TT respectively, and they were 1.1, 72.7 and 26.2 % for genotypes CC, CG and GG. Among three SNPs, only C629G showed significant association (P < 0.05) with average antler yield in Wusan Sika deer, while no SNP was significant in Dongfeng Sika deer. These preliminary results implied that the identified SNPs of MTNR1A gene might influence the antler yield in Wusan Sika deer.

Identification of microRNA-18a as a novel regulator of the insulin-like growth factor-1 in the proliferation and regeneration of deer antler.

To investigate the effect of miR-18a on the regulation of the insulin-like growth factor (IGF-1) during growth of antlers in sika deer, miRNA Chip analysis, Target Scan and real-time PCR analysis were used to identify miRNAs that bind to the 3'-UTR of IGF-1. An miR-18a mimic was transfected into antler cartilage cells and the expression levels were quantified by real-time PCR. Dual luciferase assays revealed that miR-18a binds to the 3'-UTR of the IGF-1 gene thus indicating this to be a target gene regulated by miR-18a. MTT assays and cell-cycle analyses confirmed that miR-18a significantly inhibited proliferation of cartilage cells. In contrast, transfection of miR-18a inhibitors increased proliferation. Furthermore, Western blot analysis showed that over-expression of miR-18a down-regulated IGF-1 protein levels while IGF-1 expression was increased after transfection of miR-18a inhibitors. Thus, miR-1 down-regulated IGF-1 expression thus implicating miR-18a as an important regulator of antler proliferation.