Arxius de Miscel·lània Zoològica. Volum 17 (2019) Pàgines: 209-218

mtDNA confirms the presence of Moschus leucogaster (Ruminantia, Moschidae) in Gaurishankar Conservation Area, Nepal

Shrestha, B. B., Khatiwada, J. R., Thanet, D. R.


Paraules clau

Moschus leucogaster, Anàlisi filogenètica, Àrea de Conservació de Gaurishankar, Nepal


Shrestha, B. B., Khatiwada, J. R., Thanet, D. R., 2019. mtDNA confirms the presence of Moschus leucogaster (Ruminantia, Moschidae) in Gaurishankar Conservation Area, Nepal. Arxius de Miscel·lània Zoològica, 17: 209-218, DOI:



Data de recepció:


Data d'acceptació:


Data de publicació:



mtDNA confirms the presence of Moschus leucogaster (Ruminantia, Moschidae) in Gaurishankar Conservation Area, Nepal
Musk deer (genus Moschus), an endangered mammal, is not only of great concern for its conservation, but also of great interest to understand its taxonomic and phylogenetic associations in Nepal. The aim of this study was to identify the taxonomic status of musk deer in Gaurishankar Conservation Area (GCA) using mitochondrial genomic data of cytochrome b (370 bps) through phylogenetic analysis of all the species of musk deer. The results showed that the species found in GCA is confirmed as Himalayan musk deer Moschus leucogaster further expanding its distributional range in Nepal.

Key words: Moschus leucogaster, Phylogenetic analysis, Gaurishankar Conservation Area, Nepal


El ADN mitocondrial confirma la presencia de Moschus leucogaster (Ruminantia, Moschidae) en el Área de Conservación de Gaurishankar, Nepal
El ciervo almizclero del Himalaya (género Moschus), un mamífero amenazado, no es solo objeto de gran preocupación por lo que respecta a su conservación sino también de gran interés para entender sus asociaciones taxonómicas y filogenéticas en Nepal. El objetivo de este estudio ha sido identificar el estatus taxonómico del ciervo almizclero del Himalaya en el Área de Conservación de Gaurishankar (GCA) utilizando datos genómicos mitocondriales del citocromo b (370 bps) mediante análisis filogenéticos de todas las especies de ciervo almizclero. Los resultados han mostrado que la especie hallada en GCA es el ciervo almizclero del Himalaya Moschus leucogaster, lo que amplía en mayor medida su área de distribución en Nepal.

Palabras clave: Moschus leucogaster, Análisis filogenético, Área de Conservación de Gaurishankar, Nepal


L’ADN mitocondrial confirma la presència de Moschus leucogaster (Ruminantia, Moschidae) a l’Àrea de Conservació de Gaurishankar, Nepal
El cérvol mesquer de ventre blanc (gènere Moschus), un mamífer amenaçat, no és només objecte de gran preocupació pel que fa a la conservació sinó també de gran interès per entendre les seves associacions taxonòmiques i filogenètiques al Nepal. L’objectiu d’aquest estudi ha estat identificar l’estatus taxonòmic del cérvol mesquer de ventre blanc a l’Àrea de Conservació de Gaurishankar (GCA) utilitzant dades genòmiques mitocondrials del citocrom b (370 bps) mitjançant anàlisis filogenètiques de totes les espècies de cérvol mesquer. Els resultats han mostrat que l’espècie trobada es el cérvol mesquer de ventre blanc Moschus leucogaster, la qual cosa amplia encara més la seva àrea de distribució al Nepal.

Paraules clau: Moschus leucogaster, Anàlisi filogenètica, Àrea de Conservació de Gaurishankar, Nepal


With the advent of biotechnology, genetic information can be obtained through animal’s degraded remains such as bones, dried skins, faeces and fossils (Ramón-Laca et al., 2015; Silva et al., 2015). Genetic information can be accrued from DNA contained with these animal specimens. DNA test is the most popular method for identification of species. In this system, DNA sequence of mitochrondrial cytochrome b (cyt-b) gene has been widely used for species identifications (Khatiwada et al., 2015). Cyt-b is probably the best-known mitochondrial gene with respect to function and structure of its protein product (Esposti et al., 1993). It is used as a valuable tool for the construction of evolutionary relationship among population, species, and higher taxa (Harrison, 1989; Su et al., 1999) since it contains both slowly and rapidly evolving codon positions, as well as more conservative and more variable regions or domains (Meyer and Wilson, 1990; Irwin et al., 1991; Cantatore et al., 1994; Farias et al., 2001).

