Illumina-based de novo transcriptome sequencing and analysis of Chinese forest musk deer

Analysis and comparison of blood metabolome of forest musk deer in musk secretion and non-secretion periods

Musk is an important animal product, but the musk secretion mechanism of forest musk deer (Moschus berezovskii) is still unclear. The musk synthesis process in forest musk deer is extremely complex, and many raw materials are directly or indirectly derived from forest musk deer blood. In this study, metabolomics was used to analyze the blood of forest musk deer in secretory and non-secretory phases for the first time, aim at explaining the secretion mechanism from the perspective of blood metabolism. We found that P450-related, choline-related, axonal regeneration and other pathways and related metabolites were significantly enriched during the musk secretion of forest musk deer. These pathways and metabolites related to P450 and choline in blood may have important implications for the mechanism of musk secretion in forest musk deer, because blood components were closely related to musk components and could provide raw materials for musk synthesis in musk gland cells.

Single-nucleus transcriptomics and chromatin accessibility analysis of musk gland development in Chinese forest musk deer (Moschus berezovskii)

Musk secreted by male forest musk deer (Moschus berezovskii) musk glands is an invaluable component of medicine and perfume. Musk secretion depends on musk gland maturation; however, the mechanism of its development remains elusive. Herein, using single cell multiome ATAC + gene expression coupled with several bioinformatic analyses, a dynamic transcriptional cell atlas of musk gland development was revealed, and key genes and transcription factors affecting its development were determined. Twelve cell types, including two different types of acinar cells (Clusters 0 and 10) were identified. Single-nucleus RNA and single-nucleus ATAC sequencing analyses revealed that seven core target genes associated with musk secretion (Hsd17b2, Acacb, Lss, Vapa, Aldh16a1, Aldh7a1, and Sqle) were regulated by 12 core transcription factors (FOXO1, CUX2, RORA, RUNX1, KLF6, MGA, NFIC, FOXO3, ETV5, NR3C1, HSF4, and MITF) during the development of Cluster 0 acinar cells. Kyoto Encyclopedia of Genes and Genomes enrichment showed significant changes in the pathways associated with musk secretion during acinar cell development. Gene set variation analysis also revealed that certain pathways associated with musk secretion were enriched in 6-year-old acinar cells. A gene co-expression network was constructed during acinar cell development to provide a precise understanding of the connections between transcription factors, genes, and pathways. Finally, intercellular communication analysis showed that intercellular communication is involved in musk gland development. This study provides crucial insights into the changes and key factors underlying musk gland development, which serve as valuable resources for studying musk secretion mechanisms and promoting the protection of this endangered species.

Integrated multi-omics analysis reveals insights into Chinese forest musk deer (Moschus berezovskii) genome evolution and musk synthesis

Among the artiodactyls, male animals belonging to the Family Moschidae have a unique tissue, the musk gland, with the capability of musk synthesis. However, the genetic basis of musk gland formation and musk production are still poorly understood. Here, musk gland tissues from two juvenile and three adult Chinese forest musk deer (Moschus berezovskii) were utilized to analyze genomic evolution events, evaluate mRNA profiles and investigate cell compositions. By performing genome reannotation and comparison with 11 ruminant genomes, three expanded gene families were identified in the Moschus berezovskii genome. Transcriptional analysis further indicated that the musk gland displayed a prostate-like mRNA expression pattern. Single-cell sequencing revealed that the musk gland is composed of seven distinguishable cell types. Among them, sebaceous gland cells and luminal epithelial cells play important roles in musk synthesis, while endothelial cells master the regulation of cell-to-cell communication. In conclusion, our study provides insights into musk gland formation and the musk-synthesizing process.

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.

ITRAQ-based quantitative proteomics analysis of forest musk deer with pneumonia

Pneumonia can seriously threaten the life of forest musk deer (FMD, an endangered species). To gain a comprehensive understanding of pneumonia pathogenesis in FMD, iTRAQ-based proteomics analysis was performed in diseased (Pne group) lung tissues of FMD that died of pneumonia and normal lung tissues (Ctrl group) of FMD that died from fighting against each other. Results showed that 355 proteins were differentially expressed (fold change ≥ 1.2 and adjusted P-value < 0.05) in Pne vs. Ctrl. GO/KEGG annotation and enrichment analyses showed that dysregulated proteins might play vital roles in bacterial infection and immunity. Given the close association between bacterial infection and pneumonia, 32 dysregulated proteins related to Staphylococcus aureus infection, bacterial invasion of epithelial cells, and pathogenic Escherichia coli infection were screened out. Among these 32 proteins, 13 proteins were mapped to the bovine genome. Given the close phylogenetic relationships of FMD and bovine, the protein-protein interaction networks of the above-mentioned 13 proteins were constructed by the String database. Based on the node degree analysis, 5 potential key proteins related to pneumonia-related bacterial infection in FMD were filtered out. Moreover, 85 dysregulated proteins related to the immune system process were identified given the tight connection between immune dysregulation and pneumonia pathogenesis. Additionally, 12 proteins that might function as crucial players in pneumonia-related immune response in FMD were screened out using the same experimental strategies described above. In conclusion, some vital proteins, biological processes, and pathways in pneumonia development were identified in FMD.

