Isolation and characterization of tannin-degrading bacteria from the rumen of wild Hokkaido sika deer (Cervus nippon yezoensis).

We isolated tannin-degrading bacteria from the rumen of wild Hokkaido sika deer and characterized their phylogeny and tannase activity in relation to sample sources. The condensed tannin level was higher in all deer rumen samples (n = 20) than in forage-fed cattle rumen samples (n = 6), whereas no hydrolyzable tannins were detected in any of the rumen samples. Rumen bacteria were enumerated on nonselective brain heart infusion (BHI) agar medium and then transferred onto tannic acid-containing BHI agar plates to screen for bacteria only showing growth (tannin-resistant bacteria) and those showing both growth and a clear zone (tannin-degrading bacteria). Summer samples provided only tannin-resistant bacteria, none of which showed tannin-degrading activity. Although winter samples also provided tannin-resistant bacteria, most isolates exhibited tannin-degrading activity. A total of 70 isolates exhibiting tannin-degrading activity were classified as Streptococcus bovis group based on 16S rRNA gene sequencing and further classified into two groups, either group A or group B. Group A consisted of isolates showing weak tannase activity, whereas group B included a majority of the isolates exhibiting high tannase activity. These results suggest that wild Hokkaido sika deer develop tannin-degrading Streptococcus in the rumen during winter, which allows access to woody food materials rich in tannins.

Deer antlers: the fastest growing tissue with least cancer occurrence.

Deer antlers are a bony organ solely able to acquired distinct unique attributes during evolution and all these attributes are against thus far known natural rules. One of them is as the fastest animal growing tissue (2 cm/day), they are remarkably cancer-free, despite high cell division rate. Although tumor-like nodules on the long-lived castrate antlers in some deer species do occur, but they are truly benign in nature. In this review, we tried to find the answer to this seemingly contradictory phenomenon based on the currently available information and give insights into possible clinic application. The antler growth center is located in its tip; the most intensive dividing cells are resident in the inner layer of reserve mesenchyme (RM), and these cells are more adopted to osteosarcoma rather than to normal bone tissues in gene expression profiles but acquire their energy mainly through aerobic oxidative phosphorylation pathway. To counteract propensity of neoplastic transformation, antlers evolved highly efficient apoptosis exactly in the RM, unparalleled by any known tissues; and annual wholesale cast to jettison the corps. Besides, some strong cancer suppressive genes including p53 cofactor genes and p53 regulator genes are highly positively selected by deer, which would have certainly contributed to curb tumorigenesis. Thus far, antler extracts and RM cells/exosomes have been tried on different cancer models either in vitro or in vivo, and all achieved positive results. These positive experimental results together with the anecdotal folklore that regular consumption of velvet antler is living with cancer-free would encourage us to test antlers in clinic settings.

Characterization of complete mitochondrial genome of three Horse flies of the genus Tabanus (Diptera: Tabanidae): comparative analysis.

BACKGROUND: Tabanidae (Horse-Flies or Deer-Flies) are one of the most economically important as well as medically significant haematophagous insect family within the order Diptera. Members of this group are also responsible for the mortality of substantial number of live-stock every year. Due to their pathogen transmission potential and vector competencies makes them an important insect group to study. Till now, mitochondrial genome of 18 species of tabanids were available. METHODS AND RESULTS: The complete mitogenome of three species T. diversifrons (15,809 bp), T. rubidus (15,878 bp) and T. tenens (15,872 bp) were generated by Next generation sequencing method. They consist 37 genes, with a positive AT skew and a negative GC skew. The gene order of these three species is similar to the typical gene arrangement of infra-order Tabanomorpha. Most of the tRNAs showed typical clover-leaf secondary structure except trnS1, which lacks the DHU arm. The sliding window analysis showed that the nad4L is the most conserved while atp8, and nad6 are the most variable genes. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated that all PCGs under the purifying selection. Phylogeny revealed Chrysops and Haematopota are monophyletic while species of Hybomitra are nested within the polyphyletic clade of Tabanus. T. diversifrons exhibits sister relationship with Atylotus miser. Two morphologically divergent species T. rubidus and T. tenens are found to be genetically similar and indistinguishable by mitochondrial genome. CONCLUSIONS: The hypervariable genes like atp8 and nad6 can be used as molecular markers for the identification of recently diverged lineages of family Tabanidae. Further, to address uncertainties arising from the two morphological divergent species, it is imperative to obtain data from nuclear gene markers.

