Blog

The musk chemical composition and microbiota of Chinese forest musk deer males | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  6, Article number: 18975 (2016 ) Cite this article china hpmc e

Male musk deer secrete musk from the musk gland located between their naval and genitals. Unmated male forest musk deer generate a greater amount of musk than mated males, potentially allowing them to attract a greater number of females. In this study, we used gas chromatography and mass spectrometry (GC/MS) to explore musk chemical composition of the musk pods of captive mated and unmated sexually mature Chinese forest musk deer and used next-generation sequencing to intensively survey the bacterial communities within them. Analysis of the chemical composition of the musk showed that unmated males have more muscone and cholesterol. Features of the musk16S rRNA gene showed that mated Chinese forest musk deer have both a greater Shannon diversity (p < 0.01) and a greater number of estimated operational taxonomic units than unmated ones; many bacterial genera were overrepresented in unmated Chinese forest musk deer males. Members of these genera might be involved in musk odor fermentation. PICRUSt analysis revealed that metabolic pathways such as aldosterone-regulated sodium reabsorption, metabolism of terpenoids and polyketides, flavone and flavonol biosynthesis and isoflavonoid biosynthesis were enriched in the musk of unmated Chinese forest musk deer males.

The dwarf musk deer or Chinese forest musk deer (Moschus berezovskii Flerov) is an endangered artiodactyl species native to southern and central China and northernmost Vietnam1,2. Males and females are visibly very similar; neither carries antlers, but males are armed with two abnormally long canine teeth (Figure S1) in the upper jaw that are used in fights with rivals. On the abdomen, the sexually mature males carry the eponymous musk-pod3,4, into which a musk gland (Figure S2) secretes the musk, generally described as of the color and consistency of ‘moist gingerbread’ and which appears to serve for attracting the females5,6. Its secretion has a specific musk odor and therefore the secretion’s chemical composition may be involved in chemical communication, potentially encoding information about sex and maturity7.

All animals harbor communities of microbes that profoundly affect their biology, often in beneficial ways8,9. It is also becoming obvious that symbiotic microbes can extend the host’s behavioral phenotypes in beneficial ways, including facilitating their feeding, antipredator, reproductive and communicative behaviors10,11. The fermentation hypothesis for mammalian chemical communication posits that bacteria in the scent glands of mammals generate odorous metabolites used by their hosts for communication and that variation in host chemical signals is a product of underlying variation in the bacterial communities inhabiting the scent glands12. The musk owes its odor to a preponderance of a particular ketone, known as muscone6 which is generated by the preputial gland (Figure S2). Muscone is claimed to be the sexual attractant for musk deer females13. Previous studies have shown that myocardial ketone body metabolism can be regulated by the gut microbiota during nutrient deprivation14 giving us a hint that there may be interaction between the microbiome and compounds produced by the host and potentially, musk pod microbiota and musk generation in Chinese forest musk deer. However, little is known about the composition of the microbiota of musk pod of male Chinese forest musk deer. In this study, we used next-generation sequencing to characterize and compare the bacterial communities in the musk pods of mated and unmated Chinese forest musk deer males (Moschus berezovskii Flerov).

The fresh musk production records of captive Chinese forest musk deer males aged between 3.5–4.5 years were summarized in Table S1. The largest pouch obtained from an adult musk deer was from an unmated male (UM) (31.532 g), whilst the smallest, was only 0.422 g, in a mated male (MM). We found that musk production after mating was significantly less (p < 0.001) in mated than unmated deer. These 10 samples were collected and sent for sequencing. Sequence information of each sample was summarized in Table S1. The raw sequences of this study have been deposited in the Sequence Read Archive (accession number SRP059015).

