Full-length transcriptome and microRNA sequencing reveal the specific gene-regulation network of velvet antler in sika deer with extremely different velvet antler weight

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.

PRRX1/miR-143-3p signaling regulates homeostasis of antler reserve mesenchymal cells

The molecular regulation mechanisms for maintaining the homeostasis of mesenchymal stem cells still remains poorly defined. Antler reserve mesenchymal cells (RM cells) persist through the whole rapid antler growth stage as a reserved stem cell population capable of division and differentiation, that makes the RM cells a unique model in stem cell regulation and cancer mechanism studies. Herein, we sequenced and analyzed the extracellular vesicles (EVs) of RM cells in the growth center of antler, and identified a high expression level of miR-143-3p and its target genes IGF1R, TGFβ1, BMP2, etc. The upstream positive regulatory factor PRRX1 of miR-143-3p was identified through ATAC and CUT-taq analysis, combined with dual luciferase assay. We showed that PRRX1 overexpression resulted in a decreased proliferation of RM cells and induced a higher expression of miR-143-3p. miR-143-3p enriched EVs derived from PRRX1 overexpression RM cells had an inhibitory effect on RM cells, osteosarcoma 143B cells (considered as excessive proliferation model for RM cells) and in vivo in osteosarcoma bearing mice, and the mRNA and protein levels of IGF1R were significantly reduced. We confirmed that miR-143-3p enriched EVs inhibited 3D culture induced chondrogenic differentiation of RM cells and xenogeneic antler chondrogenesis through targeting TGFβ1 and BMP2. Together, PRRX1 was identified as an activator of miR-143-3p, and higher amounts of miR-143-3p in EVs of RM cells could inhibit excessive proliferation, and help maintain the undifferentiated state of RM cells. We conclude that PRRX1/miR-143-3p signaling was a regulator of homeostasis of antler RM cells and was a potential regulator of osteosarcoma. Our findings are essential for advancing medical and biological sciences, providing new theoretical foundations and strategies for cancer treatment and tissue regeneration.

Bone metabolism associated with annual antler regeneration: a deer insight into osteoporosis reversal

Osteoporosis, a metabolic disorder, remains challenging to treat due to limited understanding of its underlying mechanism. The annual cycle of "cyclic physiological osteoporosis (CPO)" and its full reversal in male deer represents a unique natural model for studying this condition. Deer antlers, weighing up to 25 kg/pair, derive over 60% of their mineral contents from deer skeleton during mineralization. Based on the literature, we propose to divide CPO and its reversal into two phases: Phase I (approximately 115 days): from hard antler casting to the end of antler linear growth, marked by simultaneous robust antler ossification and CPO development; and Phase II (up to 165 days): from end of Phase I to the onset of antler skin shedding, characterized by complete antler mineralization and CPO reversal. This review analyzes the paradoxical occurrence of robust antler ossification and skeleton CPO within the same endocrine microenvironment during phase I; total antler mineralization and full reversal of deer skeleton CPO in phase II. Furthermore, we will discuss potential insights for osteoporosis treatment using deer materials from the period of Phase II. Our goal is to identify novel substances and therapies that could be applied in clinical setting to effectively treat osteoporosis.

MicroRNA PC-3p-2869 Regulates Antler Growth and Inhibits Proliferation and Migration of Human Osteosarcoma and Chondrosarcoma Cells by Targeting CDK8, EEF1A1, and NTN1

