Search

The annual regrowth of deer antlers provides a valuable model for studying organ regeneration in mammals. We describe a single-cell atlas of antler regrowth. The earliest-stage antler initiators were mesenchymal cells that express the paired related homeobox 1 gene ( PRRX1 + mesenchymal cells). We also identified a population of “antler blastema progenitor cells” (ABPCs) that developed from the PRRX1 + mesenchymal cells and directed the antler regeneration process. Cross-species comparisons identified ABPCs in several mammalian blastema. In vivo and in vitro ABPCs displayed strong self-renewal ability and could generate osteochondral lineage cells. Last, we observed a spatially well-structured pattern of cellular and gene expression in antler growth center during the peak growth stage, revealing the cellular mechanisms involved in rapid antler elongation.

Context Deer antlers offer a premium model for investigating the mechanisms underlying arguably the most rapid cartilage formation and remodelling system. Although the cartilage formation process in the antler has been relatively intensively studied, thus far, at a molecular level, the cartilage remodelling has essentially been untouched. Aims To construct miRNA–mRNA regulatory networks for both the cartilage formation and remodel zones in the antler tip. Methods The tissues from both the cartilage formation zone (FZ) and remodel zone (RZ) in rapid growing antlers of sika deer were sampled, profiles of both mRNA and miRNA from these samples were sequenced and analysed, miRNA–mRNA regulatory networks for these two zones were constructed, and their encoded/targeted differentially expressed genes (DEGs) were identified through bioinformatics analysis. Key results In total, 3703 DEGs in the FZ over the RZ were identified, with 1615 being upregulated and 2088 downregulated. The upregulated DEGs in the FZ were found to be mainly enriched in cell proliferation and chondrogenesis/osteogenesis, whereas those in the RZ were enriched in the formation of chondroclasts and osteoclasts. In total, 308 unique mature miRNAs were detected including 110 significantly differentially expressed miRNAs. These miRNAs are predicted to target extracellular matrix proteins, growth factors and receptors, and transcriptional factors, all related to cartilage formation and remodelling. To verify the reliability of our datasets, we successfully tested the regulatory function of one of the top 10 hub miRNAs, miR-155, in vitro. Conclusions The miRNA–mRNA regulatory networks for cartilage formation zone (FZ) in relation to cartilage remodel zone (RZ) were successfully constructed, and validated, which has laid the foundation for the identification of potent growth factors and novel regulation system in bone formation through endochondral ossification. Implications We believe that our datasets are reliable for further mining potent growth factors and novel regulation systems for rapid cartilage formation, remodelling and bone fracture repair by using this unique model, the deer antler.

The growth of antler is driven by endochondral ossification in the growth center of the apical region. Antler grows faster than cancer tissues, but it can be stably regulated and regenerated periodically. To elucidate the molecular mechanisms of how antler grows rapidly without carcinogenesis, in this study, we used RNA-seq technology to evaluate the changes of miRNA and mRNA profiles in antler at four different developmental stages, including 15, 60, 90, and 110 days. We identified a total of 55004 unigenes and 246 miRNAs of which, 10182, 13258, 10740 differentially expressed (DE) unigenes and 35, 53, 27 DE miRNAs were identified in 60-day vs. 15-day, 90-day vs. 60-day, and 110-day vs. 90-day. GO and KEGG pathway analysis indicated that DE unigenes and DE miRNA were mainly associated with chondrogenesis, osteogenesis and inhibition of oncogenesis, that were closely related to antler growth. The interaction networks of mRNA-mRNA and miRNA-mRNA related to chondrogenesis, osteogenesis and inhibition of oncogenesis of antler were constructed. The results indicated that mRNAs (COL2A1, SOX9, WWP2, FGFR1, SPARC, LOX, etc.) and miRNAs (miR-145, miR-199a-3p, miR-140, miR-199a-5p, etc.) might have key roles in chondrogenesis and osteogenesis of antler. As well as mRNA (TP53, Tpm3 and ATP1A1, etc.) and miRNA (miR-106a, miR-145, miR-1260b and miR-2898, etc.) might play important roles in inhibiting the carcinogenesis of antler. In summary, we constructed the mRNA-mRNA and miRNA-mRNA regulatory networks related to chondrogenesis, osteogenesis and inhibition of oncogenesis of antler, and identified key candidate mRNAs and miRNAs among them. Further developments and validations may provide a reference for in-depth analysis of the molecular mechanism of antler growth without carcinogenesis.