Musk was separated by diethyl ether and absolute ethyl alcohol solution treatment. All potential musk chemical compositions of Chinese forest musk deer were listed in Table S2. The table only lists the top 15 peak areas of each sample. It was found that both solvents can extract 3-Methylcyclopentadecanone (muscone). Alkanes account for about 80% of the diethyl ether extract; while ethyl alcohol extract contains more steroid, aliphatic ester and cholesterol. After rapid conventional gas chromatography and mass spectrometry (GC/MS), the abundances of the main chemical constituents of the musk samples of mated and unmated Chinese forest musk deer males were determined (Fig. 1, Table 1 and Table S2). Percentages of muscone in unmated males are greater than mated males, 22.59% versus 18.55% and 35.76% versus 26.19% for diethyl ether and ethyl alcohol extract respectively. On the contrary, percentages of 1,1-diethoxy-Ethane (acetal) in MM are greater than UM (22.8 versus 19.6%).

GC/MS chromatograms of the musk diethyl ether extraction of mated males (A) and unmated males (C), musk ethyl alcohol extraction of mated males (B) and UM (D). Only composition proportions greater than 5% were marked in all panels. The identities and relative abundances of the compounds are listed in Table 1.

In this study, we characterized the musk microbiota of 5 mated male Chinese forest musk deer and compared them with that of 5 unmated males living in the same facility in Chongqing Institute of Medicinal Plant Cultivation. We performed multiplex pyrosequencing of the V3-V4 hyper-variable regions of 16S rRNA gene. 16S rRNA gene analysis using an operational taxonomic unit (OTU) definition of sharing ≥97% nucleotide sequence identity revealed that the bacterial communities in mated and unmated musk pods were markedly different but each was dominated by fermentative anaerobes. We generated a dataset consisting of 142,179 filtered, high-quality, classifiable 16S rRNA gene sequences with an average of 15572.80 ± 1166.56 (SE) and 12863.00 ± 2230.59 (SE) per sample in mated and unmated males. More than 95% of the sequences in all samples were found to belong to the five most populated bacterial phyla (namely Firmicutes, Actinobacteria, Fusobacteria, Proteobacteria and Bacteroidetes) and one archaea phyla (Euryarchaeota) (Table 2). For bacteria, Firmicutes, Actinobacteria, Fusobacteria and Bacteroidetes were present in higher proportions in UM than in MM deers’ microbiota (32.32% versus 5.45%, 26.63% versus 5.50%, 14.65% versus 0.47% and 4.92% versus 1.72%, respectively). In contrast, Proteobacteria were more abundant in mated than in unmated males (29.44% versus 11.88%) (Table 2). For archaea, only Euryarchaeota were present in significantly higher proportions in MM than in UM deer (52.67% versus 7.67%, p < 0.01, Mann-Whitney U test). The majority of OTUs belonged to the division Firmicutes (167 ± 13.8) and Proteobacteria (351.6 ± 40.7) for UM and MM group, respectively (Table 2). The largest single OTU of UM and MM samples, is a member of the division Actinobacteria and genus Corynebacterium (26.02%) and the division Euryarchaeota and genus Methanosaeta (35.51%) respectively (Fig. 2).

Composition variation of the unmated and mated microbiota analyzed employing 16S rRNA as biomarkers.

Bars represent the relative abundance of bacterial taxa (genus and order). The top 20 abundant genus or orders are shown.

Four alpha diversity measures were calculated, including Shannon’s diversity index, ACE diversity index, observed species (observed OTUs) and Chao1 (estimated OTUs) (Fig. 3). Rarefaction curves using Shannon measure of alpha diversity were generated from musk samples obtained from both MM and UM deer (Figure S3). Each additional sample adds OTUs to the plot that were not observed in previous samples. The curve becomes asymptotic as the OTU number saturates and each sample adds an increasingly smaller number of new OTUs, indicating adequate coverage for the environment being tested. We found significant difference in Shannon diversity (p < 0.01) between MM and UM musk samples (Fig. 3A). For community richness comparison, mated musk had significantly higher number of observed and estimated (Chao1) OTUs than unmated musk (p < 0.05, Fig. 3C,D and Figure S3). No significant differences in ACE diversity indices were observed between MM and UM musk samples (Fig. 3B).

