Diversity of Medicinal Plants among Different Forest-use Types of the Pakistani Himalaya

Diversity of Medicinal Plants among Different Forest-use Types of the Pakistani Himalaya Medicinal plants collected in Himalayan forests play a vital role in the livelihoods of regional rural societies and are also increasingly recognized at the international level. However, these forests are being heavily transformed by logging. Here we ask how forest transformation influences the diversity and composition of medicinal plants in northwestern Pakistan, where we studied old-growth forests, forests degraded by logging, and regrowth forests. First, an approximate map indicating these forest types was established and then 15 study plots per forest type were randomly selected. We found a total of 59 medicinal plant species consisting of herbs and ferns, most of which occurred in the old-growth forest. Species number was lowest in forest degraded by logging and intermediate in regrowth forest. The most valuable economic species, including six Himalayan endemics, occurred almost exclusively in old-growth forest. Species composition and abundance of forest degraded by logging differed markedly from that of old-growth forest, while regrowth forest was more similar to old-growth forest. The density of medicinal plants positively correlated with tree canopy cover in old-growth forest and negatively in degraded forest, which indicates that species adapted to open conditions dominate in logged forest. Thus, old-growth forests are important as refuge for vulnerable endemics. Forest degraded by logging has the lowest diversity of relatively common medicinal plants. Forest regrowth may foster the reappearance of certain medicinal species valuable to local livelihoods and as such promote acceptance of forest expansion and medicinal plants conservation in the region.


Introduction
Medicinal plants can provide a significant source of income for rural people in developing countries, especially through the sale of wild harvested material, which contributes 15-30 % to the total income of poorer households (Hamilton 2004). The importance of medicinal plants is increasingly being recognized in the Himalayan region (Larsen et al. 2000;Olsen 2005) and particularly from ecological, social, and economic perspectives (Arnold and Pérez 2001;Negi et al. 2011). Despite their significance to rural livelihoods, the abundance and diversity of understory medicinal plant species are threatened by changes in the structural attributes of the overstory in native oldgrowth forest as a result of degradation (Gilliam 2007;Liira et al. 2007;Ticktin et al. 2007;Wyatt and Silman 2010). Changes in the structure of Himalayan forests have been reported in many studies (Kumar and Ram 2005;Nath et al. 2005;Schickhoff 1995), with most being attributed to anthropogenic activities. For example, degradation in Himalayan forests has caused rapid declines in both diversity and abundance of understory species (Mishra et al. 2004). On the other hand, in other regions of the world, it has been reported that where degraded forest are allowed to regrow, the probability of recurrence of understory species that 1 Received 1 November 2011; accepted 19 September 2012; published online 14 December 2012 had previously disappeared is greater (Barbier et al. 2008;Parrotta et al. 1997). However, such studies are scarce from Himalayan forests, so the evaluation of the extent to which forest degradation and subsequent forest regrowth brings about changes in the diversity of understory flora has not been possible.
In the mountainous region of Pakistan, local people collect around 600 medicinal plant species (approximately 10 % of the total reported plant species) to support their livelihoods (Shinwari 2010). Most of these species are understory herbs and ferns. Pakistan, with a forest cover of 4.8 % (Lubna 2001) is undergoing severe deforestation at a rate of 1.5 % per year (FAO 2005). At this rate, the current forest cover of the country may have already been reduced to approximately 4.3 %, with associated negative effects on the associated ground flora. Another possible consequence of deforestation in Pakistan is the disastrous flood of 2010, particularly in the northwest region (Rodriguez 2010), where our study was conducted. Thus, it is widely accepted that there is a need for the expansion of forest area both to enhance livelihoods of local populations and to contribute toward flood prevention measures.
The objectives of our study were 1) to compare the diversity of medicinal plants among various forest-use types, and 2) to assess the relationship between forest tree structural variables and diversity of medicinal plants. These objectives were selected with a view toward informing potential conservation and regeneration measures for medicinal plant species in Pakistan.

