Genome-wide study on genetic diversity and phylogeny of five species in the genus Cervus

Evolutionary insights on critically endangered Kashmir red deer or hangul (Cervus hanglu hanglu) through a mitogenomic lens

Background

The Kashmir red deer or Hangul (Cervus hanglu hanglu) is the only Tarim red deer species endemic to India. With a current estimated population size of fewer than 200 individuals, this critically endangered species is confined to the greater Dachigam landscape in Jammu and Kashmir. Poaching, habitat loss and fragmentation, resource competition with livestock, and small population size are the major conservation challenges for this species.

Methods

Blood sampling was conducted from two wild Hangul individuals during radio-collaring operations at Dachigam National Park, Kashmir in 2013 and 2020, respectively. Using next-generation sequencing approach, we sequenced the 16,351 bp long mitogenome of two wild-caught Hangul individuals (1 M:1 F at ~14× and ~10× coverage, respectively) from Dachigam National Park.

Results

The annotated sequences were identical with an AT-rich composition, including 13 protein-coding genes (11,354 bp), 22 tRNA genes (1,515 bp), two ribosomal genes (2,526 bp) and a non-coding control region (917 bp) in a conserved order like other red deer species. Bayesian phylogenetic reconstruction of the red deer complex revealed two major groups: the elaphoid and the wapitoid clades. Hangul formed a distinct clade with its other subspecies C. hanglu yarkandensis and is sister to the Hungarian red deer (C. elaphus hippelaphus). Divergence time analyses suggested that the Tarim deer species group separated ~1.55 Mya from their common ancestors and Hangul diverged ~0.75 Mya from closely related C. yarkandensis, corroborating with the known paleobiogeographic events related to refugia during glaciations in the Pleistocene era. This study provides baseline information on Hangul mitogenome for further research on phylogeography and other population parameters and helps in developing suitable conservation plans for this species.

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.

Antler tine homologies and cervid systematics: A review of past and present controversies with special emphasis on Elaphurus davidianus

Antlers are the most conspicuous trait of cervids and have been used in the past to establish a classification of their fossil and living representatives. Since the availability of molecular data, morphological characters have generally become less important for phylogenetic reconstructions. In recent years, however, the appreciation of morphological characters has increased, and they are now more frequently used in addition to molecular data to reconstruct the evolutionary history of cervids. A persistent challenge when using antler traits in deer systematics is finding a consensus on the homology of structures. Here, we review early and recent attempts to homologize antler structures and objections to this approach, compare and evaluate recent advances on antler homologies, and critically discuss these different views in order to offer a basis for further scientific exchange on the topic. We further present some developmental aspects of antler branching patterns and discuss their potential for reconstructing cervid systematics. The use of heterogeneous data for reconstructing phylogenies has resulted in partly conflicting hypotheses on the systematic position of certain cervid species, on which we also elaborate here. We address current discussions on the use of different molecular markers in cervid systematics and the question whether antler morphology and molecular data can provide a consistent picture on the evolutionary history of cervids. In this context, special attention is given to the antler morphology and the systematic position of the enigmatic Pere David's deer (Elaphurus davidianus).

Phylogeny and evolution of the genus Cervus (Cervidae, Mammalia) as revealed by complete mitochondrial genomes

Mitochondrial DNA (mtDNA) lineages are recognized as important components of intra- and interspecific biodiversity, and allow to reveal colonization routes and phylogeographic structure of many taxa. Among these is the genus Cervus that is widely distributed across the Holarctic. We obtained sequences of complete mitochondrial genomes from 13 Cervus taxa and included them in global phylogenetic analyses of 71 Cervinae mitogenomes. The well-resolved phylogenetic trees confirmed Cervus to be monophyletic. Molecular dating based on several fossil calibration points revealed that ca. 2.6 Mya two main mitochondrial lineages of Cervus separated in Central Asia, the Western (including C. hanglu and C. elaphus) and the Eastern (comprising C. albirostris, C. canadensis and C. nippon). We also observed convergent changes in the composition of some mitochondrial genes in C. hanglu of the Western lineage and representatives of the Eastern lineage. Several subspecies of C. nippon and C. hanglu have accumulated a large portion of deleterious substitutions in their mitochondrial protein-coding genes, probably due to drift in the wake of decreasing population size. In contrast to previous studies, we found that the relic haplogroup B of C. elaphus was sister to all other red deer lineages and that the Middle-Eastern haplogroup E shared a common ancestor with the Balkan haplogroup C. Comparison of the mtDNA phylogenetic tree with a published nuclear genome tree may imply ancient introgressions of mtDNA between different Cervus species as well as from the common ancestor of South Asian deer, Rusa timorensis and R. unicolor, to the Cervus clade.