Deer species has both scientific research and economic value, and half of these species, however, are listed as endangered animals. For the conservation purpose, we designed a novel universal deer-specific PCR primer pair based on an evolutionarily conservative coding sequence (i.e., CEP295NL gene) across some deer species. This primer pair was successfully amplified and sequenced, showing around ~ 540 bp in cervids. Validation results showed that it can be utilized to develop a reliable and simple diagnostic tool for distinguishing other closely related species, as well as possibly interspecific identification amongst cervids.

Deer antlers are the only mammalian organs capable of complete renewal. Antler renewal is a stem cell-based [antler stem cells (ASCs)] process. Maintenance and activation of the ASCs require them to be located in a specialized microenvironment (niche), and to interact with the cells resident in the niche. Based on previous experiments we found that niche of the ASCs is provided by the closely associated enveloping skin, which currently was known includes dermal papilla cells (DPCs) and epidermal cells. Antler generation/regeneration are triggered by the interactions between ASCs and the niche. In the present study, we established an in vitro co-culture system in which ASCs and DPCs, were cultured together to mimic the in vivo state. A MLEFF strategy was adopted to identify the interactive molecules from the co-culture system. In total, 128 molecules were identified and over 60% belonged to exosomes. Important biological processes that were activated by these molecules included osteoblast differentiation, angiogenesis, and the PI3K-AKT signaling pathway. In so doing, we have significantly simplified the process for identifying interactive molecules, which may be the key signals for triggering antler formation/renewal. Further study of these molecules will help us to gain insights into the mechanism of mammalian organ regeneration.

Canine distemper virus (CDV) causes a highly contagious disease in a wide range of carnivores. The hemagglutinin (H) protein of viruses shows the highest variability and plays an important role in modulation of viral antigenicity, virulence, and receptor recognition. Since 2012, canine distemper (CD) outbreaks in fur-bearing animals (minks, foxes, raccoon dogs) caused by CDV variants with I542N and Y549H substitutions in the H protein have been frequently reported in China. To characterize the molecular evolutionary dynamics and epidemiological dynamics of CDV, 235 H gene sequences of CDV wild-type strains collected from 22 countries between 1975 and 2015, including 44 strains predominant in fur-bearing animals in China, were analyzed. The phylogenetic relationships and evolutionary rates of the CDV strains were determined by Bayesian phylogenetics. The CDV strains clustered into distinct geographic genotypes, irrespective of the species of isolation. All the variant strains formed a distinct monophyletic cluster and belonged to the F sub-genotype within the Asia-1 genotype-currently the predominant sub-genotype in fur-bearing animals in China. Evolutionary analysis suggested that the variant strains originated in 2006. Furthermore, the selection pressure analysis revealed that the Y549H substitution was under positive selection pressure for adaptation toward the fur-bearing animals. The residue at position 549 also showed structural interaction with the V domain of the mink signaling lymphocyte-activation molecule (SLAM) receptor based on the homology modeling of the H-SLAM complex. Our results suggested that the Y549H substitution contributed to the molecular adaptation of CDV variants in the fur-bearing animals during the viral evolutionary phase in China.