Transmission and lesion progression of treponeme-associated hoof disease in captive elk (Cervus canadensis).

Treponeme-associated hoof disease (TAHD) is a debilitating disease of free-ranging elk (Cervus canadensis) in the northwestern U.S. While treponemes are associated with lesions, the etiology and transmissibility between elk are unknown. Our objective was to determine whether the disease can be environmentally transmitted to captive elk. Four individually housed treatment elk and 2 control elk were challenged with soil mixed with inoculum prepared from free-ranging elk hooves from TAHD-positive elk or autoclaved hooves from normal elk, respectively. The inoculum for each group was applied to the interdigital space and added to pre-existing soil in each pen. Eight challenges were conducted at 1-4-week intervals and lesion development was assessed during a 138-day challenge period that was followed by a 170-day monitoring period to document lesion progression. All treatment elk, but no control elk, developed gross and histologic lesions consistent with TAHD. Treponema phylotypes similar to those in bovine digital dermatitis in cattle were detected using 16S rRNA gene amplicon sequencing from lesions in all treatment elk, but no control elk, during the challenge period. Lesions progressed from ulcerations in the interdigital space to extensive ulceration and underrunning of the hoof capsule by 35 and 173 days following the initial inoculation, respectively. Lameness in treatment elk was correlated with lesion development (R = 0.702, p≤0.001), and activity of infected elk was reduced during the challenge (p≤0.001) and monitoring periods (p = 0.004). Body condition was significantly lower in treatment than control elk 168 days following the initial inoculation (p = 0.05) and at each individual elk's study endpoint (p = 0.006). Three of 4 treatment elk were euthanized when they reached humane endpoints, and one elk recovered. These results provide direct evidence that TAHD is a transmissible infectious disease in elk. As such, actions that reduce transmission risk can support disease management and prevention.

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.

Retinoic acid alters metalloproteinase action in red deer antler stem cells.

Metalloproteinases (MMP)s regulate developmental processes, control angiogenesis and wound healing, participate in the formation of immune receptors, and are expressed in stem cells. Retinoic acid (RA) is a potential modulator of these proteinases. The aim was to determine (1) MMPs' action in antler stem cells (ASCs) before and after differentiation into adipo-, osteo-, and chondrocytes and (2) the effect of RA on modifying MMP action in ASCs. Antler tissue from pedicle was collected approximately 40 days after antler casting, post mortem from healthy breeding five year old males (N = 7). The cells were isolated from the pedicle layer of periosteum after skin separation and cultured. The pluripotency of the ASCs was evaluated by mRNA expression for NANOG, SOX2, and OCT4. ASCs were stimulated with RA (100nM) and differentiated for 14 days. The MMP (1-3) and TIMP(1-3) (tissue inhibitor of MMPs) mRNA expression was determined in the ASCs, their concentrations in the ASCs and the medium after RA stimulation as well as profiles of mRNA expression for MMPs: 1-3 and TIMPs: 1-3 during differentiation of ASC to osteocytes, adipocytes and chondrocytes. RA increased MMP-3 and TIMP-3 mRNA expression and output (P < 0.05) and not influenced on MMP-1 and TIMP-1 mRNA expression and output in ASC (P > 0.05). Depending on differentiation of ASC to osteocytes, adipocytes or chondrocytes, MMPs`and TIMPs`expression profile fluctuates for all studied proteases and its inhibitors. The studies demand continuation considering the role of proteases in stem cells physiology and differentiation. The results may be relevant for the study of cellular processes during the cancerogenesis of tumor stem cells.