Identification and Characterization of Alternative Splicing Variants and Positive Selection Genes Related to Distinct Growth Rates of Antlers Using Comparative Transcriptome Sequencing

Simple Summary

The size of antlers varies among species; antlers of the wapiti (Cervus canadensis xanthopygus) grow much faster than those of its close relative the sika deer (Cervus nippon hortulorum) in the same growing period. This contrast provides a potential model for comparative studies for the identification of potent growth factors and unique regulatory systems. In the present study, the reference transcriptomes of the antler reserve mesenchyme (RM) tissue of wapiti and sika deer were constructed using single molecule real time sequencing data. The expression profiling, positive selection, and alternative splicing of the antler transcripts were compared, and interactive relationships and expression patterns of hub genes were identified and analysed. We identified that RNA Binding Motif Protein X-Linked (RBMX) gene was under strongly positive selection. One gene found to interact with RBMX was methyltransferase-like 3 (METTL3), an oncogene that could promote translation of cancer cell proteins. There was a contrasting relationship in expression level between RBMX and METTL3 genes in the RM tissue. We believe our study can provide a better understanding of rapid antler growth at the molecular level in particular and endochondral ossification in general.

Abstract

The molecular mechanism underlying rapid antler growth has not been elucidated. The contrast of the wapiti and sika deer antler provides a potential model for comparative studies for the identification of potent growth factors and unique regulatory systems. In the present study, reference transcriptomes of antler RM tissue of wapiti and sika deer were constructed using single molecule real time sequencing data. The expression profiling, positive selection, and alternative splicing of the antler transcripts were compared. The results showed that: a total of 44,485 reference full-length transcripts of antlers were obtained; 254 highly expressed transcripts (HETs) and 1936 differentially expressed genes (DEGs) were enriched and correlated principally with translation, endochondral ossification and ribosome; 228 genes were found to be under strong positive selection and would thus be important for the evolution of wapiti and sika deer; among the alternative splicing variants, 381 genes were annotated; and 4 genes with node degree values greater than 50 were identified through interaction network analysis. We identified a negative and a positive regulator for rapid antler growth, namely RNA Binding Motif Protein X-Linked (RBMX) and methyltransferase-like 3 (METTL3), respectively. Overall, we took advantage of this significant difference in growth rate and performed the comparative analyses of the antlers to identify key specific factors that might be candidates for the positive or negative regulation of phenomenal antler growth rate.

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

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 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 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.

Integrative Analyses of Antler Cartilage Transcriptome and Proteome of Gansu Red Deer (Cervus elaphus kansuensis) at Different Growth Stages

Simple Summary

Velvet antler is the only organ in mammals that can completely and circularly regenerate, which involves the co-development of a variety of tissues including cartilage. Thus, velvet antler can provide an ideal model for studying chondrogenesis, endochondral ossification and rapid tissue growth. However, the mechanism of rapid growth and regeneration of velvet antler is still unclear. In this study, we conducted integrated analysis of the transcriptome and proteome of antler cartilage tissues at different growth stages. The results showed that gene13546 and its coding protein rna13546 were annotated in Wnt signaling pathway. They may play roles in antler rapid growth and regeneration.

Abstract

