Ethnopharmacological studies of indigenous plants in Kel village, Neelum Valley, Azad Kashmir, Pakistan

Study area

Neelum valley is located in the north-east part of Muzaffarabad at 900–6325 m elevation above sea level. It lies between 73°–75° E longitude and 32°–35° N latitude, covering an area of 3737 km² [4]. The Kel is a lush green village, has a hill station, and a tourist spot in the upper part of the Neelum Valley, and lies between 34.8063° N 74.3460° E at an altitude of 2554 m. The climate of the study area is of temperate type where the winter is very cold (average temperature − 2.0 °C) and the summer is pleasant (average temperature 37.0 °C), and the average rainfall is 165 cm annually [18].

The region is characterized by its remoteness, long distance from urban centers, difficult mountainous terrain, and a lack of government services, including modern health care facilities. The area has poorly developed road and other infrastructure. The people of the area rely on sustainable agriculture. Main crops include corn (Zea mays L.), turnip (Brasica rapa L.), and bean (Phaseolus vulgaris L.) in an integrated system. A high proportion of local people are associated with livestock. A number of the main occupations are associated with summer tourism, including rest house managers, tour guides, shop keepers, restaurant workers, and jeep drivers. In light of these demographic changes, it is vital to document the local knowledge of medicinal plant usage in this area before such information declines or is lost completely.

Data collection

Ethnobotanical fieldwork of this quantitative study was conducted during April 2016 to October 2016 following the method of [19]. Standard ethnobotanical methods such as participant observation and open and semi-structured interviews were used to gather the information [20, 21]. All the participants were native of the study area. The purpose, method, and nature of the research were explained before the informants and prior informed consent (PIC) were strictly followed during the field survey. A sum of 20 key informants, comprising of 14 men and 6 women, were selected after the initial survey and many discussions. Demographic information is presented in Table 1 which showed that 20 key informants were selected for the interviews, and out of these, four were between the ages of 30–45, six were ranging 46–60, eight were between the age of 61–75, and two were above 76. Educational status of the informant revealed that there were illiterate (40%), primary and middle (20% each), and secondary and graduate (10% each). Occupation wise, there were housewives (20%), shopkeepers and farmers (15%), labors and teachers (10% each), and local healers (30%). Selection of the informants was made based on their popularity in the study area. They were well known in the area due to their expertise to perform as a medicinal practitioner and had sound information on medicinal plants. In addition to this, information obtained from female informants about the use of indigenous plant in different ailments was compared with the information obtained from the male informants. It was observed that female informants have more knowledge about the utilization of local plants in the preparation and administration of local drugs, which mirror their part in household administration and infection treatment with a specific end goal to keep the family healthy. Meanwhile, their role as a plant collector particularly in rough and steep mountainous tracts of the area was found to be less as compared to men and traditional healers. The same findings were made during the ethnobotanical survey by [22] in the high mountainous region of Chail valley, District Swat, Pakistan, and [23] in Batan Island, Philippines, who additionally found that the ladies assume an imperative part in the arrangement of customary medications utilizing restorative plants. Information was gathered from the informants by using the standard method of [24].

Ethical approval for this study was obtained from the headmen of the studied areas. All respondents were asked to sign a Prior Informed Consent (PIC) form after the objectives and possible consequences of the study had been explained. The PIC form was translated into the Urdu language; however, participants were not subjected to any clinical treatment.

Data preservation

Plant specimens for each species were collected, pressed, dried. Before mounted on the herbarium sheets they were sprayed with the help of a preservative 1% HgCl₂ solution. For identification of the plant specimens, authentication of data, botanical names and families of each plant specimen were confirmed with the help of herbaria comparison, taxonomic literature, manuals and Flora of Pakistan [25]. Whereas the International Plant Name Index (IPNI), Scopus, Web of Science and Google Scholar, Catalog of vascular plants of West Pakistan and Kashmir [26] were also consulted to obtain correct botanical names. Jain and Rao [27] was followed to assign voucher specimens number to each plant and mounting of plants on herbarium sheet was done by standard herbarium techniques [28]. Plants were then deposited at the Herbarium Department of Botany, University of Poonch Rawalakot with voucher numbers.

Quantitative ethnobotanical data analysis

