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

Abstract

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.

Keywords

acinar cells / Chinese forest musk deer / musk gland development / single cell multiome ATAC + gene expression

References

[1] BeyeaMM,HeslopCL, SawyezCG et al. (2007). Selective up-regulation of LXR-regulated genes ABCA1, ABCG1, and APOE in macrophages through increased endogenous synthesis of 24(S), 25-epoxycholesterol. Journal of Biological Chemistry 282,5207–5216.

[2] BrockerC,LassenN, EsteyT et al. (2010). Aldehyde dehydrogenase 7A1 (ALDH7A1) is a novel enzyme involved in cellular defense against hyperosmotic stress. Journal of Biological Chemistry 285,18452–18463.

[3] BuenrostroJD,WuBJ, LitzenburgerUM et al. (2015). Single-cell chromatin accessibility reveals principles of regulatory variation. Nature 523,486–490.

[4] BakkenTE,HodgeRD, MillerJA et al. (2018). Single-nucleus and single-cell transcriptomes compared in matched cortical cell types. PLoS ONE 13,e0209648.

[5] BuenrostroJD,CorcesMR, LareauCA et al. (2018). Integrated single-cell analysis maps the continuous regulatory landscape of human hematopoietic differentiation. Cell 173,1535–1548.

[6] ButlerA,HoffmanP, SmibertP,Papalexi E,SatijaR (2018). Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nature Biotechnology 36,411–420.

[7] CaoXH,ZhouYD (2007). Progress on anti-inflammatory effects of musk. China Pharmacy 18,1662–1665.

[8] ChristmannRB,Sampaio-Barros P,StifanoG et al. (2014). Association of interferon-and transforming growth factor β-regulated genes and macrophage activation with systemic sclerosis-related progressive lung fibrosis. Arthritis & Rheumatology 66,714–725.

[9] ChungNC,StoreyJD (2015). Statistical significance of variables driving systematic variation in high-dimensional data. Bioinformatics 31,545–554.

[10] CraigSEL,Brady-Kalnay SM (2015). Regulation of development and cancer by the R2B subfamily of RPTPs and the implications of proteolysis. Seminars in Cell & Developmental Biology 37,108–118.

[11] CharkoftakiG,ChenY, HanM et al. (2017). Transcriptomic analysis and plasma metabolomics in Aldh16a1-null mice reveals a potential role of ALDH16A1 in renal function. Chemico-Biological Interactions 276,15–22.

[12] CarterRA,Bihannic L,RosencranceC et al. (2018). A single-cell transcriptional atlas of the developing murine cerebellum. Current Biology 28,2910–2920.

[13] CantalapiedraCP,Hernández-Plaza A,LetunicI,BorkP,Huerta-Cepas J (2021). eggNOG-mapper v2: Functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Molecular Biology and Evolution 38,5825–5829.

[14] Castro-MondragonJA,Riudavets-Puig R,RauluseviciuteI et al. (2022). JASPAR 2022: The 9th release of the open-access database of transcription factor binding profiles. Nucleic Acids Research 50,D165–D173.

[15] DeerwesterS,DumaisST, FurnasGW,Landauer TK,HarshmanR (1990). Indexing by latent semantic analysis. Journal of the American Society for Information Science and Technology 41,391–407.

[16] DobinA,DavisCA, SchlesingerF et al. (2013). STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 29,15–21.

[17] FengQQ,LiuJT (2015). Progress on pharmacological activity of muscone. Food and Drug 17,212–214.

[18] GaudetP,Livstone MS,LewisSE,ThomasPD (2011). Phylogenetic-based propagation of functional annotations within the Gene Ontology consortium. Briefings in Bioinformatics 12,449–462.

[19] GalbraithMD,Andrysik Z,PandeyA et al. (2017). Cdk8 kinase activity promotes glycolysis. Cell Reports 21,1495–1506.

[20] GiordaniL,HeGJ, NegroniE et al. (2019). High-dimensional single-cell cartography reveals novel skeletal muscle-resident cell populations. Molecular Cell 74,609–621.

