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.

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.

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

BACKGROUND: Previous investigations of phylogeny in Cervus recovered many clades without whole genomic support. METHODS: In this study, the genetic diversity and phylogeny of 5 species (21 subspecies/populations from C. unicolor, C. albirostris, C. nippon, C. elaphus and C. eldii) in the genus Cervus were analyzed using reduced-representation genome sequencing. RESULTS: A total of 197,543 SNPs were identified with an average sequencing depth of 16 x. A total of 21 SNP matrices for each subspecies/population and 1 matrix for individual analysis were constructed, respectively. Nucleotide diversity and heterozygosity analysis showed that all 21 subspecies/populations had different degrees of genetic diversity. C. eldii, C. unicolor and C. albirostris showed relatively high expected and observed heterozygosity, while observed heterozygosity in C. nippon was the lowest, indicating there was a certain degree of inbreeding rate in these subspecies/populations. Phylogenetic ML tree of all Cervus based on the 21 SNP matrices showed 5 robustly supported clades that clearly separate C. eldii, C. unicolor, C. albirostris, C. elaphus and C. nippon. Within C. elaphus clade, 4 subclades were well differentiated and statistically highly supported: C. elaphus (New Zealand), C. e. yarkandensis, C. c. canadensis and the other grouping the rest of C. canadensis from China. In the C. nippon clade, 2 well-distinct subclades corresponding to C. n. aplodontus and other C. nippon populations were separated. Phylogenetic reconstruction indicated that the first evolutionary event of the genus Cervus occurred approximately 7.4 millions of years ago. The split between C. elaphus and C. nippon could be estimated at around 3.6 millions of years ago. Phylogenetic ML tree of all samples based on individual SNP matrices, together with geographic distribution, have shown that there were 3 major subclades of C. elaphus and C. canadensis in China, namely C. e. yarkandensis (distributed in Tarim Basin), C. c. macneilli/C. c. kansuensis/C. c. alashanicus (distributed in middle west of China), and C. c. songaricus/C. c. sibiricus (distributed in northwest of China). Among them, C. e. yarkandensis was molecularly the most primitive subclade, with a differentiation dating back to 0.8-2.2 Myr ago. D statistical analysis showed that there was high probability of interspecific gene exchange between C. albirostris and C. eldii, C. albirostris and C. unicolor, C. nippon and C. unicolor, and there might be 2 migration events among 5 species in the genus Cervus. CONCLUSIONS: Our results provided new insight to the genetic diversity and phylogeny of Cervus deer. In view of the current status of these populations, their conservation category will need to be reassessed.

[Present situation and analysis of deer industry standardization at home and abroad].

Deer is valuable all over the body,which is rich in nutritional value and medicinal value. Deer breeding and processing are very advanced in North America and New Zealand where many related standards have been published. The development of Chinese deer industry lack standard and normal management,neither standards' number nor coverage area formed complete frame structure. The international standards like Panax ginseng and P. notoginseng were more lacked. This paper makes a classification statistics on standardization organizations at home and abroad,foreign standards,Chinese national standards,industry standards,local standards and enterprise standards. The classes,contents,ages,implementation and promotion and demonstration area construction of standards were compared and analyzed. We found Chinese deer industry standards were deficient in coverage,uniformity,innovation,repeatability and support. And we give advises for the construction of industry quality standard system,organizational mobility and ideology of consumers,hoping to boost the standard construction and promote international competitiveness of Chinese deer industry.

Multi-Fault Diagnosis of Gearbox Based on Improved Multipoint Optimal Minimum Entropy Deconvolution.