S100A4 is a multiple-function protein highly expressed in tumor or stem cells. We found S100A4 was a novel protein partner for heat shock protein 47 (HSP47) in deer antlerogenic periosteum cells (AP cells), indicating that S100A4 could bind with HSP47. S100A4 had both calcium-dependent and calcium-independent patterns (labeled as SCd and SCi, respectively) to execute different biological activities. Homology models of HSP47, SCd and SCi were constructed. HSP47:collagen model, HSP47:collagen I-V, HSP47:SCd and HSP47:SCi complexes were built using ZDOCK software. Together with free SCd and SCi, 200 ns molecular dynamic (MD) simulations were performed to analyze binding free energies and SCi/SCd conformational changes. The energetic results showed that SCi had the strongest affinity to HSP47, and followed by collagens. SCd had little interaction with HSP47. Decomposition energy results showed that collagen model interacted with HSP47 mainly though neutral amino acids. When SCi bound with HSP47, the majority of mediated amino acids were charged. These results indicated that SCi could compete with collagen on the binding site of HSP47. Root mean square fluctuation (RMSF) values and cross-correlation matrices of principal component analysis (PCA) were calculated to evaluate the SCi/SCd structural variation during MD simulation. Both HSP47 and Ca2+ could stabilize the conformation of SCi/SCd. The loops interacting with Ca2+s and linking the two EF-hand motifs were impacted particularly. The relative moving directions of α-helices in EF-hands were distinct by the binding effect of HSP47 and Ca2+. We found that SCi may regulate the differentiation of AP cells by disturbing the interaction between HSP47 and collagen. Communicated by Ramaswamy H. Sarma.

Deer antlers are extraordinary mammalian organs that can fully regenerate annually. Antler renewal is a stem cell-based epimorphic process and antler stem (AS) cells can initiate de novo generation of antlers in postnatal mammals. However, although being called stem cells, the AS cells have not been characterized at molecular level based on the stem cell criteria. Comprehensive characterization of the AS cells would undoubtedly help to decipher the mechanism underlying the full regeneration of deer antlers, the only case of stem cell-based epimorphic regeneration in mammals. In the present study, three types of AS cells (antlerogenic periosteal cells APCs, for initial pedicle and first antler formation; pedicle periosteal cells PPC, for annual antler regeneration; and reserve mesenchyme cells RMCs, for rapid antler growth), were isolated for comprehensive molecular characterization. A horn-growth-related gene, RXFP2, was found to be expressed only in AS cells lineages but not in the facial periosteal cells (FPCs, locates geographically in the vicinity of the APCs or PPCs), suggesting the RXFP2 might be a specific marker for the AS cell lineage in deer. Our results demonstrated that AS cells expressed classic MSC markers including surface markers CD73, CD90, CD105 and Stro-1. They also expressed some of the markers including Tert, Nestin, S100A4, nucleostemin and C-Myc, suggesting that they have some attributes of the ESCs. Microinjection of male APC into deer blastocysts resulted in one female foetus (110 days gestation) recovered with obvious pedicle primordia with both male and female genotype detected in the ovary. In conclusion, the AS cells should be defined as MSCs but with partial attributes of ESCs.

Tarim red deer (Cervus elaphus yarkandensis) is the only subspecies of red deer (of 22 subspecies) from Central Asia. This species is a desert dweller of the Tarim Basin of southern Xinjiang, China, and exhibits some unique adaptations to the dry and extreme hot climate. We report here the assembly of a Tarim red deer genome employing a 10X Genomics library, termed CEY_v1. Our genome consisted of 2.6 Gb with contig N50 and scaffold N50 of 275.5 Kb and 31.7 Mb, respectively. Around 96% of the assembled sequences were anchored onto 34 chromosomes based on the published high-quality red deer genetic linkage map. More than 94% BUSCOs complete genes (including 90.5% single and 3.6% duplicated ones) were detected in the CEY_v1 and 20,653 genes were annotated. The CEY_v1 is expected to contribute to comparative analysis of genome biology, to evolutionary studies within Cervidae, and to facilitating investigation of mechanisms underlying adaptation of this species to the extreme dry and hot climate., Measurement(s)DNA • genome • sequence_assembly • sequence feature annotationTechnology Type(s)DNA sequencing • sequence assembly process • sequence annotationSample Characteristic - OrganismCervus elaphus Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12410564