Identification of a novel risk factor for chronic wasting disease (CWD) in elk: S100G single nucleotide polymorphism (SNP) of the prion protein gene (PRNP).

Prion diseases are fatal and malignant infectious encephalopathies induced by the pathogenic form of prion protein (PrPSc) originating from benign prion protein (PrPC). A previous study reported that the M132L single nucleotide polymorphism (SNP) of the prion protein gene (PRNP) is associated with susceptibility to chronic wasting disease (CWD) in elk. However, a recent meta-analysis integrated previous studies that did not find an association between the M132L SNP and susceptibility to CWD. Thus, there is controversy about the effect of M132L SNP on susceptibility to CWD. In the present study, we investigated novel risk factors for CWD in elk. We investigated genetic polymorphisms of the PRNP gene by amplicon sequencing and compared genotype, allele, and haplotype frequencies between CWD-positive and CWD-negative elk. In addition, we performed a linkage disequilibrium (LD) analysis by the Haploview version 4.2 program. Furthermore, we evaluated the 3D structure and electrostatic potential of elk prion protein (PrP) according to the S100G SNP using AlphaFold and the Swiss-PdbViewer 4.1 program. Finally, we analyzed the free energy change of elk PrP according to the S100G SNP using I-mutant 3.0 and CUPSAT. We identified 23 novel SNP of the elk PRNP gene in 248 elk. We found a strong association between PRNP SNP and susceptibility to CWD in elk. Among those SNP, S100G is the only non-synonymous SNP. We identified that S100G is predicted to change the electrostatic potential and free energy of elk PrP. To the best of our knowledge, this was the first report of a novel risk factor, the S100G SNP, for CWD.

Antler-derived microRNA PC-5p-1090 inhibits HCC cell proliferation, migration, and invasion by targeting MARCKS, SMARCAD1, and SOX9.

The capability of microRNAs (miRNAs) to regulate gene expression across species has opened new avenues for miRNA-based therapeutics. Here, we investigated the potential of PC-5p-1090 (miR-PC-1090), a miRNA found in deer antlers, to control the malignant phenotypes of hepatocellular carcinoma (HCC) cells. Using Cell Counting Kit-8 and transwell assays, we found that heterologous expression of miR-PC-1090 inhibited HCC cell proliferation, migration, and invasion. Bioinformatics analysis indicated that predicted miR-PC-1090 targets, including MARCKS, SMARCAD1, and SOX9, were significantly elevated in HCC tissues, and their high expressions were associated with poor overall survival of HCC patients. Moreover, mechanistic investigations revealed that miR-PC-1090 promoted the degradation of MARCKS and SMARCAD1 mRNAs and hindered the translation of SOX9 mRNA by recognizing their 3' untranslated regions. Subsequent loss-of-function and rescue experiments confirmed the involvement of MARCKS, SMARCAD1, and SOX9 in miR-PC-1090-suppressed HCC cell proliferation, migration, and invasion. Notably, MARCKS knockdown induced the downregulation of phosphorylated MARCKS and a corresponding upregulation of phosphorylated AKT in HCC. Conversely, miR-PC-1090 repressed MARCKS phosphorylation and effectively circumvented the activation of the PI3K/AKT pathway. Furthermore, miR-PC-1090 regulates the Wnt/ß-catenin pathway through SMARCAD1- and SOX9-mediated reduction of ß-catenin expression. Overall, our results illustrate the tumor-suppressive activity and molecular mechanism of antler-derived miR-PC-1090 in HCC cells, indicating its potential as a multiple-target agent for HCC treatment.

Galectin-1 promotes angiogenesis and chondrogenesis during antler regeneration.