Variations in alpha diversity of the musk microbiota.

(A,B) Comparisons of Shannon diversity (A), ACE diversity indices (B), the number of observed OTUs (C) and the number of estimated OTUs (Chao1) (D) between mated and unmated Chinese forest musk deer males by Mann-Whitney test. All alpha diversity metrics were calculated upon the rarified OTU subsets. In all panels, boxes represent the interquartile range (IQR) between the first and third quartiles. The lines inside boxes represent the median.

Beta-diversity comparisons were performed using the Greengenes reference tree and computing phylogeny-based UniFrac distances15 between samples within QIIME. For exploring relationships between MM and UM musk samples’ phylogeny, both principal coordinate analysis (PCoA) on UniFrac distance matrices within QIIME and phylogenetic network analysis of community similarity16 were used to reconstruct networks from UniFrac distance matrices with SplitsTree417 using the neighbour-netagglomerative method18. The relationship between the musk microbiotas of UM and MM were visualized by a dendrogram. Each branch on the tree represents one musk microbiota (Fig. 4). As expected, the neighbour-net network reconstructed from the unweighted UniFrac distance matrix on this data set still shows a clear pattern of clustering by their mating status (Fig. 4). Principal coordinate analyses (PCoA) were also used to examine the relationships between musk microbiota of MM and UM Chinese forest musk deer (Figure S4) and similar clustering patterns were observed. On the PCoA plot, each circle represents one gut microbiota. Consistent with the dendrograms, the musk microbiotas of the MM were distinct from those of UM. These results suggest that the microbiota may play a potential role in chemical composition of musk in Chinese forest musk deer.

Clustering analysis of the evolution of the musk microbiotas of the mated and unmated Chinese forest musk deer males.

Adjacent bar charts show the phylum level classification (as determined by the RDP classifier) for each of the sequences per sample.

The community structures of musk gland microbiota are significantly different between mated and unmated Chinese forest musk deer males. Members of the Corynebacterium, Fusobacterium, Jeotgalicoccus, Mannheimia and Planomicrobium are five of the several bacterial genera that were overrepresented in the UM microbiota. We did not find any evidence of the genus Mannheimia in the MM samples. The genus Mannheimia contains five species from the family Pasteurellaceae that have been recently reclassified19. The five named species within this genus are M. haemolytica, M. varigena, M. glucosida, M. ruminalis and M. granulomatis19. Several unnamed taxa are also distinct from these named species20. The overrepresentation of this genus and other genus in unmated Chinese forest musk deer males indicates that members of this genus might be involved in fermentation.

When the third-level KEGG pathways of these bacteria were examined, genes controlling metabolic pathways, including aldosterone-regulated sodium reabsorption, endocrine and other factor-regulated calcium reabsorption, xylene degradation and aprolactam degradation (Fig. 5) were enriched in unmated males. Aldosterone-regulated sodium reabsorption and endocrine and other factor-regulated calcium reabsorption belong to the excretory system pathway; xylene degradation and aprolactam degradation belong to the xenobiotics biodegradation and metabolism pathway, but the relationships of these pathways with musk odor is still unknown.

Predicted function of musk microbiota between musk of mated and unmated Chinese forest musk deer males.

The gene copy numbers of five samples within the same sample group were pooled. Log10-transformed values for each functional gene (row) were used. KEGG pathway of the third level was shown in the heatmap. Statistical significance between the two groups of the gene distribution were performed using Mann-Whitney U-test with cutoffs of p < 0.01.