Under complicated conditions, the extraction of a multi-fault in gearboxes is difficult to achieve. Due to improper selection of methods, leakage diagnosis or misdiagnosis will usually occur. Ensemble Empirical Mode Decomposition (EEMD) often causes energy leakage due to improper selection of white noise during signal decomposition. Considering that only a single fault cycle can be extracted when MOMED (Multipoint Optimal Minimum Entropy Deconvolution) is used, it is necessary to perform the sub-band processing of the compound fault signal. This paper presents an adaptive gearbox multi-fault-feature extraction method based on Improved MOMED (IMOMED). Firstly, EEMD decomposes the signal adaptively and selects the intrinsic mode functions with strong correlation with the original signal to perform FFT (Fast Fourier transform); considering the mode-mixing phenomenon of EEMD, reconstruct the intrinsic mode functions with the same timescale, and obtain several intrinsic mode functions of the same scale to improve the entropy of fault features. There is a lot of white noise in the original signal, and EEMD can improve the signal-to-noise ratio of the original signal. Finally, through the setting of different noise-reduction intervals to extract fault features through MOMED. The proposed method is compared with EEMD and VMD (Variational Mode Decomposition) to verify its feasibility.

Epidermal Growth Factor Promotes Proliferation and Migration of Follicular Outer Root Sheath Cells via Wnt/ß-Catenin Signaling.

BACKGROUND/AIMS: To investigate the effect and molecular mechanism of EGF on the growth and migration of hair follicle outer root sheath (ORS) cells. METHODS: Intact anagen hair follicles were isolated from mink skin and cultured with EGF in vitro to measure ORS daily growth. Meanwhile, purified primary ORS cells were treated or transfected with EGF, and their proliferation and migration were assessed by MTT assay and transwell assay, respectively. The signaling pathway downstream of EGF was characterized by using the Wnt/ß-catenin signaling inhibitor, XAV-939. RESULTS: EGF of 2-20 ng/ml, not higher or lower, promoted the growth of follicular ORS in vitro. EGF treatment or overexpression promoted the proliferation and migration of ORS cells. Moreover, EGF stimulation induced nuclear translocation of ß-catenin, and upregulated the expression of Wnt10b, ß-catenin, EGF receptor and SOX9. Inhibition of Wnt/ß-catenin signaling by XAV-939 significantly reduced the basal and EGF-enhanced proliferation and migration of ORS cells. In addition, a number of follicle-regulatory genes, such as Survivin, Msx2 and SGK3, were upregulated by EGF in the ORS cells, which was also inhibited by XAV-939. CONCLUSION: EGF promotes the proliferation and migration of ORS cells and modulates the expression of several follicle-regulatory genes via Wnt/ß-catenin signaling.

Transcriptome analysis of sika deer in China.

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.

Bacterial community composition and fermentation patterns in the rumen of sika deer (Cervus nippon) fed three different diets.

Sika deer (Cervus nippon) rely on microorganisms living in the rumen to convert plant materials into chemical compounds, such as volatile fatty acids (VFAs), but how the rumen bacterial community is affected by different forages and adapt to altered diets remains poorly understood. The present study used 454-pyrosequencing of bacterial 16S ribosomal RNA (rRNA) genes to examine the relationship between rumen bacterial diversity and metabolic phenotypes using three sika deer in a 3 × 3 latin square design. Three sika deer were fed oak leaves (OL), corn stover (CS), or corn silage (CI), respectively. After a 7-day feeding period, when compared to the CS and CI groups, the OL group had a lower proportion of Prevotella spp. and a higher proportion of unclassified bacteria belonging to the families Succinivibrionaceae and Paraprevotellaceae (P<0.05). Meanwhile, the concentration of isobutyrate was significantly lower (P<0.05) in the OL group than in the CS and CI groups. There was no significant change of dominant bacterial genera in the OL group after 28 days of feeding. Conversely, total volatile fatty acids (TVFAs) showed an increase after 28 days of feeding, mainly due to the increasing of acetate, propionate, and valerate (P<0.05). The interplay between bacteria and metabolism in the OL group differed from that in the CS and CI groups, especially for the interaction of TVFAs and acetate/propionate. Overall, the current study suggested that Prevotella spp. played critical roles in the fermentation of feed in the rumen of sika deer. However, the differences in interplay patterns between rumen bacterial community composition and metabolic phenotypes were altered in the native and domesticated diets indicating the changed fermentation patterns in the rumen of sika deer.

Novel amdoparvovirus infecting farmed raccoon dogs and arctic foxes.

