The first high-quality reference genome of sika deer provides insights for high-tannin adaptation

Population genomics of sika deer reveals recent speciation and genetic selective signatures during evolution and domestication

BACKGROUND

Population genomic analysis can reconstruct the phylogenetic relationship and demographic history, and identify genomic selective signatures of a species. To date, fundamental aspects of population genomic analyses, such as intraspecies taxonomy, evolutionary history, and adaptive evolution, of sika deer have not been systematically investigated. Furthermore, accumulating lines of evidences have illustrated that incorrect species delimitation will mislead conservation decisions, and even lead to irreversible mistakes in threatened species.

RESULTS

In this study, we resequenced 81 wild and 71 domesticated sika deer representing 10 main geographic populations and two farms to clarify the species delimitation, demographic and divergence histories, and adaptive evolution of this species. First, our analyses of whole genomes, Y chromosomes and mitochondrial genomes revealed substantial genetic differentiation between the continental and Japanese lineages of sika deer, representing two phylogenetically distinct species. Second, sika deer in Japan were inferred to have experienced a "divergence-mixing-isolation" evolutionary scenario. Third, we identified four candidate genes (XKR4, NPAS3, CTNNA3, and CNTNAP5) possibly involved in body size regulation of sika deer by selective sweep analysis. Furthermore, we also detected two candidate genes (NRP2 and EDIL3) that may be associated with an important economic trait (antler weight) were under selection during the process of domestication.

CONCLUSION

Population genomic analyses revealed that the continental and Japanese lineages represent distinct phylogenetic species. Moreover, our results provide insights into the genetic selection signatures related to body size differences and a valuable genomic resource for future genetic studies and genomics-informed breeding of sika deer.

DNA Metabarcoding Illuminates Seasonal Dietary Pattern and Niche Partitioning by Three Sympatric Herbivores

Diet composition is among the most critical dimensions of animal ecology, yet seasonal dietary diversity has rarely been investigated in sympatric herbivores. This study used DNA metabarcoding to conduct an analysis of seasonal variations in diet composition and trophic niches for sympatric sika deer, Reeves' muntjac, and Chinese hare in Taohongling National Nature Reserve (TNNR). The results showed that Smilax china (11.79%) was the leading food eaten by sika deer in summer, whereas dominated by Rubus spp. (36.42%) and Loropetalum chinense (25.48%) in winter; Rubus spp. accounted for the majority of Reeves' muntjac's diet throughout the year. In comparison, the Chinese hare primarily consumed Smilax china from winter to spring but changed to Poa annua (10.81%) and Setaria viridis (23.05%) in summer and fall. Compared to other seasons, significant differences (Shannon index, p < 0.05) occurred in spring and summer, showing higher diversity of food items across the three herbivorous. Nonmetric multidimensional scaling (NMDS) analysis suggested significant partitioning in the food items of Chinese hares compared to the two ruminants. Both sika deer and Reeves' muntjac occupied a wider niche breadth and dietary diversity in summer, reflecting generalised feeding habits (Sd: Ba = 0.06; Rm: Ba = 0.04) and lower in fall (Ba = 0.01) with stronger selectivity and specialization, which was consistent with the optimal foraging theory. Notably, no significant difference was indicated in seasonal niche breadth for Chinese hare (p > 0.05). The niche overlap indices were 0.989 (fall) and 0.831 (winter) between sika deer and Reeves' muntjac, indicating a higher dietary similarity and overlap. However, differences in foraging plant taxa and abundance ratios may facilitate dietary niche partitioning. The diet of herbivores reflected plant-herbivore interactions and seasonal diet differences were correlated with feeding strategies, which facilitate coexistence and reduce competition of co-occurring species in the food dimension.

Single-Cell Analysis of Molecular Mechanisms in Rapid Antler Osteogenesis During Growth and Ossification Stages

