Deer antler is the only completely regenerable organ in mammals. During the rapid growth period, the antler proliferates even faster than cancerous tissue growth. However, the proliferation and development of antler have been in a stable and controllable growth cycle. In this study, we analyzed the time series expression data of nine samples from mesenchyme layer in three male sika deer in the early period of the antler with a saddle-like appearance (30 days), the rapid growth period of the antler with two branches (60 days), and the final period of the antler with three branches (90 days). Whole Transcriptome sequencing results show that in the 30 d versus 60 d group, 1,464 genes, 85 long noncoding RNAs (lncRNAs), and 61 miRNAs were identified as differentially expressed; 1,748 genes, 138 lncRNAs, and 78 miRNAs were identified as differentially expressed in 30d versus 90d group; and 816 differentially expressed genes (DEGs), 49 differentially expressed lncRNAs (DE lncRNAs), and 24 differentially expressed miRNA (DE miRNAs) were identified in 60d versus 90d group. A total of 182 miRNA-mRNA interaction pairs and 89 miRNA-lncRNA interaction pairs were screened from DEGs, DE miRNAs, and DE lncRNAs to construct the ceRNA regulatory network (ceRNET). In summary, we identified candidate mRNAs, miRNAs and lncRNAs that regulate the development of antler tip. It may lay the foundation for further investigating the molecular mechanism of antler rapid growth and development.

Despite the scientific and medicinal importance of diploid sika deer (Cervus nippon), its genome resources are limited and haplotype-resolved chromosome-scale assembly is urgently needed. To explore mechanisms underlying the expression patterns of the allele-specific genes in antlers and the chromosome evolution in Cervidae, we report, for the first time, a high-quality haplotype-resolved chromosome-scale genome of sika deer by integrating multiple sequencing strategies, which was anchored to 32 homologous groups with a pair of sex chromosomes (XY). Several expanded genes (RET, PPP2R1A, PPP2R1B, YWHAB, YWHAZ, and RPS6) and positively selected genes (eIF4E, Wnt8A, Wnt9B, BMP4, and TP53) were identified, which could contribute to rapid antler growth without carcinogenesis. A comprehensive and systematic genome-wide analysis of allele expression patterns revealed that most alleles were functionally equivalent in regulating rapid antler growth and inhibiting oncogenesis. Comparative genomic analysis revealed that chromosome fission might occur during the divergence of sika deer and red deer (Cervus elaphus), and the olfactory sensation of sika deer might be more powerful than that of red deer. Obvious inversion regions containing olfactory receptor genes were also identified, which arose since the divergence. In conclusion, the high-quality allele-aware reference genome provides valuable resources for further illustration of the unique biological characteristics of antler, chromosome evolution, and multi-omics research of cervid animals.

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.

The most southern population of reindeer (Rangifer tarandus) inhabits northeastern China, but the migration route and origin of this population have not been confirmed. The sequences of mitochondrial DNA control regions from domestic and wild herds from Eurasia and China were analysed. The results showed that the Chinese reindeer population originated independently from north-central Russian domestic herds, belonging to a large reindeer population that was present across Beringia during the last glacial period. Some studies have reported that the Chinese reindeer population is closely related to wild forest reindeer herds in Russia. Our results, however, indicate that wild forest reindeer herds of southeastern Russia contributed little or nothing to the Chinese reindeer herd gene pool. Chinese reindeer herds have a much greater genetic similarity to more northerly distributed tundra-type herds that inhabit open areas. The present findings will be essential for future conservation planning for Chinese reindeer.

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.