[21] Guerrero-JuarezCF,Dedhia PH,JinSQ et al. (2019). Single-cell analysis reveals fibroblast heterogeneity and myeloid-derived adipocyte progenitors in murine skin wounds. Nature Communications 10,650.

[22] HeikinheimoKA,LaineMA, RitvosOV,Voutilainen RJ,HoganBL,LeivoIV (1999). Bone morphogenetic protein-6 is a marker of serous acinar cell differentiation in normal and neoplastic human salivary gland. Cancer Research 59,5815–5821.

[23] HänzelmannS,Castelo R,GuinneyJ (2013). GSVA: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 14,7.

[24] HeL,WangWX, LiLH et al. (2014). Effects of crowding and sex on fecal cortisol levels of captive forest musk deer. Biological Research 47,48.

[25] HsiaLT,AshleyN, OuaretD,Wang LM,WildingJ,BodmerWF (2016). Myofibroblasts are distinguished from activated skin fibroblasts by the expression of AOC3 and other associated markers. PNAS 113,E2162–E2171.

[26] HuP,LiuJ, ZhaoJJ et al. (2018). Single-nucleus transcriptomic survey of cell diversity and functional maturation in postnatal mammalian hearts. Genes & Development 32,1344–1357.

[27] HuH,MiaoYR, JiaLH,Yu QY,ZhangQ,GuoAY (2019). AnimalTFDB 3.0: A comprehensive resource for annotation and prediction of animal transcription factors. Nucleic Acids Research 47,D33–D38.

[28] HeS,WangLH, LiuY et al. (2020). Single-cell transcriptome profiling of an adult human cell atlas of 15 major organs. Genome Biology 21,294.

[29] HuF,YangY (2020). Application of single cell sequencing technology in adipose tissue. Chinese Journal of Cell Biology 42,2282–2288.

[30] HaoYH,HaoS, Andersen-NissenE et al. (2021). Integrated analysis of multimodal single-cell data. Cell 184,3573–3587.

[31] InagakiK,OtsukaF, SuzukiJ et al. (2006). Involvement of bone morphogenetic protein-6 in differential regulation of aldosterone production by angiotensin II and potassium in human adrenocortical cells. Endocrinology 147,2681–2689.

[32] IndumathiS,Dhanasekaran M,RajkumarJS,SudarsanamD (2013). Exploring the stem cell and non-stem cell constituents of human breast milk. Cytotechnology 65,385–393.

[33] JensenM,HansenB, De MeytsP,Schäffer L,UrsøB (2007). Activation of the insulin receptor by insulin and a synthetic peptide leads to divergent metabolic and mitogenic signaling and responses. Journal of Biological Chemistry 282,35179–35186.

[34] JieH,XuZX, GaoJ et al. (2021). Differential expression profiles of microRNAs in musk gland of unmated and mated forest musk deer (Moschus berezovskii). PeerJ 9,e12710.

[35] JinSQ,Guerrero-Juarez CF,ZhangLH et al. (2021). Inference and analysis of cell-cell communication using CellChat. Nature Communications 12,1088.

[36] KaushikVK,KavanaM, VolzJM et al. (2009). Characterization of recombinant human acetyl-CoA carboxylase-2 steady-state kinetics. Biochimica et Biophysica Acta 1794,961–967.

[37] KorsunskyI,MillardN, FanJ et al. (2019). Fast, sensitive and accurate integration of single-cell data with Harmony. Nature Methods 16,1289–1296.

[38] LuML,HuangYW, LinSX (2002). Purification, reconstitution, and steady-state kinetics of the trans-membrane 17 beta-hydroxysteroid dehydrogenase 2. Journal of Biological Chemistry 277,22123–22130.

[39] LiuQ,Zerbinatti CV,ZhangJ et al. (2007). Amyloid precursor protein regulates brain apolipoprotein E and cholesterol metabolism through lipoprotein receptor LRP1. Neuron 56,66–78.

