[Research progress on molecular genetics of forest musk deer].

Forest musk deer is one of the large-scale farming musk deer animals with the largest population at the same time. The male musk deer can secrete valuable medicines, which has high medicinal and economic value. Due to the loss of habitat and indiscriminate hunting, the numbers of wild population specie and the distribution have been drastically reduced. Therefore, in-depth understanding of the molecular genetics progress of forest musk deer will pave a way for musk deer protection and breeding. In this review, the progress associated with the molecular marker, genetic classification, artificial breeding, musk secretion and disease in past decades were reviewed, in order to provide a theoretical basis for subsequent molecular genetic researches in forest musk deer.

Identification of Potential miRNA-mRNA Regulatory Network Associated with Growth and Development of Hair Follicles in Forest Musk Deer.

In this study, sRNA libraries and mRNA libraries of HFs of FMD were constructed and sequenced using an Illumina HiSeq 2500, and the expression profiles of miRNAs and genes in the HFs of FMD were obtained at the anagen and catagen stages. In total, 565 differentially expressed unigenes (DEGs) were identified, 90 of which were upregulated and 475 of which were downregulated. In the BP category of GO enrichment, the DEGs were enriched in the processes related to HF development and differentiation, including the hair cycle regulation and processes, HF development, skin epidermis development, regulation of HF development, skin development, the Wnt signaling pathway, and the BMP signaling pathway. Through KEGG analysis it was found that DEGs were significantly enriched in pathways associated with HF development and growth. A total of 186 differentially expressed miRNAs (DEmiRNAs) were screened (p < 0.05) in the HFs of FMD at the anagen stage vs. the catagen stage, 33 of which were upregulated and 153 of which were downregulated. Through DEmiRNA-mRNA association analysis, we found DEmiRNAs and target genes that mainly play regulatory roles in HF development and growth. The enrichment analysis of DEmiRNA target genes revealed similarities with the enrichment results of DEGs associated with HF development. Notably, both sets of genes were enriched in key pathways such as the Notch signaling pathway, melanogenesis, the cAMP signaling pathway, and cGMP-PKG. To validate our findings, we selected 11 DEGs and 11 DEmiRNAs for experimental verification using RT-qPCR. The results of the experimental validation were consistent with the RNA-Seq results.

Molecular cloning and characterization of caprine calpastatin gene.

Calpastatin (CAST) is an important gene for meat quality traits in livestock and poultry. The cDNA of caprine CAST gene was amplified for the first time using RACE-PCR. Results showed the full-length cDNA of caprine CAST gene (Accession no. GU944861) was 2435 base pair (bp) and contained a 2187 bp open reading frame encoding a protein with 728 amino acid residues. Bioinformatic analysis indicated that caprine CAST cDNA was 89.8-95.4, 83.5-92.2, 72.8-81.8 and 69.8-73.5% identical to sheep, cattle, pig and human CAST cDNA. It was predicted that caprine CAST contained four conserved domains with 42 serine phosphorylation loci, 18 threonine phosphorylation loci, 1 tyrosine phosphorylation locus and 5 specific PKC phosphorylation loci. This work provided an important experimental basis for further research on the function of CAST in goat.

[Cloning and tissue expression of CAST transcript â ¡ in goat].

Calpastatin (CAST) gene is closely related with meat quality in livestock and poultry. Based on the bovine and ovine mRNA sequences, the cDNA of CAST Ⅱ gene in goat was amplified successfully for the first time by using RACE-PCR. Results showed that CAST Ⅱ of goat was 2474 bp in length with an open reading frame (ORF) 1695 bp long and encoded 564 amino acids, and there were four conserved domains and one conserved seven-peptide domain in amino acids sequences. Bioinformation analysis indicated that its secondary structures mainly were random coil and helical regions, and contained rich hydrophobic regions, certain phosphorylation sites, and protein kinase C (PKC) sites. Meanwhile, analysis of tissue expression of the gene in Tianfu meat goat demonstrated it was expressed in seven selected tissues. When the goat was of 6-month age, the highest expression was observed in longissimusdorsi, which was significantly higher than that of crureus (P<0.05) and other internal organ tissues (P<0.01).Furthermore, the expression of CAST II increased with the rise of the age and became the highest when the goat was at three-year age.

Distribution patterns of microsatellites and development of its marker in different genomic regions of forest musk deer genome based on high throughput sequencing.

Forest musk deer (Moschus berezovskii, FMD) is an endangered artiodactyl species, male FMD produce musk. We have sequenced the whole genome of FMD, completed the genomic assembly and annotation, and performed bioinformatic analyses. Our results showed that microsatellites (SSRs) displayed nonrandomly distribution in genomic regions, and SSR abundances were much higher in the intronic and intergenic regions compared to other genomic regions. Tri- and hexanucleotide perfect (P) SSRs predominated in coding regions (CDSs), whereas, tetra- and pentanucleotide P-SSRs were less abundant. Trifold P-SSRs had more GC-contents in the 5'-untranslated regions (5'UTRs) and CDSs than other genomic regions, whereas mononucleotide P-SSRs had the least GC-contents. The repeat copy numbers (RCN) of the same mono- to hexanucleotide P-SSRs had different distributions in different genomic regions. The RCN of trinucleotide P-SSRs had increased significantly in the CDSs compared to the transposable elements (TEs), intronic and intergenic regions. The analysis of coefficient of variability (CV) of P-SSRs showed that the RCN of mononucleotide P-SSRs had relative higher variation in different genomic regions, followed by the CV pattern of RCN: dinucleotide P-SSRs > trinucleotide P-SSRs > tetranucleotide P-SSRs > pentanucleotide P-SSRs > hexanucleotide P-SSRs. The CV variations of RCN of the same mono- to hexanucleotide P-SSRs were relative higher in the intron and intergenic regions, followed by that in the TEs, and the relative lower was in the 5'UTR, CDSs and 3'UTRs. 58 novel polymorphic SSR loci were detected based on genotyping DNA from 36 captive FMD and 22 SSR markers finally showed polymorphism, stability, and repetition.