The velvet antler is a unique model for cancer and regeneration research due to its periodic regeneration and rapid growth. Antler growth is mainly triggered by the growth center located in its tip, which consists of velvet skin, mesenchyme and cartilage. Among them, cartilage accounts for most of the growth center. We performed an integrative analysis of the antler cartilage transcriptome and proteome at different antler growth stages. RNA-seq results revealed 24,778 unigenes, 19,243 known protein-coding genes, and 5535 new predicted genes. Of these, 2722 were detected with differential expression patterns among 30 d, 60 d, and 90 d libraries, and 488 differentially expressed genes (DEGs) were screened at 30 d vs. 60 d and 60 d vs. 90 d but not at 30 d vs. 90 d. Proteomic data identified 1361 known proteins and 179 predicted novel proteins. Comparative analyses showed 382 differentially expressed proteins (DEPs), of which 16 had differential expression levels at 30 d vs. 60 d and 60 d vs. 90 d but not at 30 d vs. 90 d. An integrated analysis conducted for DEGs and DEPs showed that gene13546 and its coding protein protein13546 annotated in the Wnt signaling pathway may possess important bio-logical functions in rapid antler growth. This study provides in-depth characterization of candidate genes and proteins, providing further insights into the molecular mechanisms controlling antler development.

Revealing the coexistence of differentiation and communication in an endemic hare, Lepus yarkandensis (Mammalia, Leporidae) using specific-length amplified fragment sequencing

BACKGROUND

The Yarkand hare (Lepus yarkandensis Günther, 1875) is endemic to oasis and desert areas around the Tarim Basin in the Xinjiang Uyghur Autonomous Region of northwest China; however, genome-wide information for this species remains limited. Moreover, the genetic variation, genetic structure, and phylogenetic relationships of Yarkand hare from the plateau mountain regions have not been reported. Thus, we used specific-length amplified fragment sequencing (SLAF-seq) technology to evaluate the genetic diversity of 76 Yarkand hares from seven geographic populations in the northern and southwestern parts of the Tarim Basin to investigate single-nucleotide polymorphism (SNP) marker-based population differentiation and evolutionary processes. Selective sweep analysis was conducted to identify genetic differences between populations.

RESULTS

Using SLAF-seq, a total of 1,835,504 SNPs were initially obtained, of which 308,942 high-confidence SNPs were selected for further analysis. Yarkand hares exhibited a relatively high degree of genetic diversity at the SNP level. Based on pairwise FST estimates, the north and southwest groups showed a moderate level of genetic differentiation. Phylogenetic tree and population structure analyses demonstrated evident systematic phylogeographical structure patterns consistent with the geographical distribution of the hares. Hierarchical analysis of molecular variation further indicated that genetic variation was mainly observed within populations. Low to moderate genetic differentiation also occurred among populations despite a common genomic background, likely due to geographical barriers, genetic drift, and differential selection pressure of distinct environments. Nevertheless, the observed lineage-mixing pattern, as indicated by the evolutionary tree, principal component analysis, population structure, and TreeMix analyses, suggests a certain degree of gene flow between the north and southwest groups. This may be related to the migration of hares to high-altitude water sources southwest of the basin during glacial climatic oscillations, as well as river re-diffusion and oasis restoration in the basin following the glacial period. We also identified candidate genes, and their associated gene ontology terms and pathways, related to the adaptation of Yarkand hares to different environmental habitats.

CONCLUSIONS

The identified genome-wide SNPs, genetic diversity, and population structure of Yarkand hares expand our understanding of the genetic background of this endemic species and provide valuable insights into its environmental adaptation, allowing for further exploration of the underlying mechanisms.

Experimental oral transmission of chronic wasting disease to sika deer (Cervus nippon)