BACKGROUND: Deer antlers are the only known mammalian structure that undergoes full regeneration. In addition, it is peculiar because when growing, it contains vascularized cartilage. The differentiation of antler stem cells (ASCs) into chondrocytes while inducing endochondral extension of blood vessels is necessary to form antler vascularized cartilage. Therefore, antlers provide an unparalleled opportunity to investigate chondrogenesis, angiogenesis, and regenerative medicine. A study found that Galectin-1 (GAL-1), which can be used as a marker in some tumors, is highly expressed in ASCs. This intrigued us to investigate what role GAL-1 could play in antler regeneration. METHODS: We measured the expression level of GAL-1 in antler tissues and cells by immunohistochemistry, WB and QPCR. We constructed antlerogenic periosteal cells (APCs, one cell type of ASCs) with the GAL-1 gene knocked out (APCGAL-1-/-) using CRISPR-CAS9 gene editing system. The effect of GAL-1 on angiogenesis was determined by stimulating human umbilical vein endothelial cells (HUVECs) using APCGAL-1-/- conditioned medium or adding exogenous deer GAL-1 protein. The effect of APCGAL-1-/- on chondrogenic differentiation was evaluated compared with the APCs under micro-mass culture. The gene expression pattern of APCGAL-1-/- was analyzed by transcriptome sequencing. RESULTS: Immunohistochemistry revealed that GAL-1 was widely expressed in the antlerogenic periosteum (AP), pedicle periosteum (PP) and antler growth center. Western blot and qRT-PCR analysis using deer cell lines further supports this result. The proliferation, migration, and tube formation assays of human umbilical vein endothelial cells (HUVECs) showed that the proangiogenic activity of APCGAL-1-/- medium was significantly decreased (P < 0.05) compared with the APCs medium. The proangiogenic activity of deer GAL-1 protein was further confirmed by adding exogenous deer GAL-1 protein (P < 0.05). The chondrogenic differentiation ability of APCGAL-1-/- was impeded under micro-mass culture. The terms of GO and KEGG enrichment of the differentially expressed genes (DEGs) of APCGAL-1-/- showed that down-regulated expression of pathways associated with deer antler angiogenesis, osteogenesis and stem cell pluripotency, such as the PI3K-AKT signaling pathway, signaling pathways regulating pluripotency of stem cells and TGF-ß signaling pathway. CONCLUSIONS: Deer GAL-1, has strong angiogenic activity, is widely and highly expressed in deer antler. The APCs can induce angiogenesis by secreting GAL-1. The knockout of GAL-1 gene of APCs damaged its ability to induce angiogenesis and differentiate into chondrocytes. This ability is crucial to the formation of deer antler vascularized cartilage. Moreover, Deer antlers offer a unique model to explore explore how angiogenesis at high levels of GAL-1 expression can be elegantly regulated without becoming cancerous.

Single-nucleus RNA and ATAC sequencing uncovers the molecular and cellular characteristics in the musk gland of Chinese forest musk deer (Moschus berezovskii).