The primary biological significance of the musk gland secretion is to attract the females to males during the rutting season and to mark territory7. Muscone is the main active ingredient which is claimed to be the sexual attractant for musk deer females13. In this study we found that these unmated Chinese forest musk deer males secrete more musk than mated ones. By using gas chromatography and mass spectrometry we studied the chemical composition of mated and unmated males. The composition of musk from the Chinese forest musk deer preputial gland secretion coincides with data reported by Do et al. (1975)21 (see Table 1), with muscone accounting for 22.6% of secretion extract for Nepal musk deer21. It is known that acetalisation, or the organic reaction that involves the formation of an acetal, is a reversible reaction. One method of acetal formation is the nucleophilic addition of an alcohol to a ketone or an aldehyde. We found that muscone is more and acetal (1,1-diethoxy-Ethane) is less in unmated males, whilst muscone is less and acetal is more in mated males. Steroid compounds (cholestanol, cholesterol and a number of the androstane derivatives) make up a large proportion of musk secretion components, which concurs with previously reported data7. However, these previous studies only found approximately ten kinds of chemical components, while we found more than 50 components in this study. Additional studies should be carried out to understand the relationship between mating status and musk chemical compounds.

All animals are inhabited by a mixed population of microbes and contrary to popular belief, most microbes appear highly beneficial to their hosts. They are of critical importance in animal nutrition and immune defense and can serve as crucial catalysts for the effective development and functioning of host tissues. Their contribution to host behavior is also becoming increasingly appreciated. It has been hypothesized that one way they do so is by producing the components of chemical signals that animals use to communicate. Our results demonstrate that the bacterial profiles of Chinese forest musk deer diverge between mated and unmated males. The production of musk is associated with its mating status. We propose the possibility that it isn’t the bacteria that are triggering musk production, but that after mating, the deer start to produce less musk through a physiological pathway and only some bacterial strains can thrive in the new low-musk environment. We are also assuming that once the male mates, he produces little musk until the next rut, at which point he begins musk production and his microbiome changes to the unmated profile; given the anatomical position of the musk pod, is it also possible that the pod is getting “seeded” with novel bacteria from the female during the act of mating, further study is needed to confirm these hypothesis.

The microbiome plays an important role in human odour production and the number of certain microorganisms is strongly correlated with the intensity of the odour emitted22,23. Kevin et al.12, for example, showed that symbiotic bacteria underlie species-specific odors in both spotted and striped hyenas and further underlie sex and reproductive state-specific odors among spotted hyenas. The skin’s location at the interface with the outside world therefore makes it the most subject to environmental influences that will affect the microbiota24. Chinese forest musk deer males secrete a specific musk odor and therefore the secretion’s chemical composition can be involved in chemical communication7, which appears to serve for attracting females5,6. Since mated Chinese forest musk deer males and unmated males produce various amount of musk, in combination with the fact that the musk pod has an opening to contact outside, it is easy to postulate that animals may share different musk microbiota. In this study, our results showed that mated Chinese forest musk deer males and unmated males harbor different musk microbiotas. The richness of musk microbial diversity in each sample was estimated by rarefaction curves and ACE, Chao1 and Shannon indexes were calculated by MOTHUR (Fig. 2 and Figure S3). We concluded that OTU richness will increase with additional number of sequences. At least 1,732 and 3,061 OTUs were identified from UM and MM deer, respectively. Mated males showed diverse α-diversity than unmated males.

Despite some potential biases25, PICRUSt provides some insight into bacterial community functions in musk of Chinese forest musk deer. Muscone is claimed to be the sexual attractant for musk deer females13. In this study, we found many functional pathways related to ketone metabolism overrepresented in unmated Chinese forest musk deer males, for example, aldosterone-regulated sodium reabsorption, metabolism of terpenoids and polyketides, flavone and flavonol biosynthesis and isoflavonoid biosynthesis pathways. It is important to note that the PICRUSt analysis should be considered speculative because there are also significant other unrelated pathways that appear such as Parkinson’s disease etc. The relevance of application of PICRUSt to predict bacterial activities in musk of Chinese forest musk deer needs further study. Further transcriptomics and metabolomics studies are desired to confirm these discoveries and improve our understanding of the bacterial functions in the musk of Chinese forest musk deer.