[40] LangfelderP,HorvathS (2008). WGCNA: An R package for weighted correlation network analysis. BMC Bioinformatics 9,559.

[41] LiuL,XieZM (2012). Expression and function of angiopoietin-1 during mouse kidney development. Chinese Journal of Biochemical Pharmaceutics 33,430–432.

[42] LiH (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv: Genomics.

[43] LiJ,Klughammer J,FarlikM et al. (2016). Single-cell transcriptomes reveal characteristic features of human pancreatic islet cell types. EMBO Reports 17,178–187.

[44] LiJH (2020). Identification of cell diversity and specific expressed genes in chicken breast muscle using single-cell RNA sequencing (Master of Agriculture thesis). Institute of Animal Sciences Graduate School,Beijing.

[45] LiuCM,HongTT, YuL,ChenY, WangSH,Ren ZJ (2023). Single-nucleus RNA and ATAC sequencing uncovers the molecular and cellular characteristics in the musk gland of Chinese forest musk deer (Moschus berezovskii). FASEB Journal 37,e22742.

[46] MansukhaniA,Moscatelli D,TalaricoD,LevytskaV,Basilico C (1990). A murine fibroblast growth factor (FGF) receptor expressed in CHO cells is activated by basic FGF and Kaposi FGF. PNAS 87,4378–4382.

[47] MiyazakiT,LiuZJ, KawaharaA et al. (1995). Three distinct IL-2 signaling pathways mediated by bcl-2,c-myc, and lck cooperate in hematopoietic cell proliferation. Cell 81,223–231.

[48] MulliganKA,Cheyette BNR (2012). Wnt signaling in vertebrate neural development and function. Journal of Neuroimmune Pharmacology 7,774–787.

[49] MesminB,BigayJ, Moser von FilseckJ,Lacas-GervaisS,DrinG, AntonnyB (2013). A four-step cycle driven by PI(4)P hydrolysis directs sterol/PI(4)P exchange by the ER-Golgi tether OSBP. Cell 155,830–843.

[50] MorikawaM,DerynckR, MiyazonoK (2016). TGF-β and the TGF-β Family: Context-dependent roles in cell and tissue physiology. Cold Spring Harbor Perspectives in Biology 8,a021873.

[51] MaJW,WeiKK, LiuJW et al. (2020). Glycogen metabolism regulates macrophage-mediated acute inflammatory responses. Nature Communications 11,1769.

[52] Magro-LopezE,Muñoz-Fernández MÁ (2021). The role of BMP signaling in female reproductive system development and function. International Journal of Molecular Sciences 22,11927.

[53] NorrieJL,LupoMS, XuBS et al. (2019). Nucleome dynamics during retinal development. Neuron 104,512–528.

[54] OrnitzDM,XuJ, ColvinJS et al. (1996). Receptor specificity of the fibroblast growth factor family. Journal of Biological Chemistry 271,15292–15297.

[55] ParkJ,Shrestha R,QiuCX et al. (2018). Single-cell transcriptomics of the mouse kidney reveals potential cellular targets of kidney disease. Science 360,758–763.

[56] Polycarpou-SchwarzM,Groß M,MestdaghP et al. (2018). The cancer-associated microprotein CASIMO1 controls cell proliferation and interacts with squalene epoxidase modulating lipid droplet formation. Oncogene 37,4750–4768.

[57] PadyanaAK,GrossS, JinL et al. (2019). Structure and inhibition mechanism of the catalytic domain of human squalene epoxidase. Nature Communications 10,97.

[58] PervolarakisN,NguyenQH, WilliamsJ et al. (2020). Integrated single-cell transcriptomics and chromatin accessibility analysis reveals regulators of mammary epithelial cell identity. Cell Reports 33,108273.

[59] RottaR,NoackA (2011). Multilevel local search algorithms for modularity clustering. ACM Journal of Experimental Algorithmics 16,2.1–2.27.

[60] SamalB,SunY, StearnsG,Xie C,SuggsS,McNieceI (1994). Cloning and characterization of the cDNA encoding a novel human pre-B-cell colony-enhancing factor. Molecular and Cellular Biology 14,1431–1437.

[61] ShresthaMN (1998). Animal welfare in the musk deer. Applied Animal Behaviour Science 59,245–250.

[62] ShannonP,MarkielA, OzierO et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Research 13,2498–2504.

[63] St ClairRM,DumasCM, WilliamsKS et al. (2019). PKC induces release of a functional ectodomain of the guidance cue semaphorin6A. FEBS Letters 593,3015–3028.

[64] TrevinoAE,Müller F,AndersenJ et al. (2021). Chromatin and gene-regulatory dynamics of the developing human cerebral cortex at single-cell resolution. Cell 184,5053–5069.

[65] VillarJ,ZhangHB, SlutskyAS (2019). Lung repair and regeneration in ARDS: Role of PECAM1 and Wnt signaling. Chest 155,587–594.

[66] WesterlingT,Kuuluvainen E,MäkeläTP (2007). Cdk8 is essential for preimplantation mouse development. Molecular and Cellular Biology 27,6177–6182.

[67] WoodA,Houghton SL,BiswasA (2016). A comparative study of immunohistochemical myoepithelial cell markers in cutaneous benign cystic apocrine lesions. The American Journal of Dermatopathology 38,475–483.

[68] WangYL,HaCY (2018). Research progress on musk and artificial propagation technique of forest musk deer. China Journal of Chinese Materia Medica 43,3806–3810.

[69] XuZX,JieH, ChenBL et al. (2017). Illumina-based de novo transcriptome sequencing and analysis of Chinese forest musk deer. Journal of Genetics 96,1033–1040.

[70] XieT,WangYZ, DengN et al. (2018). Single-cell deconvolution of fibroblast heterogeneity in mouse pulmonary fibrosis. Cell Reports 22,3625–3640.

[71] YangQS,MengXX, XiaL,FengZJ (2003). Conservation status and causes of decline of musk deer (Moschus spp.) in China. Biological Conservation 109,333–342.

[72] YoungMD,Mitchell TJ,Vieira BragaFA et al. (2018). Single-cell transcriptomes from human kidneys reveal the cellular identity of renal tumors. Science 361,594–599.

[73] YangJM,PengGF, ShuF et al. (2021). Characteristics of steroidogenesis-related factors in the musk gland of Chinese forest musk deer (Moschus berezovskii). Journal of Steroid Biochemistry and Molecular Biology 212,105916.

[74] ZouboulisCC,XiaLQ, DetmarM et al. (1991). Culture of human sebocytes and markers of sebocytic differentiation in vitro. Skin Pharmacology and Physiology 4,74–83.

[75] ZhaoLA,LiYH, WangBX (2009). Artificial propagation technology of forest musk deer. Forestry of Shaanxi 2009 (5),45.

[76] ZhouCY (2012). Summary of forest musk deer breeding management experience. China Animal Healthcare 14,63–65.

[77] ZhaoL,ChenXJ, ZhuJ et al. (2015). Lanosterol reverses protein aggregation in cataracts. Nature 523,607–611.

[78] ZengWH,JiangS, KongXD et al. (2016). Single-nucleus RNA-seq of differentiating human myoblasts reveals the extent of fate heterogeneity. Nucleic Acids Research 44,e158.

[79] ZhouWQ,JiZC, FangWX,Ji HK (2019). Global prediction of chromatin accessibility using small-cell-number and single-cell RNA-seq. Nucleic Acids Research 47,e121.

[80] ZhengCL,JieH, ZhouL et al. (2021). Tissue structure and developmental characteristics analysis of sachet of Moschus berezouskii. Heilongjiang Animal Science and Veterinary Medicine 2021 (8),126–129,154–156.

RIGHTS & PERMISSIONS

2024 The Authors. Integrative Zoology published by International Society of Zoological Sciences, Institute of Zoology/Chinese Academy of Sciences and John Wiley & Sons Australia, Ltd.