MicroRNAs (miRNAs) play a crucial role in maintaining the balance between the rapid growth and suppression of tumorigenesis during antler regeneration. This study investigated the role of a novel miRNA, PC-3p-2869 (miR-PC-2869), in antler growth and its therapeutic potential in human osteosarcoma and chondrosarcoma. Stem-loop RT-qPCR showed that miR-PC-2869 was expressed extensively in diverse layers of antler tissues. Overexpression of miR-PC-2869 suppressed the proliferation and migration of antler cartilage cells. Similarly, heterologous expression of miR-PC-2869 reduced the proliferation, colony formation, and migration of osteosarcoma cell line MG63 and U2OS and chondrosarcoma cell line SW1353. Moreover, 18 functional target genes of miR-PC-2869 in humans were identified based on the screening of the reporter library. Among them, 15 target genes, including CDK8, EEF1A1, and NTN1, possess conserved miR-PC-2869-binding sites between humans and red deer (Cervus elaphus). In line with this, miR-PC-2869 overexpression decreased the expression levels of CDK8, EEF1A1, and NTN1 in MG63, SW1353, and antler cartilage cells. As expected, the knockdown of CDK8, EEF1A1, or NTN1 inhibited the proliferation and migration of MG63, SW1353, and antler cartilage cells, demonstrating similar suppressive effects as miR-PC-2869 overexpression. Furthermore, we observed that CDK8, EEF1A1, and NTN1 mediated the regulation of c-myc and cyclin D1 by miR-PC-2869 in MG63, SW1353, and antler cartilage cells. Overall, our work uncovered the cellular functions and underlying molecular mechanism of antler-derived miR-PC-2869, highlighting its potential as a therapeutic candidate for bone cancer.

MiRNA Profiling and Its Potential Roles in Rapid Growth of Velvet Antler in Gansu Red Deer (Cervus elaphus kansuensis)

A significant variety of cell growth factors are involved in the regulation of antler growth, and the fast proliferation and differentiation of various tissue cells occur during the yearly regeneration of deer antlers. The unique development process of velvet antlers has potential application value in many fields of biomedical research. Among them, the nature of cartilage tissue and the rapid growth and development process make deer antler a model for studying cartilage tissue development or rapid repair of damage. However, the molecular mechanisms underlying the rapid growth of antlers are still not well studied. MicroRNAs are ubiquitous in animals and have a wide range of biological functions. In this study, we used high-throughput sequencing technology to analyze the miRNA expression patterns of antler growth centers at three distinct growth phases, 30, 60, and 90 days following the abscission of the antler base, in order to determine the regulatory function of miRNA on the rapid growth of antlers. Then, we identified the miRNAs that were differentially expressed at various growth stages and annotated the functions of their target genes. The results showed that 4319, 4640, and 4520 miRNAs were found in antler growth centers during the three growth periods. To further identify the essential miRNAs that could regulate fast antler development, five differentially expressed miRNAs (DEMs) were screened, and the functions of their target genes were annotated. The results of KEGG pathway annotation revealed that the target genes of the five DEMs were significantly annotated to the "Wnt signaling pathway", "PI3K-Akt signaling pathway", "MAPK signaling pathway", and "TGF-β signaling pathway", which were associated with the rapid growth of velvet antlers. Therefore, the five chosen miRNAs, particularly ppy-miR-1, mmu-miR-200b-3p, and novel miR-94, may play crucial roles in rapid antler growth in summer.

Current Situation and Utilization of Velvet Deer Germplasm Resources in China

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.

lncRNA Sequencing of Antler Mesenchymal Tissue Revealed that the Regulatory Network of Antler Cell Proliferation and Differentiation

Antlers have been widely studied due to their unique physiological characteristics, such as rapid growth, periodic shedding and regeneration. However, little is known about how antler growth is regulated by long non-coding RNA (lncRNA). The aim of the present study was to identify the lncRNAs expression profile and explore the function of lncRNAs during the antler growth. Herein, RNA-sequencing technology (RNA-seq) was performed on the three growth periods (early developmental period: EP, middle developmental period: MP, later developmental period: LP) of male sika deer (Cervus nippon) antler, 16 differentially expressed lncRNAs (DE lncRNAs) and 11 DE lncRNAs were identified in EP vs MP and MP vs LP related to cell proliferation and cell differentiation, respectively. Finally, lncRNAs-mRNAs co-expression networks were constructed based on the identified DE lncRNAs and their potential trans-target genes. The result reveals that lncRNAs may play diverse roles in different periods of antler growth. It provides a novel perspective for revealing the molecular mechanism of antler growth.

