Conserving an Endangered Canid: Assessing Distribution, Habitat Protection, and Connectivity for the Dhole (Cuon alpinus) in Bhutan

Most canids face population declines and range contractions worldwide. Although the dhole (Cuon alpinus) is widely distributed across 10 countries in South and Southeast Asia, limited studies exist on this species. Despite its globally “Endangered” status and ecological role as an apex predator, assessments on its distribution are limited to a few landscapes and countries. This explains the lack of a dhole-specific species conservation plan in most range countries, including Bhutan where no current population estimate exists. The species has also recovered from a country-wide poisoning campaign in the 1970s and 80s. In this study, we determine the dhole's distribution pattern and assess the protection and connectivity of dhole habitat in Bhutan. We anticipated dholes to be extant within their habitat well-represented in protected areas (PAs) and biological corridors (BCs). We used 721 georeferenced dhole occurrence records and eight environmental variables in MaxEnt software to model potential dhole distribution and habitat suitability. The model output was overlaid on the spatial layers of PAs and BCs to assess habitat protection and connectivity. As anticipated, we found the dhole widely distributed in all districts, PAs, and BCs in Bhutan. Dholes were recorded at the highest elevation range limit of 4,980 m above sea level, which overlapped with the “Vulnerable” snow leopard (Panthera uncia). Our model identified 72% (27,634 km2) of the country as suitable areas for dholes, of which, 31% (11,899 km2) was highly suitable and 41% (15,735 km2) was moderately suitable. Contrary to our expectation, PAs and BCs encompassed only 29% (8,046 km2) and 12% (3,185 km2) of suitable areas for dholes, respectively. A vast majority of the areas we deemed suitable for dholes currently remain unprotected, thus making dholes more vulnerable to human persecution and local extermination. We recommend adjusting PA boundaries to fully encompass suitable dhole habitat, and also advocate improved livestock husbandry to reduce dhole related livestock predation and minimize conflict, thereby ensuring its long-term survival in Bhutan.