Musk deer (genus Moschus, Linnaeus, 1758) are widely distributed in the sub-alpine and alpine vegetation (2,500 to 4,500 m) of the Himalayan region of Nepal (Kattel, 1992). All the species of musk deer belongs to order Cetartiodactyla of family-Moschidae (IUCN, 2013). Although many morphological studies have been conducted on the taxonomy of this group, there are still controversies regarding the number of its species and sub-species and the phylogenetic relationship among them (Groves et al., 1995). Seven species within the genus Moschus are recognized in the world. They include Anhui musk deer (M. anhuiensis Wang et al., 1982 ), Kashmir musk deer (M. cupreus Grubb, 1982), Siberian musk deer (M. moschiferus Linnaeus, 1758 ), black musk deer (M. fuscus Li, 1981), Himalayan musk deer (M. leucogaster Hodgson, 1839), forest musk deer (M. berezovskii Flerov, 1929) and Alpine musk deer (M. chrysogaster Hodgson et al., 1839 ) (IUCN, 2013). Three species of musk deer ie M. chrysogaster, M. leucogaster and M. fuscus are said to be found in Nepal (Jnawali et al., 2011; Timmins and Duckworth, 2015; Wang and Harris, 2015; Harris, 2016) but their exact distribution and status are still dubious, and they all were considered as one species: M. chrysogaster before their taxonomic separation (Jnawali et al., 2011). Although all the three species of musk deer found in Nepal are categorized as Endangered by IUCN (IUCN, 2013), only M. chrysogaster is enlisted as protected mammal by National Park and Wildlife Conservation Act 1973 due to lack of documentation on the confirmation of species of musk deer (Jnawali et al., 2011). There are only few studies on taxonomic and phylogenetic analysis of musk deer in the context of Nepal (Singh et al., 2019). The Gaurishankar Conservation Area (GCA) is the newly established conservation area where molecular studies confirming the species presence are crucial for its conservation and management in the future. In this study, we aimed to determine the taxonomic status of musk deer in the GCA using mitochondrial genomic data of cytochrome b through phylogenetic analysis of all the species of musk deer.

Material and methods

Study area

GCA is located in central Nepal encompassing Ramechhap, Dolakha and Sindhupalchok districts and has an area of 2,179 km2 (fig. 1). It was declared as the ‘Conservation Area’ in January 2010, and entrusted its management for twenty years to National Trust for Nature Conservation (NTNC) in July 2010 (DNPWC, 2011). It is located in the high mountain physiographic region of Nepal which consists of 35.38 % forestland, 9.76 % shrubland and 8.79 % grassland respectively. It has 16 major vegetation types and has a great faunal diversity that includes 34 species of mammals, 16 species of fishes, 10 species of amphibians, 8 species of lizards, 14 species of snakes and 235 species of birds (DNPWC, 2013). Besides musk deer (Moschus spp.), endangered species found in the conservation area are snow leopard (Panthera uncia), clouded leopard (Neofelis lupus), leopard cat (Felis benghalensis), red panda (Ailurus fulgens), wolf (Canis lupus), and chinese pangolin (Manis pentadactyla) (DNPWC, 2013; Shrestha and Meng, 2014).

Fig. 1. Map showing the study area: A, location of Gaurishankar Conservation Area in Nepal; B, location of Lapche of Dolakha District where the samples were collected; C, location of the 39 sampling sites in the Lapche area of GCA. (The triangles indicate the fresh latrines from where samples were collected.) Fig. 1. Mapa que muestra el área de estudio: A, localización del Área de Conservación de Gaurishankar en Nepal; B, localizavión de Lapche, en el distrito de Dolakha donde se recolectaron las muestras; C, localización de los 39 puntos de muestreo en el área de Lapche de GCA. (Los triángulos indican los puntos donde se recogieron muestras fecales.)


Our study included the surrounding areas of Risan Gumbo Himal from Hum danda to Gumbo danda at Lapche area of Lamabagar VDC of Dolakha district which encompasses an elevation range of 3.500 to 4.200 m and lies between 28º 6′ 7” and 28º 7′ 3” N latitude and 86º 9′ 59” and 86º 10′ 52” E longitude. It has four types of vegetation: Betula forest mostly dominated by Betula utilis; mixed forest having a mixed species of Betula utilis, Abies spectabilis, Sorbus spp., Rhododendron campanulatum, Salix spp. and Juniperus indica; rhododendron forest mostly dominated by Rhododendron campanulatum; and alpine scrub mostly dominated by shrubby rhododendron species i.e., Rhododendron lepidotum, Rhododendron ciliatum and Rhododendron anthopogan.