Comprehensive transcriptome analysis of sika deer antler using PacBio and Illumina sequencing

Antler is the fastest growing and ossifying tissue in animals and it is a valuable model for cartilage/bone development. To understand the molecular mechanisms of chondrogenesis and osteogenesis of antlers, the PacBio Sequel II and Illumina sequencing technology were combined and used to investigate the mRNA expression profiles in antler tip, middle, and base at six different developmental stages, i.e., at 15th, 25th, 45th, 65th, 100th and 130th growth days. Consequently, we identified 24,856 genes (FPKM > 0.1), including 8778 novel genes. Besides, principal component analysis (PCA) revealed a significant separation between the growth stage (25th, 45th and 65th days) and ossification stage (100th and 130th days). COL2A1 gene was significantly abundant in the growth stage, whereas S100A7, S100A12, S100A8, and WFDC18 genes were abundant at the ossification stage. Subsequently screened to 14,765 significantly differentially expressed genes (DEGs), WGCNA and GO functional enrichment analyses revealed that genes related to cell division and chondrocyte differentiation were up-regulated, whereas those with steroid hormone-mediated signaling pathways were down-regulated at ossification stages. Additionally, 25 tumor suppressor genes and 11 oncogenes were identified and were predicted to interact with p53. Co-regulation of tumor suppressor genes and oncogenes is responsible for the special growth pattern of antlers. Together, we constructed the most complete sika deer antler transcriptome database so far. The database provides data support for subsequent studies on the molecular mechanism of sika deer antler chondrogenesis and osteogenesis.

RNA sequencing-based identification of microRNAs in the antler cartilage of Gansu red deer (Cervus elaphus kansuensis)

BACKGROUND

The velvet antler is a complex mammalian bone organ with unique biological characteristics, such as regeneration. The rapid growth stage (RGS) is a special period in the regeneration process of velvet antler.

METHODS

To elucidate the functions of microRNAs (miRNAs) at the RGS of antler development in Gansu red deer (Cervus elaphus kansuensis), we used RNA sequencing (RNA-seq) to analyze miRNA expression profiles in cartilage tissues of deer antler tips at three different growth stages.

RESULTS

The RNA-seq results revealed 1,073 known and 204 novel miRNAs, including 1,207, 1,242, and 1,204 from 30-, 60-, and 90-d antler cartilage tissues, respectively. To identify key miRNAs controlling rapid antler growth, we predicted target genes of screened 25 differentially expressed miRNAs (DEMs) and specifically expressed miRNAs (SEMs) in 60 d and annotated their functions. The KEGG results revealed that target genes of 25 DEMs and 30 SEMs were highly classified in the "Metabolic pathways", "Pathways in cancer", "Proteoglycans in cancer" and "PI3K-Akt signaling pathway". In addition, a novel miRNA (CM008039.1_315920), highly enriched in "NF-kappa B signaling pathway", may need further study.

CONCLUSIONS

The miRNAs identified in our study are potentially important in rapid antler growth. Our findings provide new insights to help elucidate the miRNA-mediated regulatory mechanisms involved during velvet antler development in C. elaphus kansuensis.

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

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.

Peptide Biomarkers Discovery for Seven Species of Deer Antler Using LC-MS/MS and Label-Free Approach

Deer antler is a globally widely used precious natural medicine and the material of deer horn gelatin. However, identification of deer antler species based on traditional approaches are problematic because of their similarity in appearance and physical-chemical properties. In this study, we performed a comprehensive antler peptidome analysis using a label-free approach: nano LC-Orbitrap MS was applied to discover peptide biomarkers in deer adult beta-globin (HBB_A), and HPLC-Triple Quadrupole MS was used to verify their specificity. Nineteen peptide biomarkers were found, on which foundation a strategy for antlers and a strategy for antler mixtures such as flakes or powder are provided to identify seven species of deer antler including Eurasian elk (Alces alces), reindeer (Rangifer tarandus), white-tailed deer (Odocoileus viginianus), white-lipped deer (Przewalskium albirostris), fallow deer (Dama dama), sika deer (Cervus nippon), and red deer (Cervus elaphus) simultaneously. It is worth noting that our search found that the HBB_A gene of sika deer, red deer, and North American wapiti (Cervus canadensis) in China may have undergone severe genetic drifts.