Ethnobotanical data was collected using various quantitative indices including use value (UV), relative frequency of citation (RFC), informant consensus factor (Fic), fidelity level (FL), data matrix ranking (DMR), priority ranking (PR), and Jaccard index (JI). Association between indices was tested using correlation analysis.

UVs

The use value is a quantitative measure of the relative importance of the species. It is used to record the most important plant species in the study area based upon the number of uses cited by the number of people. It was calculated following the standard protocol of [29].\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ \mathrm{UV}=\genfrac{}{}{0pt}{}{\left(\Sigma \mathrm{Ui}\right)}{N} $$\end{document}UV=ΣUiNwhere Ui is the number of use reports cited by each informant for a given species and N is the total number of informants.

RFC

The relative frequency of citation was calculated to assess the incidence of one particular plant species used for the treatment of particular disease or disease category. It was calculated by using following formula as described by [30].

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ \mathrm{RFC}=\frac{\mathrm{FC}}{N} $$\end{document}RFC=FCNwhere FC is the number of informants reporting particular uses of a species and N is the total number of informants.

Fic

Informant consensus factor was used to check the similarity on the informant’s information for each use category and also to check the authenticity of the work by using the following formula:\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ \mathrm{Fic}=\frac{N_{ur}-{N}_t}{N_{ur}-1} $$\end{document}Fic=Nur−NtNur−1where N_ur is the number of use reports in each category, N_t is the number of species used. The factor provides a range of 0 to 1 [31].

FL

The fidelity level (FL) was determined for the most frequently used category by calculating the percentage of informants claiming the use of a certain plant for the same major purpose [28].\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ \mathrm{FL}=\frac{\mathrm{Np}}{N}\times 100 $$\end{document}FL=NpN×100where Np is the number of informants that claim a use of a plant species for a particular use and N is the number of informants that use the plants.

DMR

Data matrix ranking was conducted using eight multi-use plants commonly reported by informants following protocol as described by [32]. Preferentially selected 10 informants were asked to give value to each trait on the basis of benefits obtained from each plant.

PR

Threats to plants of the study area were assessed by obtaining the data from six selected key informants following [20, 28]. This information is beneficial to determine the leading threats to native flora and also helps to suggest the necessary appropriate conservation measures.

JI

To find out the resemblance of indigenous data among various ethnic groups, this study was compared with already published work from the surrounding areas by using the Jaccard index [33].\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ \mathrm{JI}=\frac{c\times 100}{\left(a+b\right)-c} $$\end{document}JI=c×100a+b−c

In this equation, a is the number of species of our area, b is the number of species of the neighboring area, and c is the number of common species in both areas.

Pearson correlation

To find the correlation between RFC and UV values, Pearson correlation was carried out using SPSS version 16.00 Statistical software. The R² was determined to assess species variation in terms of used value.

Documentation of medicinal plants

The present study was conducted in Kel village and its allied areas in the upper parts of Neelum Valley, which is far and remote area of Azad Jammu and Kashmir (Fig. 1). Each species recorded in this study was used in curing different diseases by the local people of the area. Table 1 presents the demographic data of informants while Table 2 provided detailed information of 50 plant species belonging to 33 families used against 13 disease categories. The present work is an effort to enumerate the ethnobotanical data based on which extensively used, high-rated medicinal plant can be selected for searching bioactive compounds to treat ailments. There are only few studies from the area of Kashmir and adjoining areas [1, 2, 18, 34–37]. A total of 50 species were documented to manage miscellaneous human and livestock ailments in the study area. The exercise to utilize the medicinal plants for the welfare of the native people is being used by many rural communities in Pakistan [21, 38].

Fig. 1

Map of the studied area with locations within the upper part of District Neelum

Demographic data

Demographic data showed that there were 70% male and 30% female informants. On the basis of age, the informants were classified into four major groups, i.e., informants of 30–45 years (20%), 46–60 years (30%), 61–75 years (40%), and 76 years and above (10%). During the interview, it was observed that indigenous knowledge regarding the use of medicinal plants was more prevailing in uneducated individuals, i.e., (40%) and similar information was declining with increase in the education level of people with primary and middle level education (20%) and secondary level and graduate (10%). This might be due to the reason that illiterate people have long direct contact with medicinal plants and are well familiar to the usage of these plants. Whereas, highly educated people have little interest in traditional medicine due to their high exposure to the modern education and hence are not involved in learning and practicing ethnobotanical knowledge. Similar results were recorded by [1] in Pakistan, [39] in Ethopia, and [40] in Uganda.