Pellet sample collection

We opportunistically collected 39 fecal pellet samples (faeces) between the period of 2013 and 2014 from Lapche area of GCA (fig. 1) and preserved in 95 % ethanol for further molecular analysis. For each pellet sample, the sampling date, the location and geographical coordinates were recorded. Those sample pellets were used for laboratory analysis for species identification.

DNA extraction, PCR amplification and sequencing

A dry pellet from each sample was cut into smaller pieces using sterile scissors and tweezers and further processed for DNA extraction as recommended in protocol by Qiagen QIAamp DNA Stool kit (Qiagen, Valencia, CA, USA). Polymerase chain reaction (PCR) was used to amplify the mitochondrial genes cyt-b. The primers and PCR conditions were used as described by Pan et al. (2015). PCR products were visualized in 1.5 % agarose gel and positive PCR products were sequenced following bi-directional sequencing from ABI 3100 automated sequencer.

Sequence analysis

All available nucleotide sequences of the cyt-b gene of Moschus species were downloaded from the NCBI GenBank database and also data source made available by Pan et al. (2015) and Singh et al. (2019) (table 1). The sequences of Alces alces americana, Ovis aries and Tragulus kanchil were also downloaded to be used as outgroups. Among the 39 pellet samples, only 10 pellet samples yielded complete sequence. All nucleotide sequences were assembled by SeqMan and visually checked to determine the accuracy of the variables site identified by the program. All the sequences were then aligned with ClustalW in BIOEDIT Version 7.1.9 (Thompson et al., 1994) using the default settings. All of newly determined sequences were deposited in GenBank under accession numbers (MN720942-MN720951). Phylogenetic analysis was conducted using Maximum likelihood (ML). The maximum likelihood analysis was conducted in MEGA7 with 1,000 bootstraps (Kumar et al., 2016).

Table 1. GenBank accession number (N) of specimens used in the phylogenetic analysis. Tabla 1. Número de entreda en GenBank (N) de los especímenes utilizados en el análisis filogenético.


Results and discussion

Phylogenetic relationship of musk deer from GCA

The aligned dataset of cyt-b sequence contained 370 bps including 69 variable sites and 51 parsimony informative sites (excluding outgroups). The phylogenetic relationships strongly support genus Moschus as a monophyletic clade with higher bootstrap supports (bootstrap = 99). Molecular data analysis suggested that the population of musk deer in Lapche, GCA, were genetically similar with M. leucogaster from western Nepal and Qinhai, Tibet, China and were nested together in a ML tree (fig. 2). The uncorrected genetic divergence of the cyt-b gene sequences among the M. leucogaster population of Lapche, GCA, Manang, Nepal and Tibet, China were ranges from 0.00 to 0.3% respectively (table 2). Whereas, relatively low genetic divergences between M. leucogaster and its closest relatives M. chrysogaster and M. fuscus were 1.4 % and 1.7 % respectively.

Fig. 2. Maximum likelihood tree estimated based on mtDNA cyt–b sequences. Values on branches of the tree show the bootstrap support value for maximum likelihood. A. americana, O. aries and T. kanchil were selected as outgroups. Sample names correspond to those given in table 1. Fig. 2. Árbol de máxima probabilidad estimada basado en secuencias de ADN mitocondrial (cyt–b). Los valores de las ramas del árbol muestran el valor de soporte bootstrap para la máxima probabilidad. A. americana, O. aries y T. kanchil se seleccionaron como exogrupos. Los nombres de las muestras corresponden a los indicados en la tabla 1.


Table 2. Genetic uncorrected p–distance of the cyt–b sequences of the genus Moschus used in this study. Only two samples were used for this analysis as there was no genetic variation among 10 samples from Lapche. Tabla 2. Distancia genética no corregida (p–distance) de las secuencias cyt–b del género Moschus utilizadas en este estudio. Utilizamos únicamente dos muestras para este análisis dado que no existieron variaciones genéticas entre las 10 muestras de Lapche.