The Chinese forest musk deer (FMD; Moschus berezovskii) is an endangered artiodactyl mammal. Musk secreted by the musk gland of male has extremely high economic and medicinal value. However, the molecular and cellular characteristics of the musk gland have not been studied. Here, we investigated the diversity and transcriptional composition of musk gland cell types and the effect of cell type-specific chromatin accessibility on gene expression using single-nucleus RNA sequencing (snRNA-seq) and single-nucleus ATAC sequencing (snATAC-seq) association analysis. Based on uniform manifold approximation and projection (UMAP) analysis, we identified 13 cell types from the musk gland, which included two different acinar cells (cluster 0 and cluster 10). Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that many pathways related to musk secretion were enriched in acinar cells. Our analysis also revealed acinar cell core transcription factors and core target genes, and further constructed acinar cell-specific regulatory networks. In cluster 0, 11 core target genes (Nedd4l, Adcy9, Akr1c1, Vapb, Me1, Acsl1, Acss3, Srd5a1, Scnn1a, Acadm, and Nceh1) possibly related to musk secretion were regulated by 24 core transcription factors (SP3, NFIC, NR6A1, EHF, RUNX1, TFAP2A, RREB1, GRHL2, NFIB, ELF1, MAX, KLF5, REL, HES1, POU2F3, TFDP1, NR2C1, ATF7, MEIS1, NR4A2, NFIA, PBX1, ZNF652, and NFKB1). In cluster 10, four core target genes (Akr1c1, Pcca, Atp1b1, and Sgk1) possibly related to musk secretion were regulated by 10 core transcription factors (BARX2, EHF, PBX1, RUNX1, NFIB, FOXP1, KLF3, KLF6, ETV6, and NR3C2). Moreover, the credibility of snRNA-seq and snATAC-seq data was verified by fluorescence in situ hybridization and immunohistochemistry. Finally, cell communication analysis demonstrated that the two types of acinar cells mainly have communications in musk secretion-related processes. In conclusion, we provided important insights and invaluable resources for the molecular and cellular characteristics of the musk gland, which will lay a foundation for the study of musk secretion mechanism in the future.

MicroRNA let-7f-5p regulates PI3K/AKT/COX2 signaling pathway in bacteria-induced pulmonary fibrosis via targeting of PIK3CA in forest musk deer.

Background: Recent studies have characterized that microRNA (miRNA) is a suitable candidate for the study of bleomycin/LPS-induced pulmonary fibrosis, but the knowledge on miRNA in bacteria-induced pulmonary fibrosis (BIPF) is limited. Forest musk deer (Moschus berezovskii, FMD) is an important endangered species that has been seriously affected by BIPF. We sought to determine whether miRNA exist that modulates the pathogenesis of BIPF in FMD. Methods: High-throughput sequencing and RT-qPCR were used to determine the differentially expressed miRNAs (DEmiRNAs) in the blood of BIPF FMD. The DEmiRNAs were further detected in the blood and lung of BIPF model rat by RT-qPCR, and the targeting relationship between candidate miRNA and its potential target gene was verified by dual-luciferase reporter activity assay. Furthermore, the function of the candidate miRNA was verified in the FMD lung fibroblast cells (FMD-C1). Results: Here we found that five dead FMD were suffered from BIPF, and six circulating miRNAs (miR-30g, let-7f-5p, miR-27-3p, miR-25-3p, miR-9-5p and miR-652) were differentially expressed in the blood of the BIPF FMD. Of these, let-7f-5p showed reproducibly lower level in the blood and lung of the BIPF model rat, and the expression levels of PI3K/AKT/COX2 signaling pathway genes (PIK3CA, PDK1, Akt1, IKBKA, NF-κB1 and COX2) were increased in the lung of BIPF model rats, suggesting that there is a potential correlation between BIPF and the PI3K/AKT/COX2 signaling pathway. Notably, using bioinformatic prediction and experimental verification, we demonstrated that let-7f-5p is conserved across mammals, and the seed sequence of let-7f-5p displays perfect complementarity with the 3' UTR of PIK3CA gene and the expression of the PIK3CA gene was regulated by let-7f-5p. In order to determine the regulatory relationship between let-7f-5p and the PI3K/AKT/COX2 signaling pathway in FMD, we successfully cultured FMD-C1, and found that let-7f-5p could act as a negative regulator for the PI3K/Akt/COX2 signaling pathway in FMD-C1. Collectively, this study not only provided a study strategy for non-invasive research in pulmonary disease in rare animals, but also laid a foundation for further research in BIPF.

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.

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 (HBBA), 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 HBBA gene of sika deer, red deer, and North American wapiti (Cervus canadensis) in China may have undergone severe genetic drifts.

First report of a strong association between genetic polymorphisms of the prion protein gene (PRNP) and susceptibility to chronic wasting disease in sika deer (Cervus nippon).

