Differential expression profiles of microRNAs in musk gland of unmated and mated forest musk deer (Moschus berezovskii)

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.