Sika deer are known to prefer oak leaves, which are rich in tannins and toxic to most mammals; however, the genetic mechanisms underlying their unique ability to adapt to living in the jungle are still unclear. In identifying the mechanism responsible for the tolerance of a highly toxic diet, we have made a major advancement by explaining the genome of sika deer. We generated the first high-quality, chromosome-level genome assembly of sika deer and measured the correlation between tannin intake and RNA expression in 15 tissues through 180 experiments. Comparative genome analyses showed that the UGT and CYP gene families are functionally involved in the adaptation of sika deer to high-tannin food, especially the expansion of the UGT family 2 subfamily B of UGT genes. The first chromosome-level assembly and genetic characterization of the tolerance to a highly toxic diet suggest that the sika deer genome may serve as an essential resource for understanding evolutionary events and tannin adaptation. Our study provides a paradigm of comparative expressive genomics that can be applied to the study of unique biological features in non-model animals.

Cervidae represents a family that is not only rich in species diversity but also exhibits a wide range of karyotypes. The controversies regarding the phylogeny and classification of Cervidae still persist. The flourishing development of the genomic era has made it possible to address these issues at the genomic level. Here, the genomes of nine species were used to explore the phylogeny and chromosomal evolutionary events of Cervidae. By conducting whole-genome comparisons, we identified single-copy orthologous genes across the nine species and constructed a phylogenetic tree based on the single-copy orthologous genes sequences, providing new insights into the phylogeny of Cervidae, particularly the phylogenetic relationship among sika deer, red deer, wapiti and Tarim red deer. Gene family analysis revealed contractions in the olfactory receptor gene family and expansions in the histone gene family across eight Cervidae species. Furthermore, synteny analysis was used to explore the chromosomal evolutionary events of Cervidae species, revealing six chromosomal fissions during the evolutionary process from Bovidae to Cervidae. Notably, specific chromosomal fusion events were found in four species of Cervus, and a unique chromosomal fusion event was identified in Muntiacus reevesi. Our study further completed the phylogenetic relationship within the Cervidae and demonstrated the feasibility of inferring species phylogeny at the whole-genome level. Additionally, our findings on gene family evolution and the chromosomal evolutionary events in eight Cervidae species lay a foundation for comprehensive research of the evolution of Cervidae.

Velvet deer are not only a representative special economic animal but also an important part of livestock. With the increasing awareness of international competition for germplasm resources in China, more and more attention has been paid to the protection and utilization of germplasm resources. However, there is poor understanding about velvet deer resources. Therefore, we are providing a comprehensive introduction of Chinese velvet deer germplasm resources from the aspects of ecological distribution, domestication and breeding.

Although Cervus elaphus (Linnaeus, 1758) has been well studied, the subspecific taxonomy of Cervus canadensis populations in Qinghai and Gansu, China, is still controversial, and the mitochondrial characteristics of Cervus elaphus (Linnaeus, 1758) remain incompletely understood. We assembled 89 mitogenomes of C. canadensis from five geographical populations across Qinghai and Gansu. Phylogenetic analysis confirmed that the 89 individuals are taxonomically classified as C. c. kansuensis. Nucleotide compositions showed a higher abundance of adenine and cytosine compared to guanine and thymine in both complete mitogenomes and mitochondrial PCGs. Codon usage analysis revealed a strong preference towards A-ending codons (68.04% of over-represented codons, RSCU > 1.6) in mitochondrial PCGs, with systemic avoidance of G-ending codons (53.30% of unused codons, RSCU = 0). The CAMs of 13 PCGs are reported for the first time. Furthermore, the ENC plot showed that the codon usage of all PCGs was biased except for gene ATP8. The PR2 bias plot showed that gene ND6 exhibited bias towards T3 and G3, whereas the other genes preferred A3 and C3. Both the ENC-plot and PR2 bias plot suggested that natural selection played an important role in the forces driving codon usage bias in mitochondrial PCGs. Our results demonstrate the subspecific status of C. canadensis distributed in Qinghai and Gansu as C. c. kansuensis, and provide insights into the mitochondrial characteristics of C. c. kansuensis. The mitogenome sequences assembled in this study provide valuable data for further understanding of the Cervus elaphus (Linnaeus, 1758) mitogenome.

Sika deer are known to prefer oak leaves, which are rich in tannins and toxic to most mammals; however, the genetic mechanisms underlying their unique ability to adapt to living in the jungle are still unclear. In identifying the mechanism responsible for the tolerance of a highly toxic diet, we have made a major advancement in the elucidation of the genomics of sika deer. We generated the first high-quality, chromosome-level genome assembly of sika deer and measured the correlation between tannin intake and RNA expression in 15 tissues through 180 experiments. Comparative genome analyses showed that the UGT and CYP gene families are functionally involved in the adaptation of sika deer to high-tannin food, especially the expansion of UGT genes in a subfamily. The first chromosome-level assembly and genetic characterization of the tolerance toa highly toxic diet suggest that the sika deer genome will serve as an essential resource for understanding evolutionary events and tannin adaptation. Our study provides a paradigm of comparative expressive genomics that can be applied to the study of unique biological features in non-model animals.

Both the productive efficiency and physical health of sika deer are strongly linked to their body conformation phenotypes. Breeding sika deer with excellent growth traits using molecular breeding technologies has become essential. Phenotypic data of 12 body conformation traits were measured for 613 sika deer across three age groups, including body traits and cephalic traits. Genetic typing was performed using the Sika Deer 100K SNP Liquid Chip, and genetic parameters were estimated through animal models to obtain the heritability and genetic correlation of traits within each age group. Single-trait and multi-trait GWASs were conducted, using GEMMA software to identify gene variants significantly associated with sika deer body conformation traits. Most of the 12 body conformation traits exhibited moderate to high heritability. The single-trait GWAS identified 49 SNPs (38 candidate genes), while multi-trait GWAS detected 134 SNPs (80 candidate genes), including 114 novel loci. A total of 163 SNPs and 196 candidate genes were identified, with 17 genes detected by both methods. These genes may be involved in pain perception, cell cycle regulation, immune response, and protein ubiquitination. Two significant KEGG pathways were enriched: steroid hormone biosynthesis and drug metabolism–cytochrome P450. Collectively, these detected loci and genes may serve as potential genetic resources for marker-assisted breeding, contributing to subsequent genetic improvement of body conformation in sika deer.