Antlers, as the only fully regenerable bone tissue in mammals, serve as an exceptional model for investigating bone growth, mineralization, articular cartilage repair, and the pathophysiology of osteoporosis. Nevertheless, the exact molecular mechanisms governing osteogenesis, particularly the dynamic cellular interactions and signaling pathways coordinating these processes, remain poorly characterized. This study used single-cell RNA sequencing (scRNA-seq) on the 10× Genomics Chromium platform, combined with bulk-RNA sequencing results, to comprehensively analyze molecular regulatory mechanisms in rapid antler osteogenesis. The results showed that eight cell types were identified in sika deer antler during the growth and ossification stages: mesenchymal, chondrocyte, osteoblast, pericyte, endothelial, monocyte/macrophage, osteoclast, and NK cells. Chondrocytes were predominantly found during the growth stage, while osteoblasts were more abundant during the ossification stage. Mesenchymal cells were subclassified into three subcategories: MSC_1 (VCAN and SFRP2), MSC_2 (TOP2A, MKI67), and MSC_3 (LYVE1 and TNN). MSC_3 was predominantly present during the growth stage. During the growth stage, MSC_1 and MSC_2 upregulated genes related to vasculature development (COL8A1, NRP1) and cell differentiation (PTN, SFRP2). During the ossification stage, these subcategories upregulated genes involved in the positive regulation of p53 class mediator signal transduction (RPL37, RPL23, RPS20, and RPL26), osteoblast differentiation (SPP1, IBSP, BGLAP), and proton-motive ATP synthesis (NDUFA7, NDUFB3, NDUFA3, NDUFB1). Endothelial cells were categorized into five subpopulations: Enc_1 (SPARCL1, VWF), Enc_2 (MCM5), Enc_3 (ASPM, MKI67), Enc_4 (SAT1, CXCL12), and Enc_5 (ZFHX4, COL6A3). Combined scRNA-seq and bulk RNA-seq analysis revealed that the ossification stage's upregulation genes included osteoclast- and endothelial cell-specific genes, while the growth stage's upregulation genes were mainly linked to collagen organization, osteoblast differentiation, mitotic cell cycle, and chondrocyte differentiation. Overall, this study offers a detailed single-cell analysis of gene expression patterns in antlers during the growth and ossification stages, providing insights into the molecular mechanisms driving rapid osteogenesis.

Draft genome of the endangered visayan spotted deer (Rusa alfredi), a Philippine endemic species

The Visayan Spotted Deer (VSD), or Rusa alfredi, is an endangered and endemic species in the Philippines. Despite its status, genomic information on R. alfredi, and the genus Rusa in general, is missing. This study presents the first draft genome assembly of the VSD using the Illumina short-read sequencing technology. The resulting RusAlf_1.1 assembly has a 2.52 Gb total length, with a contig N50 of 46 Kb and scaffold N50 size of 75 Mb. The assembly has a BUSCO complete score of 95.5%, demonstrating the genome's completeness, and includes the annotation of 24,531 genes. Our phylogenetic analysis based on single-copy orthologs revealed a close evolutionary relationship between R. alfredi and the genus Cervus. RusAlf_1.1 represents a significant advancement in our understanding of the VSD. It opens opportunities for further research in population genetics and evolutionary biology, potentially contributing to more effective conservation and management strategies for this endangered species.

BRCA1 is involved in sustaining rapid antler growth possibly via balancing of the p53/endoplasmic reticulum stress signaling pathway

BACKGROUND

Regeneration is the preferred approach to restore the structure and function after tissue damage. Rapid proliferation of cells over the site of damage is integral to the process of regeneration. However, even subtle mutations in proliferating cells may cause detrimental effects by eliciting abnormal differentiation. Interestingly deer antlers, arguably the fastest regenerating mammalian tissue, have not been reported, thus far, to grow malignant tumors. They provide a mammalian model to understand the possible mechanism by which rapid regeneration is achieved while avoiding the development of malignancies. Antler regeneration is based on the proliferation and differentiation of antler stem cells (AnSCs).

RESULTS

We identified 39 hub genes which may function in regulating the balance between rapid proliferation and genomic stability in the AnSCs during antler regeneration. Among these 39 genes, the tumor suppressor gene, BRCA1, was found to be more sensitive to DNA damage in the AnSCs compared to that in the deer somatic cells, and BRCA1 deletion in the AnSCs via CRISPR/Cas9 resulted in significantly higher levels of DNA damage. Lack of BRCA1 promoted cell apoptosis and cell senescence and inhibited cell proliferation and cell self-renewal. RNA-seq results showed that in the absence of BRCA1, the p53 signaling pathway was significantly up-regulated. Associated with this change, the cell apoptosis and cell senescence-relevant-genes, CDKN1A, CDKN2A and Fas were over expressed, but the expression of cell-cycle-progression-related genes was inhibited. In addition, BRCA1 expression levels were found to be more sensitive to endoplasmic reticulum stress (ERS) in the AnSCs compared to the somatic cells. Deletion of BRCA1 gene aggravated ERS and ERS-induced cell apoptosis.