In summary, we characterized the musk microbiotas of the mated and unmated Chinese forest musk deer males and found that, according to the 16S rRNA based community structure analysis, the musk microbiota of these deer diverged based on their different mating condition. We also identified bacterial taxa differentially represented between the mated and unmated males. Further studies are needed to examine differential microbial roles in the hosts’ biochemical pathways and physiology.

The behavioral observations were carried out with 10 healthy Chinese forest musk deer males, kept in open-air yards, between 15th Nov 2013 and 15th Dec 2013, which coincided with the rutting season. Each male, along with 4 healthy females (3.5–4.5 years old with the similar body size and never mated before) was housed in individual breeding yards (about 500 m2), each yard have one big play ground and five small living rooms (about 8 m2 each) (Figure S5). Mating behavior was monitored by HIKVISION video cameras (Hangzhou, China), artificial observation by farmers and the pregnancy of any females. After behavioral monitoring was complete, the farmers separated males and females with the least disturbance, as the females are easily frightened in the early stages of pregnancy. In order to prevent fighting among males, they were housed separately in 10 living rooms (12 m2) and there is a playground in front of those rooms (about 60 m2). Females were left in the original breeding yards. If we observed a male mating with females, they were classified as mated males. Meanwhile, if we did not observe any mating behavior or pregnancy they were classified as unmated males.

The mating season of the musk deer is from November to January26. The fresh musk specimens from 10 individuals (5 mated and 5 unmated) were obtained, weighed and recorded in early October (2014) (Table S1). They were all born in captivity and had not mated before. Their age range from 4.5 to 5.5 years old (Table S1). All samples of the musk pouches used for bacterial surveys were obtained directly from the musk gland scent pods of Chinese forest musk deer in the Chongqing Institute of Medicinal Plant Cultivation (Chongqing, China, altitude: 678m).

Musk deer were anesthetized with Telazol (4.5 mg/kg) delivered from a CO2 -powered darting rifle, especially while the musk deer was alone and lying down. A special spoon was used for musk collection from the pod. Musk samples were placed in sterile cryogenic vials, stored in liquid nitrogen and transported to Sichuan Agricultural University, where they remained frozen at −80 °C until their bacterial profiles were determined. All animal experimentation were approved by the Institutional Animal Care and Use Committee of Sichuan Agricultural University under permit number DKY- S20143204. The methods were carried out in accordance with the approved guidelines.

After microbial analysis, the rest of fresh musk was dried by exposure to phosphorus pentoxide (P2O5). We weighed the musk every 30 minutes until its weight ceased to change. Because a very small amount of musk remained for the mated group after bacterial DNA extraction and musk dehydration so we pooled three samples together (taking 0.07 gram of each sample, total 0.21 gram). After blending, each 0.1 gram of musk was dissolved in 2.5 ml absolute ethyl alcohol (purity > 99.7%) or 2.5 ml diethyl ether (purity > 99.5%), then extracted by ultrasonic for two hours. The dissolved samples were centrifuged at 13000g for 5 minutes and the liquid supernatant was used for chemical components analysis. The same method was applied for unmated (sample UM1, UM2 and UM3) males.

A Shimadzu GCMS-QP 2010 plus gas chromatograph (Shanghai, China), coupled to a Shimadzu QP 5000 mass spectrometer which was equipped with a split/split less injector and a DB5-MS column [(30 m × 0.25 mm i.d., 0.25 mm film thickness), (Agilent J&W, Agilent Technologies, CA, USA)] was used for the chromatographic analysis of the targeted compounds. The mass spectral data base NIST08.LIB was used to analyze the data. The confidence coefficient of the data above 80% was adopted.

The injector and mass transfer-line temperatures were set at 280 °C and 290 °C, respectively. The column temperature was programmed at 60 °C for 5 min, increased to 150 °C at a rate of 10 °C/min and held for 5 min, then ramped at 5 °C/min to 280 °C (held for 5 min). The total run time was 58.5 min. The column pressure was set to 49.5 kPa and the column flow rate was 0.9 ml/min. The linear velocity was 34.7 cm/sec. The flux purge was 3.0 ml/min. The solvent cut time was set to 3.0 min. Sample injection was performed by autoinjection. The injection volume was 1.0 μL. The split ratio of diethyl ether extract and ethyl alcohol extract was set to 5:1 and 10:1, respectively. Helium (purity ≥ 99.99%) was used as the carrier gas at a flow rate of 1.0 ml/min. Selective ion monitoring (SIM) mode was adopted for the determination of the analytes.