STUDY AREA
The study was carried out in the Khyber Pakhtunkhwa (KPK) province, the former Northwest Frontier Province of Pakistan, which constitutes 40 % of the country's forested area (Lubna 2001). Forests cover 17 % of the province (WWF-P 2004) and stretch across the Himalayas, Hindu Kush, and Karakurum mountain ranges. Covering an area of 3,312 ha, the Ayubia National Park (ANP) is one of 24 national parks in Pakistan and is situated between 33º-01' and 34º-38' N latitude and 73º-20' and 73º-30' E longitude (Fig. 1b). The study covered the ANP and its surrounding forests, which together constitute 8,978 ha in area. The study area ranged from 1,220 m to 2,865 m in altitude and lies within reach of the monsoon, resulting in a mean annual precipitation of 1,500 mm. There is snow cover between November and March and the mean annual temperature is 12°C (WWF-P 2004). The natural vegetation in the study area is Himalayan moist forest. The dominant tree species are Abies pindrow Royle, Cedrus deodara G. Don, Pinus roxburghii Sargent, Pinus wallichiana A.B. Jacks, Quercus incana Roxb, and Taxus wallichiana Zacc (Shafiq 2003).
Approximately 50,000 people currently live in twelve villages around the ANP. Fuelwood consumption of each household in the study area has been estimated at 11,000 kg/year dry weight (Aumeeruddy et al. 2004;Hussain 2003). Annually, 13 tons of fresh fodder (grasses and herbaceous species) are collected per household, of which 34 % are medicinal plants (Rabia and Khan 2004). The collection of fuelwood and fodder are mostly the responsibility of women. In a region similar to the study area, around 47 % of children are involved in the collection, 55 % of women in processing, and 70 % of men in the marketing of medicinal plants. About 10 % of the local population is employed full-time in the collection and marketing of plants with approximately 30 % being employed part-time (Ahmad 2003). Thus these extractive activities represent locally important economic activities in the region.

FOREST ADMINISTRATION HISTORY, FOREST-USE TYPES, AND PLOT SELECTION
The administration and management approach of the forests dates back to 1857 following British colonization, when they were designated as "reserved" forests and "Guzara" forests. Reserved forests are defined as forests in which all human disturbance is prohibited, whereas Guzara forests are intended to meet domestic requirements (forest products) of local people. A part of the reserved forest was declared a national park in 1984 (Ayubia National Park, ANP), and the management was handed over to the KPK wildlife department. Guzara forests and reserve forests outside the ANP are administered by the KPK forest department (Farooque 2002).
For the present study we selected three forestuse types: old-growth forest, forest degraded by logging, and regrowth forest. Old-growth forest refers to forest with little or no human disturbance, and is characterized by an abundance of large diameter trees and a mainly closed forest canopy. Old-growth forest can mainly be found in the ANP and in surrounding reserve forests and should consequently be afforded the highest level of protection from resource exploitation.
Forests degraded by logging are considerably affected by selective tree felling; many stumps of timber species such as Cedrus deodara and Pinus wallichiana and preferred fuelwood species such as Quercus spp. can be found. They occur in both reserve and Guzara forests and have a relatively open canopy. In addition, this forest-use type also supported livestock grazing, fodder collection, and the gathering of valuable medicinal plants species before the introduction of protection measures, which have restricted such anthropogenic activities since 2008 with the support of the community and World Wide Fund for Nature Pakistan.
Regrowth forests are areas naturally regenerating after heavy logging (for domestic and commercial purposes) of old-growth forests. They consist mainly of narrow-diameter trees with a high stem density. They generally occur within reserve forests and the ANP. The existence of some remnant trees identifies these forest-use types as regrowth forests, as opposed to secondary forest identified in a study by Putz and Redford (2010). Such forests have been free from logging and other anthropogenic activities for approximately 50 years, and, like oldgrowth forest, this forest-use type is also afforded some statutory protection.
In addition to the above forest-use types, there are some agroforests and derived woodlands in the Guzara forests (Adnan and Hölscher 2011) from which local communities living around the ANP meet most of their requirement for fodder, fuelwood, and medicinal plants. Despite statutory protection, some human interference for the collection of such forest resources was observed in certain parts of the ANP and reserve forests (Aumeeruddy et al. 2004). This is due to their proximity to the community, complexities in the land tenure laws, and the actual status of what is supposed to be Guzara forest, reserved forests, or the ANP, which is sometimes unclear to local people.
Data for the study were developed in three stages. First, from July to October 2008, a reconnaissance survey was conducted for the identification of selected forest-use types and the development of a GIS based map (Fig. 1b). Fifteen study plots were randomly selected within a sample area of 1,406 ha within each forest-use type and within an altitudinal range of between 2,200 and 2,400 masl. Random sample points were selected using the ILWIS software (Version 3.4: ITC 2007) which generates random values, one for the X axis and one for the Y axis, with the point of intersection representing the sample point. Each random point represented the center of each sample plot. Next, field data on forest structure were collected from the selected plots between June and September 2009. This was followed by a brief survey of the diversity of medicinal plants between August and September 2010.

