Comparative Ecological Perspectives of Two Ancient Lineages of Gray Wolves: Woolly Wolf (Canis lupus chanco) and Indian Wolf (Canis lupus pallipes)

Geographical isolation can often lead to speciation, and two disconnected populations of the same species living in drastically different bioclimatic regions provide an opportunity to understand the process of speciation. The Woolly wolf is found in the cold-arid, Trans-Himalayan landscape, while the Indian wolf inhabits the semi-arid grasslands of Central India. Both the lineages of wolves from India have generated scientific debate on their taxonomic status in recent years. In this study, we collected data and reviewed published literature to document the ecological and behavioral differences between the Woolly wolf and the Indian wolf. Most studies have used genetic data; hence we discuss variation in spatial ecology, habitat preferences, vocalization, diet diversity and cranial measurements of these two subspecies. The spatial ecology of two lineages was compared from the data on three Woolly and ten Indian wolves tagged with GPS collars. The telemetry data shows that there has been no difference in the day-night movement of Woolly wolves, whereas Indian wolves show significant high displacement during the night. The BBMM method indicated that Woolly wolf home ranges were three times larger than the Indian wolf. The Woolly wolf diet is comprised of 20 different types of food items, whereas the Indian wolf diet consists of 17 types. The Woolly and Indian wolf largely depend upon domestic prey base, i.e., 48.44 and 40.34%, respectively. We found no differences in the howling parameters of these subspecies. Moreover, the Woolly wolf skull was significantly longer and broader than the Indian wolf. Wolves of India are ancient and diverged from the main clade about 200,000–1,000,000 years ago. Their genetic and ecological evolution in different bioclimatic zones has resulted in considerable differences as distinct subspecies. The present study is a step in understanding ecological differences between two important, genetically unique subspecies of wolves.