CONCLUSIONS

Our results revealed that BRCA1 is involved in sustaining rapid antler growth possibly via promotion of DNA damage repair that acts to maintain genome stability while protecting cells from p53/ERS-induced cell death. Understanding the mechanisms underlying the role played by BRCA1 in the process of antler regeneration is of great significance not only for regenerative medicine, but also for the understanding of cancer development.

High-quality sika deer omics data and integrative analysis reveal genic and cellular regulation of antler regeneration

The antler is the only organ that can fully regenerate annually in mammals. However, the regulatory pattern and mechanism of gene expression and cell differentiation during this process remain largely unknown. Here, we obtain comprehensive assembly and gene annotation of the sika deer (Cervus nippon) genome. We construct, together with large-scale chromatin accessibility and gene expression data, gene regulatory networks involved in antler regeneration, identifying four transcription factors, MYC, KLF4, NFE2L2, and JDP2, with high regulatory activity across the whole regeneration process. Comparative studies and luciferase reporter assay suggest the MYC expression driven by a cervid-specific regulatory element might be important for antler regenerative ability. We further develop a model called combinatorial TF Oriented Program (cTOP), which integrates single-cell data with bulk regulatory networks and find PRDM1, FOSL1, BACH1, and NFATC1 as potential pivotal factors in antler stem cell activation and osteogenic differentiation. Additionally, we uncover interactions within and between cell programs and pathways during the regeneration process. These findings provide insights into the gene and cell regulatory mechanisms of antler regeneration, particularly in stem cell activation and differentiation.

Parallel gene expansions drive rapid dietary adaptation in herbivorous woodrats

How mammalian herbivores evolve to feed on chemically defended plants remains poorly understood. In this study, we investigated the adaptation of two species of woodrats (Neotoma lepida and N. bryanti) to creosote bush (Larrea tridentata), a toxic shrub that expanded across the southwestern United States after the Last Glacial Maximum. We found that creosote-adapted woodrats have elevated gene dosage across multiple biotransformation enzyme families. These duplication events occurred independently across species and substantially increase expression of biotransformation genes, especially within the glucuronidation pathway. We propose that increased gene dosage resulting from duplication is an important mechanism by which animals initially adapt to novel environmental pressures.

Insights from homozygous signatures of cervus nippon revealed genetic architecture for components of fitness

This study investigates the genomic landscape of Sika deer populations, emphasizing the detection and characterization of runs of homozygosity (ROH) and their contribution towards components of fitness. Using 85,001 high-confidence SNPs, the investigation into ROH distribution unveiled nuanced patterns of autozygosity across individuals especially in 2 out of the 8 farms, exhibiting elevated ROH levels and mean genome coverage under ROH segments. The prevalence of shorter ROH segments (0.5-4 Mb) suggests historical relatedness and potential selective pressures within these populations. Intriguingly, despite observed variations in ROH profiles, the overall genomic inbreeding coefficient (FROH) remained relatively low across all farms, indicating a discernible degree of genetic exchange and effective mitigation of inbreeding within the studied Sika deer populations. Consensus ROH (cROH) were found to harbor genes for important functions viz., EGFLAM gene which is involved in the vision function of the eye, SKP2 gene which regulates cell cycle, CAPSL involved in adipogenesis, SPEF2 which is essential for sperm flagellar assembly, DCLK3 involved in the heat stress. This first ever study on ROH in Sika deer, to shed light on the adaptive role of genes in these homozygous regions. The insights garnered from this study have broader implications in the management of genetic diversity in this vulnerable species.

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.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Comparative Genomics Reveal Phylogenetic Relationship and Chromosomal Evolutionary Events of Eight Cervidae Species

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.

Analysis of Whole-Genome as a Novel Strategy for Animal Species Identification