Characteristics of steroidogenesis-related factors in the musk gland of Chinese forest musk deer (Moschus berezovskii)

Musk secreted by Chinese forest musk deer (FMD; Moschus berezovskii) is a highly valuable ingredient in the fields of perfumery and medicine, and the main factor affecting the production of musk is the androgen level of male FMD. To clarify whether the musk gland of FMD can synthesize androgen, we compared and analyzed the expression patterns of steroid hormone biosynthesis-related genes in the musk gland and testis of FMD by RNA-seq and RT-qPCR. We obtained 33,308 and 38,602 unigenes from the musk gland and testis, respectively, and 26,780 co-expressed unigenes. Analysis of co-expressed genes revealed that 12,647 genes were annotated to 11,484 Gene Ontology terms and 10,941 genes were annotated to 6120 pathways, including several pathways important in metabolic and synthetic activity. Next, 21 steroid hormone synthesis-related genes were screened from the transcriptome of the musk gland of 4-month-old FMD. The expression levels of three key genes of steroid hormone biosynthesis (CYP11A1, CYP17A1, and HSD3B) in the musk gland differed from their expression levels in the testis based on RT-qPCR. Furthermore, immunohistochemistry indicated that CYP11A1, CYP17A1, and HSD3B were localized in the glandular tubular columnar cells of the musk gland. These results suggested that the musk gland of male FMD has the potential to locally synthesize steroid hormone and thus plays a critically important role in musk secretion.

The development of genome-wide single nucleotide polymorphisms in blue wildebeest using the DArTseq platform

Blue wildebeest (Connochaetes taurinus taurinus) are economically important antelope that are widely utilised in the South African wildlife industry. However, very few genomic resources are available for blue wildebeest that can assist in breeding management and facilitate research. This study aimed to develop a set of genome-wide single nucleotide polymorphism (SNP) markers for blue wildebeest. The DArTseq genotyping platform, commonly used in polyploid plant species, was selected for SNP discovery. A limited number of published articles have described the use of the DArTseq platform in animals and, therefore, this study also provided a unique opportunity to assess the performance of the DArTseq platform in an animal species. A total of 20,563 SNPs, each located within a 69 bp sequence, were generated. The developed SNP markers had a high average scoring reproducibility (>99%) and a low percentage missing data (~9.21%) compared to other reduced representation sequencing approaches that have been used in animal studies. Furthermore, the number of candidate SNPs per nucleotide position decreased towards the 3' end of sequence reads, and the ratio of transitions (Ts) to transversions (Tv) remained similar for each read position. These observations indicate that there was no read position bias, such as the identification of false SNPs due to low sequencing quality, towards the tail-end of sequencing reads. The DArTseq platform was also successful in identifying a large number of informative SNPs with desirable polymorphism parameters such as a high minor allele frequency (MAF). The Bos taurus genome was used for the in silico mapping of the marker sequences and a total of 6020 (29.28%) sequences were successfully mapped against the bovine genome. The marker sequences mapped to all of the bovine chromosomes establishing the genome-wide distribution of the SNPs. Moreover, the high observed Ts:Tv ratio (2.84:1) indicate that the DArTseq platform targeted gene-rich regions of the blue wildebeest genome. Finally, functional annotation of the marker sequences revealed a wide range of different putative functions indicating that these SNP markers can be useful in functional gene studies. The DArTseq platform, therefore, represents a high-throughput, robust and cost-effective genotyping platform, which may find adoption in several other African antelope and animal species.

Whole-genome sequencing of wild Siberian musk deer (Moschus moschiferus) provides insights into its genetic features

Background

Siberian musk deer, one of the seven species, is distributed in coniferous forests of Asia. Worldwide, the population size of Siberian musk deer is threatened by severe illegal poaching for commercially valuable musk and meat, habitat losses, and forest fire. At present, this species is categorized as Vulnerable on the IUCN Red List. However, the genetic information of Siberian musk deer is largely unexplored.

Results

Here, we produced 3.10 Gb draft assembly of wild Siberian musk deer with a contig N50 of 29,145 bp and a scaffold N50 of 7,955,248 bp. We annotated 19,363 protein-coding genes and estimated 44.44% of the genome to be repetitive. Our phylogenetic analysis reveals that wild Siberian musk deer is closer to Bovidae than to Cervidae. Comparative analyses showed that the genetic features of Siberian musk deer adapted in cold and high-altitude environments. We sequenced two additional genomes of Siberian musk deer constructed demographic history indicated that changes in effective population size corresponded with recent glacial epochs. Finally, we identified several candidate genes that may play a role in the musk secretion based on transcriptome analysis.