MEDICINAL PLANTS AND USES
Medicinal plants include plants that are traditionally used for the treatment of various human diseases and in ethnoveterinary medicine (EVM). Species identification in the field was facilitated by local collectors of medicinal plants and by a taxonomist from the World Wide Fund for Nature Pakistan (WWF-P). Plant samples were collected from the field, and herbarium vouchers were stored at WWF-Pakistan's Nathiagali office. All medicinal plant vouchers were corrected for their botanical names, family names, and publication authors according to the software Index Kewensis version 2.0 (1997). From the available literature and collector information, local names and uses were documented (Shafiq 2003;Sher and Khan 2006), along with reproductive strategies and light requirements [Refer to Appendix in Electronic Supplementary Material (ESM)] (Fern 1997).
Most local names of medicinal plants were recorded in two languages: Hindko (spoken in the study area) and Pashto (spoken in most parts of northwest Pakistan). Medicinal plants are used whole or in parts-leaves, whole roots or tubers, stems, flowers and seeds, or fruits. Almost all medicinal plant species also have a secondary use such as fodder, vegetables, insecticides, fruits, narcotics, cosmetics, handicrafts, dyes, tea, and others. Medicinally, these plants are used to treat a number of conditions ranging from fever to cardiac problems [Refer to Appendix in Electronic Supplementary Material (ESM)] (Gillani et al. 2010;Ibrar 2003).

FOREST TREE INVENTORY AND MEDICINAL PLANTS ASSESSMENT
Plot design consisted of tree inventory plots of 20×20 m (400 m 2 ) for data collection of the stand structural variables (Table 1). Within each inventory plot, a sample plot strip of 20×2 m (40 m 2 ) ( Fig. 1c) was subdivided into 10 quadrates of 4 m 2 , and plant density, species richness, and species frequency of medicinal plants were recorded therein. The density of each plant species was recorded by counting the number of individuals; the frequency of each medicinal plant was recorded as the percentage of plots that included the species. The estimations of density and frequency of medicinals were carried out following the protocol outlined by Curtis and McIntosh (1951). Species richness, Shannon-Wiener diversity index H / , and Shannon evenness of medicinal plants (Magurran 2004) per plot strip were calculated. The Jaccard abundancebased similarity index (Chao et al. 2005) and the Sørensen similarity index based on species presence-absence (Magurran 2004) were applied to quantify medicinal species turnover within forestuse types. The Jaccard index was also applied to quantify similarity of given forest-use type with old-growth forest on the basis of medicinal plant abundance. Higher values of these indices correspond to a higher similarity in species composition, lower beta diversity, and lower species turnover (Chao et al. 2005).
Indicator species analysis (Dufrêne and Legendre 1997) was used to calculate indicator values for each medicinal plant species in each forest-use type as the product of the species' mean abundance in a particular forest-use type relative to other forest-use types and the proportion of plots in that forest-use type where it was present.
The indicator values represent the degree to which species vary among forest-use types in both relative frequency and relative abundance. The prices of medicinal plant species were taken from the local market in Abbottabad city in 2010.