Total bacteria DNA was extracted from each 0.4 g musk samples by using MO BIO UltraClean fecal DNA kit (Carlsbad, CA, USA). Bacterial 16S rRNA genes in musk were PCR amplified using a set of two broadly conserved, degenerate primers targeting the V3–V4 variable regions of the 16S gene (341F: 5′–CCC TAC ACG ACG CTC TTC CGA TCT NXX XXX XXC CTA CGG GNG GCW GCA G–3′; 805R: 5′–GTG ACT GGA GTT CCT TGG CAC CCG AGA ATT CCA GAC TAC HVG GGT ATC TAA TCC–3′), the forward primer containing a 7-bp error-correcting barcode unique to each sample.

PCR reaction was conducted using a Mastercycler Gradient Thermal Cycler (Eppendorf) with the following condition: 94 °C for 3 min (initial dissociation), followed by 30 cycles of [94 °C for 30 s (denaturation), 60 °C for 45 s (annealing) and 72 °C for 1 min (extending)], followed by a final extension of 8 minutes at 72 °C. Negative non template controls were run for each barcoded primer pair to test for reagent contamination. Each sample was amplified and then pooled together for analysis. The PCR products were analyzed by gel electrophoresis and purified by using a PCR Purification kit (Sangon, Shanghai, China). Pyrosequencing was performed on an Illumina MiSeq 2 × 250 platform.

QIIME27 pipeline was mostly used for data quality controls and analyses. Raw sequence reads were filtered to meet minimum and maximum length of 200 bp and 450 bp with no ambiguous base calls. Sequences that contained more than one ambiguous base call (N) and did not also have a complete barcode and primer at one end were eliminated. A read was discarded if it was identified as a putative chimera by UCHIME28. The sequences that passed the above procedure were denoised in order to correct for potential sequencing errors at 99% level using UCLUST28. Operational Taxonomic Unit (OTUs) were defined using UCLUST v1.1.579, which have been proven to generate satisfactory and comparable numbers of OTU29. The degree of similarity between sequences was defined as ≥97% as a threshold to obtain OTU identity at the species level. RDP classifier (version 2.10) software was used30 to classify the sequences according to the taxonomy proposed by Garrity et al.31, maintained at the Ribosomal Database Project (RDP 10 database, Update18). Richness, Shannon Index, ACE index, Chao1 and Coverage were included in Alpha diversity analysis by using MOTHUR32. Mann-Whitney U test was used for significance test. Weighted unifrac and unweighted unifrac metric distances, each with 10× subsampling, were calculated to obtain beta diversity index, which used phylogenetic tree to compute phylogenetic relatedness of the bacterial community between samples. These results were viewed by Principal Coordinate Analysis (PCoA)15. A two-sided Student’s t-test was used for significance test of beta diversity difference between mated and unmated musk deer. Metagenome functional content from 16S rRNA was predicted using PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States)33. The R package ‘pheatmap’ was used for data analysis and plotting.

How to cite this article: Li, D. et al. The musk chemical composition and microbiota of Chinese forest musk deer males. Sci. Rep. 6, 18975; doi: 10.1038/srep18975 (2016).

He, L. et al. Effects of crowding and sex on fecal cortisol levels of captive forest musk deer. Biological research 47, 48, 10.1186/0717-6287-47-48 (2014).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Wang , YH & Moschus berezovskii , RB The IUCN Red List of Threatened Species .Version 2014.3.Downloaded on 30 December 2014 (2008), http://www.iucnredlist.org/details/13894/0