INTRODUCTION
With 13 genera encompassing 37 species in 81% of countries, canids are widespread across diverse habitats and environments in all continents except Antarctica (Fleming et al., 2017). As apex predators, large canids can influence prey populations and cause trophic cascades when their populations fluctuate (Newsome and Ripple, 2015;Fleming et al., 2017). Globally, canid distributions overlap with human-modified landscapes (Srivathsa et al., 2019b) to pose conservation challenges (Lamb et al., 2020), resulting in population declines caused by habitat loss, prey decrease, human persecution, disease, and overexploitation (Ripple et al., 2014). The most notable examples include the extermination of the Mexican wolf (Canis lupus baileyi) from its natural range (Brown, 1983), a 93% range contraction of the African wild dog (Lycaon pictus; Woodroffe and Sillero-Zubiri, 2020), eradication of gray wolves (Canis lupus) from most of the United States and Europe (Mech, 1995), and extinction of the Falkland Islands wolf (Dusicyon australis; Sillero-Zubiri, 2015).
The dhole (Cuon alpinus, Pallas 1811), or Asiatic wild dog, is one of the most widely distributed members of the 10 canid species described from Asia (Din et al., 2013). It is an apex social carnivore that preys mostly on ungulates (Kamler et al., 2012 in forested areas across most of South and Southeast Asia and parts of China (Srivathsa et al., 2014;Kamler et al., 2015). Although largely restricted to protected areas (Kamler et al., 2015), dholes have also been recorded in unprotected secondary forests, multi-use forest fragments, and agro-forest plantations adjacent to protected areas (Srivathsa et al., 2014). Despite being shy and elusive with infrequent contacts with humans , this canid has been extirpated from ∼82% of its original range through human persecution and habitat loss (Wolf and Ripple, 2017). It is currently listed as globally "Endangered" by the IUCN based on an estimated population of 4,500-10,500 individuals comprising <2,500 adults, with India housing substantial populations concentrated to the south of the Ganges River in the Western Ghats and central forested regions (Kamler et al., 2015).
Although endangered, the dhole has received less conservation attention than other charismatic carnivores (Widodo et al., 2020). Much of the focus on dholes is related to their depredation on livestock (e.g., Katel et al., 2015;Srivathsa et al., 2020). While basic information from camera-trapping (presence and relative abundance) has been reported for sitespecific dhole populations in India (Datta et al., 2008;Bashir et al., 2014), Myanmar (Rao et al., 2005), Peninsular Malaysia (Kawanishi and Sunquist, 2008), and Laos (Johnson et al., 2006), there remains a lack of country-specific consolidated distribution data and range maps for dholes. This hinders the assessment of their population status (Karanth et al., 2009) for both in-country and regional conservation planning (Srivathsa et al., 2014;Punjabi et al., 2017), because knowledge of suitable sites where species occur and survive can aid in conservation planning (Papeş and Gaubert, 2007).
Several modeling studies on dhole distribution and occupancy have been carried out at varying scales. At local scales, parkwide potential dhole distribution modeling was carried out by Namgyal and Thinley (2017) in Bhutan's Jigme Dorji National Park and by Rahman et al. (2018) in Indonesia's Ujung Kulon National Park, whereas Singh et al. (2020) recently reported on dhole occupancy in India's Dampa Tiger Reserve. At the broader landscape scale, Srivathsa et al. (2014), Punjabi et al. (2017), and Srivathsa et al. (2019a) modeled dhole occupancy across the Western Ghats of Karnataka, India. Similarly, Srivathsa et al. (2019b) also assessed occupancy by dholes across the Pench-Kanha Landscape in Madya Pradesh, India. Recently, Widodo et al. (2020) modeled potential dhole distribution across the Rimbang Baling and Tesso Nilo landscapes in Sumatra. At the country level, Jenks et al. (2012) used maximum entropy modeling to predict potential dhole distribution in Thailand based on dhole occurrence data from 15 protected areas while Karanth et al. (2009) used occupancy modeling to predict areas of dhole occurrence in India.
In Bhutan, dholes are apex predators similar to tigers (Panthera tigris; Thinley et al., 2018) and snow leopards (Panthera uncia; (Leki and Shrestha, 2018). There is no current population estimate for dholes in Bhutan. They were, however, almost extirpated from the country in the 1970s and 80s by mass poisoning campaigns due to blames over persistent livestock depredation (Wang and Macdonald, 2006;Thinley et al., 2011;Namgyal and Thinley, 2017). Because dholes are known to control populations of wild pig (Sus scrofa), the principal crop-raiding species in Bhutan (Wangchuk, 2004;Thinley et al., 2018), it is believed that wild pig populations in Bhutan substantially increased and intensified crop damage after the mass extermination of dholes (Wangchuk, 2004). Despite the dhole population recovering and re-establishing itself in Bhutan from the late 1990s with some probable recolonization from the neighboring Indian states of Assam and West Bengal (Wangchuk, 2004), little is known on its current distribution in Bhutan (Namgyal and Thinley, 2017). The dhole is still not listed as a protected species in Schedule I of Bhutan's Forests and Nature Conservation Act despite being globally endangered (Namgyal and Thinley, 2017). Therefore, it is important to ascertain dhole distribution in Bhutan at the landscape level to promote efficient research and planning decisions (Guisan et al., 2013). The only previous attempt to discern dhole distribution in Bhutan was by Wangchuk (2004) who interviewed 67 field forestry staff members and surveyed residents in 18 villages across seven dzongkhags (districts) of Gasa, Paro, Punakha, Thimphu, Trongsa, Wangduephodrang, and Zhemgang (Figure 1). However, no distribution map was produced in addition to the anecdotes of localities where dholes were present or absent. As such, there is no information on how well dhole habitats are protected within Bhutan's protected area network.
An adequate knowledge of distribution, including associated environmental and anthropogenic variables, further enables appropriate modeling to predict additional areas where dholes are likely to occur, both within and outside protected areas, as recently demonstrated by Thinley et al. (2021) for the tiger in Bhutan. Because dholes are prominently linked to livestock predation in Bhutan (Thinley et al., 2011;Katel et al., 2015;Rajaratnam et al., 2016), determining and modeling their distribution further enables an identification of potential human-wildlife conflict hotspots to prioritize mitigation efforts . Juxtaposing conflict hotspots and habitat protection against the current and predicted dhole distribution can provide the spatial framework to develop an appropriate dhole conservation plan for Bhutan.
In this study, we investigated dhole distribution based on current presence records across Bhutan and assessed the protection and connectivity of its habitat. Specifically, we modeled potential dhole distribution and suitable habitat coverage in Bhutan's protected areas (PAs) and interconnecting biological corridors (BCs). Based on Bhutan's strong conservation policy and high proportion of land dedicated to nature conservation, we anticipated dholes to be extant in their habitats well-represented in PAs and BCs.

Study Area
Bhutan (Figure 1) is one of the least-populated countries of Asia with only 735,553 people (NSB, 2017). Situated in the eastern Himalayas (Figure 1 inset) and administratively divided into 20 dzongkhags (district), Bhutan encompasses an area of 38,394 km 2 (NSB, 2018) and is well-known for its rich biodiversity. Approximately 71% of the country is forested (FRMD, 2016) and encompasses 11,248 species of plants and animals, including 129 mammal species (NBC, 2017). This mammalian community includes four wild canid species, namely dhole, Tibetan wolf (Canis lupus chanco), golden jackal (Canis aureus), and red fox (Vulpes vulpes) . The country's Constitution mandates the government to maintain a minimum of 60% forest cover for eternity. As such, more than half of the country (51%) is designated as a protected area network (Figure 1) comprising protected areas and interconnecting biological corridors (Dorji et al., 2019). The country's topography is mountainous and highly rugged with a pronounced elevation range from 97 m above sea level (a.s.l.) in the southern foothills to 7,750 m a.s.l. in the greater Himalayas near the Chinese border (Tshering et al., 2020).

Modeling Potential Dhole Distribution or Habitat
We modeled potential dhole distribution and habitat suitability in Bhutan using MaxEnt program version 3.4.3 (Phillips et al., 2020) which uses a widely employed maximum entropy method of modeling species distribution based on presence-only data (Phillips et al., 2006). In MaxEnt, georeferenced occurrence points of a target species are associated with environmental variables to yield a spatial layer representing the most widespread probability of its presence, given constraints imposed by these environmental layers (Elith et al., 2011). The resulting layer also constitutes a habitat suitability layer of a species and depicts its realized niche, which is a subset of its fundamental niche (Phillips et al., 2006). MaxEnt modeling was chosen because our dataset solely comprised dhole presence points, thereby, precluding the use of a boosted regression tree based on additional absence data (Yu et al., 2020) and species occupancy modeling based on repeated surveys (MacKenzie et al., 2006).
We utilized a database of 721 georeferenced dhole occurrence points in Bhutan from 2014 to 2019 based on cumulative pooling of dhole records from: (a) a 2014-2015 nationwide camera trapping survey on tigers (Thinley et al., 2015) up to 4,500 m a.s.l.; (b) a 2015-2016 nationwide camera trapping survey on snow leopards  between 3,500 and 5,500 m a.s.l.; (c) camera trap and sign surveys across protected areas between 2014 and 2019; (d) SMART (Spatial Monitoring and Reporting Tool) patrolling reports submitted by wildlife personnel between 2015 and 2019; and (e) camera trapping of wildlife in catchment studies  in Wangchu and Kholongchhu sub-basins between 2018 and 2019.
Following Thinley et al. (2021), we selected eight environmental (geophysical and anthropogenic) variables deemed to influence dhole occurrence based on its ecology (Johnsingh and Acharya, 2013), particularly prey selectivity (Wangchuk, 2004;Wang and Macdonald, 2009;Thinley et al., 2011) and habitat use (Aryal et al., 2015;Namgyal and Thinley, 2017). All environmental variables were processed in ArcMap version 10.7.1 where spatial layers were classified into 10 categories based on suitability to dholes and standardized with a spatial resolution of 30×30 m cell corresponding to the spatial resolution of the Digital Elevation Model (DEM), Transverse Mercator projection, DRUKREF 03 coordinate system, and geographic extent of Bhutan. Dhole points co-occurring within the same 30×30 m cell were omitted to minimize errors from spatial autocorrelation (Kanagaraj et al., 2011).
Categorization of spatial layers for environmental variables was based on Thinley et al. (2021) as follows: 1) Elevation: reclassified from Bhutan's Digitial Elevation Model (Jarvis et al., 2006) into 10 classes such that lower elevations were ranked more suitable than higher elevations ( Table 1); 2) Slope: extracted from the DEM such that lower slope classes were ranked more suitable than higher classes ( Table 1); 3) Prey richness: prey densities are ecological determinants of carnivore densities (Karanth et al., 2004) and dholes generally prey on medium to large ungulates (Wang and Macdonald, 2009;Srivathsa et al., 2020). We, therefore, merged the distribution of nine potential wild prey species (Wangchuk, 2004;Wang and Macdonald, 2009;Thinley et al., 2011)-wild pig, sambar (Rusa unicolor), muntjac (Muntiacus muntjac), Himalayan goral (Naemorhedus goral), serow (Capricornis sumatraensis), hog deer (Axis porcinus), alpine musk deer (Moschus chrysogaster), gaur (Bos gaurus), and Bhutan takin (Budorcas taxicolor whitei). Spatial layers were obtained from the Field Guide to Mammals of Bhutan , rasterized, and reclassified with higher suitability values assigned to areas with higher prey richness ( Table 1); 4) Distance from rivers and streams: rasterized from the Drainage Map of Bhutan based on Euclidean distance; areas closer to rivers and streams were assigned higher suitability values than those farther away ( Table 1);  (2017), and Thinley et al. (2021), in our generated MaxEnt model we used the default setting of 500 iterations with a convergence threshold of 0.00001, a regularization multiplier of 1, and a maximum background of 10,000 background pseudo-absence points with 50% random tests and 10 replicates. We assessed model performance using AUC (Area under the Receiver Operating Curve) whereby values ≤0.5 indicate very poor fit, >0.5 indicate good fit, and equal to 1 indicate perfect fit (Fielding and Bell, 1997). The resulting MaxEnt probability surface was exported to ArcMap and reclassified into the probabilities of dhole occurrence ("highly probable, " "moderately probable, " and "not probable"). These probabilities correspond to the suitability surfaces ("highly suitable, " "moderately suitable, " and "unsuitable") because a relatively high number of dhole occurrence points collected across a wide range of environmental conditions optimally reflects the dhole's fundamental niche.

Assessing Dhole Protection and Habitat Connectivity
We overlaid protected area (PA) and biological corridor (BC) layers on the dhole suitability layer to assess how much of the habitat that we deemed suitable for dholes is encompassed within PAs (assessing dhole protection) and corridors (assessing dhole habitat connectivity between protected areas). Using the Field Calculator tool in ArcMap, we computed areas (km 2 ) for each of the dhole suitability classes falling within and outside PAs and BCs.

Dhole Distribution
Dholes were distributed throughout Bhutan across all 20 districts and in all PAs and BCs (Figure 2)  Himalayas (Figure 1). They were also present in almost all habitat types ranging from sub-tropical forests in the lowlands to alpine meadows in the uplands. The highest concentration of dhole occurrence records was observed in JDNP whereas the least was recorded in Wangchuck Centennial National Park (WCNP) (Figure 2).

Suitable Areas for Dholes
The MaxEnt model (AUC = 0.73 for training data; 0.72 for test data) predicted that suitable area for dholes covered 72% (27,634 km 2 ) of Bhutan (Figure 3; Table 3). This comprised 31% (11,899 km 2 ) of highly suitable area and 41% (15,735 km 2 ) of moderately suitable area (Figure 3; Table 3). Among the eight environmental variables, model prediction (or gain) was maximally influenced by slope (26.4%), followed by distance from human settlement (24.3%), elevation (16.5%), and land cover (13.1%; Table 4). Distance from water bodies ( Table 4) contributed least (1%) to model prediction. Overall, suitable area for dholes coincided with flat and moderately flat area situated further away from human settlements in forested areas below 5,000 m a.s.l. The remaining 28% (10,760 km 2 ) of modeled area was unsuitable for dholes (Figure 3; Table 3) and overlapped with steep areas closer to human settlements and roads, and with areas above 5,000 m a.s.l. which were either too cold or permanently covered with snow and glaciers.

Dhole Distribution
We mapped the first-ever nation-wide distribution of dholes in Bhutan (Figure 2), which also constituted the first study of its kind amongst dhole range countries in the Himalayan Mountains. Dholes were present in high elevation alpine meadows and screes, and their distribution overlapped with "Vulnerable" snow leopards which are known to range between 3,400 and 5,186 m a.s.l. . Dholes are reported to be sympatric with snow leopards in the arid region of the Altun Mountains in western China (Xue et al., 2015). We documented the highest elevation occurrence (4,980 m a.s.l.) for dholes in Bhutan, which also constitutes the highest across its entire distribution range, thus representing an uppermost range limit for dholes. Whether this range extreme is attributed to climate change or prey availability, including livestock such as yak (Bos grunniens), needs further investigation. To date, high altitude livestock predation in Bhutan has been mainly attributed to snow leopards (Rajaratnam et al., 2016). Further investigation is required on a potential dietary overlap with the snow leopard (Namgyal and Thinley, 2017) to ascertain the competitive impact of dholes on this iconic flagship carnivore species for the montane Himalayas.
As anticipated, we also discovered that the dhole has the widest distribution amongst large carnivores in Bhutan. It was previously not reported in the three eastern districts of Trashigang, Samdrupjongkhar and Pemagatashel (Wangchuk, 2004), but is now confirmed to be present in all 20 districts of Bhutan. Its wide distribution also indicates the possibility of frequent conflicts with humans due to livestock depredation (Namgyal and Thinley, 2017) given that dholes are principal livestock predators in western Bhutan (Katel et al., 2015;Tshering and Thinley, 2017) and neighboring Arunachal Pradesh in India (Lyngdoh et al., 2014). Increased prevalence of livestock predation by dholes is a distinct possibility in Jigme Dorji National Park (Thinley et al., 2011) which harbored the highest concentrations of dhole occurrence records. Livestock loss presents a significant socioeconomic setback for rural farmers in Bhutan's predominantly agrarian society because the loss of yak results in sizable income loss to upland pastoralists, while the loss of cattle compromises agricultural production and nutrition for lowland agro-pastoralists (Sangay and Vernes, 2008;Rajaratnam et al., 2016).
In relation to agricultural crop loss, wild pigs are responsible for more crop damage compared to any other species in Bhutan. Such losses have major negative impacts on subsistence farmers throughout Bhutan. In fact, the crop damage caused by wild pigs is so severe that wild pigs are considered a national pest with significant funding allocated by the government to control this species. Consequently, the controlling effect by dhole on wild pig populations (Wangchuk, 2004) suggests that increasing the number and expanding the distribution of dholes in Bhutan may lead to lowered numbers of wild pig, which ultimately, would benefit local farmers.
Slope was the top contributor to dhole distribution and habitat suitability (Table 4), indicating the influence of Himalayan rugged topography on dhole distribution in Bhutan. Aryal    (2015) also determined slope to be the most significant factor influencing dhole presence and occupancy in Nepal's Dhorpatan Hunting Reserve, whereby dholes used gently sloped land more than steep areas. This is likely related to the cursorial hunting strategy of dholes which is more effective on flatter terrain and gentle slopes compared to more rugged terrain or areas with steeper slopes (Kamler et al., 2012). In contrast to Jenks et al. (2012), contribution of prey richness to modeled potential dhole distribution in our study was minimal and much lower than that of land cover. This is possibly because forested habitats throughout Bhutan naturally harbor a high richness of potential prey species . We do, however, acknowledge the greater role of prey density in influencing dhole distribution. This variable was excluded due to the lack of its spatial layer, which constitutes a limitation to our study.

Dhole Conservation Mismatch
We assessed whether protected areas and biological corridors are adequate for conserving dhole habitats. Although protected areas in Bhutan were designated for multiple species, we found them inadequate for protecting and connecting wide-ranging species like the dhole, similar to that observed for the tiger (Thinley et al., 2021). Contrary to expectations, only 29% of total suitable habitat for dholes was available in protected areas, reflecting an increased vulnerability of the species to further human persecution and habitat loss through livestock depredation and land-use changes outside protected areas. Similarly, only 30% of potential habitat for dholes in Thailand was encompassed within protected areas (Jenks et al., 2012) while 41% of areas occupied by dhole were inside protected wildlife reserves in the Western Ghats landscape within Karnataka, India (Srivathsa et al., 2014). As such, the current level of landscape protection may be insufficient to support functional dhole meta-populations (Bargelt et al., 2020) across Bhutan as equally demonstrated by the meager 33-35% of suitable areas captured by Bhutan's protected areas for the wide ranging tiger (Thinley et al., 2021).

Dhole as a Potential Umbrella Species
Dholes are also estimated to require five times more land area than large-bodied carnivores such as tigers (Kamler et al., 2012), mainly because of the social structure of populations living in exclusive territories (Johnsingh, 1982), unlike solitary tigers with overlapping territories (Carter et al., 2015). Large space requirements in conjunction with a hypercarnivorous diet (Van Valkenburgh, 1991) make dholes more vulnerable to extirpation, as evidenced by their disappearance from more reserves than tigers (Woodroffe and Ginsberg, 1998). Most reserves in Asia are typically focused on conserving umbrella species such as tigers (e.g., Wikramanayake et al., 2011) and greater onehorned rhinoceros (Rhinoceros unicornis) in Nepal (e.g., Aryal et al., 2017) to highlight conservation efforts. As exemplified above, umbrella species are typically large-bodied animals which require large areas to ensure species persistence. They are good surrogates for overall biodiversity but are more sensitive than other species to human activities, ecosystem changes, and habitat destruction whilst having the largest land requirements and most stringent ecological needs (Woodroffe and Ginsberg, 2005).
Recently, Kaszta et al. (2020) advocated habitat prioritization for the medium-sized clouded leopard (Neofelis nebulosa) as a good indicator and focal species to address limited resources for immediate biodiversity conservation actions to conserve forest ecosystems and forest-dependent biodiversity in Southeast Asia. Based on this premise and the results of our study, we suggest that the dhole, given its unusually high requirements for space and prey, could constitute a more effective umbrella species in Asia and a driver for future designation of protected areas and/or expansion of current reserves with dhole metapopulations.

Dhole as a Keystone Species
In most terrestrial ecosystems, large carnivores have been identified as "keystone species" (Terborgh et al., 1999;Ripple et al., 2014) based on the premise that keystone species exert disproportionately larger influence on an ecosystem relative to their abundance (Power et al., 1996). Although the top-down effects of dholes are loosely described (see Thinley et al., 2018), previous research showed that dholes killed more ungulates than sympatric leopards (Panthera pardus) and tigers (Venkataraman, 1999), indicating that dholes have a greater impact on ungulate numbers relative to other large Asian carnivore species. As such, there may be a parallel to pack-living wolves in North America which cause trophic cascades primarily via predation on large ungulates which, in turn, affects vegetation growth patterns across the landscape (Ripple et al., 2001;Beschta and Ripple, 2009). Further repercussions from the decline of dholes are exemplified in Bhutan, whereby wild pig populations thrived and increased crop depredation after dholes were exterminated from many areas (Wangchuk, 2004). The dhole may, therefore, be a top keystone carnivore in Bhutan, and probably across its range. Consequently, their presence may have a greater impact on biodiversity in Asia, compared to other large carnivores. Thus, if protected areas in Asia are to preserve entire ecosystem functions, adequate protection and conservation of dholes should be equally considered.

Conservation Implications
Based on our study, Bhutan is a stronghold for dholes in the eastern Himalayas due to their widespread distribution and availability of large tracts of suitable habitat. In order to ensure long-term survival of the dhole in Bhutan, we recommend replicating the recommendations of Thinley et al. (2021) for tigers in terms of readjusting PA boundaries to encompass prime dhole habitats and extending dhole conservation efforts outside the PAs. However, local people exhibit resentment to dholes and their conservation due to persistent dhole-related livestock predation (Katel et al., 2015). This is principally driven by socioeconomic losses because pastoral communities experiencing human-carnivore conflict tend to have low income with low tolerance to carnivores and their conservation (Ahmad et al., 2016). As such, there is likelihood for retaliation against dholes reminiscent of historic poisoning efforts against the species in the 1970-80s (Wangchuk, 2004). Therefore, dhole conservation efforts both within and outside the PAs need to incorporate efforts to improve livestock husbandry like grazing livestock in and around villages, including stall-feeding and cooperative herding of livestock in forests during the day (Katel et al., 2015). Tshering and Thinley (2017) further recommended stock improvement, fodder development, pasture development, and livestock insurance schemes to reduce livestock predation by dholes, while Sangay and Vernes (2008) advocated non-grazing of livestock in depredation hotspots. Dhole-specific livestock insurance schemes are also feasible to alleviate socio-economic loss from livestock predation by dholes. Local people should also be educated about the positive impacts that dholes have on ecosystems, including suppression of major crop-destroying species such as wild pigs. Based on this control of crop depredators, rural farmers might be inclined to adopt livestock protection measures against dholes as well. This will alleviate livestock losses in addition to a reduction in crop damage. We further advocate the listing of the dhole in Schedule I of Bhutan's Forests and Nature Conservation Act, urging for increased legal protection of this globally endangered canid. Lastly, we recommend other dhole range countries to conduct a similar study on the role of protected areas in conserving dhole populations to ensure species viability across its distribution area.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

ETHICS STATEMENT
The study was reviewed and approved by Department of Forests and Park Services, Bhutan.

AUTHOR CONTRIBUTIONS
PT: conceptualization, data collection, data curation, formal analysis, and writing-original draft. RR: conceptualization, validation, writing-review, and editing. JK: validation, writingreview, and editing. CW: data collection, data curation, formal analysis, and writing-review. All authors contributed to the article and approved the submitted version.

FUNDING
This study was funded by the Royal Government of Bhutan.

ACKNOWLEDGMENTS
This study was conducted as part of the annual performance agreement (APA) of the Department of Forests and Park Services (DoFPS) that was duly approved by the Royal Government of Bhutan (reference number UWICER/ADM-18/2019/Camp Thimphu dated 10/01/2019). We thank Lobzang Dorji, the Director of DoFPS, for administrative approval and support. We also thank relevant wildlife officials and Chief Forest Officers of all territorial forest divisions and protected areas for sharing dhole occurrence data.