Family contribution

The most dominant family with respect to number of species used was Asteraceae (eight species) followed by Lamiaceae (five species) and Fabaceae (three species). Berberidaceae, Dryopteridaceae, Plantaginaceae, and Rosaceae were represented by two species while all other families shared one species each (Fig. 2).

Fig. 2

Number of species represented by each species

The high number of species presented by these families reflects that the study area has rich biodiversity and the inhabitants have sound knowledge of the usage of these plants in locally prepared medicines. These widely distributed families in different parts of the Kel village and its surrounding areas and their extensive use might be credited to their wide range of bioactive compounds [40]. Asteraceae is largely described to have high number of bioactive compounds therefore contributing to the high use for medicinal purposes [41, 42].

Plant part(s) used and mode of administration

Leaves (30%) and the whole plant (23%) were the most frequently used plant parts as shown in Fig. 3 followed by roots (14%), bark (11%), seeds (8%), fruits (6%) aerial parts and stem (3%) and flowers (2%). Leaves are extensively used in herbal medicines due to the presence of active secondary metabolites [43, 44], and many communities around the world also utilized leaves for the preparation of different herbal products [1, 2]. High proportion of leaves and roots might due to easy availability, more easily extractable phytochemicals, crude drugs, and other mixtures which have proven as valuable regarding phytotherapy [45].

Fig. 3

Plant parts used by traditional healers

Medicinal plants employed in herbal medicines were into nine classes (Fig. 4). The most common methods of preparation of herbal remedies was decoction (20%), extract (17%), infusion (15%), powder and juice (13% each), paste (11%), poultice (7%), and oil (8%). The same results were documented by [2] who reported that most plants were utilized in decoctions, followed by juice and powder. Decoctions are frequently reported as the major forms of preparation in ethnopharmacological studies because they are easy to make [2]. The powder is prepared by crushing plant parts after drying them in shade, and the paste is prepared by grinding fresh or dried plant parts with oil or water [1]. Extraction is also reported as one of the major forms of ethnobotanical practices due to the fact that during extraction, large quantities of components are released which can be used to cure ailment rapidly [46].

Fig. 4

Mode of preparation of medicinal plants

Use values

A numerical method for data inquiry that measures types of uses associated with particular species and relative importance of species in a family (Vendruscolo and Mentz 2006). In our study shown in Table 2, the plant with the highest UVs was Berberis lycium (2.05) followed by Impatiens glandulifera (1.95), Artemisia scoparia (1.75), Ageratum conozoides (1.75), and Achillea millefolium (1.7). Other important species with respect to UVs were Artemisia vulgaris (1.65), Taraxacum officinale (1.45), Inula grandiflora (1.35), Asparagus adsendens (1.25), Cynoglossum lanceolatum (1.2), and Hedera nepalensis (1.1).

High use values of a species show its consistent use for the treatment of various diseases, and high usage reported by the number of informants shows that the plant is well recognized by the local inhabitants and is well utilized in ethomedicinal purposes [1]. The five top-ranked medicinal plants based on UVs were the mostly used plants for various ailments. Among these, B. lyceum was used for jaundice and diarrhea prepared in the form of paste and poultice whereas, I. glandulifera was commonly used for cancer and skin allergies prepared as extract. Infusion of A. scoparia was used to treat hepatitis, jaundice, and stomach disorders. Hot infusion of A. conozoides prepared with leaves or whole plant was used to cure leprosy, fevers, diarrhea, and rheumatism. Leaves of A. millefolium are masticated to cure toothache and juice prepared from leaves is used in earache. Tea is made from the plant as remedy to cure TB, abdominal pain, and fever. It has been observed that those plant which are used over and again are more likely to be biological active [47]. On the other hand, plant with low UVs and RFC values are not essentially insignificant. The low value may reflect the unawareness of people of that area and lack of understanding of proper use of medicinal flora [48].

Relative frequency of citation