Survival crises stalk many animals, especially endangered and rare animals. Accurate species identification plays a pivotal role in animal resource conservation. In this study, we developed an animal species identification method called Analysis of whole-GEnome (AGE), which identifies species by finding species-specific sequences through bioinformatics analysis of the whole genome and subsequently recognizing these sequences using experimental technologies. To clearly demonstrate the AGE method, Cervus nippon, a well-known endangered species, and a closely related species, Cervus elaphus, were set as model species, without and with published genomes, respectively. By analyzing the whole genomes of C. nippon and C. elaphus, which were obtained through next-generation sequencing and online databases, we built specific sequence databases containing 7,670,140 and 570,981 sequences, respectively. Then, the species specificities of the sequences were confirmed experimentally using Sanger sequencing and the CRISPR-Cas12a system. Moreover, for 11 fresh animal samples and 35 commercially available products, our results were in complete agreement with those of other authoritative identification methods, demonstrating AGE's precision and potential application. Notably, AGE found a mixture in the 35 commercially available products and successfully identified it. This study broadens the horizons of species identification using the whole genome and sheds light on the potential of AGE for conserving animal resources.

Haplotype-resolved Genome of Sika Deer Reveals Allele-specific Gene Expression and Chromosome Evolution

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.

Diet composition and interspecific niche of Taohongling Sika deer (Cervus nippon kopschi) and its sympatric Reeve’s muntjac (Muntiacus reevesi) and Chinese hare (Lepus sinensis) in winter (Animalia, Mammalia)

Species co-existence depends on how organisms utilize their environment and resources. Little is known about the winter diet composition and sympatric co-existence of South China sika deer and its companion species in Taohongling. In this study, high-throughput sequencing and metabarcoding trnL were used to study the diet composition and interspecific relationship including sika deer, Reeve’s muntjac, and Chinese hare. Our results show that 203 genera in 90 families are contained in the diet of sika deer, 203 genera in 95 families for Reeve’s muntjac, and 163 genera in 75 families for Chinese hare. Sika deer fed on Rubus chingii, Loropetalum chinense, and Eurya japonica in winter, accounting for 75.30%; Reeve’s muntjac consumed mainly R. chingii, E. japonica, and Euonymus grandiflorus, accounting for 68.80%, and Chinese hare mainly fed on R. chingii, Smilax china, and Rhus chinensis, accounting for 41.98%. The Shannon index showed no significant difference between groups (p &gt; 0.05). The NMDS analysis found considerable overlap among three species. Sika deer and Reeve’s muntjac consumed similar forage plants but varied greatly in Chinese hare, which occupied the widest choice in winter, resulting in higher diet breadth and increased dietary divergence, thereby reducing competition and facilitating coexistence. The diet niche overlap index among them, as represented by Pianka’s index, ranging from 0.62 between sika deer and Chinese hare to 0.83 between sika deer and Reeve’s muntjac, which indicated a more similar niche and potential competition in closely related species. Our findings provide a new diet perspective of three herbivores, leading to a more comprehensive understanding of resource partitioning and species coexistence.

Duplication, Loss, and Evolutionary Features of Specific UDP-Glucuronosyltransferase Genes in Carnivora (Mammalia, Laurasiatheria)

UDP-glucuronosyltransferases (UGTs) are one of the most important enzymes for xenobiotic metabolism or detoxification. Through duplication and loss of genes, mammals evolved the species-specific variety of UGT isoforms. Among mammals, Carnivora is one of the orders that includes various carnivorous species, yet there is huge variation of food habitat. Recently, lower activity of UGT1A and 2B were shown in Felidae and pinnipeds, suggesting evolutional loss of these isoforms. However, comprehensive analysis for genetic or evolutional features are still missing. This study was conducted to reveal evolutional history of UGTs in Carnivoran species. We found specific gene expansion of UGT1As in Canidae, brown bear and black bear. We also found similar genetic duplication in UGT2Bs in Canidae, and some Mustelidae and Ursidae. In addition, we discovered contraction or complete loss of UGT1A7-12 in phocids, some otariids, felids, and some Mustelids. These studies indicate that even closely related species have completely different evolution of UGTs and further imply the difficulty of extrapolation of the pharmacokinetics and toxicokinetic result of experimental animals into wildlife carnivorans.