A new amdoparvovirus, named raccoon dog and fox amdoparvovirus (RFAV), was identified in farmed sick raccoon dogs and arctic foxes. Phylogenetic analyses showed that RFAV belongs to a new species within the genus Amdoparvovirus of the family Parvoviridae. An RFAV strain was isolated in Crandell feline kidney cell culture.

The First High-quality Reference Genome of Sika Deer Provides Insights into 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 genome of sika deer. We generated the first high-quality, chromosome-level genome assembly of sika deer and measured the correlation between tannin intake and RNA expression in 15 tissues through 180 experiments. Comparative genome analyses showed that the UGT and CYP gene families are functionally involved in the adaptation of sika deer to high-tannin food, especially the expansion of the UGT family 2 subfamily B of UGT genes. The first chromosome-level assembly and genetic characterization of the tolerance to a highly toxic diet suggest that the sika deer genome may serve as an essential resource for understanding evolutionary events and tannin adaptation. Our study provides a paradigm of comparative expressive genomics that can be applied to the study of unique biological features in non-model animals.

Rabies virus expressing dendritic cell-activating molecules enhances the innate and adaptive immune response to vaccination.

Our previous studies indicated that recruitment and/or activation of dendritic cells (DCs) is important in enhancing the protective immune responses against rabies virus (RABV) (L. Zhao, H. Toriumi, H. Wang, Y. Kuang, X. Guo, K. Morimoto, and Z. F. Fu, J. Virol. 84:9642-9648). To address the importance of DC activation for RABV vaccine efficacy, the genes for several DC recruitment and/or activation molecules, e.g., granulocyte-macrophage colony-stimulating factor (GM-CSF), macrophage-derived chemokine (MDC), and macrophage inflammatory protein 1α (MIP-1α), were individually cloned into RABV. The ability of these recombinant viruses to activate DCs was determined in vitro and in vivo. Infection of mouse bone marrow-derived DCs with each of the recombinant viruses resulted in DC activation, as shown by increased surface expression of CD11c and CD86 as well as an increased level of alpha interferon (IFN-α) production compared to levels observed after infection with the parent virus. Intramuscular infection of mice with each of the viruses recruited and/or activated more DCs and B cells in the periphery than infection with the parent virus, leading to the production of higher levels of virus-neutralizing antibodies. Furthermore, a single immunization with recombinant RABV expressing GM-CSF or MDC protected significantly more mice against intracerebral challenge with virulent RABV than did immunization with the parental virus. Yet, these viruses did not show more virulence than the parent virus, since direct intracerebral inoculation with each virus at up to 1 × 10(7) fluorescent focus units each did not induce any overt clinic symptom, such as abnormal behavior, or any neurological signs. Together, these data indicate that recombinant RABVs expressing these molecules activate/recruit DCs and enhance protective immune responses.

Mapping the morphogenetic potential of antler fields through deleting and transplanting subregions of antlerogenic periosteum in sika deer (Cervus nippon).

Morphogenetic fields are a localised and regionally regulated group of cells capable of responding to signals leading to the development of organs. In this study, we sought to determine if antlers develop from such a field. We divided antler fields into four subregions: anterior, posterior, medial and lateral. The antlerogenic periosteum (AP) in each subregion (half of the AP) was deleted and then transplanted into an ectopic site. Antlers form from the cells exclusively residing in the AP, which is located in an antler field. The morphogenetic potential of each subregion was assessed by the antler growth from both the defective field and the transplantation site. The results showed that when the AP anterior half was absent, the fields formed antlers missing the first tine, whereas when the anterior half was present, the ectopic sites regenerated antlers containing the first tine. When the medial half was deleted, the fields could only grow spike antlers, and when the medial half was present, the ectopic sites developed branched antlers. In contrast, the antler fields were able to compensate the defects caused by ablation of the posterior or the lateral half to form relatively normal antlers; and the ectopic sites containing these grafted halves only formed spike antlers. Therefore, antler morphogenetic information was primarily held in the AP anterior-medial halves. This study substantiates the presence of morphogenetic fields in regulating the distinct pattern of antler growth, and demonstrates that antler development is a useful model for the study of morphogenetic fields.

Comparison of miRNA and mRNA Expression in Sika Deer Testes With Age.