Conclusions

Here, we present a high-quality draft genome of wild Siberian musk deer, which will provide a valuable genetic resource for further investigations of this economically important musk deer.

Comparative Genomics Reveals the Genetic Mechanisms of Musk Secretion and Adaptive Immunity in Chinese Forest Musk Deer

The Chinese forest musk deer (Moschus berezovskii; FMD) is an artiodactyl mammal and is both economically valuable and highly endangered. To investigate the genetic mechanisms of musk secretion and adaptive immunity in FMD, we compared its genome to nine other artiodactyl genomes. Comparative genomics demonstrated that eight positively selected genes (PSGs) in FMD were annotated in three KEGG pathways that were related to metabolic and synthetic activity of musk, similar to previous transcriptome studies. Functional enrichment analysis indicated that many PSGs were involved in the regulation of immune system processes, implying important reorganization of the immune system in FMD. FMD-specific missense mutations were found in two PSGs (MHC class II antigen DRA and ADA) that were classified as deleterious by PolyPhen-2, possibly contributing to immune adaptation to infectious diseases. Functional assessment showed that the FMD-specific mutation enhanced the ADA activity, which was likely to strengthen the immune defense against pathogenic invasion. Single nucleotide polymorphism-based inference showed the recent demographic trajectory for FMD. Our data and findings provide valuable genomic resources not only for studying the genetic mechanisms of musk secretion and adaptive immunity, but also for facilitating more effective management of the captive breeding programs for this endangered species.

microRNA and Other Small RNA Sequence Profiling across Six Tissues of Chinese Forest Musk Deer (Moschus berezovskii)

The Chinese forest musk deer (Moschus berezovskii) is an economically important species distributed throughout southwest China and northern Vietnam. Occurrence and development of disease are aggravated by inbreeding and genetic diversity declines in captive musk deer populations. Deep transcriptomics investigation may provide a promising way to improve genetic health of captive and wild FMD population. MicroRNAs (miRNAs), which regulate gene expression by targeting and suppressing of mRNAs, play an important role in physiology and organism development control. In this study, RNA-seq technology was adopted to characterize the miRNA transcriptome signature among six tissues (heart, liver, spleen, lung, kidney, and muscle) in Chinese forest musk deer at two years of age. Deep sequencing generated a total of 103,261,451 (~87.87%) good quality small RNA reads; of them 6,622,520 were unique across all six tissues. A total of 2890 miRNAs were identified, among them 1129 were found to be expressed in all tissues. Moreover, coexpression of 20 miRNAs (>2000RPM) in all six tissues and top five highly expressed miRNAs in each tissue implied the crucial and particular function of them in FMD physiological processes. Our findings of forest musk deer miRNAs supplement the database of transcriptome information for this species and conduce to our understanding of forest musk deer biology.

The draft genome sequence of forest musk deer (Moschus berezovskii)

Background

The forest musk deer, Moschus berezovskii, is one of seven musk deer (Moschus spp.) and is distributed in Southwest China. Akin to other musk deer, the forest musk deer has been traditionally and is currently hunted for its musk (i.e., global perfume industry). Considerable hunting pressure and habitat loss have caused significant population declines. Consequently, the Chinese government commenced captive breeding programs for musk harvesting in the 1950s. However, the prevalence of fatal diseases is considerably restricting population increases. Disease severity and extent are exacerbated by inbreeding and genetic diversity declines in captive musk deer populations. It is essential that knowledge of captive and wild forest musk deer populations' immune system and genome be gained in order to improve their physical and genetic health. We have thus sequenced the whole genome of the forest musk deer, completed the genomic assembly and annotation, and performed preliminary bioinformatic analyses.

Findings

A total of 407 Gb raw reads from whole-genome sequencing were generated using the Illumina HiSeq 4000 platform. The final genome assembly is around 2.72 Gb, with a contig N50 length of 22.6 kb and a scaffold N50 length of 2.85 Mb. We identified 24,352 genes and found that 42.05% of the genome is composed of repetitive elements. We also detected 1,236 olfactory receptor genes. The genome-wide phylogenetic tree indicated that the forest musk deer was within the order Artiodactyla, and it appeared as the sister clade of four members of Bovidae. In total, 576 genes were under positive selection in the forest musk deer lineage.

Conclusions

We provide the first genome sequence and gene annotation for the forest musk deer. The availability of these resources will be very useful for the conservation and captive breeding of this endangered and economically important species and for reconstructing the evolutionary history of the order Artiodactyla.