Antler regeneration, a stem cell-based epimorphic process, has potential as a valuable model for regenerative medicine. A pool of antler stem cells (ASCs) for antler development is located in the antlerogenic periosteum (AP). However, whether this ASC pool is homogenous or heterogeneous has not been fully evaluated. In this study, we produced a comprehensive transcriptome dataset at the single-cell level for the ASCs based on the 10x Genomics platform (scRNA-seq). A total of 4,565 ASCs were sequenced and classified into a large cell cluster, indicating that the ASCs resident in the AP are likely to be a homogeneous population. The scRNA-seq data revealed that tumor-related genes were highly expressed in these homogeneous ASCs: i.e. TIMP1, TMSB10, LGALS1, FTH1, VIM, LOC110126017 and S100A4. Results of screening for stem cell markers suggest that the ASCs may be considered as a special type of stem cell between embryonic (CD9) and adult (CD29, CD90, NPM1 and VIM) stem cells. Our results provide the first comprehensive transcriptome analysis at the single-cell level for the ASCs, and identified only one major cell type resident in the AP and some key stem cell genes, which may hold the key to why antlers, the unique mammalian organ, can fully regenerate once lost.

Sika deer (Cervus Nippon) and wapiti (Cervus elaphus) are closely related species and their hybridization can result in significant allele-shift of their gene pool. Additive genetic effects and putative heterotic effects of their hybridization on growth performance could confer considerable economic advantage in deer farming. Here, we used double-digest restriction site-associated DNA sequencing technology (ddRAD-seq) and detected ∼320,000 genome-wide SNPs from 30 captive individuals: 7 sika deer, 6 wapiti and 17 F1 hybrids (reciprocal cross). By screening observed heterozygosity of each SNP across four taxonomic groups, we report for the first time a resource of 2,015 putative diagnostic SNP markers (species-specific SNPs for sika deer and wapiti), which can be used to design tools for assessing or monitoring the degree of hybridization between sika deer and wapiti. These ddRAD-seq data and SNP datasets are also valuable resources for genome-wide studies, including trait discovery for breeders of domestic deer.

Annual antler renewal is a stem cell–based epimorphic process driven by antler stem cells (ASCs) resident in antlerogenic periosteum (AP). Antlerogenic periosteal cells express a high level of S100A4, a metastasis-associated protein, which intrigued us to explore what role S100A4 could play in antler regeneration. The present study set out to investigate expression and effects of S100A4 in the ASCs and their progeny. The results showed that not only did cells from the AP express a high level of S100A4, but also the pedicle periosteum and the antler growth center. In the antler growth center, we found S100A4-positive cells were specifically located in blood vessel walls and in vascularized areas. In vitro, recombinant deer S100A4 protein stimulated the proliferation of the AP cells, promoted proliferation, migration and tube formation of human vascular endothelial cells, and enhanced migration of Hela cells, but not AP cells. These findings demonstrated that S100A4 in the ASCs may play a significant role in stimulating angiogenesis, proliferation, but not motility, of ASCs. Deer antlers offer a unique model to explore how rapid cell proliferation with a high level of S100A4 expression is elegantly regulated without becoming cancerous.

Background: Cervids have evolved very successful means for survival and thriving to adapt to various climates and environments. One of these successful means might be the effective and efficient way of communication. To support this notion, cervids are well equipped with a variety of skin glands that distribute in different body regions. However, studies relevant to adaptive evolution in cervids, particularly on olfactory reception at the molecular level, have thus far not been reported. Objective: To provide valuable insights into molecular evidence for the adaptive evolution of olfactory-related gene in cervids. Methods: Based on recently sequenced genomes of cervids and closely-related-species, we performed comparative genomic analysis at genome level using bioinformatics tools. Results: Tree topology strongly supported that Bovidae was the sister group of Moschidae and both formed a branch that was then clustered with Cervidae. Expansion of heavy chain genes of the dynein family and 51 rapidly evolving genes could be associated with adaptation of cilia that serve as sensory organelles and act as cellular antennae. Based on the branch-site model test along the deer branch spanning 7–21 mammalian species, 14 deer olfactory receptor genes were found to be undergoing positive selection pressure and 89 positive selection sites (probability > 60%) had amino acid substitutions unique to deer. Conclusion: This study, for the first time, provides significant molecular evidence for adaption of olfactory-related genes of cervids according to their olfactory behavior.