Prion diseases are incurable neurodegenerative disorders caused by proteinase K-resistant prion protein (PrPSc ) derived from normal prion protein (PrPC ) encoded by the prion protein gene (PRNP). Although the cervid PRNP gene plays a pivotal role in the pathological mechanism of chronic wasting disease (CWD), there is no existing association analysis between susceptibility to CWD and genetic polymorphisms of the PRNP gene in sika deer. We investigated genetic polymorphisms of the PRNP gene using amplicon sequencing in sika deer. In addition, to identify a genetic susceptibility factor, we compared the genotype, allele and haplotype frequencies of the PRNP gene between CWD-positive and CWD-negative sika deer. Furthermore, to assess the effect of the genetic polymorphisms on sika deer prion protein (PrP), we performed in silico analysis using PolyPhen-2, PROVEAN and AMYCO. Finally, we analysed the tertiary structure and electrostatic potential of sika deer PrP based on single nucleotide polymorphisms (SNPs) using the SWISS-MODEL and Swiss-PdbViewer programs. We found a total of 24 SNPs of the PRNP gene, including 22 novel SNPs (10 synonymous SNPs and 12 nonsynonymous SNPs), in sika deer. Among the nonsynonymous SNPs, we found a strong association of susceptibility to CWD with c.56G > A (Ser19Asn). In addition, we found that c.56G > A (Ser19Asn), c.296A > T (His99Leu) and c.560T > A (Val187Asp) were predicted to have damaging effects on sika deer PrP. Furthermore, we observed significant alterations in the electrostatic potential of sika deer PrP by genetic polymorphisms of the 187Asp allele. To the best of our knowledge, this was the first association study between genetic polymorphisms of the PRNP gene and susceptibility to CWD in sika deer.

Genomic architecture of phenotypic extremes in a wild cervid.

Identifying the genes underlying fitness-related traits such as body size and male ornamentation can provide tools for conservation and management and are often subject to various selective pressures. Here we performed high-depth whole genome re-sequencing of pools of individuals representing the phenotypic extremes for antler and body size in white-tailed deer (Odocoileus virginianus). Samples were selected from a tissue repository containing phenotypic data for 4,466 male white-tailed deer from Anticosti Island, Quebec, with four pools representing the extreme phenotypes for antler and body size after controlling for age. Our results revealed a largely homogenous population but detected highly divergent windows between pools for both traits, with the mean allele frequency difference of 14% for and 13% for antler and body SNPs in outlier windows, respectively. Genes in outlier antler windows were enriched for pathways associated with cell death and protein metabolism and some of the most differentiated windows included genes associated with oncogenic pathways and reproduction, processes consistent with antler evolution and growth. Genes associated with body size were more nuanced, suggestive of a highly complex trait. Overall, this study revealed the complex genomic make-up of both antler morphology and body size in free-ranging white-tailed deer and identified target loci for additional analyses.

Reference gene screening for analyzing gene expression in the heart, liver, spleen, lung and kidney of forest musk deer.

The screening of reference genes for real-time quantitative PCR (qPCR) in forest musk deer (FMD) tissue is of great significance to the basic research on FMD. However, there are few reports on the stability analysis of FMD reference genes so far. In this study, We used qPCR to detect the expression levels of 11 reference gene candidates (18S rRNA, beta-actin [ACTB], glyceraldehyde-3-phosphate dehydrogenase [GAPDH], TATA box-binding protein [TBP], hypoxanthine phosphoribosyltransferase 1 [HPRT1], tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta polypeptide [YWHAZ], hydroxymethylbilane synthase [HMBS], eukaryotic translation elongation factor 1 alpha 1 [EEF1A1], succinate dehydrogenase complex flavoprotein subunit A [SDHA], peptidylprolyl isomerase B [PPIB], and ubiquitin C [UBC]) in heart, liver, spleen, lung and kidney of FMD. After removing 18S rRNA on account of its high expression level, geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to evaluate the expression stability of the remaining genes in the five organs, and further comprehensive ranking was calculated by RefFinder. According to the results, the selected reference genes with the most stable expression in the heart of FMD are SDHA and YWHAZ, while in the liver are ACTB and SDHA; in the spleen and lung are YWHAZ and HPRT1; in the kidney are YWHAZ and PPIB. The use of common reference genes in all five organs is not recommended. The analyses showed that tissue is an important variability factor in genes expression stability. Meanwhile, the result can be used as a reference for the selection of reference genes for qPCR in further study.