The pellet samples used in this study confirmed the presence of Himalayan musk deer M. leucogaster in GCA. Although three species of musk deer: Moschus fuscus, Moschus chrysogaster and Moschus leucogaster are said to be distributed in Nepal (Timmins and Duckworth, 2015; Wang and Harris, 2015; Harris, 2016), most of studies regarding distribution, habitat ecology, latrines and its associated plant composition and diversity, gastro intestinal parasites of musk deer conducted in central and eastern Nepal have considered musk deer species as M. chrysogaster (Aryal et al., 2010; Aryal and Subedi, 2011; Subedi et al., 2012; Shrestha and Moe, 2015; Achhami et al., 2016). The potential misidentification of musk deer species is due to its secretive behaviour and similar morphological characteristics (Groves et al., 1995; Su et al., 2000; Guha et al., 2007). Musk deer is shy and nocturnal hence difficult to detect in daytime (Green, 1986). Additionally, if they are seen in daytime they hide in shrub understory and even if they are encountered in forest and open area, they are visible only for a few seconds (Singh et al., 2019). Furthermore, it is not easy to identify the species of musk deer by considering their physical appearance and morphological characteristics (Groves et al., 1995; Singh et al., 2019). Singh et al. (2019) confirmed the presence of Himalayan musk deer M. leucogaster in Manang and Kaski district. Based on their study, we can predict that M. leucogaster distribution may extend up to central and eastern parts of Nepal. GCA lies in the central part of Nepal close to Kaski and Manang district. This study confirms the presence of Moschus leucogaster in central Nepal and further supports the distribution of this species may encompass the Central and Eastern Nepal including Langtang National Park, Makalu Barun National Park, Sagarmatha National Park and Kanchanjunga Conservation Area and adjoining areas of Central and Eastern Nepal.


This study was financially supported by Mohamed bin Zayed Species Conservation Fund, UAE (Project no. 13257839) and Idea Wild. We would like to thank Department of National Parks and Wildlife Conservation, Nepal and National Trust for Nature conservation for granting research permission in GCA. Dr. Janak Raj Khatiwada is supported by Chinese Academy of Sciences President’s International Fellowship Initiative (2018PB0016) postdoctoral fellowship.