IGF1R and LOX Modules Are Related to Antler Growth Rate Revealed by Integrated Analyses of Genomics and Transcriptomics

Deer antlers are organs of bone and have an extremely rapid growth rate. Thus far, the molecular mechanism underlying rapid antler growth has not been properly elucidated, and key genes driving this growth rate have not been fully identified. In this study, based on the newly assembled high-quality sika deer genome, we conducted an integrated analysis of genome-wide association analysis (GWAS) and weighted gene co-expression network analysis (WGCNA) using genome resequencing data from our previous GWAS, with weight and transcriptome sequencing data of faster- vs. slower-growing antlers of sika deer. The expressions of key genes were verified using Fragments Per Kilobase of transcript per Million fragments mapped (FPKM) in different tissue zones of the antler growth center, different types of sika deer tissues and antler tissues collected from faster and slower growth rates. The results show that a total of 49 genes related to antler growth rate were identified, and most of those genes were enriched in the IGF1R and LOX modules. The gene regulation network of antler growth rate through the IGF1R pathway was constructed. In conclusion, the integration of GWAS and WGCNA analyses had great advantages in identifying regulatory genes of complex antler growth traits over using singular methods individually, and we believe that our findings in the present study can provide further insight into unveiling the mechanism underlying extraordinary fast antler growth rate in particular, as well as the regulatory mechanism of rapid tissue proliferation in general.

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

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.

Development of 16 novel EST-SSR markers for species identification and cross-genus amplification in sambar, sika, and red deer

Deer genera around the globe are threatened by anthropogenic interference. The translocation of alien species and their subsequent genetic introgression into indigenous deer populations is particularly harmful to the species of greatest conservation concern. Products derived from deer, including venison and antler velvet, are also at risk of fraudulent labeling. The current molecular markers used to genetically identify deer species were developed from genome sequences and have limited applicability for cross-species amplification. The absence of efficacious diagnostic techniques for identifying deer species has hampered conservation and wildlife crime investigation efforts. Expressed sequence tag-simple sequence repeat (EST-SSR) markers are reliable tools for individual and species identification, especially in terms of cross-species genotyping. We conducted transcriptome sequencing of sambar (Rusa unicolor) antler velvet and acquired 11,190 EST-SSRs from 65,074 newly assembled unigenes. We identified a total of 55 unambiguous amplicons in sambar (n = 45), which were selected as markers to evaluate cross-species genotyping in sika deer (Cervus nippon, n = 30) and red deer (Cervus elaphus, n = 46), resulting in cross-species amplification rates of 94.5% and 89.1%, respectively. Based on polymorphic information content (>0.25) and genotyping fidelity, we selected 16 of these EST-SSRs for species identification. This marker set revealed significant genetic differentiation based on the fixation index and genetic distance values. Principal coordinate analysis and STRUCTURE analysis revealed distinct clusters of species and clearly identified red-sika hybrids. These markers showed applicability across different genera and proved suitable for identification and phylogenetic analyses across deer species.

Dynamic mRNA and miRNA expression of the head during early development in bighead carp (Hypophthalmichthys nobilis)

Background

Head of fish species, an exquisitely complex anatomical system, is important not only for studying fish evolution and development, but also for economic values. Currently, although some studies have been made on fish growth and body shapes, very limited information is available on the molecular mechanism of head development.

Results