Embryo implantation in the mink follows the pattern of many carnivores, in that preimplantation embryo diapause occurs in every gestation. Details of the gene expression and regulatory networks that terminate embryo diapause remain poorly understood. Illumina RNA-Seq was used to analyze global gene expression changes in the mink uterus during embryo diapause and activation leading to implantation. More than 50 million high quality reads were generated, and assembled into 170,984 unigenes. A total of 1684 differential expressed genes (DEGs) in uteri with blastocysts in diapause were compared to the activated embryo group (p < 0.05). Among these transcripts, 1527 were annotated as known genes, including 963 up-regulated and 564 down-regulated genes. The gene ontology terms for the observed DEGs, included cellular communication, phosphatase activity, extracellular matrix and G-protein couple receptor activity. The KEGG pathways, including PI3K-Akt signaling pathway, focal adhesion and extracellular matrix (ECM)-receptor interactions were the most enriched. A protein-protein interaction (PPI) network was constructed, and hub nodes such as VEGFA, EGF, AKT, IGF1, PIK3C and CCND1 with high degrees of connectivity represent gene clusters expected to play an important role in embryo activation. These results provide novel information for understanding the molecular mechanisms of maternal regulation of embryo activation in mink.

Background With the unprecedented rapid growth rate (up to 2.75 cm/day), velvet antler is an invaluable model for the identification of potent growth factors and signaling networks for extremely fast growing tissues, mainly cartilage. Antler growth center (AGC) locates in its tip and consists of five tissue layers: reserve mesenchyme (RM), precartilage (PC), transition zone (TZ), cartilage (CA) and mineralized cartilage (MC). The aim of this study was to investigate the transcription dynamics in the AGC using RNA-seq technology. Results Five tissue layers in the AGC were collected from three 3-year-old male sika deer using our previously reported sampling method (morphologically distinguishable). After sequencing (15 samples; triplicates/tissue layer), we assembled a reference transcriptome de novo and used RNA-seq to measure gene expression profiles across these five layers. Nine differentially expressed genes (DEGs) were selected from our data and subsequently verified using qRT-PCR. The results showed a high consistency with the RNA-seq results (R2 = 0.80). Nine modules were constructed based on co-expression network analysis, and these modules contained 370 hub genes. These genes were found to be mainly involved in mesenchymal progenitor cell proliferation, chondrogenesis, osteogenesis and angiogenesis. Combination of our own results with the previously published reports, we found that Wnt signaling likely plays a key role not only in stimulating the antler stem cells or their immediate progeny, but also in promoting chondrogenesis and osteogenesis during antler development. Conclusion We have successfully assembled a reference transcriptome, generated gene expression profiling across the five tissue layers in the AGC, and identified nine co-expressed modules that contain 370 hub genes and genes predorminantly expressed in and highly relevant to each tissue layer. We believe our findings have laid the foundation for the identification of novel genes for rapid proliferation and chondrogenic differentiation of antler cells.

Nanog plays a crucial role in the maintenance of stem cell pluripotency. Annual full regeneration of deer antlers has been shown to be a stem cell-based process, and antler stem cells (ASCs) reportedly express Nanog. In the present study, we found that Nanog RNA expressed by ASCs was a pseudogene (Nanog-ps). The coding sequence of Nanog-ps was 93.1% homologous to that of bovine Nanog, but with two missing nucleotides after position 391. Deletion of the two nucleotides in Nanog-ps resulted in a frame-shift mutation, suggesting that Nanog-ps would not encode a normal Nanog protein. Overexpression of Nanog-ps failed to affect downstream genes of Nanog or to enhance cell proliferation in the ASCs. However, this pseudogene was transcribed in the ASCs and encoded a nuclear protein; the expression levels of Nanog-ps were also related to the degree of stemness in antler cells. Here, we reported this pseudogene, because it could serve as a useful marker for identifying ASCs and evaluating the degree of their stemness.

Sika deer is of great commercial value because their antlers are used in tonics and alternative medicine and their meat is healthy and delicious. The goal of this study was to generate transcript sequences from sika deer for functional genomic analyses and to identify the transcripts that demonstrate tissue-specific, age-dependent differential expression patterns. These sequences could enhance our understanding of the molecular mechanisms underlying sika deer growth and development. In the present study, we performed de novo transcriptome assembly and profiling analysis across ten tissue types and four developmental stages (juvenile, adolescent, adult, and aged) of sika deer, using Illumina paired-end tag (PET) sequencing technology. A total of 1,752,253 contigs with an average length of 799 bp were generated, from which 1,348,618 unigenes with an average length of 590 bp were defined. Approximately 33.2 % of these (447,931 unigenes) were then annotated in public protein databases. Many sika deer tissue-specific, age-dependent unigenes were identified. The testes have the largest number of tissue-enriched unigenes, and some of them were prone to develop new functions for other tissues. Additionally, our transcriptome revealed that the juvenile-adolescent transition was the most complex and important stage of the sika deer life cycle. The present work represents the first multiple tissue transcriptome analysis of sika deer across four developmental stages. The generated data not only provide a functional genomics resource for future biological research on sika deer but also guide the selection and manipulation of genes controlling growth and development.