STATISTICAL ANALYSIS
The Mann-Whitney test was applied to test differences in the mean values of medicinal plant variables among given forest-use types. For a given forest-use type, a Spearman correlation was performed between tree stand structural variables and medicinal plant variables. Across forest-use types, a detrended correspondence analysis (DCA) was used to identify the variables of forest tree stand structure most probably related to the densities of medicinal plants and to differentiate medicinal plant species in relation to their correlation with the axes. Thus, nine forest treestand structural variables (basal area, tree canopy cover, stem density, height, leaf area index, Shannon evenness, Shannon index, litter thickness, litter cover) were standardized to ensure equal weight was attributed to each of the variables. Data on the density of medicinal plants were logarithmically transformed as DCA based on the assumption of normally distributed data.
Indicator values of medicinal plant species among forest-use types were tested for statistical significance using a Monte Carlo test at 1,000 permutations; significance was indicated by "*".

COMPARISON OF MEDICINAL PLANTS' DIVERSITY AND ABUNDANCE AS A WHOLE BETWEEN FOREST-USE TYPES
A total number of 59 medicinal plant species (herbs and ferns) occurred on the study plots, of which 41 were perennial, thirteen were annual, and five were biennial. Twenty-nine species were used for their leaves and 20 species for their roots and tubers. In addition to medicinal uses, 24 plants were used as fodder and eleven as vegetables. Most of the medicinal plant species (28) were used to treat various intestinal diseases, followed by 21 for stomach diseases [Refer to Appendix in Electronic Supplementary Material (ESM)].
The highest number of the medicinal plant species (55) occurred in old-growth forest, 40 species grew in forest degraded by logging, and 49 species grew under regrowth forest. The mean species richness of medicinal plants (18.1 per 40 m 2 plot) in old-growth forest plots was almost two times higher than in forest degraded by logging (Table 2). For medicinal plants, presence/ absence in old-growth forest was observed with significantly higher Sørensen similarity index within the plots (0.42) in comparison to forest degraded by logging (0.29). For medicinal plant abundance, a significantly higher Jaccard similarity index with old-growth forest was recorded for regrowth forest (0.34) than for forest degraded by logging (0.20) (Fig. 2).
The mean density of all medicinal plants (20 m -2 ) in regrowth forest plots was greater than that in forest degraded by logging (13.7 m -2 ). Similarly, the mean Shannon diversity index H \ for regrowth forest (2.3) was also significantly higher than forest degraded by logging (1.8). The mean Shannon evenness of medicinal plants under both regrowth forest (0.84) and forest degraded by logging (0.80) was significantly higher than in old-growth forest (0.77). Oldgrowth forest and regrowth forest showed no significant differences for mean species' density, richness, or Shannon diversity index H \ (Table 2). Twenty medicinal species were identified as significant indicator species for given forest types. Fourteen were indicator species of old-growth forest, including species endemic to the Himalayan region such as Aconitum heterophyllum, Arum venosum, Coptis tecta, Impatiens bicolor, Paenonia emodi, and Podophyllum emodi. Four species were indicator species of regrowth forest, and only two were indicator species of forest degraded by logging. The local market prices of these medicinal plant species ranged between 0.2 and 9.0 USD kg -1 (Table 3).

RELATIONSHIPS BETWEEN TREE STAND STRUCTURE AND MEDICINAL PLANTS
Tree canopy cover in old-growth forest showed a significantly positive correlation with density of medicinal plants (r 00.69, n 015, p < 0.01). However, tree canopy cover was significantly negatively correlated with Shannon evenness of medicinal plants under old-growth forest (r0−0.58, p<0.05). In contrast, tree canopy cover under forest degraded by logging showed a significantly negative correlation with density of medicinal plants (r0 −0.87, p<0.01). The Shannon H \ of trees showed a significantly positive correlation with Shannon H \ of all medicinal plants under each forest-use type ( Table 4).
As indicated by the DCA on the 45 plots of three forest-use types, the axis-1 (eigenvalue 0 0.36) of forest stand structural data correlated   aparine, and Trifolium repens, which implies that these species decrease with increasing tree canopy cover. The densities of 33 medicinal plants did not show any correlation with axis-1, such as Achillea millefolium, Agrimonia eupatoria and Gentiana kurroo, and they were canopy-cover indifferent species (Fig. 3).