Relative frequency citation was calculated to determine the most common occurring medicinal plants used in various disorders as reported by the local informants. In this study, the RFC value ranges from 0.1 to 0.75 as shown in Table 2. The plant with highest RFC was Berberis lycium (0.75). Other significant plants with higher RFC were Cynoglossum lanceolatum (0.65) and Impatiens glandulifera and Achillea millefolium (0.60 each), Ageratum conozoides, and Swertia petiolate (0.55 each).

Plants with high RFC value are very famous among the local people of the area. These plants could be used in research baseline for subsequent assessment of phytochemical profiling, and in future drug discoveries [49]. Such type of plants should be exposed to further pharmacological investigations to appraise the growth of commercial yields [50]. Due to high anthropogenic pressure, such plant species should be given priority for sustainable conservation [2].

Informant consensus factor

The Fic value ranged from 0.71 to 1.00, calculated on the base of use reports against each disease category as shown in Table 3, inferring a high consensus value among participants. The ailments were classified into 13 different categories. Disease categories urinary system, cardiac diseases, baldness, and abortion and miscarriage showed maximum (100%) consensus, since the informants approved of using only one plant species for each category. The most cited plants under these categories Cynoglossum lanceolatum, Prunella vulgaris, Urtica dioica and Viburnum cotinifolium. Other important disease categories were liver disorders and cancerous diseases with Fic values of 0.94 each, and the most cited plants for these categories were Artemisia vulgaris and Taxus wallichiana.

The high ICF values recorded in the Kel village reflected the high dependability of native people on local flora especially for urinary, cardiac, and miscarriage diseases, whereas, low Fic values for ear diseases indicate less consistency of informer’s knowledge [51]. More often a high value of Fic is linked with few species used for one ailment group and has high use value [52], whereas a low value is linked with high number of species presenting high level of disagreement among informants to treat a specific ailment class.

Ailment classification with high Fic values prevailing in the territory may be because of the reason that native people are utilizing water sources on daily basis without having any hygienic or precautionary measurements. Disease of bronchitis, cough, flue, and fever are dominating in the mountainous areas. This is due to the reason that this area has cold and temperate type of climate, and indigenous people sustain their lives by working in forest and fields [21].

Fidelity level

Plants were being used alone or as combined for treat different diseases. Table 3 showed the fidelity level (FL) values for species used for specific ailment as reported by the informants. Plants with higher FL value were more frequently used than those with lower Fl value. In our study, the species having the highest FL value were Berberis lyceum (95%) used in jaundice, hepatitis, typhoid, fever, and tuberculosis disorders, followed by Dioscorea bulbifera, Impatiens glandulifera (90%) used in rheumatism, joint pain, and Artemisia vulgaris (90%) used in liver disorders. Cynoglossum lanceolatum was the plant with the lowest FL value (25%) being used to treat urinary disorders. Plants that have higher number of FL value are considered as model plants that can be employed in further ethnopharmacological research [53].

Direct matrix ranking

Direct matrix ranking (DMR) for tree species showed that Olea ferruginea had the most multipurpose use other than medicinal uses, construction, and cash income (56), firewood (28), fodder (20), fruit and food (20), and hedge and fence (12) (Table 4). Other tree species with high DMR value were Quercus incana and Jaglan regia. The lowest DMR value was recorded for Ailanthus excelsa (Table 5). These result indicated that native, poor people of the area are highly reliant on the tree species for construction, household furniture, fencing, and roof and wall thatching. Similar outcomes were reported by [54] in the East Zone and [55] in the North Zone of Ethopia. The multiple roles of some plants other than their medicinal significance are used in eastern medicine; nevertheless, they play an important role in the routine life of the local people [21]. Livestock raring overgrazing were common activities in the study area that is posing pressure on the flora of the region. Such activities can make the soil more compact resulting in the low seed growth and germination [56].

Overuse of plant resources for construction, fuel wood, and agricultural purposes made them vulnerable in the study area; thus, an urgent conservation action is needed to save the fast eroding tree flora. Necessary steps should be taken by the forest and other concerned departments, and a range land management plan should be chalk out in order to recover the local vegetation in the study area.

PR