In this study, RNA sequencing (RNA–Seq) and small RNA sequencing (sRNA–Seq) technologies were used to conduct integrated analysis for the head of bighead carp at different development stages, including 1, 3, 5, 15 and 30 Dph (days post hatch). By RNA-Seq data, 26 pathways related to growth and bone formation were identified as the main physiological processes during early development. Coupling this to sRNA–Seq data, we picked out six key pathways that may be responsible for head development, namely ECM receptor interaction, TNF signaling pathway, osteoclast differentiation, PI3K–Akt signaling pathway, Neuroactive ligand–receptor interaction and Jak–STAT signaling pathway. Totally, 114 important candidate genes from the six pathways were obtained. Then we found the top 20 key genes according to the degree value by cytohubba, which regulated cell growth, skeletal formation and blood homeostasis, such as pik3ca, pik3r1, egfr, vegfa, igf1 and itga2b. Finally, we also acquired 19 key miRNAs playing multiple roles in the perfection of various tissues in the head (such as brain, eye and mouth) and mineralization of head bone system, such as let–7e, miR–142a–5p, miR–144–3p, miR–23a–3p and miR–223.

Conclusions

Results of this study will be informative for genetic mechanisms of head development and also provide potential candidate targets for the interaction regulation during early growth in bighead carp.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08387-x.

Comprehensive transcriptome characterization of Grus japonensis using PacBio SMRT and Illumina sequencing

The red-crowned crane (Grus japonensis) is an endangered species distributed across southeast Russia, northeast China, Korea, and Japan. Here, we sequenced for the first time the full-length unreferenced transcriptome of red-crowned crane mixed samples using a PacBio Sequel platform. A total of 359,136 circular consensus sequences (CCS) were obtained via clustering to remove redundancy. A total of 303,544 full-length non-chimeric sequences were identified by judging whether CCS contained 5' and 3' adapters, and the poly(A) tail. Eight samples were sequenced using Illumina, and PacBio sequencing data were corrected according to the collected Illumina data to obtain more accurate full-length transcripts. A total of 4,100 long non-coding RNAs, 13,115 simple sequences repeat loci and 29 transcription factor families were identified. The expression of lncRNAs and TFs in pancreas was lowest comparing with other tissues. Many enriched immune-related transmission pathways (MHC and IL receptors) were identified in the spleen. This study will contribute to a better understanding of the gene structure and post-transcriptional regulatory network, and provide references for future studies on red-crowned cranes.

SWATH-MS Quantitative Proteomic Analysis of Deer Antler from Two Regenerating and Mineralizing Sections

Simple Summary

Deer antler is a unique and astonishing case of annual regeneration in mammalians. Several studies have pointed out the potential for use of velvet antler extract as a nutraceutical supplement, among others, because of its anti-cancer activity. The study of antler regeneration and growth allow us to identify the main proteins and regulatory pathways involved in cell differentiation and regeneration. For this purpose, two sections of antlers (tips and middle sections) using ribs as controls were analyzed from a proteomic point of view. A total of 259 proteins mainly associated with antioxidant mechanisms and Wnt signalling pathways could be responsible for deer antler regeneration and these proteins may be linked to human health benefits. Further studies should be focused on discovering which proteins from velvet antler extracts are associated with these beneficial effects.

Abstract

Antlers are the only organ in the mammalian body that regenerates each year. They can reach growth rates of 1–3 cm/day in length and create more than 20 cm²/day of skin in the antler tips (their growth centers). Previous proteomic studies regarding antlers have focused on antler growth centers (tips) compared to the standard bone to detect the proteins involved in tissue growth. However, proteins of cell differentiation and regeneration will be more accurately detected considering more growing tissues. Thus, we set out to compare proteins expressed in antler tips (the highest metabolism rate and cell differentiation) vs. middle sections (moderate cell growth involving bone calcification), using ribs as controls. Samples were obtained in mid-June with antlers’ phenology corresponding to the middle of their growth period. Quantitative proteomic analysis identified 259 differentially abundant proteins mainly associated with antioxidant metabolic mechanisms, protein formation and Wnt signalling pathway, meanwhile, the mid antler section was linked to blood proteins. The high metabolic rate and subsequent risk of oxidative stress also seem to have resulted in strong antioxidant mechanisms. These results suggest that redox regulation of proteins is a key factor in the model of deer antler regeneration.