Discussion
Multi-use medicinal plants are an important part of the understory flora and play an essential role both socioeconomically and ecologically in the region. The present study indicates variation in the diversity of understory medicinal plant species among studied forest-use types. In particular, oldgrowth forests had a higher diversity of medicinal plants compared to forests degraded by logging. Species density, species richness, Shannon diversity, and similarity of medicinal plants were also higher in old-growth forest than in forest degraded by logging. A higher density of medicinal plants was found in relation to higher tree canopy cover, while a higher diversity of medicinal plants was related to higher tree diversity in old-growth forest. These findings support other studies on Himalayan forests that indicate a higher diversity and abundance of understory plant species in undisturbed compared to disturbed native forest (Bhuyan et al. 2003;Sundriyal and Sharma 1996;Uniyal et al. 2010). The effect of a more open canopy may explain the lower values of medicinal plant variables under forest degraded by logging in comparison to oldgrowth forest. Disturbed overstory has been proposed as one of the reasons for declines in diversity of plant species in the Khasi hills of Meghalaya, India (Mishra et al. 2004). In contrast, some studies showed lower diversity of medicinal species under undisturbed, abundant overstory in comparison to disturbed and more open overstories (Aguilar-Støen and Moe 2007;Nath et al. 2005;Voeks 1996Voeks , 2004. However, a comparison of such studies with our study is complicated by variation in a number of prevailing factors between the study areas. In our study, lower tree canopy cover was related to a higher density of medicinal plants under forest degraded by logging, most probably due to the presence of certain species adapted to sunny conditions. In addition, the past land-use history of forest degraded by logging in our study revealed overcollection of medicinal plants that may have been associated with a lower abundance and diversity of medicinal plants. Regrowth forest was found to have a more diverse medicinal flora than forest degraded by logging. Studies have confirmed that if a forest is allowed to regenerate and is protected, the understory species that had once disappeared from it also regenerate (Islam et al. 2001;Keersmaeker et al. 2011;Parrotta et al. 1997). Forest regrowth alters conditions for the ground vegetation, such as soil fertility, light, temperature, and moisture (Barbier et al. 2008;Carnevale and Montagnini 2002;Shankar et al. 1998), all of which affect competition dynamics at ground level. In our study, higher medicinal plant density and diversity were related to higher tree diversity under forest regrowth. Moreover, regrowth and old-growth forests exhibited only minor differences in medicinal plant density, richness, and diversity. Similar results were found by Harrelson and Matlack (2006) in their study of Mesophytic old and young forest stands in southeastern Ohio, U.S.A., suggesting that regrowth forest supports the recovery of medicinal understory plants to a level similar to that found in old-growth forest when it was protected and allowed to naturally regrow.
Medicinal plants under forest degraded by logging and regrowth forest were more evenly distributed than those under old-growth forest. Low evenness of species under old-growth forest might have resulted from high variations in the densities of medicinal plants that are adapted to various growing conditions. This is evident from the results of this study, wherein higher tree canopy cover resulted in low evenness of medicinal plants under old-growth forest. A study on Quercus leucotrichophora forest in the Garhwal Himalaya likewise showed lower evenness of herb species under undisturbed forest than under highly-disturbed forest (Uniyal et al. 2010). And D'Amato et al. (2009) found lower herbaceous species evenness under old-growth forest of Tsuga canadensis due to lower species richness when compared to second-growth forest.
Across forest-use types, higher tree canopy cover was associated with the increased abundance of certain medicinal plants. This included the most valuable species such as Aconitum heterophyllum, Podophyllum emodi, Valeriana jatamansi, and Viola canescens, which are also adapted to deep-shade conditions. The density of such medicinal plant species was found to be higher under old-growth forest in comparison to other forest-use types. This influence might be due to the natural effect of crown cover on the amount of available light for ground vegetation (Shanley and Luz 2003). The abundant overstory provides habitat that promotes shade-adapted understory species while suppressing sun-adapted species (Balandier et al. 2006;Sharma et al. 2009).
The old-growth forests harbor many species that are endemic to the Himalayan region. Studies have confirmed that nearly half of the plant population of the Himalayan region is endemic (Kala 2004;Samant and Dhar 1997). For example, Aconitum heterophyllum, an endemic species of the northwestern Himalayas (Kala 2004) that is critically endangered (Beigh et al. 2006), was only encountered in old-growth forest and was absent from forest degraded by logging and regrowth forests. Other species that were only found in old-growth forest, such as Arum venosum, Coptis tecta, and Impatiens bicolor, are all endemic to the Himalayan region and vulnerable (Ali 2011;WWF-P 2004). The reasons for the occurrence of these plants only in old-growth forest may be manifold, but one is that they are quite likely adapted to shady conditions. Moreover, regrowth forest and forest degraded by logging might lack the ability to provide a proper environment for the reappearance of these species due to past anthropogenic factors. Other species, such as Paeonia emodi and Podophyllum emodi, occurred only in old-growth and regrowth forests. Although distributed throughout the Himalayas, these species have been reported as highly vulnerable in the study area by WWF-Pakistan. Comparably shadier conditions than that of forest degraded by logging and about four decades of limited human interference may explain the reappearance of such species under regrowth forest. Some species were only found in old-growth forest and forest degraded by logging, such as Gentiana kurroo. Their presence in forest degraded by logging rather than in regrowth forest might be due to the closer proximity to old-growth forest.
In our study, low tree canopy cover appeared to favor species such as Galium aparine, which are more commonly adapted to sunny conditions. The density of these species under regrowth forest and forest degraded by logging was higher than under old-growth forest, perhaps reflecting the more open canopies in the former forest-use types. Some of these species only occurred in forest degraded by logging and regrowth forest, such as Chrysanthemum leucanthemum. Medicinal plant species were also observed to be indifferent to canopy cover, such as Achillea millefolium and Agrimonia eupatoria. These species appear to tolerate partially shady conditions (Fern 1997), although other variables may be relevant as well.
The medicinal plant species recorded in this study are valuable in terms of their traditional usage as medicines at the local level. Moreover, about one-third represent direct economic resources, as indicated by market prices equal to or greater than 0.2 USD kg -1 . Species such as Aconitum heterophyllum, Valeriana jatamansi, and Viola canescens, which occurred in higher abundance and frequency in old-growth forest, are traded along an established chain to the city markets and beyond. Demand for such medicinal plants in Pakistan is increasing at present 15 % annually (Shinwari 2010). For example, the total annual consumption of Aconitum heterophyllum increased sharply from 740 kg in the 1980s (Williams and Ahmad 1999) to around ten tons in 2003, or roughly a thirteen-fold increase over twenty years. In the ANP, the extraction of A. heterophyllum has increased from 64 kg in 2003 (Ibrar 2003) to 500 kg in 2007, or an eight-fold increase in four years. The absence of such species in regrowth forest and forest degraded by logging may well encourage increased harvest in oldgrowth forest.
In conclusion, old-growth forest in the Pakistani Himalaya sustains a highly diverse medicinal flora and plays an important role as refuge for vulnerable endemics. Forest degraded by logging has the lowest diversity of medicinal plants and, to a large extent, harbors more ubiquitous medicinal species. If regrowth is allowed in these degraded areas, this may foster the reappearance of certain medicinal species valuable to local livelihoods and thereby promote local acceptance of forest expansion and medicinal plant conservation in the region.