Biodiversity Science ›› 2019, Vol. 27 ›› Issue (9): 960-969.doi: 10.17520/biods.2019101

• Special Feature: Snow Leopards: Survey, Research and Conservation Status in China • Previous Article     Next Article

Metabarcoding diet analysis of snow leopards (Panthera uncia) in Wolong National Nature Reserve, Sichuan Province

Qi Lu1, 2, Qiang Hu3, Xiaogang Shi3, Senlong Jin3, Sheng Li1, 2, Meng Yao1, 2, *()   

  1. 1 School of Life Sciences, Peking University, Beijing 100871
    2 Institute of Ecology, Peking University, Beijing 100871
    3 Wolong National Nature Reserve, Wenchuan, Sichuan 623004
  • Received:2019-03-26 Accepted:2019-06-06 Online:2019-09-25
  • Meng Yao

As the apex predator of plateau ecosystems in Central Asia and the Qinghai-Tibet Plateau, the snow leopard (Panthera uncia) plays an essential role in maintaining food-web structure and ecosystem stability. Learning the diet composition and dynamics of the snow leopard is important for understanding its role in ecosystem functioning and interspecific interactions. Previous diet analyses of the snow leopard have been based mainly on morphological identification of food debris in the feces, though the accuracy of this practice has been broadly debated. The Qionglai Mountains are located at the southeast edge of the snow leopard range, harboring a small and relatively isolated population of snow leopards that are barely studied. Using non-invasive sampling, we collected 38 putative snow leopard fecal samples in the Wolong National Nature Reserve in the Qionglai Mountains. To identify the fecal origin, we extracted the fecal DNA and amplified the mitochondrial DNA 16S rRNA gene fragment. Twenty-two fecal samples were identified as originating from snow leopards. Subsequently, vertebrate universal primers and a snow leopard-specific blocking oligo were used to amplify the food components in the fecal DNA, and then high-throughput sequencing was performed to analyze the diet composition of snow leopards. The blue sheep (Pseudois nayaur) was detected in 67% of the samples and was found to be the main staple food of snow leopards’ diet. The domestic yak (Bos grunniens) appeared in 33% of the fecal samples, also accounting for a high proportion of the snow leopard diet. In addition, pikas (Ochotona spp.) and birds were found in a small number of fecal samples. Therefore, wild prey was found to be the main food source for snow leopards in Wolong. However, livestock (yak) also accounted for a relatively large proportion of their diet.

Key words: snow leopard (Panthera uncia), non-invasive sampling, high-throughput sequencing, barcoding, diet analysis, Qionglai Mountains

Table 1

Sequences of the primers used in the study"

引物 Primer 序列 Sequence (5′-3′) 产物长度 Product length 退火温度 Annealing temperature 参考文献 Reference
16S F: GAGAAGACCCTATGGAGC ~380 bp 50℃ Xiong et al, 2016
12SV5 F: TAGAACAGGCTCCTCTAG ~140 bp 60℃ Riaz et al, 2011
- - Shehzad et al, 2012a

Fig. 1

The locations of snow leopard feces collected in Wolong National Nature Reserve (NNR), Sichuan Province. Snow leopard range map is from the IUCN Red List of Threatened Species."

Table 3

A comparison of diet component of snow leopards from different regions"

Study area
No. of
Species identification
Diet analysis
Prey taxa (FO) (%)
Wild ungulates
Unidentified items
Baltistan, Pakistan
49 DNA 形态学
24.5 6.1 4.2 69.3 59.2 32.6 Anwar et al, 2011
South Gobi,
81 DNA DNA 79.0 - 1.2 19.8 - - Shehzad et al, 2012a
Manang, Nepal
213 形态学
51.6 53.5 1.4 17.8 - - Oli et al, 1994
PinValley, Nepal
51 形态学
56.9 3.9 - 23.6 25.5 5.9 Bagchi & Mishra,
Kibber, Nepal
44 形态学
29.6 6.8 15.9 38.5 27.3 19.5 Bagchi & Mishra,
Murghab, Tajikistan
36 DNA 形态学
78.0 45.0 - - - - Kachel et al, 2016
Madiyan, Tajikistan
20 DNA 形态学
50.0 55.0 - 2.0 - - Kachel et al, 2016
Ladakh, India
173 形态学
35.2 25.5 4.6 22.1 60.7 - Chundawat &
Rawat, 1994
Qinghai, China
206 形态学
37.6 44.3 - 6.3 8.0 1.4 Schaller, 1988
Qinghai, China
307 形态学
36.9 58.9 2.3 9.8 2.9 34.2 Liu et al, 2003
Qinghai, China
250 形态学
67.6 17.6 - 16.0 18.0 6.4 Li, 2012
Xinjiang, China
16 DNA 形态学
50.0 112.5 81.3 156.3 25.0 - Wang et al, 2014
Sichuan, China
18 DNA DNA 66.7 11.1 11.1 33.3 - - 本研究 This study

Table 2

Prey species detected in the molecular diet analysis of snow leopard fecal samples (N = 18) in Wolong National Nature Reserve"

Prey taxa
No. of occurrence
FO (%)
RRA (%)
Pseudois nayaur
12 66.7 57.1
Bos grunniens
6 33.3 27.4
高原鼠兔Ochotona curzoniae 1 5.5 0.5
红耳鼠兔Ochotona erythrotis 1 5.5 5.6
鹡鸰科Motacillidae 1 5.5 5.6
雉科 Phasianidae 1 5.5 3.8
总计 Total 100
1 Akrim F, Mahmood T, Max T, Nadeem MS, Qasim S, Andleeb S ( 2018) Assessment of bias in morphological identification of carnivore scats confirmed with molecular scatology in north-eastern Himalayan region of Pakistan. PeerJ, 6, e5262.
2 Alexander JS, Gopalaswamy AM, Shi K, Riordan P ( 2015) Face value: Towards robust estimates of snow leopard densities. PLoS ONE, 10, e0134815.
3 Anwar MB, Jackson R, Nadeem MS, Janecka JE, Hussain S, Beg MA, Muhammad G, Qayyum M ( 2011) Food habits of the snow leopard Panthera uncia (Schreber, 1775) in Baltistan, Northern Pakistan. European Journal of Wildlife Research, 57, 1077-1083.
4 Bagchi S, Mishra C ( 2006) Living with large carnivores: Predation on livestock by the snow leopard (Uncia uncia). Journal of Zoology, 268, 217-224.
5 Bocci A, Lovari S, Khan MZ, Mori E ( 2017) Sympatric snow leopards and Tibetan wolves: Coexistence of large carnivores with human-driven potential competition. European Journal of Wildlife Research, 63(6), 92.
6 Boyer F, Mercier C, Bonin A, Bras YL, Taberlet P, Coissac E ( 2016) OBITOOLS: A unix-inspired software package for DNA metabarcoding. Molecular Ecology Resources, 16, 176-182.
7 Chen LD, Liu XH, Fu BJ ( 1999) Evaluation on giant panda habitat fragmentation in Wolong Nature Reserve. Acta Ecologica Sinica, 19, 291-297. (in Chinese with English abstract)
[ 陈利顶, 刘雪华, 傅伯杰 ( 1999) 卧龙自然保护区大熊猫生境破碎化研究. 生态学报, 19, 291-297.]
8 Chetri M, Odden M, Sharma K, Flagstad O, Wegge P ( 2019) Estimating snow leopard density using fecal DNA in a large landscape in north-central Nepal. Global Ecology and Conservation, 17, e00548.
9 Chundawat RS, Rawat GS ( 1994) Food habits of the snow leopard in Ladakh. In: International Snow Leopard Trust and Northwest Plateau Institute of Biology (eds Fox JL, Jizeng D), pp.127-132. International Snow Leopard Trust, Seattle.
10 Coissac E ( 2012) OligoTag: A program for designing sets of tags for next-generation sequencing of multiplexed samples. Methods in Molecular Biology, 888, 13-31.
11 Deagle BE, Thomas AC, Mcinnes JC, Clarke LJ, Vesterinen EJ, Clare EL, Kartzinel TR, Eveson JP ( 2019) Counting with DNA in metabarcoding studies: How should we convert sequence reads to dietary data? Molecular Ecology, 28, 391-406.
12 Janecka JE, Jackson R, Zhang YG, Li DQ, Munkhtsog B, Buckley BV, Murphy WJ ( 2008) Population monitoring of snow leopards using noninvasive collection of scat samples: A pilot study. Animal Conservation, 11, 401-411.
13 Jumabay-Uulu K, Wegge P, Mishra C, Sharma K ( 2014) Large carnivores and low diversity of optimal prey: A comparison of the diets of snow leopards Panthera uncia and wolves Canis lupus in Sarychat-Ertash Reserve in Kyrgyzstan. Oryx, 48, 529-535.
14 Kachel SM, McCarthy KP, McCarthy TM, Oshurmamadov N ( 2016) Investigating the potential impact of trophy hunting of wild ungulates on snow leopard Panthera uncia conservation in Tajikistan. Oryx, 51, 597-604.
15 Li J ( 2012) Ecology and Conservation Strategy of Snow Leopard (Panthera uncia) in Sanjiangyuan Area on the Tibetan Plateau.PhD dissertation, Peking University, Beijing. (in Chinese with English abstract)
[ 李娟 ( 2012) 青藏高原三江源地区雪豹(Panthera uncia)的生态学研究及保护.博士学位论文, 北京大学, 北京.]
16 Li J, Yin H, Wang DJ, Jiagong Z, Lu Z ( 2013 a) Human-snow leopard conflicts in the Sanjiangyuan Region of the Tibetan Plateau. Biological Conservation, 166, 118-123.
17 Li J, Schaller GB, McCarthy TM, Wang DJ, Jiagong Z, Cai P, Basang L, Lu Z ( 2013 b) A communal sign post of snow leopards (Panthera uncia) and other species on the Tibetan Plateau, China. International Journal of Biodiversity, 2013, 1-8.
18 Li J, Xiao LY, Lu Z ( 2016) Challenges of snow leopard conservation in China. Science China Life Sciences, 59, 637-639.
19 Li S, Wang DJ, Lu Z, McShea WJ ( 2010) Cats living with pandas: The status of wild felids within giant panda range, China. Cat News, 52, 20-23.
20 Liu CG, Zheng SW, Ren JR ( 2003) Research on foods and foods source about snow leopard (Panthera uncia). Journal of Shaanxi Normal University (Natural Science Edition), 31, 154-159. (in Chinese with English abstract)
[ 刘楚光, 郑生武, 任军让 ( 2003) 雪豹的食性与食源调查研究. 陕西师范大学学报(自然科学版), 31, 154-159.]
21 Liu G, Ning Y, Xia XF, Gong MH ( 2018) The application of high-throughput sequencing technologies to wildlife diet analysis. Acta Ecologica Sinica, 38, 3347-3356. (in Chinese with English abstract)
[ 刘刚, 宁宇, 夏晓飞, 龚明昊 ( 2018) 高通量测序技术在野生动物食性分析中的应用. 生态学报, 38, 3347-3356.]
22 Lovari S, Boesi R, Minder I, Mucci N, Randi E, Dematteis A ( 2009) Restoring a keystone predator may endanger a prey species in a human-altered ecosystem: The return of the snow leopard to Sagarmatha National Park. Animal Conservation, 12, 559-570.
23 Ma M, Bariushaa M, Xu F, Mardan T, Yin SJ, Wei SD ( 2005) Markings as indicator of snow leopard in field survey, in Xinjiang. Chinese Journal of Zoology, 40(4), 34-39. (in Chinese with English abstract)
[ 马鸣, Bariushaa M, 徐峰, 买尔旦, 殷守敬, 魏顺德 ( 2005) 新疆雪豹调查中的痕迹分析. 动物学杂志, 40(4), 34-39.]
24 Mallon DP, Jackson RM ( 2017) A downlist is not a demotion: Red List status and reality. Oryx, 51, 605-609.
25 McCarthy T, Mallon D, Sanderson EW, Zahler P, Fisher K ( 2016) What is a snow leopard? Biogeography and status overview. In: Snow Leopards (eds McCarthy T, Mallon D), pp. 23-42 . Academic Press, London.
26 Messmer TA ( 2000) The emergence of human-wildlife conflict management: Turning challenges into opportunities. International Biodeterioration and Biodegradation, 45, 97-102.
27 Myers N, Mittermeier RA, Mittermeier CG, Fonseca GA, Kent J ( 2000) Biodiversity hotspots for conservation priorities. Nature, 403, 853-858.
28 Oli MK, Taylor IR, Rogers ME ( 1994) Snow leopard Panthera uncia predation of livestock: An assessment of local perceptions in the Annapurna Conservation Area, Nepal. Biological Conservation, 68, 63-68.
29 Pei NC, Chen BF ( 2013) DNA barcoding of life: A classification of uses according to function and scale after ten years of development. Biodiversity Science, 21, 616-627. (in Chinese with English abstract)
[ 裴男才, 陈步峰 ( 2013) 生物DNA条形码: 十年发展历程、研究尺度和功能. 生物多样性, 21, 616-627.]
30 Pompanon F, Deagle BE, Symondson WO, Brown DS, Jarman SN, Taberlet P ( 2012) Who is eating what: Diet assessment using next generation sequencing. Molecular Ecology, 21, 1931-1950.
31 Riaz T, Shehzad W, Viari A, Pompanon F, Taberlet P, Coissac E ( 2011) ecoPrimers: Inference of new DNA barcode markers from whole genome sequence analysis. Nucleic Acids Research, 39, e145.
32 Riordan P, Shi K ( 2016) China: Current state of snow leopard conservation in China. In: Snow Leopards (eds McCarthy T, Mallon D), pp. 523-531. Academic Press, London.
33 Ripple WJ, Estes JA, Beschta RL, Wilmers CC, Ritchie EG, Hebblewhite M, Berger J, Elmhagen B, Letnic M, Nelson MP, Schmitz OJ, Smith DW, Wallach AD, Wirsing AJ ( 2014) Status and ecological effects of the world’s largest carnivores. Science, 343, 1241484.
34 Schaller GB, Ren JR, Qiu MJ ( 1988) Status of the snow leopard Panthera uncia in Qinghai and Gansu provinces, China. Biological Conservation, 45, 179-194.
35 Shehzad W, McCarthy TM, Pompanon F, Purevjav L, Coissac E, Riaz T, Taberlet P ( 2012 a) Prey preference of snow leopard (Panthera uncia) in South Gobi, Mongolia. PLoS ONE, 7, e32104.
36 Shehzad W, Riaz T, Nawaz MA, Miquel C, Poillot C, Shah SA, Pompanon F, Coissac E, Taberlet P ( 2012 b) Carnivore diet analysis based on next-generation sequencing: Application to the leopard cat (Prionailurus bengalensis) in Pakistan. Molecular Ecology, 21, 1951-1965.
37 Shi XG, Hu Q, Li JQ, Tang Z, Yang J, Li WJ, Shen XL, Li S ( 2017) Camera-trapping surveys of the mammal and bird diversity in Wolong National Nature Reserve, Sichuan Province. Biodiversity Science, 25, 1131-1136. (in Chinese with English abstract)
[ 施小刚, 胡强, 李佳琦, 唐卓, 杨建, 李文静, 申小莉, 李晟 ( 2017) 利用红外相机调查四川卧龙国家级自然保护区鸟兽多样性. 生物多样性, 25, 1131-1136.]
38 Shrestha B, Aihartza J, Kindlmann P ( 2018) Diet and prey selection by snow leopards in the Nepalese Himalayas. PLoS ONE, 13, e0206310.
39 Taberlet P, Coissac E, Pompanon F, Brochmann C, Willerslev E ( 2012) Towards next-generation biodiversity assessment using DNA metabarcoding. Molecular Ecology, 21, 2045-2050.
40 Tang Z, Yang J, Liu XH, Wang PY, Li ZY ( 2017) Research on snow leopards (Panthera uncia) using camera-trapping in Wolong National Nature Reserve, China. Biodiversity Science, 25, 62-70. (in Chinese with English abstract)
[ 唐卓, 杨建, 刘雪华, 王鹏彦, 李周园 ( 2017) 基于红外相机技术对四川卧龙国家级自然保护区雪豹(Panthera uncia)的研究. 生物多样性, 25, 62-70.]
41 Tseng ZJ, Wang X, Slater GJ, Takeuchi GT, Li Q, Liu J, Xie G ( 2014) Himalayan fossils of the oldest known pantherine establish ancient origin of big cats. Proceedings of the Royal Society B: Biological Sciences, 281, 20132686.
42 Valentini A, Pompanon F, Taberlet P ( 2009) DNA barcoding for ecologists. Trends in Ecology & Evolution, 24, 110-117.
43 Vallet D, Petit EJ, Gatti S, Levréro F, Ménard N ( 2007) A new 2CTAB/PCI method improves DNA amplification success from faeces of Mediterranean (Barbary macaques) and tropical (lowland gorillas) primates. Conservation Genetics, 9, 677-680.
44 Wang J, Laguardia A, Damerell PJ, Riordan P, Shi K ( 2014) Dietary overlap of snow leopard and other carnivores in the Pamirs of Northwestern China. Chinese Science Bulletin, 59, 3162-3168.
45 Xiong MY, Wang DJ, Bu HL, Shao XN, Zhang D, Li S, Wang RJ, Yao M ( 2017) Molecular dietary analysis of two sympatric felids in the mountains of Southwest China biodiversity hotspot and conservation implications. Scientific Reports, 7, 41909.
46 Xiong MY, Shao XN, Long Y, Bu HL, Zhang D, Wang DJ, Li S, Wang RJ, Yao M ( 2016) Molecular analysis of vertebrates and plants in scats of leopard cats (Prionailurus bengalensis) in southwest China. Journal of Mammalogy, 97, 1054-1064.
47 Zhou YY, Feng JC, Duo HR, Yang HL, Li J, Li DQ, Zhang YG ( 2014) Population survey and genetic diversity of snow leopards Panthera uncia as revealed by fecal DNA. Acta Theriologica Sinica, 34, 138-148. (in Chinese with English abstract)
[ 周芸芸, 冯金朝, 朵海瑞, 杨海龙, 李娟, 李迪强, 张于光 ( 2014) 基于粪便DNA的青藏高原雪豹种群调查和遗传多样性分析. 兽类学报, 34, 138-148.]
[1] Gui-Feng GAO. (2020) Techniques and methods of microbiomics and their applications . Chin J Plant Ecol, 44(生态技术与方法专辑): 0-0.
[2] Ben-Feng HAN Xin Zhou Xue Zhang. (2020) Verification of virus identity and host association using genomics technologies . Biodiv Sci, 28(5): 0-0.
[3] Shao Xinning, Song Dazhao, Huang Qiaowen, Li Sheng, Yao Meng. (2019) Fast surveys and molecular diet analysis of carnivores based on fecal DNA and metabarcoding . Biodiv Sci, 27(5): 543-556.
[4] Li Hanxi, Huang Xuena, Li Shiguo, Zhan Aibin. (2019) Environmental DNA (eDNA)-metabarcoding-based early monitoring and warning for invasive species in aquatic ecosystems . Biodiv Sci, 27(5): 491-504.
[5] Liu Shanlin. (2019) DNA barcoding and emerging reference construction and data analysis technologies . Biodiv Sci, 27(5): 526-533.
[6] Li Meng, Wei Tingting, Shi Boyang, Hao Xiyang, Xu Haigen, Sun Hongying. (2019) Biodiversity monitoring of freshwater benthic macroinvertebrates using environmental DNA . Biodiv Sci, 27(5): 480-490.
[7] Zhang Xue, Li Xing’an, Su Qinzhi, Cao Qina, Li Chenyi, Niu Qingsheng, Zheng Hao. (2019) A curated 16S rRNA reference database for the classification of honeybee and bumblebee gut microbiota . Biodiv Sci, 27(5): 557-566.
[8] Chen Zhixiang, Yao Xueying, Stephen R. Downie, Wang Qizhi. (2019) Assembling and analysis of Sanicula orthacantha chloroplast genome . Biodiv Sci, 27(4): 366-372.
[9] Hu Jianlin,Liu Zhifang,Ci Xiuqin,Li Jie. (2019) Use of DNA Barcoding in Identifying Tropical Trees from Dipterocarpaceae . Chin Bull Bot, 54(3): 350-359.
[10] JIANG Yu-Ling, CHEN Xu-Hui, MIAO Qing, QU Bo. (2019) Difference in fungal communities between in roots and in root-associated soil of nine orchids in Liaoning, China . Chin J Plant Ecol, 43(12): 1079-1090.
[11] Dandan Lang,Min Tang,Xin Zhou. (2018) Qualitative and quantitative molecular construction of plant-pollinator network: Application and prospective . Biodiv Sci, 26(5): 445-456.
[12] Hou Qinxi, Ci Xiuqin, Liu Zhifang, Xu Wumei, Li Jie. (2018) Assessment of the evolutionary history of Lauraceae in Xishuangbanna National Nature Reserve using DNA barcoding . Biodiv Sci, 26(3): 217-228.
[13] Xiaojuan Deng, Jianli Liu, Xingfu Yan, Peigui Liu. (2018) Community composition of bacteria associated with ascocarps of Tuber indicum using traditional culture method and Roche 454 high-throughput sequencing . Biodiv Sci, 26(12): 1318-1324.
[14] Qiaoli Li, Na Yan, Qiong Song, Junzhan Guo. (2018) Complete Chloroplast Genome Sequence and Characteristics Analysis of Morus multicaulis . Chin Bull Bot, 53(1): 94-103.
[15] Jinfeng Hao, Xiaohong Zhang, Yusong Wang, Jinlin Liu, Yongchao Zhi, Xinjiang Li. (2017) Diversity investigation and application of DNA barcoding of Acridoidea from Baiyangdian Wetland . Biodiv Sci, 25(4): 409-417.
Full text