Integrated analysis of miRNA and mRNA transcriptomic reveals antler growth regulatory network

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.

Bioactive Molecular Discovery Using Deer Antlers as a Model of Mammalian Regeneration

The ability to activate and regulate stem cells during wound healing and tissue regeneration is a promising field that is resulting in innovative approaches in the field of regenerative medicine. The regenerative capacity of invertebrates has been well documented; however, in mammals, stem cells that drive organ regeneration are rare. Deer antlers are the only known mammalian structure that can annually regenerate to produce a tissue containing dermis, blood vessels, nerves, cartilage, and bone. The neural crest derived stem cells that drive this process result in antlers growing at up to 2 cm/day. Deer antlers thus provide superior attributes compared to lower-order animal models, when investigating the regulation of stem cell-based regeneration. Antler stem cells can therefore be used as a model to investigate the bioactive molecules, biological processes, and pathways involved in the maintenance of a stem cell niche, and their activation and differentiation during organ formation. This review examines stem cell-based regeneration with a focus on deer antlers, a neural crest stem cell-based mammalian regenerative structure. It then discusses the omics technical platforms highlighting the proteomics approaches used for investigating the molecular mechanisms underlying stem cell regulation in antler tissues.

Antiquity and fundamental processes of the antler cycle in Cervidae (Mammalia)

The origins of the regenerative nature of antlers, being branched and deciduous apophyseal appendages of frontal bones of cervid artiodactyls, have long been associated with permanent evolutionary precursors. In this study, we provide novel insight into growth modes of evolutionary early antlers. We analysed a total of 34 early antlers affiliated to ten species, including the oldest known, dating from the early and middle Miocene (approx. 18 to 12 million years old) of Europe. Our findings provide empirical data from the fossil record to demonstrate that growth patterns and a regular cycle of necrosis, abscission and regeneration are consistent with data from modern antlers. The diverse histological analyses indicate that primary processes and mechanisms of the modern antler cycle were not gradually acquired during evolution, but were fundamental from the earliest record of antler evolution and, hence, explanations why deer shed antlers have to be rooted in basic histogenetic mechanisms. The previous interpretation that proximal circular protuberances, burrs, are the categorical traits for ephemerality is refuted.

Comparative Metabolomics Study Revealed Difference in Central Carbon Metabolism between Sika Deer and Red Deer Antler

The antler regeneration has been well studied for the past two decades and adopted in the regenerative medicine model for studying on developmental biology. Despite our growing knowledge of functional molecules regulating antler regeneration, we still do not know whether antler from different deer species possess the exact same mechanism or not. Our previous comparative study between sika deer and red deer suggests that the metabolic pathways between them are profoundly different based on protein level. Therefore, the metabolomic technology is used to identify and quantify the metabolites in antler samples, providing interesting insights into differential metabolite profile of antlers between sika deer and red deer. The distinct metabolic characteristics of sika deer compared to red deer provide an opportunity to explain why the red deer antler with a larger size. The enrichment analysis of differential metabolites showed that three pathways including glycine and serine metabolism, methionine metabolism, and pterine biosynthesis had a significant difference between two antler groups.

Development of novel EST microsatellite markers for genetic diversity analysis and correlation analysis of velvet antler growth characteristics in Sika deer

BACKGROUND

Sika deer is one of the most popular and valued animals in China. However, few studies have been conducted on the microsatellite of Sika deer, which has hampered the progress of genetic selection breeding. To develop and characterize a set of microsatellites for Sika deer which provide helpful information for protection of Sika deer natural resources and effectively increase the yield and quantity of velvet antler.

RESULTS