To elucidate the complex physiological process of testis development and spermatogenesis in Sika deer, this study evaluated the changes of miRNA and mRNA profiles in the four developmental stages of testis in the juvenile (1-year-old), adolescence (3-year-old), adult (5-year-old), and aged (10-year-old) stages. The results showed that a total of 198 mature, 66 novel miRNAs, and 23,558 differentially expressed (DE) unigenes were obtained; 14,918 (8,413 up and 6,505 down), 4,988 (2,453 up and 2,535 down), and 5,681 (2,929 up and 2,752 down) DE unigenes, as well as 88 (43 up and 45 down), 102 (44 up and 58 down), and 54 (18 up and 36 down) DE miRNAs were identified in 3- vs. 1-, 5- vs. 3-, and 10- vs. 5-year-old testes, respectively. By integrating miRNA and mRNA expression profiles, we predicted 10,790 mRNA-mRNA and 69,883 miRNA-mRNA interaction sites. The target genes were enriched by GO and KEGG pathways to obtain DE mRNA (IGF1R, ALKBH5, Piwil, HIF1A, BRDT, etc.) and DE miRNA (miR-140, miR-145, miR-7, miR-26a, etc.), which play an important role in testis development and spermatogenesis. The data show that DE miRNAs could regulate testis developmental and spermatogenesis through signaling pathways, including the MAPK signaling pathway, p53 signaling pathway, PI3K-Akt signaling pathway, Hippo signaling pathway, etc. miR-140 was confirmed to directly target mutant IGF1R-3'UTR by the Luciferase reporter assays. This study provides a useful resource for future studies on the role of miRNA regulation in testis development and spermatogenesis.

The combination of SMRT sequencing and Illumina sequencing highlights organ-specific and age-specific expression patterns of miRNAs in Sika Deer.

Due to the lack of high-quality Sika Deer (Cervus nippon) transcriptome and sRNAome across multiple organs or development stages, it is impossible to comprehensively analyze the mRNA and miRNA regulatory networks related to growth, development and immunity response. In this study, we used single molecule-real time sequencing (SMRT-seq) and Illumina sequencing methods to generate transcriptome and sRNAome from ten tissues and four age groups of Sika Deer to help us understand molecular characteristics and global miRNA expression profiles. The results showed that a total of 240,846 consensus transcripts were generated with an average length of 2,784 bp. 4,329 Transcription factors (TFs), 109,000 Simple Sequence Repeats (SSRs) and 18,987 Long non-coding RNAs (LncRNAs) were identified. Meanwhile, 306 known miRNAs and 143 novel miRNAs were obtained. A large number of miRNAs showed organ-specific and age-specific differential expression patterns. In particular, we found that the organ-specific miRNAs were enriched in the brain, some of which shared only between the brain and adrenal. These miRNAs were involved in maintaining specific functions within the brain and adrenal. By constructing miRNA96mRNA interaction networks associated with Sika Deer immunity, we found that miRNAs (miR-148a, miR-26a, miR-214, let-7b, etc.) and mRNAs (CD6, TRIM38, C3, CD163, etc.) might play an important role in the immune response of Sika Deer spleen. Together, our study generated an improved transcript annotation for Sika Deer by SMRT-seq and revealed the role of miRNA in regulating the growth, development and immunity response of Sika Deer.