Green, M. The distribution, status and conservation of the Himalayan musk deer Moschus chrysogaster. Biological Conservation 35, 347–375, 10.1016/0006-3207(86)90094-7 (1986).

Green, M. J. Scent-marking in the Himalayan musk deer (Moschus chrysogaster). Journal of Zoology 1, 721–737 (1987).

Feng, W., You, Y., Yong, H., Li, G. & Gu, D. Histological observation of Chinese forest musk deer (Moschus berezovskii Flerov) musk glands (in Chinese). The Zoological Journal, 13–15 (1981).

Hawkins, T. H. Musk and the musk deer. Nature 166, 262 (1950).

Article  CAS  ADS  PubMed  Google Scholar 

Sokolov, V. E., Kagan, M. Z., Vasilieva, V. S., Prihodko, V. I. & Zinkevich, E. P. Musk deer (Moschus moschiferus): Reinvestigation of main lipid components from preputial gland secretion. Journal of chemical ecology 13, 71–83, 10.1007/BF01020352 (1987).

Article  CAS  PubMed  Google Scholar 

Gilbert, S. F., Sapp, J. & Tauber, A. I. A symbiotic view of life: We have never been individuals. The Quarterly review of biology 87, 325–341 (2012).

McFall-Ngai, M. et al. Animals in a bacterial world, a new imperative for the life sciences. Proceedings of the National Academy of Sciences 110, 3229–3236 (2013).

Article  CAS  ADS  Google Scholar 

McCune, B., Grace, J. B. & Urban, D. L. Analysis of ecological communities. Vol. 28 (MjM software design Gleneden Beach, OR, 2002).

Lee, Z. M.-P., Bussema, C. & Schmidt, T. M. rrnDB: documenting the number of rRNA and tRNA genes in bacteria and archaea. Nucleic acids research 37, D489–D493 (2009).

Article  CAS  PubMed  Google Scholar 

Theis, K. R. et al. Symbiotic bacteria appear to mediate hyena social odors. Proceedings of the National Academy of Sciences 110, 19832–19837 (2013).

Article  CAS  ADS  Google Scholar 

Theimer, E. T. Fragrance Chemistry: The Science of the Sense of Smell. Ch. 12, 433 (Mookherjee, B. D. & Wilson, R. A.) The chemistry and fragrance of natural musk compounds (2012).

Crawford, P. A. et al. Regulation of myocardial ketone body metabolism by the gut microbiota during nutrient deprivation. Proceedings of the National Academy of Sciences 106, 11276–11281 (2009).

Article  CAS  ADS  Google Scholar 

Mount, D. W. Using the Basic Local Alignment Search Tool (BLAST). CSH protocols 2007, pdb top17, 10.1101/pdb.top17 (2007).

Parks, D. H. & Beiko, R. G. Measuring community similarity with phylogenetic networks. Mol Biol Evol 29, 3947–3958, 10.1093/molbev/mss200 (2012).

Article  CAS  PubMed  Google Scholar 

Huson, D. H. & Bryant, D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 23, 254–267, 10.1093/molbev/msj030 (2006).

Article  CAS  PubMed  Google Scholar 

Bryant, D. & Moulton, V. Neighbor-net: an agglomerative method for the construction of phylogenetic networks. Mol Biol Evol 21, 255–265, 10.1093/molbev/msh018 (2004).

Article  CAS  PubMed  Google Scholar 

Angen, O., Mutters, R., Caugant, D. A., Olsen, J. E. & Bisgaard, M. Taxonomic relationships of the [Pasteurella] haemolytica complex as evaluated by DNA-DNA hybridizations and 16S rRNA sequencing with proposal of Mannheimia haemolytica gen. nov., comb. nov., Mannheimia granulomatis comb. nov., Mannheimia glucosida sp. nov., Mannheimia ruminalis sp. nov. and Mannheimia varigena sp. nov. International journal of systematic bacteriology 49 Pt 1, 67–86 (1999).

Article  CAS  PubMed  Google Scholar 