Achhami, B., Sharma, H. P., Bam, A. B., 2016. Gastro intestinal parasites of musk deer (Moschus chrysogaster Hodson, 1839) in Langtang National Park, Nepal. Journal of Institute of science and Technology, 21: 71-75.
Aryal, A., Raubenheimer, D., Subedi, S., Kattle, B., 2010. Spatial habitat overlap and Habitat preference of Himalayan musk deer (Moschus chrysogaster) in Sagarmatha (Mt. Everest) National park, Nepal. Current Research Journal of Biological Science, 2: 217-225.
Aryal, A., Subedi, A., 2011.The conservation and potential habitat of the Himalayan musk deer, Moschus chrysogaster in the protected ares of Nepal. International Journal of Conservation Science, 2: 127-141.
Cantatore, P., Roberti, M., Pesole, G., Ludovico, A., Milella, F., Gadaleta, M. N., Saccone, C., 1994. Evolutionary analysis of cytochrome b sequences in some Perciformes: evidence for a slower rate of evolution than in mammals. Journal of Molecular Evolution, 39: 589–597.
DNPWC, 2011. Anunal Report, Department of National Parks and Wildlife Conservation, Babarmal, Kathmandu, Nepal.
– 2013. Protected Areas of Nepal [in Nepali], Department of National Parks and Wildlife Conservation, Babarmal, Kathmandu, Nepal.
Esposti, D. M., Vries, S. De., Crimi, M., Ghelli, A., Patarnello, T., Meyer, A., 1993. Mitochondrial cytochrome b: evolution and structure of the protein. Biochimica et Biophysica Acta, 1143: 243–271.
Farias, I. P., Orti, G., Sampaio, I., Schneider, H., Meyer, A., 2001. The cytochrome b genes a phylogenetic marker: The limits of resolution for analyzing relationships among Cichlid fishes. Journal of Molecular Evolution, 53: 89-103.
Green, M. J., 1986. The distribution, status and conservation of the Himalayan musk deer Moschus chrysogaster. Biological conservation, 35: 347–375.
Groves, C. P., Yingxiang, W., Grubb, P., 1995. Taxonomy of musk-deer, genus Moschus (Moschidae, Mammalia). Acta Theriologica Sinica, 15: 181–197.
Guha, S., Goyal, S. P., Kashyap, V. K., 2007. Molecular phylogeny of musk deer: a genomic view with mitochondrial 16S rRNA and cytochrome b gene. Molecular Phylogenetics and Evolution, 42: 585–597.
Harrison, R. G., 1989. Animal mitochondrial DNA as genetic marker in population and evolutionary biology. Trends in Ecology and Evolution, 4: 6–11.
Harris, R., 2016. Moschus chrysogaster. The IUCN Red List of Threatened Species 2016: e.T13895A61977139,–1.RLTS.T13895A61977139.en [Accessed on 27 October 2019].
Irwin, D. M., Kocher, T. D., Wilson, A. C., 1991. Evolution of cytochrome b gene in mammals. Molecular Biology and Evolution, 2:13–34.
IUCN, 2013. IUCN Red List of Threatened Species. version 2013.2, [Accessed on 13 July 2014].
Jnwali, S. R., Baral, H. S., Lee, S., Acharya, K. P., Upadhyay, G. P., Pandey, M., Shrestha, R., Joshi, D., Laminchhane, B. R., Griffiths, J., Khatiwada, A. P., Subedi, N., Amin, R. (Compliers), 2011. The Status of Nepal Mammals: The National Red List Series. Department of National Parks and Wildlife Conservation, Kathmandu, Nepal.
Kattel, B., 1992. Ecology of the Himalayan musk deer in Sagarmatha National Park, Nepal. PhD thesis, Colorado State University.
Khatiwada, J. R., Wang, B., Ghimire, B. S., Vasudevan, K., Paudel, S., Jiang, J. P., 2015. A New Species of the Genus Tylototriton (Amphibia: Urodela: Salamandridae) from Eastern Himalaya. Asian Herpetological Research, 6: 245−256.
Kumar, S., Stecher, G., Tamura, K., 2016. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Molecular Biology and Evolution, 33: 1870-1874.
Meyer, A., Wilson, C., 1990. Origin of tetrapods inferred from their mitochondrial DNA affiliation to lungfish. Journal of Molecular Evolution, 31: 359–364.
Pan, T., Wang, H., Hu, C., Sun, Z., Zhu, X., Meng, T., Meng, X., Zhang, B., 2015. Species delimitation in the genus Moschus (Ruminantia: Moschidae) and its high-plateau origin. Plos One, 10: e0134183.
Ramón-Laca, A., Soriano, L., Gleeson, D., Godoy, J. A., 2015. A simple and effective method for obtaining mammal DNA from faeces. Wildlife Biology, 21: 195-204.
Shrestha, B. B., Meng, X., 2014. Spring habitat preferance, association and threats of himalayan musk deer (Moschus leucogaster) in Gaurishankar conservation area, Nepal. International Journal of Conservation Science, 5: 539-540.
Shrestha, B. B., Moe, S. R., 2015. Plant diversity and composition associated with Himalayan musk deer latrine sites. Zoology and Ecology, 25: 295-304.
Silva, T. L., Godinho, R., Castro, D., Abáigar, T., Brito, J. C., Alves, P. C., 2015. Genetic identification of endangered North African ungulates using noninvasive sampling. Molecular Ecology Resources, 15:652-661.
Singh, P. B., Khatiwada, J. R., Saud, P., Jiang, Z., 2019. mtDNA analysis confirms the endangered Kashmir musk deer extends its range to Nepal. Scientific Reports, 9: 4895.
Subedi, A., Aryal, A., Koirala, R. K., Timilsina, Y. P., Meng, X., McKenzie, F., 2012. Habitat ecology of Himalayan musk deer (Moschus chrysogster) in Manaslu Conseration area, Nepal. International Journal of Zoological Research, 8: 81-89.
Su, B., Wang, Y., Lan, H., Wang, W., Zhang, Y., 1999. Phylogenetic study of complete cytochrome b genes in Musk deer (Genus Moschus) using museum samples. Molecular Phylogenetic and Evolution, 12: 241-249.
Su, B., Wang, Y., Wang, Q., 2000. Mitochondrial DNA sequences imply Anhui musk deer a valid species in genus Moschus. Zoological Research, 22: 169–173.
Thompson, J. D., Higgins, D. G., Gibson, T. J., 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Research, 22: 4673–4680.
Timmins, R. J., Duckworth, J. W., 2015. Moschus leucogaster. The IUCN Red List of Threatened Species 2015: e.T13901A61977764,–2.RLTS.T13901A61977764.en [Accessed on 27 October 2019].
Wang, Y., Harris, R., 2015. Moschus fuscus. The IUCN Red List of Threatened Species 2015: e.T13896A61977357,–4.RLTS.T13896A61977357.en [Accessed on 27 October 2019].


Contingut ressenyat a: