Food Preference of Rhesus Monkey ( Macaca mulatta) in the Margalla Hills National Park, Islamabad, Pakistan

a


O n l i n e F i r s t A r t i c l e
In non-human primates, Rhesus monkey shows the widest range of geographical and ecological distribution in the world. It is found naturally in Afghanistan, India, Pakistan, China, Laos, Bangladesh, Bhutan, Burma, Nepal, Thailand, and Vietnam. This species is found in a variety of habitats throughout its range, including urban areas (Fleagle, 1988). The Rhesus monkey is one of the two primate species found in Pakistan, the other one is Grey Langur (Semnopithecus entellus) distributed in district Mansehra, Pallas Kohistan, and Azad Jammu and Kashmir (AJ K) (Roberts, 1997). In Pakistan, the distribution of rhesus monkey is restricted to mountainous areas having forest cover, typically associated with Himalayan moist temperate forest. It extends in the northwest from Kafir valleys of southern Chitral, southwards though Dir and eastwards through Swat/Kohistan and Hazara district (Roberts, 1997). The subspecies Macaca mulatta villosa occurs in lower Kaghan valley around Paras and Shogran and also in Neelam valley of AJK. The same sub species occurs throughout the Murree hills. Individuals from the outer foothill of Himalayan range such as Margalla Hills appear to be much smaller in size (Roberts, 1997). No doubt, it is much more widespread in distribution than langur but Rhesus monkey has become very rare both in Swat Kohistan and Chitral due to persecution by local villagers (Roberts, 1997). Currently, Rhesus monkey is considered Near Threatened in Pakistan (Sheikh and Molur, 2005). Globally it has got the status of 'LC' (Least Concern) (Timmins et al., 2015). In all habitat types, feeding and resting are the major activities of the rhesus macaques' day and they spend the rest of their time traveling, grooming, playing, and other activities (Seth and Seth, 1986). Rhesus monkeys dominate in the tropical, subtropical, and temperate forests below 3,000 m a.s.l. all over Nepal. Assamese monkeys (Macaca assamensis) were patchily distributed along rivers in the tropical and subtropical areas. Both species principally utilized forests parapatrically. Discontinuous distribution of Assamese monkeys probably appeared because of the expansion of rhesus monkey distribution in the mid-and late-Pleistocene (Wada, 2005). They are synanthropic, thriving in human-altered environments, including urban of some communities. This adaptive characteristic as evolutionary strategy has made Rhesus macaques the most widely distributed and successful primates in the world (Hasan et al., 2013). Habitat use and positional behavior by sympatric Rhesus macaques and Aasamees macaques were interspecific in limestone habitat at Nonggang Nature Reserve, southwestern Guangxi, China. These differences in positional behavior and habitat use are linked to differences in limb length and body size in the two species, and may be explained by the spatial distribution of preferred foods and structure of the forest in different areas of the limestone hills (Huang et al., 2015). Rhesus monkeys are diurnal in feeding activity and forage to a considerable extent on the ground. Studies in India indicate that they usually move regularly along a rough circuit within a particular territory, feeding as they go and sleeping in a fresh place each night (Roberts, 1997). It is observed that if they are undisturbed and food supplies are plentiful they will feed in the same area for several consecutive days, returning at night to the same group of trees. In the spring and early summer, they eat a lot of grasses and forbs and have been observed greedily pulling up handfuls (Roberts, 1997).
Rhesus monkey is largely vegetarian. Its diet includes leaves, flowers, fruits, berries and seeds of many species of plants, grass, grains and algae from ponds. It also eats insects, spiders and is known to eat small birds, lizards, or similar small animals. Macaque eats fruits, berries, leaves, flowers, seeds and bark from over 70 species of plants (Lindburg, 1977). Insects (termites, grass-hoppers, ants and beetles) and occasionally honey combs of wild bees are eaten when available. Sometimes during the sunny part of the monsoon season a group spends the whole day catching and eating grasshoppers (Lindburg, 1977). There is seasonal variability in the consumption of more important food plants while foraging, the monkeys quickly fill their cheek pouches with food, especially small fruits and berries, and then sit in a safe place and in an alert position while continuing to chew (Lindburg, 1977).
No study has been done on foraging and food composition of Rhesus monkey in the past in Margalla Hills National Park (MHNP), Islamabad. The present study was conducted on important aspect of Rhesus monkey i.e. food composition and preference in Margalla Hills National Park Islamabad. It provided basic data on food composition and preference of Rhesus monkey in various seasons, which will be helpful in habitat management and improvement of this species in the park.

Study sites
The study was conducted in the MHNP, covering an area of 14,786 hectares which runs along the northeast border of the federal capital city of Islamabad (Fig. 1). The Margalla Hills are one of the western most extensions of Indo-Himalayan ecosystem. The hills represent a contact zone with the arid Irano-Saharan ecosystem, which extends southwesterly. The MHNP hosts 250 species of birds, 38 O n l i n e

F i r s t A r t i c l e
of mammals, at least 13 taxa of reptiles and numerous taxa of other animals (GOP, 1992). Topography is rugged and elevation ranges from 456 m to 1580 m. The general aspect is southerly and terrain is interspersed with both large and small valleys. The rocks have been observed to date back to Jurassic and Triassic ages. Soils are dark, with high mineral contents and are capable of supporting good tree growth despite being shallow (Anwar and Chapman, 2000). The climate is sub-tropical semi-arid. The region lies in the monsoon belt and experiences two rainy seasons. Winter season lasts from January to March and summer season from July to September with average rainfall 1900 mm. There have been occasional incidents of light snowfall in severe winters. Average minimum and maximum temperature is -3.9 o C and 46.6 o C, respectively (PMD, 2018) (https://nwfc.pmd.gov.pk/new/monthlyreports.php). The flora of the park is mainly subtropical evergreen scrub forest on the lower slopes and subtropical chir pine forest at higher elevations. The vegetation has been classified into five major phytocommunities on the basis of physiognomy, floristic composition and dominance of the vegetation. These include: Olea-Acacia, Acacia-Carissa, Olea-Carissa, Myrisine-Dodonea and Pinus-Quercus community (Anwar and Chapman, 2000).
There are over 30 settlements in and immediately around the periphery of the National Park. The total population in the settlements in national park is close to 70,000. The largest settlement is in Nurpur with a population of just over 15,000 persons, followed by Chauntra with 12,000, and Kot Hathial, Gokena, Talhar and Shah Allah Ditta with about 9,000 persons each. Major threats to the national park are encroachments, woodcutting and grazing, wild fires, alien invasive plant species, quarries and increasing number of visitors (IWMB, 2019).

Field observations for food preference
As a first step, a reconnaissance survey of study area was conducted to know the potential feeding sites of the monkeys and to divide the area into five study sites (Kuwait-hostel, E-7 hill side road, wildlife view point, Daman-e-koh and Lohe dunde) on the basis of distribution of monkey groups and different food resources. Five monkey groups were named as group A, group B, group C, group D and group E according to their potential feeding sites, Kuwait-hostel, E-7 hill side road, wildlife view point, Daman-e-koh and Lohe dunde, respectively. The second step involved the detailed analytical study on food composition of rhesus monkey using various techniques given as follows: Food composition of the species based on field observation was studied by continuous-recording sampling method and instantaneous and scan sampling by employing two observation methods (Altmann, 1974). Food remnants method was used in combination with fecal analysis and visual observation. It is useful when the animal was not directly visible but the researcher was close enough to identify the species on which it was feeding. When the animal had moved on, the feeding location was investigated and collected food remnants were characterized (Ullrey et al., 2003).

Sample collection
A total of 150 fecal pellets were collected from MHNP in winter and summer season. None of the fecal pellets from which samples were taken was judged to be more than 24 h old. Fecal analysis technique was used by collecting the samples of fecal matter of the species from various sites of study area. Both macroscopic and microscopic techniques were used to find out the indigestible hard parts such as seeds and fibrous material of ingested food (Ullrey et al., 2003).

Slide preparation
In different seasons, plant species from different study area were collected. Reference slides were prepared following William (1962) and Ward (1970). The required vegetative parts of the plants were obtained and dried. These fresh specimens/dried tissues are soaked in plant soaking solutions (distilled water, ethyl alcohol, and

F i r s t A r t i c l e
glycerin (1:1:1)) for a night then washed with tap water for about 10-20 min each specimen of plant tissue was ground in virtis homogenizer with distilled water. These contents were poured in microsieve. This microsieve was composed of 6cm long hollow cylinder having 0.05mm mesh of stainless steel wire that is fitted with a rubber stopper at one end of cylinder in such a way that it could be left filled with 1% sodium hypochlorite for clearing the specimen and was kept soaked in a sodium hypochlorite solution of 5% chlorax and 4 parts of distilled water (1:4) for 20-30 min. To neutralize the basic effect of sodium hypochlorite equal amount of dilute acetic acid was added to the tissues, were placed in mordant solution for 15-30 min, and then this distilled water was dripped into the sieve to remove any basic residues. The contents were placed in hematoxylin stain for 10-15 min then washed with tap water. On a clean slide a drop of Apathy of mounting medium (100cc distilled water and 100g gum Arabic) was placed. The stained plant material was mixed with this mounting medium with a wet camel brush and the material was uniformly spread over 22x40 mm of slide. Two drops of mounting medium were added to the plant material and were covered with glass cover of 22x40 mm and pressed tightly by a peril eraser for uniform contact of glass cover and slide. Labeling of slides was done for identification and was left at room temperature overnight for fastening of material on these slides.

Slide interpretation
The segments in each fecal slide were identified on the basis of shape, size and relationship of short cells, structure of cell wall of long cells, stomatal patterns along with guard cells, nature and arrangement of epidermal hair. Two observers, trained in the procedures of Holechek and Gross (1982), analyzed both diet and fecal samples using AmScope B490 binocular compound microscope. Samples were analyzed at 100X, although 200X magnification was sometimes used for higher resolution (Holechek and Valdez, 1985).

Food composition
Dietary plant species found in the fecal sample were confirmed after a detailed examination of all cell features by making comparison with the reference plant microphotographic key. The relative frequency of a plant species in the fecal samples was determined and shown as a relative importance value (RIV) (Jnawali, 1995).
The relative importance value was the total number of segments identified for a particular dietary plant species divided by the total number of all the counts made in the sample, multiplied with 100 (Jnawali, 1995). The data on provisioning feeding and scavenging were collected directly by observations and questionnaires were used to interview the visiting tourists. The field study techniques of focal animal sampling and all animals sampling (Altmann, 1974) were used to find out feeding bout, feeding rate and feeding behavior.

Statistical analysis
The data obtained on foraging and food composition were subjected to statistical analysis using Estimate S (8.20) and one-way ANOVA test (Steel et al., 1997).

Food composition based on field observations
Food composition based on flora Rhesus monkeys exploited a wide range of food varieties, but they concentrated on small number of selected plant species. They consumed 30 plant species in the study area (Table I). Their diet included leaves, fruit, flowers, buds, seeds, petioles and wood extract. Fruits were the most important energy rich components of diet but leaves, seeds and flowers were important too. Fruits generally contain relatively large quantity of simple sugars and are readily useable source of energy. Leaves were the major parts of plant consumed overall which probably satisfied certain nutritional requirements of the Rhesus monkey. They prefer young leaves of plants over mature ones which contain fairly high concentration of calcium (Biddulp, 1959) while young leaves have high percentage of raw proteins (Struhsaker and Oates, 1975).
Food composition determined on the basis of field observations showed seasonal variations linked to current vegetation structure of the locality and abundance of seasonal plants in study area. Groups A, B and C shared almost the same feeding ecology due to the same vegetative structure in those localities.   In feeding range of each particular group, food composition also varied seasonally. Group A, group B and group C consumed 15, 11 and 13 plant species, respectively. In general, young leaves, buds and mature fruits were more commonly consumed than flowers and mature leaves. This is due to nutritional requirements in terms of higher protein contents and lower cellulose levels found in these plant parts compared to mature leaves (Hladik, 1977). Su and Lee (2001) reported that food composition of Formosan monkeys (Macaca cyclopis) consisted of 56 plant species including fruits, stems, leaves, flowers and shoots. There was seasonal variation in the proportion of each food type in the diet of the macaques. Fruits comprised the highest proportion (53.8%), flowers (7.32%), leaves (14.92%), stems (11.76%), shoots (2.44%) and insects (9.76%) in the diet of Rhesus monkey. The proportion of fruits was higher mainly in summer and they mostly consumed the fruits of two plant species Trema orientalis and Ficus irisana. their proportion was high during the winter. Insects were mainly consumed during summer (July and August) due to insect abundance during rainy (July-August) season and their feeding reduced in September and October due to shortage of insects and severe interspecific competition with insectivorous birds. Hence, the macaques switched on to plant matter in the fall and winter seasons. Food composition of red howler monkey (Allouatta seniculus) by field observations showed that they fed on 195 plant species from 45 families and major feeding categories were mature fruits (21.5%), young leaves (54%) and flowers (12.6%). Other feeding categories were old leaves, unripe fruits, bark and moss and termitarium soil. The monkeys were less selective in feeding than other Howler monkey groups Alouatta spp. because 19 plant species accounted only for one percent of their food composition and showed seasonal diet variations (Julliot and Sabatier, 1993).

Provisioning
Rhesus monkeys consumed over a dozen kinds of provisioned food items in the study area. Major food items ( Fig. 2A) included junk foods (slanties, lays, biscuits etc.) (20.4%) and white bread slices (13.6%) followed by Chapatties (3.4%) and fruits > (4%). Consumption rates of various food items varied and depended on size and texture of provisioned food items (Fig. 2B). Present study showed that there is no pronounced variation in food items offered by the tourists to Rhesus monkeys in the park except in seasonal fruits which tourists carried with them. Ciani and Chiarelli (1988) studied qualitative and quantitative parameters of foraging in macaques. Assessment of ecological preferences and individual foraging strategies in food gathering was done in the forest of Simla, India. Their studies showed that monkeys had short bout duration indicating high consumption rate for provisioned and town food resources compared to free ranging monkeys. Foraging of different age and sex classes show pronounced discriminations. Adult males spent much less time for feeding than other classes and showed a preference for urban food resources. Contrary to adult males, lactating mothers in adult females were sluggish and giving extra time to foraging. Sub-adults showed the intermediate characteristics of foraging compared to adult females and males. Infants showed high dependency on lactation and low foraging to town food resources. The quality and exploitation of food patches greatly influenced the foraging in the study area. These monkeys showed competition in foraging; particularly males showed agnostic feeding for access to nutritious town food resources.

Scavenging
A total of 19 types of food items were found in scavenging food composition of the monkeys in the study

Food Preference of Rhesus Monkey (Macaca mulatta)
sites (Table II) which were qualitatively analyzed as major and minor food components based on the availability in the study sites. Of these food items, 15 were major while remaining 4 types were minor food items. The plant based and cooked expired food stuff were major food items while tetra packs, canned food material was minor in the scavenging food composition (Table II). The solid waste bins from CDA were the major source of the food in the human surrounding/ intervened sites (Kuwait Hostel, E-7 sector, Wildlife View point, Daman-e-Koh) in the study area.

Food composition based on fecal analysis Mechanical fecal analysis
The mechanical fecal analysis of 150 fecal pellets (Fig. 3) of Rhesus monkey showed presence of indigestible hard parts including seeds, seed coats, blades of forbs and soil components (Table III). The seeds and seed coats of different plant species were identified as they slightly changed in color, texture and shape during digestion (Fig.  4). The seeds of 7 plant species were found including guava (Psidium guajava) (8.67%), oranges (Citrus sinensis) (5.34%), pomegranate (Punica granatum) (2.67%), tomatoes (Solanum lycopersicum) (4%), pears (Pyrus pashia) (2%), Kinjal (Cassia fistula) (12%) and Ficus bengalensis, Ficus carica and Ficus religiosa (7.34%). The seed coats found in fecal matter were pieces of ground nuts (10.67%), tomatoes (4%), Kinjal fruits (3%) and grams (9%). The results showed occurrence of high percent frequency of seed and seed coats consumed by monkeys (Table III). Fruits and seeds not observed during the field observations were also found as major components of their diet.
Blades of grasses and forbs found in the analyzed fecal matter included Elusine indica (32%), Asparagus officinalis (2.6%) and Cynodon dactylon (27.67%). Pebbles (18%) were the only soil components found in fecal matter but their presence in the feces did not support the geophagic (soil eating) habit of Rhesus monkeys in the study area as no field observation regarding this matter was recorded by the researchers.

Micro-histological fecal analysis
The results of micro-histological fecal analysis ( Supplementary Fig. 1) showed significant differences in relative importance of preferred food plant species in feces of Rhesus monkey. The percent occurrence of the plant species in fecal samples revealed diet selective pressure in habitat of the monkeys of each feeding range in the study area and showed that during the spring and summer seasons monkeys exploited diverse food resources and all feeding categories (Table IV). The presence of cells and cell fragments (Microhistolgical reference fecal slide plates) of Carissa opaca 136 (90.67%), Cassia fistula 29 (90.67%), Cynodon dactylon 54 (36%), Dalbergia sissoo 113 (75.34%) and Ficus bengalensis 32 (21.34%) in the slides showed their relative importance as dominant plant species in the seasonal diets.

Plant species identified from fecal samples during autumn and winter seasons
Rhesus monkey showed restricted foraging preference during autumn and winter season (Storr, 1961). Relatively small proportion of plant species in fecal samples showed that during autumn and winter seasons the monkeys exploited restricted food resources and selective feeding categories such as mainly leaves (Table IV). The presence of the cells and cell fragments of Carissa opaca 119 (79%), Cassia fistula 24 (16%), Cynodon dactylon 34 (22.67%), Dalbergia sissoo 113 (75.34%) and Ficus bengalensis 44 (29.34%) in the slides showed their relative importance in the seasonal diets (Table IV)

Estimation of species richness
A comparison of species observed mean (runs) in fecal samples with other species richness estimators such as ACE mean, ICE Mean, Chao 1 Mean, Chao 2 Mean and Jack 1 Mean using Estimate S.8.2.0 was made for data analysis (Fig. 5). Total plant species observed in the micro-histological fecal analysis were 29 as compared to 30 observed in field observations consumed by Rhesus monkey in study area. ICE Mean shows the highest possible food plant species in the study area would be 41 but Jack 1 mean confirmed the total observed species were 29 in more than 500 fecal samples. The microscope technique can only be a useful tool if observers have a 90% or more recognition level of the plant species being examined, and if they become careful not to overestimate species with stellate trichomes or hairs. Identification of such species should be based only on those trichomes attached to recognizable epidermal tissues and cell pattern and/ or stomata pattern on peridermal tissues (Alipayo et al., 1992). The overall food composition of Rhesus monkey in the study area consists of 83% plants material,

O n l i n e F i r s t A r t i c l e
S. Aslam et al.
14% provisioned material and 03% scavenged diet (Fig.  6). Su and Lee (2001) while studying feeding habits of Formosan rock macaque (Macaca cyclopis) through fecal analysis and field observations found that monkeys were mainly frugivorous and fed on 51 plants species and insects of five orders but differences were found in the results of both of these methods. They analyzed 101 fecal samples and found the frequencies of occurrence of fruits (96%), plant bodies (98%), flowers (3%), and animal matter (82.3% ). They found 46 species of fruits by seeds in the feces and remains of partially digested leaves and stem of one plant species. Lindburg (1977) reported the food composition of Macaca mulatta in Siwalik forest in north India based on fecal analysis which showed 65-70 percent frugivorous food habit. Su and Lee (2001) reported the relative importance of fruits (46.2%), plant bodies (47.6%), and animal matter (6.2%) respectively in fecal samples of Formosan rock macaque. Flowers were found in traces and relative importance in fecal samples varied seasonally (Sanders et al., 1980). Fruits and animal matter increased during spring and summer while plants in autumn and winter.

Food resources preferences
The food resources exploited by different age/sex classes were significantly different as the 'p' value is greater than the significance level (Table V). The food resources exploited by different age/sex classes were calculated on the basis of total number of 5 min interval scan samples for each food resource and the percentage value were estimated (Table VI). Food resource preference was mainly based on availability and relative abundance of the food resources in the study area (Hanya et al., 2011). Adult females (young leaves 52%, mature leaves 5%, fruits 16%, seeds 15% buds 9% and flowers 3%) and juveniles (young leaves 52%, mature leaves 4%, fruits 13%, seeds 11%, buds 14% and flowers 6%) were more selective in feeding and would climb to the highest canopies of Carissa opaca and other food plant species to gather young buds. It could be possible that adult females (Fig. 7A) and juveniles (Fig. 7B) preferred young leaves/ buds and were using this resource as an additional source of protein. Lactating females were found to consume seeds. Seeds are extremely high in protein and fatty acids (Heller et al., 2002), and thus, seeds likely represent adult females' means of obtaining adequate levels of these nutritional components.  Agetsuma and Nakagawa (2007) studied the foraging behaviour of Japanese free ranging monkey (Macaca fuscata) in two different climatic habitats: a cool temperate habitat (Kinkazan Island) and a warm temperate habitat (Yakushima Island). Both habitat types differ significantly in diet composition and daily activity patterns. Time spent feeding on Kinkazan Island was 1.7 times that on Yakushima Island. They concluded that two factors were responsible for these: (1) the energy required for thermoregulation of monkeys on Kinkazan Island is greater than that on Yakushima Island; and (2) the food quality, which affects the intake speed of available energy, is lower on Kinkazan Island. However, monkeys in both habitats increased their moving time and decreased their feeding time when they fed on foods of relatively high quality.
Analysis of food resources preference of five age/ sex classes of Rhesus monkeys revealed that it mainly O n l i n e

F i r s t A r t i c l e
varied due to nutritional requirement and physiological condition of monkeys (Kumar and Solanki, 2003). Adult females and juveniles were more selective in feeding and would climb to the highest canopies of food plant species to gather young buds, males were rarely engaged in this behaviour. The food resources preference, social dominance and spatial distribution of particular food plant species were factors which gave rise to feeding agnostics in the monkeys. The analysis of food resources preference revealed that it also mainly varied showing higher food resource preference particularly during the mating season as did the other age/ sex classes of Rhesus monkey (Hanya and Chapman, 2013). Adult males showed low preference for young leaves over mature leaves but showed almost equal preference for the fruits as other groups (Fig. 7C). Ciani (1986) reported that males showed high inter-troop feeding aggressions in the town where they were competing for highly valuable food resources. The sub-adults showed almost similar food resources preference as the juveniles. They showed higher preference for young leaves over mature leaves but fruit and seeds were given nearly the same preference as did the adult males (Fig. 7D). Rhesus monkeys preferred succulent foods over the non-succulents. Most of the plants consumed by Rhesus monkeys were succulents as these plants store excess water in their roots, leaves and shoots (Sajeva and Costanzo, 1997), providing the macaques an excellent source of water and additional nutrients. Adult males, juvenile males and, to some extent, juvenile females engaged vigorously in feeding activity than adult females, making it necessary to fulfill water needs through vegetation during foraging bouts and thus likely that they would more frequently exploit this resource.

CONCLUSIONS
Rhesus monkey is an omnivore species in the Margalla Hills National Park. It actively selected feeding sites having a human influence in the form of a tourist resort, residential areas or garbage dumps which helped them in acquiring food easily. The major natural plant source dietary items in their diet were Carissa opaca, Cynodon dactylon, Cassia fistula, Dalbergia sissoo, Elusine indica and Ficus bengalensis. Rhesus monkey showed preference for fruits on leaves in their diet. Rhesus monkey showed competition in the provisioned feeding. Provision of food (14%) by tourists accounted an important source in the food composition of Rhesus monkey in the study area. Rhesus monkey regularly visited garbage dumps daily and scavenged on selected food remnants, as a part of their food. Presence of water bodies in their foraging ranges was an important factor but scarcity of water in their feeding ranges was also noted.

RECOMMENDATIONS
i. Extraction/ cutting of natural vegetation in some areas of Rhesus monkey habitat in the national park must be controlled; this is the leading cause of habitat as well as natural foods destruction in the study area. ii. Protection and management of plant species consumed by Rhesus monkey should be done on priority basis by concerned authorities. This will give increased food availability and reduce competition for them. iii. There is immediate need to build up water sources such as ponds for drinking and bathing purposes for Rhesus monkey by concerned authorities. iv. Provisioning feeding by tourists should be minimized or stopped as this poses serious risks for transmission of contagious diseases from humans to monkeys and thus threatens their survival. v. Serious and immediate steps should be taken for the regular collection of garbage in the study areas which attract Rhesus monkey for scavenging that may lead to outbreak of diseases in monkeys and thus put their survival at risk. vi. There is need to conduct research on feeding ecology, food preference based on nutritional requirements of Rhesus monkey in MHNP.

O n l i n e F i r s t A r t i c l e ACKNOWLEDGMENT
The authors acknowledge the Islamabad Wild-life Management Board (IWMB) for granting permission for data collection in Margalla Hills National Park (MHNP) Islamabad. Last but not least, we highly appreciate Dr. Maqsood Anwar (late), Department of Wildlife Management at Arid Agriculture University for initial support in conceiving the idea and planning of this study.

Funding
This study was self sponsored by the researchers. No funding was provided by any organization for this study.

Ethical statement
No animals and plants were harmed in the data collection from the field.

Supplementary material
There is supplementary material associated with this article. Access the material online at: https://dx.doi. org/10.17582/journal.pjz/20210724200757

Statement of conflict of interest
The authors have declared no conflict of interest.