We conducted a transcriptome survey of Sika deer using next-generation sequencing technology. One hundred eighty-two thousand two hundred ninety-five microsatellite markers were identified in the transcriptome, 170 of 200 loci were successfully amplified across panels of 140 individuals from Shuangyang Sika deer population. And 29 loci were found to be obvious polymorphic. Number of alleles is from 3 to 14. The expected heterozygosity ranged from 0.3087 to 0.7644. The observed heterozygosity ranged from 0 to 0.7698. The polymorphism information content values of those microsatellites varied ranged from 0.2602 to 0.7507. The marker-trait association was tested for 6 important and kernel characteristics of two-branched velvet antler in Shuangyang Sika deer through one-way analysis of variance. The results showed that marker-trait associations were identified with 8 different markers, especially M009 and M027.

CONCLUSIONS

This study not only provided a large scale of microsatellites which were valuable for future genetic mapping and trait association in Sika deer, but also offers available information for molecular breeding in Sika deer.

Whole Transcriptome Analysis of Mesenchyme Tissue in Sika Deer Antler Revealed the CeRNAs Regulatory Network Associated With Antler Development

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

Reclassification of velvet antler portions following transcriptomic analysis

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.

Full-length transcriptome sequences of Agropyron cristatum facilitate the prediction of putative genes for thousand-grain weight in a wheat-A. cristatum translocation line

Background

Agropyron cristatum (L.) Gaertn. (2n = 4x = 28; genomes PPPP) is a wild relative of common wheat (Triticum aestivum L.) and provides many desirable genetic resources for wheat improvement. However, there is still a lack of reference genome and transcriptome information for A. cristatum, which severely impedes functional and molecular breeding studies.

Results

Single-molecule long-read sequencing technology from Pacific Biosciences (PacBio) was used to sequence full-length cDNA from a mixture of leaves, roots, stems and caryopses and constructed the first full-length transcriptome dataset of A. cristatum, which comprised 44,372 transcripts. As expected, the PacBio transcripts were generally longer and more complete than the transcripts assembled via the Illumina sequencing platform in previous studies. By analyzing RNA-Seq data, we identified tissue-enriched transcripts and assessed their GO term enrichment; the results indicated that tissue-enriched transcripts were enriched for particular molecular functions that varied by tissue. We identified 3398 novel and 1352 A. cristatum-specific transcripts compared with the wheat gene model set. To better apply this A. cristatum transcriptome, the A. cristatum transcripts were integrated with the wheat genome as a reference sequence to try to identify candidate A. cristatum transcripts associated with thousand-grain weight in a wheat-A. cristatum translocation line, Pubing 3035.

Conclusions

Full-length transcriptome sequences were used in our study. The present study not only provides comprehensive transcriptomic insights and information for A. cristatum but also proposes a new method for exploring the functional genes of wheat relatives under a wheat genetic background. The sequence data have been deposited in the NCBI under BioProject accession number PRJNA534411.

Comparative proteomics analysis reveals the difference during antler regeneration stage between red deer and sika deer

Deer antler, as the only mammalian regenerative appendage, provides an optimal model to study regenerative medicine. Antler harvested from red deer or sika deer were mainly study objects used to disclose the mechanism underlying antler regeneration over past decades. A previous study used proteomic technology to reveal the signaling pathways of antler stem cell derived from red deer. Moreover, transcriptome of antler tip from sika deer provide us with the essential genes, which regulated antler development and regeneration. However, antler comparison between red deer and sika deer has not been well studied. In our current study, proteomics were employed to analyze the biological difference of antler regeneration between sika deer and red deer. The proteomics profile was completed by searching the UniProt database, and differentially expressed proteins were identified by bioinformatic software. Thirty-six proteins were highly expressed in red deer antler, while 144 proteins were abundant in sika deer. GO and KEGG analysis revealed that differentially expressed proteins participated in the regulation of several pathways including oxidative phosphorylation, ribosome, extracellular matrix interaction, and PI3K-Akt pathway.