Comprehensive transcriptome analysis of sika deer antler using PacBio and Illumina sequencing

Antler is the fastest growing and ossifying tissue in animals and it is a valuable model for cartilage/bone development. To understand the molecular mechanisms of chondrogenesis and osteogenesis of antlers, the PacBio Sequel II and Illumina sequencing technology were combined and used to investigate the mRNA expression profiles in antler tip, middle, and base at six different developmental stages, i.e., at 15th, 25th, 45th, 65th, 100th and 130th growth days. Consequently, we identified 24,856 genes (FPKM > 0.1), including 8778 novel genes. Besides, principal component analysis (PCA) revealed a significant separation between the growth stage (25th, 45th and 65th days) and ossification stage (100th and 130th days). COL2A1 gene was significantly abundant in the growth stage, whereas S100A7, S100A12, S100A8, and WFDC18 genes were abundant at the ossification stage. Subsequently screened to 14,765 significantly differentially expressed genes (DEGs), WGCNA and GO functional enrichment analyses revealed that genes related to cell division and chondrocyte differentiation were up-regulated, whereas those with steroid hormone-mediated signaling pathways were down-regulated at ossification stages. Additionally, 25 tumor suppressor genes and 11 oncogenes were identified and were predicted to interact with p53. Co-regulation of tumor suppressor genes and oncogenes is responsible for the special growth pattern of antlers. Together, we constructed the most complete sika deer antler transcriptome database so far. The database provides data support for subsequent studies on the molecular mechanism of sika deer antler chondrogenesis and osteogenesis.