Angen, O., Ahrens, P. & Bisgaard, M. Phenotypic and genotypic characterization of Mannheimia (Pasteurella) haemolytica-like strains isolated from diseased animals in Denmark. Veterinary microbiology 84, 103–114 (2002).

Article  CAS  PubMed  Google Scholar 

Do, J., Kitatsuji, E. & Yoshii, E. Study on the components of musk. I. Ether soluble components. Chemical and Pharmaceutical Bulletin 23, 629–635 (1975).

Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J Mol Biol 215, 403–410, 10.1016/S0022-2836(05)80360-2 (1990).

Article  CAS  PubMed  Google Scholar 

Jackman, P. & Noble, W. Normal axillary skin in various populations. Clinical and experimental dermatology 8, 259–268 (1983).

Article  CAS  PubMed  Google Scholar 

Sanford, J. A. & Gallo, R. L. Functions of the skin microbiota in health and disease. Seminars in immunology 25, 370–377, 10.1016/j.smim.2013.09.005 (2013).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Zeng, B. et al. The bacterial communities associated with fecal types and body weight of rex rabbits. Scientific reports 5, 10.1038/srep09342 (2015).

Yan, Y. A preliminary study on the factors affecting musk production in Siberian musk deer (in Chinese). Journal of Chinese Medicinal Materials 2, 11–13 (1985).

Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nature methods 7, 335–336 (2010).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200, 10.1093/bioinformatics/btr381 (2011).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Matsuda, F., Tsugawa, H. & Fukusaki, E. Method for assessing the statistical significance of mass spectral similarities using basic local alignment search tool statistics. Analytical chemistry 85, 8291–8297, 10.1021/ac401564v (2013).

Article  CAS  PubMed  Google Scholar 

Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and environmental microbiology 73, 5261–5267 (2007).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Garrity, G. M. et al. Introduction to the Taxonomic Outline of Bacteria and Archaea (TOBA) Release 7.7. The Taxonomic Outline of Bacteria and Archaea 7, 1–5 (2007).

Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and environmental microbiology 75, 7537–7541 (2009).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Langille, M. G. et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nature biotechnology 31, 814–821 (2013).

Article  CAS  PubMed  PubMed Central  Google Scholar 

This work was supported by the Program from Sichuan Agricultural University (02920400); National Natural Science Foundation of China (NSFC31160432) and Sichuan Provincial Department of Science and Technology Program (2015JQO023).

Li Diyan, Chen Binlong and Zhang Long contributed equally to this work.

Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, 61130, P.R. China

Diyan Li, Binlong Chen, Long Zhang, Uma Gaur, Tianyuan Ma, Nan Wu, Zhongxian Xu, Huailiang Xu, Yongfang Yao, Ting Lian, Xiaolan Fan, Deying Yang, Mingyao Yang & Qing Zhu

Laboratory of Medicinal Animal, Chongqing Institute of Medicinal Plant Cultivation, Nanchuan, 408435, Chongqing, P.R. China

Hang Jie & Guijun Zhao

Department of Anthropology and California National Primate Research Center, University of California, Davis, California

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

Conceived and designed the experiments: D.L., Q.Z. and J.S.T., Performed the experiments: D.L., L.Z., Q.Z., B.C., T.L., X.F., D.Y. and U.G., Contributed reagents/materials/analysis tools: B.C., N.W., H.J., T.M., M.Y., G.Z., H.X. and Y.Y., Wrote the paper: D.L., J.S.T. and Q.Z.

The authors declare no competing financial interests.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Li, D., Chen, B., Zhang, L. et al. The musk chemical composition and microbiota of Chinese forest musk deer males. Sci Rep 6, 18975 (2016). https://doi.org/10.1038/srep18975

DOI: https://doi.org/10.1038/srep18975

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Monthly magazines for chemistry - Chemical Monthly (2021)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

kaimaoxing hec uses Sign up for the Nature Briefing: Microbiology newsletter — what matters in microbiology research, free to your inbox weekly.