Integrated histopathological, lipidomic, and metabolomic profiles reveal mink is a useful animal model to mimic the pathogenicity of severe COVID-19 patients.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is transmitted on mink farms between minks and humans in many countries. However, the systemic pathological features of SARS-CoV-2-infected minks are mostly unknown. Here, we demonstrated that minks were largely permissive to SARS-CoV-2, characterized by severe and diffuse alveolar damage, and lasted at least 14 days post inoculation (dpi). We first reported that infected minks displayed multiple organ-system lesions accompanied by an increased inflammatory response and widespread viral distribution in the cardiovascular, hepatobiliary, urinary, endocrine, digestive, and immune systems. The viral protein partially co-localized with activated Mac-2+ macrophages throughout the body. Moreover, we first found that the alterations in lipids and metabolites were correlated with the histological lesions in infected minks, especially at 6 dpi, and were similar to that of patients with severe and fatal COVID-19. Particularly, altered metabolic pathways, abnormal digestion, and absorption of vitamins, lipids, cholesterol, steroids, amino acids, and proteins, consistent with hepatic dysfunction, highlight metabolic and immune dysregulation. Enriched kynurenine in infected minks contributed to significant activation of the kynurenine pathway and was related to macrophage activation. Melatonin, which has significant anti-inflammatory and immunomodulating effects, was significantly downregulated at 6 dpi and displayed potential as a targeted medicine. Our data first illustrate systematic analyses of infected minks to recapitulate those observations in severe and fetal COVID-19 patients, delineating a useful animal model to mimic SARS-CoV-2-induced systematic and severe pathophysiological features and provide a reliable tool for the development of effective and targeted treatment strategies, vaccine research, and potential biomarkers.

Antler development was inhibited or stimulated by cryosurgery to periosteum or skin in a central antlerogenic region respectively.

Antler development is triggered by interactions between antler stem cells resident in the antlerogenic periosteum (AP) and the niche cells in the upper portion of overlying skin mediated by diffusible molecules. These interactive cell populations are interposed by the lower portion of the skin and the subcutaneous loose connective tissue (SLCT). It is known that mechanical deletion of just the central AP (having an area equivalent to the size of a pedicle base) by cutting through the skin and SLCT effectively stimulates the marginal AP to initiate antler development. This study was designed to investigate whether the SLCT layer plays a role in antler development by acting as a physical barrier. The results showed that the marginal AP failed to give rise to an antler after the central AP was cryosurgically destroyed with the preservation of the collagen structure of the SLCT. Furthermore, antler development was significantly advanced when the collagen structures of the skin and SLCT layers were substantially attenuated by repeated sprays with liquid nitrogen while keeping the central AP intact. Therefore, we conclude that the interposing SLCT layer acts as a physical barrier between antler stem cells and the niche cell types, and that timing of antler development is primarily controlled by the permeability of the SLCT layer to the putative interactive diffusible molecules.

Stem cells responsible for deer antler regeneration are unable to recapitulate the process of first antler development-revealed through intradermal and subcutaneous tissue transplantation.

Antlers offer a unique model for the study of whether regeneration recapitulates development in a mammalian organ. Research, to date, supports the full recapitulation in antler, but a recent report that subcutaneously transplanted (ST) pedicle periosteum (PP) failed to induce that ectopic antler formation could argue against recapitulation, as antlerogenic periosteum (AP) can readily do so. However, it was not clear in that study whether the result was caused by inability of the PP to interact with the skin or owing to failure to create the required close contact to it. This study was designed to clarify this uncertainty by adopting intradermal transplantation (IT) to achieve the required close contact without the need for significant mass expansion. The results showed that IT of 1/8 of the original AP mass or more was sufficient for antler induction, whereas ST of 1/4-AP or less could not do so within 2 years. The minimum amount of AP required for antler induction using the IT approach was somewhere between 1/8 and 1/12-AP (<30 mg). The results further demonstrated that IT of 62-84 mg PP failed to induce ectopic antler formation, even if the PP had fused with the surrounding skin. Because this mass of PP was 2-3 times the minimum amount of AP required for antler induction, we conclude that PP does not recapitulate AP in induction of ectopic antler development. It is likely that PP has been restricted for antler regeneration and lost the potential to initiate antler development.

Identification of key genes associated with spermatogenesis arrest in fox hybrids using weighted gene co-expression network analysis.