The priority ranking data based on the degree of threats to plants was conducted using six key informants (Table 5). In PR, informants ranked timber mafia/exports (21.2%) as a leading threat followed by urbanization (19.35%) and construction (18.71%). Other commonly prevailing threats were agriculture extension (13.55%) and fuel and fodder (10.97%). Illegal export of high-rated species was for illegitimate trading of commercial woods [57]. Several research have been carried out on the aboriginal uses of medicinal plants in Pakistan [15, 34, 58]. Hussain et al. [59] stated that the number of endangered species is increasing day by day because of overgrazing, habitat destruction, and over exploitation of medicinal plants without any conservation plan.

Jaccard index

Ethnobotanical knowledge differs greatly between the indigenous communities because of the difference in cultural and social behaviors. Such kind of comparative analysis among communities is very useful in the sense that they can expose the significant wisdom of information among communities resulting in the findings of new drugs [60].

The results described in this our research were associated with the outcomes of 15 previous studies conducted in different surrounding areas of Pakistan, India, Nepal, and Kashmir (Table 6). The data showed that across 50 species of plants, similarity index of the data ranges from 30.77 to 1.24 whereas the variation fraction varies from 25.64 to 0.62. The maximum level of similarity was found with studies by [4] (95.65), [61] (13.04), [62] (10.57), and [7] (8.26).

The highest degree of similarity index was found in the study of [4], which revealed the same ethnic values and the same type of vegetation and geography of both areas. Furthermore, cultural exchange could have been occurred in the past between the indigenous communities resulting in the similarity in the ethnobotanical outcomes of both areas [2]. Berginia ciliata, Juglan regia, Xanthoxylum armatum, Taxus wallichiana, and Ageratum conozoides from various districts of Central and Western Nepal [63–65] were reported with similar uses. This reflected the similar climatic conditions of the adjacent areas [66] and also high adaptability of these plants to grow in different ecological zones at various elevations [67].

The minimum JI value was recorded for the work conducted by [63] (1.49) in Makawanpur District, Central Nepal [68] and in Hangu, North KPK, Pakistan, and [22] (1.96) in District Bhimber A.K, Pakistan. This difference in ethnobotanical knowledge might be because of the presence of some ecological barrier resulted in geographic isolation of species [1], and the diversities in vegetation and habitats [2]. Beside habitat isolation and difference in vegetation type, researchers have also found that ethnobotanical knowledge vary with age, sex, education level, and origin of the informant [69]. It can be concluded from these studies that geographical isolation among communities has great impact on change in vegetation type and transformation of cultural knowledge, and this might be a cause of the loss of ethnobotanical knowledge.

Pearson correlation analysis

Pearson analysis method was used to measure the strength and direction between UV and RFC variables. Both UV and RFC variables (0.500) had positive relationship of species uses by the number of informants (Table 7, Fig. 5). The coefficient of determination value (R² = 0.25) showed 25% of the difference in use value of species with respect to their citations. This indicated that this work has significant contribution in the documentation of ethnobotanical information on traditional plant uses [22, 57].

Fig. 5

Association between use value and relative frequency of citation

Novelty and future impact

The present study was compared with the previous studies conducted in different parts of the Himalayan territory and other areas as shown in Table 6, to find the novelty index. Medicinal uses of plant species reported in neighboring areas showed more resemblance compared to those documented in other areas. The data collected from the study area reveal considerable difference in plant parts used, mode of herbal preparation, and its utilization as reported from other regions.

Some of the newly documented medicinal uses and relevant plant species include Hedera nepalensis (ulcer), Inula grandiflora (liver pain), Jurinea dolomiaea (bone fracture), Saussurea lanceolata (typhoid), Impatiens glandulifera (joint pain), Betula utilis (Jaundice), Podophyllum emodi (cancer), Dryopteris ramosa (ulcer), Dryopteris stewartii (tuberculosis), Quercus ballota (dysentery), Swertia petiolata (liver pain), Fumaria officinalis (skin allergies), Plectranthus rugosus (skin allergies and diarrhea), Prunella vulgaris (heart diseases), Ailanthus excelsa (fever), Bistorta amplexicaulis (tonic), Rubus ellipticus (wound healing), and Gallium boreole (skin problems). The plant species with new medicinal uses and high RPL value could be studied further to screen bioactive compounds and their pharmacological activities to introduce novel drugs.