INTRODUCTION
Divergent selection can occur among populations with varying environmental conditions, habitats, climates, resource types and phylogeographic patterns (Schluter, 2000;Dieckmann et al., 2004). The phylogeographic patterns isolate populations that inhabit same geographic range but are unable to interact due to the presence of natural barriers, unsuitable habitat or a combination of these factors that could hinder movement (Avise et al., 1987). The long isolation of closely related populations could eventually lead to speciation by accumulating genetic changes over time. Studying and tracking the differences in species ecology, climatic conditions, and reproductive isolation between populations is essential to understand the formation of incipient species and plan appropriate conservation strategies (Moritz, 1994;Rundell and Price, 2009;Sobel et al., 2010). The mechanism of geographic isolation bringing about divergent evolution is well known in the birds, e.g., Darwin's finches in the subfamily Geospizinae (Sulloway, 1982) and five subspecies of the masked yellowthroat (Geothlypis aequinoctialis) in Central and South America (Curson et al., 1994). Among mammals too, four subspecies of Cheetah (Acinonyx jubatus) (Kitchener et al., 2017) and three subspecies of snow leopard (Panthera uncia) (Janecka et al., 2018) have evolved in different geographical areas. While speciation is an outcome of long isolation, shorter time-scale could set the populations of the same species on different evolutionary paths. The wolves in British Columbia showed strong genetic differentiation between adjacent populations, explained by habitat discontinuity between the coastal and inland regions (Muñoz-Fuentes et al., 2009).
The gray wolf (Canis lupus), one of the most widely distributed mammals and most studied species across the world (Mech and Boitoni, 2003), has had several populations worldwide that evolved in isolation. The wolves are considered one of the most resilient carnivores that have adapted and succeeded in a wide range of habitats from tundra to the deserts (Mech and Boitoni, 2003). Currently, at least forty-four wolf subspecies are described based on the variations in morphological features, behavioral aspects, and geographical distribution (Wozencraft, 2005). However, the status of many of these subspecies is uncertain and contested among scientists (Busch, 2018), making gray wolf taxonomy highly debatable (Shrotriya et al., 2012;Pilot et al., 2014;Ersmark et al., 2016). Several taxonomic resolutions are pending validation and revision at the species as well as subspecies level. Amid a muddled account of the number of species and sub-species, studies suggested that wolves from the Indian regions could form the oldest lineages (Aggarwal et al., 2003(Aggarwal et al., , 2007Sharma et al., 2004;Shrotriya et al., 2012;Werhahn et al., 2017Werhahn et al., , 2018Hennelly et al., 2021).
There are two major wolf populations in India, i.e., Woolly wolf (C. l. chanco) (also known as Tibetan or Himalayan wolf) and Indian wolf (C. l. pallipes). These wolves are supposed to have diverged from the main clade about 800,000 and 400,000 years ago, respectively (Sharma et al., 2004;Aggarwal et al., 2007;Shrotriya et al., 2012;Werhahn et al., 2017;Joshi et al., 2020). However, a recent study suggested a contradicting possibility that the Indian wolf could be basal to the Woolly wolf (Hennelly et al., 2021). The study stated that Indian and Woolly/Tibetan wolves shared the most common ancestors with Holarctic wolves 200,000 and 500,000 years, respectively (Hennelly et al., 2021). These differences in the year of divergence are mainly due to the selection of genes and the method used for the analysis. Isolated evolution of the Woolly wolf corresponds with rapid uplift of the Tibetan plateau and associated habitat modifications (Sun and Liu, 2000), which may have endowed them with the genetic adaptation for the cold and arid landscape of the Trans-Himalayas (Werhahn et al., 2018). The Indian wolf is believed to have evolved during the drier period of the Pleistocene and adapted to survive in the arid zones (Sharma et al., 2004;Singh and Kumara, 2006). After the extinction of Cheetah from India, it became the top carnivore species of the Indian open plains (e.g., semi-arid grasslands, scrublands, and grazing lands). Some of the studies proposed that the wolves of India were not subspecies but qualified as separate species viz. Canis himalayensis and Canis indica (Aggarwal et al., 2007). The debate regarding their taxonomic status is unsettled yet (Alvares et al., 2019). Nonetheless, both wolves could be presumed two distinct lineages (Pocock, 1941) and are geographically isolated, having allopatric distribution in India (Aggarwal et al., 2003). The Woolly wolf is found in the cold-arid zone of the upper Trans-Himalayas of India, covering the state of Himachal Pradesh and in two union territories, Ladakh and Jammu and Kashmir, with sightings recorded from Uttarakhand and Sikkim (Bhattacharya and Sathyakumar, 2010;Habib et al., 2013;Choudhury, 2015). The population estimates for both the lineages are not very accurate. As per the last available estimate, only 350 individuals of the Woolly wolf were found (Fox and Chundawat, 1995) which is a rough estimate for the Woolly wolf as it was assessed for the Ladakh and Spiti regions of India covering an area of 60,000 km 2 . The Indian wolf is distributed across 13 states of India, primarily in the semi-arid zone with the latest population estimate of 1200-1800 packs . Both the wolves inhabit open and grassland habitats and survive primarily outside protected areas. However, the Woolly wolf is functioning in far less anthropogenic pressure compared to the Indian wolf. Wolves in India are majorly found in non-protected areas with low natural prey populations (Habib, 2007;Maurya et al., 2011;Shrotriya et al., 2015;Habib et al., 2021b). Therefore, in most of their range in India, wolves primarily subsist on domestic livestock (Shahi, 1982;Mishra, 1997;Jethva and Jhala, 2004;Habib, 2007;Werhahn et al., 2019;Lyngdoh et al., 2020;Habib et al., 2021b) and are severely persecuted as a consequence. Although both these wolves are protected as Schedule I species under the Indian Wildlife (Protection) Act 1972. In spite of protection, conserving a species that resides in a human-dominated landscape and comes in direct conflict is a challenge.
The two lineages of wolves from India are unique and require focused conservation measures. Since both the lineages have remained isolated for a long period and functioning in two different bioclimatic zones, biological, ecological, and behavioral differences are expected to arise in these two populations. Understanding the differences is important because both subspecies are subject to management actions as they survive in a human-dominated landscape in India. In this study, we provide a comparative ecological perspective on spatial ecology, vocalization, food habits and cranial morphometry of these two lineages. The current study will help plan and implement better conservation and management actions. Moreover, information incorporating genetic analyses, morphology, and behavior using a comparative approach will help understanding key differences between the two lineages.

MATERIALS AND METHODS
We used primary as well as secondary data to elucidate differences in movement, space use, habitat use, howling characteristics, cranial morphology and diet.
Habitat Preferences, Movement, and Space Use Using Primary Data A total of 14 individuals (11 Indian wolves and three Woolly wolves) were captured, radio-collared, and monitored from 2015 to 2021 to understand the movement, space use, and habitat preferences of two wolf lineages ( Table 1). The Woolly wolves were captured in Himachal Pradesh and Indian wolves in Maharashtra, states of India. These wolves were captured using soft leg-hold traps. The traps (n = 25) were placed in a circular setup (radius = 70-80 cm), placed ∼20 cm apart from each other and tied with each other using rope (Habib, 2007). The wolf gland lure No. 100 (Stanley Hawbaker and Sons, Fort London, Pennsylvania) was used to attract the wolves toward traps. We also used fresh scats collected from different areas considering the scats were of different packs and placed around the traps to attract the wolf. The traps were connected to GSM-based device called "MinkPolice." This device operates with a magnetic switch and triggers a message on the registered email and phone numbers. Trapped wolves were held using double threaded nylon hockey net. The captured wolves were immobilized by injecting Ketamine-Xylazine (dosages based on the weight of the captured wolf) intramuscularly on their hind leg (Habib, 2007;Habib et al., 2021a). The immobilized wolves were fitted with 10 iridium, 3 Globalstar and 1 VHF collars.

Diet Diversity, Howling Parameters and Cranial Measurements Using Secondary Data
Data on the diet of both the wolf lineages were collected from published sources using specific keyword searches in Google scholar such as Tibetan wolf, Tibetan wolf diet, Himalayan wolf diet, Woolly wolf diet, Canis lupus chanco, Indian wolf diet, Indian gray wolf diet, Canis lupus pallipes. We found seven studies (two from India, three from Nepal, and two from Pakistan) to understand the Woolly wolf. For the Indian wolf diet, we found seven exclusive diet studies from Maharashtra, Gujarat, and Bihar States of India.
To understand the food habits of the Woolly and Indian wolves, we collected published data from different sources (detail provided in Supplementary Table 1). Different studies presented results in various forms, such as absolute frequency and relative frequency. To reduce the effect of study-specific variability, we first calculated the number of food items in the number of scats (Supplementary Tables 2, 3) and then calculated the relative frequency of occurrence (%RFO).
So far, there are two published studies on the howling parameter of the wolves from the Indian sub-continent (Hennelly et al., 2017;Sadhukhan et al., 2019). Hennelly et al. (2017) compared the vocalization of multiple wolf populations across the world, including the Woolly wolf and the Indian wolf and Sadhukhan et al. (2019) characterized the vocalization of the Indian wolves exclusively. Both the studies used nine parameters, such as mean frequency, maximum frequency, minimum frequency, frequency range, and end frequency, to understand the howling behavior of these two lineages. In these studies, the acoustic data were collected in the Trans-Himalayan region for Woolly wolf and Maharashtra for the Indian wolf.
The data on cranial morphometry of 14 Woolly and 20 Indian wolves were obtained from three studies, Allen (1938), Pocock (1941), and Srinivas and Jhala (2021). Data on total length (maximum length of skull from tip of rostrum to the nuchal crest); condyle basal length (distance from posterior projection of the occipital condyles to the anterior edge of pre-maxillary bones); zygomatic width (distance between outermost points of zygomatic arches), mandibular length (the distance from the jaw bone between the lower incisor of the anterior border of condyle); interorbital width (least distance between anterior orbits); postorbital width (least distance between posterior orbits) and maxillary width (distance between the central fosse of the right and left first maxillary molars), PM4 (Pre-molar 4) and M1 (molar 1) were compared to understand the differences between the two subspecies (Supplementary Figure 1). The details of studies used for comparative ecological perspectives are provided in Supplementary Table 1.

Data Analyses
The movement patterns of wolves were assessed using parameters such as daily average displacement (straight-line distance between two consecutive fixes) and average displacement during day and night. The varying inter-fix intervals were made uniform by post-processing all the data into an hourly data format (Abrahms, 2015;Habib et al., 2021a). Average displacement during the day and night were calculated by classifying locations. Average displacements of both the lineages were compared using the Wilcoxon rank sum test. Home ranges were calculated using minimum convex polygon (MCP) and Brownian bridge movement model (BBMM) (Horne et al., 2007). To understand the habitat preferences of Woolly and Indian wolves the landuse data for the states of Himachal Pradesh and Maharashtra were acquired from Bhuvan's open-source website (NRSA, 2016). 1 Habitat use and preference of the Woolly and Indian wolf lineages were analyzed using Manly's resource selection function (Manly et al., 2002). The design-III study framework of habitat use, where the animals are identified individually and both the use and the availability are measured for each one, was applied (Thomas and Taylor, 1990). The data on "use" was calculated from the number of GPS fixes and "availability" of land-use categories was calculated within the 100% MCP home range. All the measurements of movement parameters and analyses were carried out using ArcMET tool (Wall, 2014) in ArcGIS 10.6 and package "adehabitatLT, " version 0.3.25 (Calenge, 2015a) and animal movement tool (amt, version 0.1.3) (Signer et al., 2017) in program R 4.0.3 (R Core Team, 2021). Habitat use analysis was performed using the package "adehabitatHS, " version 0.3.15 in the program R, version 4.0.3 (Calenge, 2015b;R Core Team, 2021). For the comparison of cranial measurement of Woolly and Indian wolf we used two sampled unequal variance T-test and hedge's g test to evaluate the effect size. We combined the food items consumed in all the studies to represent the overall diet of wolves in the landscape and calculated the relative frequency of occurrence (RO) of each food item. We tabulated the RO of each item (number of occurrences of each prey item in scats/total number of scats × 100). We categorized the food items into wild and domestic prey and also in the body-weight classes (0-10 kg small prey; 11-70 kg medium-sized prey; >70 kg large prey) to compare the food habits. The details of studies used to comprehend the feeding habits of the two subspecies is provided in Supplementary Tables 2, 3.

Spatial Ecological Perspectives Based on the Movement Characteristics and Home Ranges of Woolly and Indian Wolves
A total of 47,240 fixes across 13 individuals (five males and eight females) of the wolves were analyzed ( Table 1) to examine daily movement, movement during day and night time, space use, and habitat use. The average hourly displacement of the Indian wolf (527.68 ± 111.52 m/hr) was significantly (W = 29, p = 0.01) higher than the Woolly wolf (338.54 ± 34.82 m/hr).
The mean hourly displacement during the day for Indian wolf (393.71 ± 172.43 m/hr) and Woolly wolf (316.70 ± 18.7 m/hr) was not significantly different (W = 12; p = 0.6). However, the mean hourly displacement during nighttime was 53% higher for the Indian wolf than the Woolly wolf and the result was significantly different (W = 28, p = 0.02) (Figure 1).

Habitat Ecology Based on the Land Use Preferences of Woolly and Indian Wolves
Resource selection function revealed a preferential selection for landuse categories by the Woolly wolf (Khi2L = 77.26, df = 6, p < 0.001) and the Indian wolf lineages (Khi2L = 520.3487, df = 35, p < 0.001). While riversides, marshes (valley habitat) (Wi = 6.44 ± 3.86SE) and scrub forests (Wi = 2.58 ± 3.39SE) were preferred by the Woolly wolf whereas the Indian wolf preferred grassland areas (Wi = 2.86 ± 0.37SE) and plantations (Wi = 2.71 ± 1.25SE). Builtup and agriculture areas were very less used in proportion to its availability by both the wolves ( Figure 4A). In addition, the Woolly wolf and the Indian wolf avoided snow cover and waterbodies, respectively ( Figure 4B).  Hennelly et al. (2017) suggested that the Indian wolf exhibited higher mean frequencies (593.5 Hz ± 211.4SE), wider frequency ranges (197.1 Hz ± 137.4SE), and longer duration (2.71 s ± 2.11SE) based on 117 howl recordings. On the other hand, 301 howls of the Woolly wolf characteristically had a lower mean frequency (428.5 Hz ± 125.7SE), shorter duration (2.56 s ± 1.68SE), and unmodulated frequency variation (101.8 Hz ± 107.1SE) ( Table 2). In contrast to Hennelly et al. (2017), the Indian wolf howl parameters reported by Sadhukhan et al. (2019) were broader in range and overlapped with the howling parameters of the Woolly wolf (

Comparative Food Habits of Woolly and Indian Wolf
We found seven studies exclusive on the diet of Woolly wolf covering the Himalayan region in India, Pakistan and Nepal. There were 20 different items (Supplementary Table 2) reported in Woolly wolf diet from 869 scats with 124.14 (± 190.94 SD) scats per study. The relative frequency of occurrences (RO) in scats for wild prey, domestic prey and vegetative matter were 41.93, 48.44, and 6.46%, respectively. The RO in scats of Woolly wolves for large, medium and small-sized prey were 32.32, 33.87, and 24.18%, respectively ( Figure 5A). About 70% of the Woolly diet consisted of goat, marmot, blue sheep, and pika, with a major contribution of cattle.
Another seven studies exclusively on the diet of Indian wolves from Rajasthan, Bihar and Maharashtra states of India revealed 17 different prey items from 6,877 scats with 982.42 (± 1363.67 SD) scats per study (Supplementary Table 3). The RO in scats for wild prey, domestic prey and vegetative matter was found 49.32, 40.34, and 9.31%, respectively. The Indian wolf diet primarily consisted of 67.95% of mediumsized prey followed by small (13.65%) and large-sized prey (7.15%) ( Figure 5B). 75% of food items consisted of blackbuck, goat, sheep and hare in the Indian wolf diet with the major contribution of livestock.

Spatial Ecological Differences Based on Home Range and Core Area Size in Woolly and Indian Wolf
The home range of Woolly wolf (95% MCP) was found eight times larger than the Indian wolf home range. Moreover, the extensive use area within the home range (50% MCP or core area) of the Woolly wolf was 17 times larger than the core area of the Indian wolf home range. The difference between the home ranges of these two lineages narrowed down in home range estimation by the advanced method of BBMM. The 95% isopleth of BBMM and 50% isopleth of BBMM of Woolly wolf were found only three times larger than those of the Indian wolf. The Woolly wolf resides in high-altitude rugged terrain areas with harsh climatic conditions. They usually used gentle slopes and valleys to move from one place to another (Lyngdoh, 2020), while MCP included cliff areas that the wolves could not use due to inaccessibility (Kenward et al., 2001;Silva et al., 2018). Different methods provided different information about home range size (Swihart and Slade, 1985;Gese et al., 1990;Silva et al., 2018). However, the MCP is widely used and can be useful to compare the home ranges of wolves from different studies. The BBMM provides reliable home range estimation, especially for central-place foragers and territorial carnivore species (Börger et al., 2008). This method is also reliable when the animals used several core areas or important travel corridors (Bullard, 1999;Horne et al., 2007;Kie et al., 2010). Moreover, we also compared the home ranges and core areas of Indian wolf (MCP) with the previous studies and found that the home ranges in the previous studies were smaller (148.49 km 2 ) (Jethva et al., 1997;Habib, 2007). The increase in the Indian wolf home ranges may probably be due to increased human pressure, habitat fragmentation, and prey depletion.
Large home ranges and core areas of the Woolly wolf could result from sparsely distributed prey species (Shrotriya et al., 2015) and suitable habitat continuity (Belongie, 2008). Wolves in mountainous regions selected low elevations, moderate slope and southwest aspect, probably to avoid snow cover and prey availability (Whittington et al., 2005). The home range for Mongolian wolves, which are closely related to the Woolly wolf, were among the largest home ranges (26,619 km 2 ) of the wolves worldwide (Kaczensky et al., 2008). Kaczensky et al. (2008) argue that long-range seasonal shift in the used area could dramatically increase the overall home-ranges. Therefore, the reason of large home range of Woolly wolf may due to the seasonal variation of space use. In contrast, the Indian wolf is surviving in patchy, fragmented habitat within high human-dominated landscape. In such landscapes where wild prey is almost absent, factors such as the type and distribution of food resources, human interference, and topography may have played a significant role in determining home-range size. The human disturbance could also be the reason for having multiple core areas in Indian wolf, as they use small undisturbed patches for resting, mating, and rearing pups (Habib, 2007).

Spatial Ecological Differences Based on Movement Characteristics of Woolly and Indian Wolf
The mean displacement (m/h) of the Indian wolves was 35.2% higher than the Woolly wolves' daily per hour movement. The daytime and nighttime displacement of the Indian wolves was 19.5 and 46.0% higher than the displacement of the Woolly wolf, respectively. The wolf is a widely distributed carnivore and its survival tactics involve a fine-tuned adaptation to local conditions. Their movement decisions are often the best functional compromise between finding food and avoiding humans (Ciucci et al., 1997). Some studies suggested that the wolf was primarily nocturnal in some areas, which allowed it to visit and move in intensive human-use areas without confrontation (Ciucci et al., 1997;Theurkauf et al., 2007). The Woolly wolf showed no difference between the day and night movement. In contrast, the Indian wolf movement aligned with the previous studies as their nighttime movement was almost twice the daytime movement. The no difference in the day night movement pattern of the Woolly wolf agreed with the findings of Vilà et al. (1995), who found that the wolf in Spain tended toward bimodal activity where nocturnality developed to cope up with human presence. Some studies suggested that wolves were most active during the night in summer (Fancy and Ballard, 1995) and during the daytime in winter (Mech, 1970).
As discussed above, the Woolly wolf home range was much larger than that of the Indian wolf. The functioning of the Woolly wolf in large home ranges but having less daily movement might be due to their use of valleys which makes the home range more linear, less compact, and larger. The higher movement of the Indian wolf in our study is in line with Daan (1981), who suggested a temporal shift in movement due to human disturbance. The other reason might be habitat fragmentation as the Indian wolves occupied fragmented areas, and the Woolly wolf inhibited continuous habitat. The results showed that the Indian wolf holds multiple core areas (Figure 2). Therefore, they need to travel more to find a secure area for resting and mating, as also supported by Ciucci et al., 1997. stated that the wolf activity mostly involved daily round-trip traveling from retreat areas to the core area of home ranges for some sort of social activity, which was also observed in previous studies (Carbyn, 1974).

Habitat Use and Preferences of Woolly and Indian Wolf
Wolves are known as grassland and openland species (Pocock, 1941;Mech, 1970). While the Indian wolves clearly showed a preference toward the grasslands and used the grassland patches more than its availability, the Woolly wolves did not show much preference toward grassland and openland (included in the barren land category) and used it according to the availability. This is due to the high availability of openland in continuity in the Trans-Himalayan region. The Woolly wolf showed an affinity to the riverside and marshes (valley habitat), which is energetically economic travel routes along with the refuge areas. Moreover, these areas are often also more sheltered from the harsh conditions inclusing cold wind. Avoidance of waterbodies by the Indian wolf could be better explained by the association of water presence with the proximity of humans in the peninsular Indian landscape. Both the wolves showed avoidance toward human-influenced areas such as agriculture and builtup more strongly so in the case of the Indian wolf. Paquet et al. (1996) found that mountain wolves did not have many options to avoid humans due to the selection of valley habitats which is often used by humans. The Indian wolf also faces higher human pressure than the Trans-Himalayan population because of the higher human population density in plains than in the Trans-Himalayan region (Mishra et al., 2009).

Differences in Howling Parameters of Woolly and Indian Wolf
Based on published literature, we compared the magnitude and pattern of howl acoustic structure to evaluate whether the longrange vocalizations showed acoustic differences between two wolf lineages in India. It is well understood that howl characteristics vary among different wolf subspecies (Kershenbaum et al., 2016). Due to the unique features of high amplitude and low frequency, a howl can travel for six kilometers or more and can be used to identify individuals for population estimation (Sadhukhan et al., 2021). Hennelly et al. (2017), compared the differences in the howls of the Himalayan wolf (aka Woolly wolf), Indian wolf, North African wolf (now African golden jackal; Sarabia et al., 2021), and Holarctic clades. The study found that small wolf subspecies, North African, Indian, and Israeli wolf had higher mean frequencies than the large wolf subspecies. Hence the body size could affect acoustic parameters.
In contrast, the study by Sadhukhan et al. (2019) reported lower mean frequencies for the Indian wolf. The difference between the Woolly wolf and the Indian wolf howls in Hennelly et al. (2017) could be an artifact of sub-sampling bias as Hennelly et al., 2017used 117 howl recordings. In contrast, Sadhukhan et al., 2019 Although both the studies sampled in the state of Maharashtra, yet their sampling years and locations varied. Both the studies sampled only a section of the Indian wolf population rendering the Indian wolf with a larger variation in mean frequencies overall.
Further, earlier studies have shown the similarities between the howls of various canid species and subspecies. Canis rufus showed similar howl type as coyotes Canis latrans, while European wolf (C. l. lupus) and Iberian wolf (C. l. signatus) howls were similar in signatures (Kershenbaum et al., 2016). The Indian wolf, Mackenzie Valley wolf (C. l. occidentalis) and Mexican wolf (C. l. baileyi) also showed similar howling signatures (Kershenbaum et al., 2016). The similarities in acoustic signatures between the two lineages may advocate the evolutionary history playing a role in howling behavior (Kershenbaum et al., 2016). We also second that canid howling is not an arbitrary signal but possesses species-specific information, reflecting that the adaptive processes of isolation or habitat features play a key role in howling behavior (Hennelly et al., 2017). Habitat and temporal variability might be the reason for different outcomes from the two studies. More studies are recommended to understand the differences in the howling parameters of these two lineages.

Differences in the Food Preferences of Woolly and Indian Wolf
Wolves are pack hunters and are known to feed on a variety of different food items. They choose their prey based on availability, abundance, pack stability, season, and habitat accessibility in the human-dominated landscapes (Imbert et al., 2016;Lyngdoh et al., 2020). Our study revealed that the diet of Woolly wolf from the Himalayan region consisted of 20 different food items from small birds, reptiles to large mammals and domestic animals such as cattle and yak. A review study on the Woolly wolf dietary spectrum across their global distribution, including Tibet and China, reported 39 different food items . The energy requirement of large carnivores makes them prone to conflict with humans as they need a large prey base (Carbone et al., 1999;Mech and Boitoni, 2003). The Woolly wolf also consumed a sufficient proportion of the large-sized prey (32.32%) and medium-sized prey (33.87%) with a considerable quantity of small prey (24.18%) consumption in their diet (Figure 5). Various studies across the wolf distribution range confirm that large prey forms the major part of the wolf diet (Imbert et al., 2016;Mengulluoglu et al., 2019;Petridou et al., 2019;Sin et al., 2019). The dependency on small prey might be because of the scarcity of the large and medium-sized animals to avoid interactions with humans or to gain and fulfill energy requirements in the harsh climatic condition. The Woolly wolf is also heavily dependent upon the domestic prey items in the Himalayan region. Livestock (Yak, Dzo cow, Goat, and Sheep) were the most consumed mammals in the Woolly wolf diet compared to sparsely distributed wild prey. Hence, it is not surprising that Woolly wolf showed a marginal preference for domestic prey than wild prey. Many studies from the Trans-Himalayas accounted that low abundance of wild prey and poor management of livestock resulted in depredation by snow leopard and Woolly wolf (Jackson and Ahlborn, 1984;Mishra, 1997;Namgail et al., 2007;Anwar et al., 2012;Suryawanshi et al., 2013;Chetri et al., 2017).
Despite the Woolly wolf 's significant domestic animal consumption, they are being tolerated by the locals in some regions (Rangarajan, 2001;Lyngdoh et al., 2020). At the same time, they are persecuted in many areas for the same reason (Mishra, 1997). Apart from domestic animal consumption, major wild prey species were Asiatic ibex (Capra sibirica), Urial (Ovis vignei), Tibetan wild ass (Equus kiang), Tibetan argali (Ovis ammon hodgsoni) and Blue sheep (Pseudois navaur) (Anwar et al., 2012;Subba, 2012;Ahmed et al., 2017;Bocci et al., 2017;Chetri et al., 2017;Werhahn et al., 2019;Lyngdoh et al., 2020;Habib et al., 2021b). However, their combined contribution was as low as 16.77% of the total food items. It could be related to the small population size and sparse distribution of the ungulate species in the region (Shrotriya et al., 2015).
The Indian wolf also primarily preyed upon domestic animals, followed by wild prey items. The Indian wolf consumed a considerable amount of vegetative material (fruits and plant items, etc.), which was absent in the Woolly wolf diet. The absence of vegetative matter in the diet of the Woolly wolf may be due to the scarcity of wild fruiting plants in the Trans-Himalayan region or missed reporting in the studies due to difficulties in identification in the studies we reviewed in this study. The Indian wolf primarily preyed upon medium-sized mammals relating to the availability of blackbuck, chinkara, and especially livestock such as sheep and goat. The Indian wolf-bearing states (Rajasthan, Andhra Pradesh, Maharashtra, Telangana) have the highest goat and sheep population (20th Livestock Census: All India Report, 2019). The average pack size of the Indian wolf varied from 1.5 to 4.7 individuals in the breeding and non-breeding season (Kumar, 1998). The preference toward medium-sized prey could be due to the smaller pack size and body size of the Indian wolf. Small packs would require less food and find it difficult to prey upon large prey species. The studies used to understand the food habits were from different protection regimes, such as the protected area of Velvadar Blackbuck Sanctuary (Jhala, 1991) and the human-dominated landscape of Maharashtra (Habib, 2007). The diet of the wolves from Velvadar Blackbuck Sanctuary was dominantly comprised of wild prey base (91.8%), whereas Maharashtra wolves' diet was dominated by domestic prey (47.8%). This shows that the wolf may prefer wild ungulates over domestic prey, depending upon the availability. Consequences of human-wolf conflict due to low prey abundance and unavailability can hinder conservation measures. Therefore, it is essential to address prey restoration and livestock security to reduce conflict and achieve better conservation management for wolves in the Himalayas and in the plains of India.

Cranial Morphometric Differences Between Woolly and Indian Wolf
The Woolly wolf exists at 3,900-5,600 m elevation across the lowoxygen region of the Himalayas Werhahn et al., 2018). Survival in such a landscape requires the Woolly wolf to face metabolic challenges such as severe oxidative stress and increased metabolic rates (Beall, 2007;Hassanin et al., 2009). Several studies have conducted a genetic analysis of the Woolly wolf, identified the genes facilitating their adaptation to cope with hypoxia (Zhang et al., 2014;Werhahn et al., 2018;Wang et al., 2020). However, their morphological adaptations against hypoxic conditions have not been paid much attention yet. Butaric and Klocke (2018) studied the adaptation of upper respiratory structures to hypoxic and cold dry air in humans occupying high and low altitudes. These adaptations in skull structures help in increased uptake and air conditioning processes. The Woolly wolf exhibited longer (total length) and broader (post and interorbital width) skulls than the Indian wolf. Hence, their skull size and structure could be an adaptation to meet respiratory demand of more oxygen (Butaric and Klocke, 2018). However, the wolves are known to conform Bergmann's rule, that is the wolves in the northern latitude are generally larger in body size (Meiri et al., 2007).

Taxonomic Dilemma
The systematics of wolves from the Indian subcontinent is less studied, remains controversial and confusing regarding their taxonomic status as sub-species and species. Indian wolves gathered the interest of the scientific community for their unique evolutionary history and uncertain taxonomic status. Although the Indian wolf was first described by Hodgson (1847), a consensus on its nomenclature is yet to be reached after 175 years with several attempts made to clear their taxonomy (Aggarwal et al., 2003;Sharma et al., 2004;Werhahn et al., 2017;Alvares et al., 2019;Joshi et al., 2020;Wang et al., 2020Wang et al., , 2021. Hodgson (1847) described the Himalayan wolf as separate species, C. laniger, noting its appearances. Later, Blanford (1888) rejected Hodgson's proposal and combined C. laniger with C. lupus, and elevated the Indian wolf to C. pallipes. However, after 50 years, Pocock (1941) described both the taxa as subspecies of C. lupus. In Pocock's scheme, C. pallipes became C. l. pallipes, and C. laniger merged with C. l. chanco. Later on, more studies suggested the uniqueness of these two taxa and identified them as the oldest lineage of wolves but did not reach any conclusion (Sharma et al., 2004;Aggarwal et al., 2007;Werhahn et al., 2017). A recently conducted study of the Woolly wolf by Joshi et al., 2020 found no evidence for C. l. chanco to be distinct species. They suggested the acceptance of Woolly wolf as C. l. chanco and not as C. langier or C. himalayensis. Their findings were additionally supported by Wang et al. (2020), who found that wolves across the Himalayas and Tibetan plateau are closely related. The Woolly wolf adapted to survive in a low oxygen environment and the Indian wolf represents two of the most endangered wolf populations (Joshi et al., 2020;Hennelly et al., 2021;Wang et al., 2021). Therefore, the wolves from India have been identified as unique and qualified as an important population and proposed as Evolutionary Significant Units (ESUs) due to their distinct evolution and adaptation (Hennelly et al., 2021).

CONCLUSION
Identifying the ecological and behavioral differences in closely related species provides understanding in the evolutionary process of speciation and helps identify a species or subspecies (Arnegard et al., 2010;Ramasindrazana et al., 2011). The literature clearly highlights the uniqueness of these two wolves from India (Sharma et al., 2004;Aggarwal et al., 2007;Shrotriya et al., 2012;Joshi et al., 2020). The genetic study strongly suggested that both the subspecies were distinct and the most ancient lineages. These lineages did not show any genetic admixture or geographic overlap with other wolves from rest of the world (Aggarwal et al., 2003(Aggarwal et al., , 2007Sharma et al., 2004). Geographical isolation and differential habitat selection of closely related adjacent wolf populations is the leading mechanism of the evolution of genetically and ecologically different subspecies (Leonard, 2014), for example, Mexican wolf, North American wolf, Italian, Iberian and Scandinavian wolf (Wayne et al., 1991;Vilà et al., 1999Vilà et al., , 2003Lucchini et al., 2004). The two geographically non-overlapping wolf subspecies from India also showed genetic as well as ecological differences exhibiting their evolutionary divergence. In this study, we found a clear difference in the spatial ecology of both the wolves.
Wolves are considered a typical grassland species in Asia and Europe (Mech, 1970). Both the lineages in India primarily choose grassland or openlands. The Woolly wolf lives in a landscape with vast openlands, hence the habitat type did not show up in preferential analysis. Plantations by the Indian wolf and marshes/riversides and shrubs by the Woolly wolf were preferred as potential refuge. The wolves are known to live in variety of habitats and are considered ecological tolerant animals (Jedrzejewski et al., 2004). The habitat use by wolves can be influenced by many factors, e.g., habitat type, prey availability and anthropogenic pressure (Meriggi et al., 1996;Ciucci and Boitoni, 1998;Mech and Boitoni, 2003). We postulate that the ecological differences between Woolly wolf and Indian wolf could be mainly due to their functioning in entirely different habitats. The home range of the Woolly wolf was significantly larger than that of the Indian wolf. The differences in home range sizes and movement patterns could be because of the Trans-Himalayan region having low prey base, less anthropogenic pressure and connected suitable patches compared to the Peninsular India having disturbed landscape with patchily distributed suitable habitats. However, to understand the ranging pattern and habitat use of Woolly wolf, we used only three individual data collared from Spiti region of trans-Himalayas. More data from different regions could help us better understanding of wolf ecology as the space use may vary with different regions based on prey base and habitat/landscape characteristics. There was no significant difference in the daytime and nighttime movement of the Woolly wolf, whereas the Indian wolf traversed more during the nighttime. Variation in the level of human disturbance could be the reason for Indian wolf being more nocturnal than Woolly wolf.
Both the lineages had significant differences in their cranial measurement, and we found that the Woolly wolf has longer and broader skulls than the Indian wolf. Moreover, the data presented in two studies on the vocalization characteristics were conflicting and did not show a clear picture of the difference. The canid howling may vary species-wise or depend upon the environment, still it possesses species-specific information which may reflect adaptive and or neutral processes of isolation (Kershenbaum et al., 2016). The studies used in this review showed contrasting results. Hennelly et al. (2017) suggested a clear difference but the results of Sadhukhan et al., 2019 suggested that howling parameters of both the lineages overlapped. The differences in their feeding ecology occurred primarily due to the availability of landscape-specific prey species. Both the wolves of India depended on the livestock for more than 50% of their diet. The wolves across the world have been reported to feed on a wide variety of food items from animal matter to vegetative matter. Their main prey in most of the areas are large and mediumsized prey depending upon the availability (Jhala, 1991;Meriggi et al., 1996;Jethva and Jhala, 2004;Chavez and Gese, 2005;Stahler and Smith, 2006;Habib, 2007;Hosseini-Zavarei et al., 2013;Newsome et al., 2016). Nevertheless, since these populations are geographically isolated, genetic, morphological data and ecological requirements suggest apparent distinctiveness. The Indian and Woolly/Tibetan wolves shared the most common ancestors with Holarctic wolves 0.2 mya (0.17-0.3 mya) and 0.5 mya (0.38-0.64 mya), respectively (Hennelly et al., 2021). These two lineages diverged and become incipient species with genetic distinctiveness a long time ago. Their geographic isolation is persistent and would facilitate the behavioral and ecological changes to intensify over time. The wolves of Asia are paid less academic attention compared to their counterparts in Europe and America. This study sheds light on the ecological and behavioral differences of the two oldest wolf lineages of the world found in India. We further suggest detailed morphological analysis and further studies should be conducted to understand the in-depth differences in ecological requirements of the subspecies.

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

ETHICS STATEMENT
The animal study was reviewed and approved by the Wildlife Institute of India, Animal Ethics Committee.

AUTHOR CONTRIBUTIONS
SK collected data, performed analysis, reviewed all the potential manuscript and reports that were included in the meta-analysis, and drafted the manuscript. SL and SoS collected data. ShS guided in analysis. BH, ShS, SG, and SL reviewed the manuscript. BH conceptualized and acquired funding of the study. All authors contributed to the article and approved the final version of the manuscript.

ACKNOWLEDGMENTS
We are thankful to the Ministry of Environment, Forest and Climate Change (MoEFCC), Himachal Pradesh and Maharashtra Forest Department for the permissions to collar animals and conduct research. Department of Science and Technology (DST) Govt. of India, MBZ Species Conservation Grant, WII Grant in Aid and MFD are duly acknowledged for funding provided to carry out the research. We are also grateful to the Field Directors, DFOs, RFOs, and especially our field assistants for providing help to track collared animal. We thank the Director, Dean, and Research Coordinator of the Wildlife Institute of India for supporting the study. We thank Zehidul Hussain and Hussain Saifi Reshamwala for their valuable comments and suggestions to improve the manuscript.