Table S3. Summary of 13,203 gene sequences matched with non-redundancy proteins of three closely related species using BLASTX. (DOCX 14 kb)

Table S2. Statistics of transcript assembly in the pipeline. (DOCX 16 kb)

Figure S2. De novo assembly of the reference transcriptome. A) Length distribution of coding and noncoding transcripts across the length bins. Triangles and dots represent the percentages of coding and noncoding transcripts in the length bins, respectively. B) Expression levels displayed differences of 3–4 orders of magnitude. For each tissue layer, transcripts with FPKM ≥0.5 are all replicates. C) In all cases, the coding transcripts (40–50%) are more abundant in the upper ranks of the distribution (Q4) than the non-coding sequences (20–30%). (PDF 468 kb)

MicroRNAs (miRNAs) can effectively regulate gene expression at the post-transcriptional level and play a critical role in tissue growth, development and regeneration. Our previous studies showed that antler regeneration is a stem cell-based process and antler stem cells reside in the periosteum of a pedicle, the permanent bony protuberance, from which antler regeneration takes place. Antlers are the only mammalian organ that can fully regenerate and hence provide a unique opportunity to identify miRNAs that are involved in organ regeneration. In the present study, we used next generation sequencing technology sequenced miRNAs of the stem cells derived from either the potentiated or the dormant pedicle periosteum. A population of both conserved and 20 deer-specific miRNAs was identified. These conserved miRNAs were derived from 453 homologous hairpin precursors across 88 animal species, and were further grouped into 167 miRNA families. Among them, the miR-296 is embryonic stem cell-specific. The potentiation process resulted in the significant regulation (>±2 Fold, q value

BackgroundReindeer (Rangifer tarandus) is the only fully domesticated species in the Cervidae family, and is the only cervid with a circumpolar distribution. Unlike all other cervids, female reindeer regularly grow cranial appendages (antlers, the defining characteristics of cervids), as well as males. Moreover, reindeer milk contains more protein and less lactose than bovids’ milk. A high quality reference genome of this specie will assist efforts to elucidate these and other important features in the reindeer.FindingsWe obtained 723.2 Gb (Gigabase) of raw reads by an Illumina Hiseq 4000 platform, and a 2.64 Gb final assembly, representing 95.7% of the estimated genome (2.76 Gb according to k-mer analysis), including 92.6% of expected genes according to BUSCO analysis. The contig N50 and scaffold N50 sizes were 89.7 kilo base (kb) and 0.94 mega base (Mb), respectively. We annotated 21,555 protein-coding genes and 1.07 Gb of repetitive sequences by de novo and homology-based prediction. Homology-based searches detected 159 rRNA, 547 miRNA, 1,339 snRNA and 863 tRNA sequences in the genome of R. tarandus. The divergence time between R. tarandus, and ancestors of Bos taurus and Capra hircus, is estimated to be 29.55 million years ago (Mya).ConclusionsOur results provide the first high-quality reference genome for the reindeer, and a valuable resource for studying evolution, domestication and other unusual characteristics of the reindeer.

Context Commercially, velvet antlers along the longitudinal axis are divided into four portions, namely, wax-like (WL), blood-colour (BC), honeycomb-like (HL) and bone (B) slices from the top to the base. However, there is no evidence at a molecular level showing the accuracy of this classification. Aims The aim of the present study was to take transcriptional approach to assess the accuracy of the traditional classification for these four portions of velvet antler, and to link the expressed mRNAs of each portion with possible functions by using bioinformatics analysis. Methods Three sticks of three-branched velvet antlers of sika deer were harvested from three anaesthetised 4-year-old sika deer. On the basis of the traditional methods used commercially, the velvet antler sticks were divided into the four portions of WL, BC, HL and B. Transcriptome sequencing was performed using Illumina HiSeq × Ten at BGI (Shenzheng, China). Key results In total, 5647 genes were obtained from the four portions. Spearman correlation analysis grouped these four portions into three clusters (WL, BC, HL+B). C-means analysis further confirmed a similar trend, indicating the accuracy of the new classification based on transcriptome analysis. Further functional analysis showed that highly expressed genes in WL, BC and HL+B were mainly related to cell cycle, cartilage development, and bone development respectively. Conclusions Four-portion classification based on traditional methods should be replaced by three-portion classification based on the mRNA expression levels. Implications We believe that this new classification can contribute to velvet antler industry, providing more accuracy in the use of velvet antlers as pharmaceuticals.

Context Immune system has been claimed as the ‘main switch’ of tissue or organ regeneration. Among immune cells, macrophages stand out as important modulators in mutiple regeneration models, such as planarian, axolotl, mammalian hair and liver. As a unique model for mammals, deer antler is considered to ideal for studying complete mammalian organ regeneration. Studies have found that antler regeneration is a stem cell-based process and antler stem cells locate in the pedicle periosteum (PP). Although the regulatory roles of the immune system in other regeneration models have been extensively studied, they remain unstudied in antler regeneration. Aims To explore the possible role of macrophages in the PP cells (PPCs). Methods We treated PPCs with a macrophage-conditioned medium (MCM) and detected effects of MCM on proliferation, migration and apoptosis of the PPCs, and identified differentially expressed genes by using the RNA-seq technique. Key results We found that MCM enhanced proliferation rate and migration rate significantly and stimulated apoptosis of the PPCs. Using the RNA-seq technique, we identified 112 differentially expressed genes in the PPCs (38 downregulated and 74 upregulated) after the MCM treatment. Furthermore, gene-ontology annotation analyses showed that the upregulated genes were mainly involved in cell adhesion, chemotaxis, wound healing, growth factor-stimulated responses, and bone formation, and the downregulated genes were involved in regulation of biosynthesis. Conclusions MCM had a great influence on the antler stem cells, and macrophages might regulate antler regeneration through altering the microenvironment and gene-expression profiles of the PPCs. Implications We believe that the results of the present study would facilitate the discovery of the roles of immune system in antler stem cells and, thus, mammalian organ regeneration in general.

Deer antlers offer a unique model to study organ regeneration in mammals. Antler regeneration relies on the pedicle periosteum (PP) cells and is triggered by a decrease in circulating testosterone (T). The molecular mechanism for antler regeneration is however, unclear. Label-free liquid chromatography-mass spectrometry (LC-MS/MS) was used to identify differentially-expressed proteins (DEPs) in the regeneration-potentiated PP (under low T environment) over the non-regeneration-potentiated PP (under high T environment). Out of total 273 DEPs, 189 were significantly up-regulated and 84 were down-regulated from these comparisons: after castration vs before castration, natural T vs before castration, and exogenous T vs before castration. We focused on the analysis only of those DEPs that were present in fully permissive environment to antler regeneration (low T). Nine transduction pathways were identified through the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, including the estrogen signaling pathway. A total of 639 gene ontology terms were found to be significantly enriched in regeneration-potentiated PP (low T) from the DEPs. Reliability of the label free LC-MS/MS was determined by qRT-PCR to estimate the expression level of selected genes. The results suggest that up-regulated heat shock proteins (HSP90AB1, HSP90B1), peptidyl-prolyl cis-trans isomerase 4 (FKBP4), mitogen-activated protein kinase 3 (MAPK3) and calreticulin (CALR) and down-regulated SHC-transforming protein 1 (SHC1), heat shock protein family A member 1A (HSPA1A) and proto-oncogene tyrosine-protein kinase (SRC) may be associated directly or indirectly with antler regeneration. Further studies are required to investigate the roles of these proteins in regeneration using appropriate in vivo models.

Embryo implantation in the mink follows the pattern of many carnivores, in that preimplantation embryo diapause occurs in every gestation. Details of the gene expression and regulatory networks that terminate embryo diapause remain poorly understood. Illumina RNA-Seq was used to analyze global gene expression changes in the mink uterus during embryo diapause and activation leading to implantation. More than 50 million high quality reads were generated, and assembled into 170,984 unigenes. A total of 1684 differential expressed genes (DEGs) in uteri with blastocysts in diapause were compared to the activated embryo group (p &lt, 0.05). Among these transcripts, 1527 were annotated as known genes, including 963 up-regulated and 564 down-regulated genes. The gene ontology terms for the observed DEGs, included cellular communication, phosphatase activity, extracellular matrix and G-protein couple receptor activity. The KEGG pathways, including PI3K-Akt signaling pathway, focal adhesion and extracellular matrix (ECM)-receptor interactions were the most enriched. A protein-protein interaction (PPI) network was constructed, and hub nodes such as VEGFA, EGF, AKT, IGF1, PIK3C and CCND1 with high degrees of connectivity represent gene clusters expected to play an important role in embryo activation. These results provide novel information for understanding the molecular mechanisms of maternal regulation of embryo activation in mink.

Tandem repeat units are only detected in the left domain of the mitochondrial DNA control region in sika deer. Previous studies showed that Japanese sika deer have more tandem repeat units than its cousins from the Asian continent and Taiwan, which often have only three repeat units. To determine the origin and evolution of these additional repeat units in Japanese sika deer, we obtained the sequence of repeat units from an expanded dataset of the control region from all sika deer lineages. The functional constraint is inferred to act on the first repeat unit because this repeat has the least sequence divergence in comparison to the other units. Based on slipped-strand mispairing mechanisms, the illegitimate elongation model could account for the addition or deletion of these additional repeat units in the Japanese sika deer population. We also report that these additional repeat units could be occurring in the internal positions of tandem repeat regions, possibly via coupling with a homogenization mechanism within and among these lineages. Moreover, the increased number of repeat units in the Japanese sika deer population could reflect a balance between mutation and selection, as well as genetic drift.

As the only known organ that can completely regenerate in mammals, deer antler is of real significance in the field of regenerative medicine. Recent studies have shown that the regenerative capacity of the antlers comes from the pedicle periosteum and the cells resident in the periosteum possess the attributes of stem cells. Currently, the molecular mechanism of antler regeneration remains unclear. In the present study, we compared the potentiated and dormant antler stem cells using isobaric tags for the relative and absolute quantification (iTRAQ) labeling of the peptides, coupled with two-dimensional liquid chromatography-tandem mass spectrometry (LC-MS/MS) to compare the proteome profiles. Proteins were identified by searching against the NCBI nr database and our own Cervine transcriptome database, and bioinformatics analysis was conducted to identify the differentially expressed proteins. Based on this searching strategy, we identified 169 differentially expressed proteins in total, consisting of 70 up- and 99 down-regulated in the potentiated vs. dormant antler stem cells. Reliability of the iTRAQ was confirmed via quantitative real-time polymerase chain reaction (qRT-PCR) to measure the expression of selected genes. We identified transduction pathways through the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, such as HIF-1 and PI3K-AKT signaling pathways that play important roles in regulating the regeneration of antlers. In summary, the initiation stage of antler regeneration, a process from dormant to potentiated states in antler stem cells, is regulated by multiple proteins and complicated signal networks.

Complete mitochondrial genome of Tianshan wapiti, Cervus elaphus songaricus, is 16,419 bp in length and contains 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and 1 control region. The phylogenetic trees were reconstructed with the concatenated nucleotide sequences of the 13 protein-coding genes using maximum parsimony (MP) and Bayesian inference (BI) methods. MP and BI phylogenetic trees here showed an identical tree topology. The monopoly of red deer, wapiti and sika deer was well supported, and wapiti was found to share a closer relationship with sika deer. Tianshan wapiti shared a closer relationship with xanthopygus than yarkandensis. Rusa unicolor and Rucervus eldi were given a basal phylogenetic position. Our phylogenetic analysis provided a robust phylogenetic resolution spanning the entire evolutionary relationship of the subfamily Cervinae.