Chronic wasting disease (CWD) affects a broad array of cervid species and continues to be detected in an expanding geographic range. Initially introduced into the Republic of Korea through the importation of CWD-infected elk (Cervus canadensis), additional cases of CWD were subsequently detected in farmed Korean elk and sika deer (Cervus nippon). Wild and farmed sika deer are found in many regions of Asia, North America, and Europe, although natural transmission to this species has not been detected outside of the Republic of Korea. In this study, the oral transmission of CWD to sika deer was investigated using material from CWD-affected elk. Pathological prion (PrP^CWD) immunoreactivity was detected in oropharyngeal lymphoid tissues of one sika deer at 3.9 months post-inoculation (mpi) and was more widely distributed in a second sika deer examined at 10.9 mpi. The remaining four sika deer progressed to clinical disease between 21 and 24 mpi. Analysis of PrP^CWD tissue distribution in clinical sika deer revealed widespread deposition in central and peripheral nervous systems, lymphoreticular tissues, and the gastrointestinal tract. Prion protein gene (PRNP) sequences of these sika deer were identical and consistent with those reported in natural sika deer populations. These findings demonstrate the efficient oral transmission of CWD from elk to sika deer.

Whole-genome sequencing of Tarim red deer (Cervus elaphus yarkandensis) reveals demographic history and adaptations to an arid-desert environment

Background

The initiation of desert conditions in the Tarim Basin in China since the late Miocene has led to the significant genetic structuring of local organisms. Tarim Red Deer (Cervus elaphus yarkandensis, TRD) have adapted to the harsh environmental conditions in this basin, including high solar radiation and temperature, aridity, and poor nutritional conditions. However, the underlying genetic basis of this adaptation is poorly understood.

Results

We sequenced the whole genomes of 13 TRD individuals, conducted comparative genomic analyses, and estimated demographic fluctuation. The ∂a∂i model estimated that the TRD and Tule elk (Cervus canadensis nannodes) populations diverged approximately 0.98 Mya. Analyses revealed a substantial influence of the Earth’s climate on the effective population size of TRD, associated with glacial advances and retreat, and human activities likely underlie a recent serious decline in population. A marked bottleneck may have profoundly affected the genetic diversity of TRD populations. We detected a set of candidate genes, pathways, and GO categories related to oxidative stress, water reabsorption, immune regulation, energy metabolism, eye protection, heat stress, respiratory system adaptation, prevention of high blood pressure, and DNA damage and repair that may directly or indirectly be involved in the adaptation of TRD to an arid-desert environment.

Conclusions

Our analyses highlight the role of historical global climates in the population dynamics of TRD. In light of ongoing global warming and the increasing incidence of droughts, our study offers insights into the genomic adaptations of animals, especially TRD, to extreme arid-desert environments and provides a valuable resource for future research on conservation design and biological adaptations to environmental change.

The systematics of the Cervidae: a total evidence approach

Systematic relationships of cervids have been controversial for decades. Despite new input from molecular systematics, consensus could only be partially reached. The initial, gross (sub) classification based on morphology and comparative anatomy was mostly supported by molecular data. The rich fossil record of cervids has never been extensively tested in phylogenetic frameworks concerning potential systematic relationships of fossil cervids to extant cervids. The aim of this work was to investigate the systematic relationships of extant and fossil cervids using molecular and morphological characters and make implications about their evolutionary history based on the phylogenetic reconstructions. To achieve these objectives, molecular data were compiled consisting of five nuclear markers and the complete mitochondrial genome of 50 extant and one fossil cervids. Several analyses using different data partitions, taxon sampling, partitioning schemes, and optimality criteria were undertaken. In addition, the most extensive morphological character matrix for such a broad cervid taxon sampling was compiled including 168 cranial and dental characters of 41 extant and 29 fossil cervids. The morphological and molecular data were analysed in a combined approach and other comprehensive phylogenetic reconstructions. The results showed that most Miocene cervids were more closely related to each other than to any other cervids. They were often positioned between the outgroup and all other cervids or as the sister taxon to Muntiacini. Two Miocene cervids were frequently placed within Muntiacini. Plio- and Pleistocene cervids could often be affiliated to Cervini, Odocoileini or Capreolini. The phylogenetic analyses provide new insights into the evolutionary history of cervids. Several fossil cervids could be successfully related to living representatives, confirming previously assumed affiliations based on comparative morphology and introducing new hypotheses. New systematic relationships were observed, some uncertainties persisted and resolving systematics within certain taxa remained challenging.