Identification of management units in threatened populations of Arnica mallotopus Makino (Asteraceae) using novel EST-SSR markers.

Arnica mallotopus is a perennial herb endemic to the snowy regions of Japan. At the southern edge of its distribution, in Kyoto Prefecture, overgrazing by sika deer and decreased snowfall have resulted in the rapid decline of A. mallotopus populations. Therefore, there is an urgent need for a conservation genetic analysis of the remaining local populations. In this study, we first developed 13 EST-SSR markers to evaluate genetic variation in A. mallotopus. The average number of alleles per locus was 5.33. Genetic analysis using these markers showed that the investigated samples were classified into two groups corresponding to landscape structure. One group isolated from a tributary of the Yura River showed a strong population bottleneck signal, likely resulting from founder effects and subsequent drifts. On the other hand, the genetic diversity of the second group in the main distribution along the Yura River was higher and less inbred. Overall, our assessment suggested recognizing the two genetic groups as management units in conservation programs for the threatened populations.

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.

Novel insights into Dhh signaling in antler chondrocyte proliferation and differentiation: Involvement of Foxa.

The desert hedgehog (Dhh) is crucial for spermatogenesis and Leydig cell differentiation, but little is known regarding its physiological function in cartilage. In this study, Dhh mRNA was abundant in antler chondrocytes, where it advanced cell proliferation concomitant with accelerated transition from the G1 to the S phase and induced elevation of the hypertrophic chondrocyte markers, Col X and Runx2. Silencing of Ptch1 resulted in appreciable Smo accumulation and enhanced rDhh stimulation of Smo, whose impediment by cyclopamine obscured the proliferative function of Dhh and alleviated its guidance of chondrocyte differentiation. Further analysis evidenced the noteworthy positive action of Smo in the bridging between Dhh and Gli transcription factors. Obstruction of Gli1 by GANT58 caused the failed stimulation of Col X and Runx2 by rDhh. Analogously, siRNA against Gli1-3 hindered chondrocyte differentiation in the context of rDhh. Simultaneously, Gli transcription factors mediated the regulation of Dhh on Foxa1, Foxa2, and Foxa3, whose knockdown impaired chondrocyte differentiation. Attenuation of Foxa antagonized the augmentation of Col X and Runx2 generated by rDhh. Collectively, Dhh signaling through its target Foxa appears to induce antler chondrocyte proliferation and differentiation.

Insight into the epigenetic landscape of a currently endogenizing gammaretrovirus in mule deer (Odocoileus hemionus).

Endogenous retroviruses (ERVs) constitute a significant part of vertebrate genomes. They originated from past retroviral infections and some of them retain transcriptional activity. The key mechanism avoiding uncontrolled ERV transcription is DNA methylation-mediated epigenetic silencing. Despite numerous studies describing the involvement of ERV activity in cellular processes, epigenetic regulation of ERVs is still poorly understood. We previously described a cervid endogenous retrovirus (CrERV) in the mule deer genome. This virus exhibits massive insertional polymorphism, suggesting recent activity. Here we employed NGS-based strategy to determine the methylation pattern of CrERV integrations in four mule deer. Besides the vast majority of methylated integrations, we identified a tiny fraction of demethylated proviral copies. These copies represent evolutionary older integrations located near gene promoters. In general, our work is a first attempt to characterize the epigenetic landscape of insertionally polymorphic ERV on a whole-genome scale and offers insight into its interactions with a host.