The silver fox and the blue fox represent different genera, but produce viable offspring. Although these hybrids show obvious heterosis, they are completely sterile due to spermatogenic arrest at the early stages of spermatogenesis, especially mitosis and meiosis I; the hybrids produce few spermatogonia and primary spermatocytes, and no secondary spermatocytes. Although the mechanisms of spermatogenic arrest have been well investigated, transcriptomic differences between hybrid and the pure-species testes have not clarified. In the present study, we used RNA sequencing (RNA-Seq) to generate testicular transcriptomic profiles for silver foxes, blue foxes, and reciprocal hybrids during the pre-breeding period and the breeding season. In total, 1,344,022 transcripts (≥200 bp) were generated; 1,057,724 genes were obtained; and 33,423 genes were shown to have fragments per kilobase of transcript per million mapped reads (FPKM) > 0.3. To identify the hub genes associated with spermatogenesis arrest, weighted gene co-expression network analysis (WGCNA) was used. Nine modules were explored. Genes in only a single module were consistently downregulated in the hybrids as compared to the pure species; these genes were significantly associated with fox hybrid male infertility. Six of the genes in this module (CATSPERD, DMRTC2, RNF17, NME5, SPEF2, SPINK2) also play key roles in mitosis and meiosis during spermatogenesis. Therefore, these six genes might be associated with fox hybrid male infertility.

Altered miRNA and mRNA Expression in Sika Deer Skeletal Muscle with Age.

Studies of the gene and miRNA expression profiles associated with the postnatal late growth, development, and aging of skeletal muscle are lacking in sika deer. To understand the molecular mechanisms of the growth and development of sika deer skeletal muscle, we used de novo RNA sequencing (RNA-seq) and microRNA sequencing (miRNA-seq) analyses to determine the differentially expressed (DE) unigenes and miRNAs from skeletal muscle tissues at 1, 3, 5, and 10 years in sika deer. A total of 51,716 unigenes, 171 known miRNAs, and 60 novel miRNAs were identified based on four mRNA and small RNA libraries. A total of 2,044 unigenes and 11 miRNAs were differentially expressed between adolescence and juvenile sika deer, 1,946 unigenes and 4 miRNAs were differentially expressed between adult and adolescent sika deer, and 2,209 unigenes and 1 miRNAs were differentially expressed between aged and adult sika deer. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that DE unigenes and miRNA were mainly related to energy and substance metabolism, processes that are closely associate with the growth, development, and aging of skeletal muscle. We also constructed mRNA-mRNA and miRNA-mRNA interaction networks related to the growth, development, and aging of skeletal muscle. The results show that mRNA (Myh1, Myh2, Myh7, ACTN3, etc.) and miRNAs (miR-133a, miR-133c, miR-192, miR-151-3p, etc.) may play important roles in muscle growth and development, and mRNA (WWP1, DEK, UCP3, FUS, etc.) and miRNAs (miR-17-5p, miR-378b, miR-199a-5p, miR-7, etc.) may have key roles in muscle aging. In this study, we determined the dynamic miRNA and unigenes transcriptome in muscle tissue for the first time in sika deer. The age-dependent miRNAs and unigenes identified will offer insights into the molecular mechanism underlying muscle development, growth, and maintenance and will also provide valuable information for sika deer genetic breeding.

Development of a nude mouse model for the study of antlerogenesis--mechanism of tissue interactions and ossification pathway.

In a previous study (Li et al., 2001. J Exp Zool 290:18-30) a nude mouse model was established to investigate deer antler development. In that study we found nude mice could support the singularly implanted antlerogenic periosteum (AP) to form pedicle-like, but not antler-like, bony protuberances. To further develop the model and at the same time to use the updated model for the investigation of antler formation, three experiments were carried out in this study. The results showed that (1) antler-like protuberances were successfully induced on the nude mouse heads via subcutaneous co-transplantation of AP and deer skin, and subsequent exposure through wounding of the deer xenografts; (2) deer skin epidermis and its attached half thickness of dermis were sufficient to interact with the AP, and the interactions were capable of transforming adult scalp skin into velvet; (3) the putative initial inductive molecules were primarily derived from the AP cellular layer, rather than fibrous layer; (4) initiation of the ossification center in the avascular cartilage of each mouse "antler" took place via metaplasia, rather than classical endochondral ossification. Further research is required to identify means for effective stimulation of calcification of the "mouse antlers" in order to create the opportunity to investigate antler regeneration using the nude mouse model. Overall, the nude mouse model, once further developed, has the potential to become a powerful tool to study underlying mechanism of antlerogenesis and organogenesis/regeneration in general.