Genomic analysis of wolves from Pakistan clarifies boundaries among three divergent wolf lineages.

Among the three main divergent lineages of gray wolf (Canis lupus), the Holarctic lineage is the most widespread and best-studied, particularly in North America and Europe. Less is known about Tibetan (also called Himalayan) and Indian wolf lineages in southern Asia, especially in areas surrounding Pakistan where all three lineages are thought to meet. Given the endangered status of the Indian wolf in neighboring India and unclear southwestern boundary of the Tibetan wolf range, we conducted mitochondrial and genome-wide sequencing of wolves from Pakistan and Kyrgyzstan. Sequences of the mitochondrial D-loop region of 81 wolves from Pakistan indicated contact zones between Holarctic and Indian lineages across the northern and western mountains of Pakistan. Reduced-representation genome sequencing of 8 wolves indicated an east-to-west cline of Indian to Holarctic ancestry, consistent with a contact zone between these two lineages in Pakistan. The western boundary of the Tibetan lineage corresponded to the Ladakh region of India's Himalayas with a narrow zone of admixture spanning this boundary from the Karakoram Mountains of northern Pakistan into Ladakh, India. Our results highlight the conservation significance of Pakistan's wolf populations, especially the remaining populations in Sindh and Southern Punjab that represent the highly endangered Indian lineage.

Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids.

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.1 Although dogs of variable size are found in the archeological record,2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.5 Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFß signaling, and skeletal formation.6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.8 The identification of a functional mutation associated with IGF1 has thus far proven elusive.6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.

Spatial Distribution and Dietary Niche Breadth of Leopard Cats (Prionailurus bengalensis) Inhabiting Margalla Hills National Park, Pakistan.

The leopard cat (Prionailurus bengalensis) is distributed throughout the Himalayan foothills in Pakistan and occurs in moist temperate and dry coniferous forests. However, the cat species is categorized as "Data Deficient" in Pakistan. In the current study we aimed to investigate the leopard cat distribution and dietary niche in the Margalla Hills National Park, the lowest part of the Himalayan foothills in Pakistan. We recorded direct (field observations, camera trapping) and indirect signs (scats) of the species on 23 trails /tracks in the study area. The leopard cat was recorded at 13 different sampling sites in the park, with an altitudinal range between 664-1441 m asl. The diet composition of the species was investigated through scat analysis, with the species identity of the scats confirmed through the mitochondrial cytB region. The species' diet comprised both animal and plant matter. The animal-based diet (in terms of frequency of occurrence) contained wild (51.75%) as well as domestic prey (7.69%), and plant species (31.47%). Wild prey included small mammals (rodents, two mongoose species, Asian palm squirrel, and Cape hare), birds, insects and snails. The domestic prey were poultry birds, sheep, goats and dogs. Consumption of wild prey was higher in summer (n = 31 scat), whereas intake of domestic prey was higher in winter (n = 37 scats). The dietary niche breadth was wider 14.84 in summer than winter 10.67. A chi-square test showed a significant difference in seasonal dietary intake of the leopard cat. The study concludes that the leopard cat feeds mainly on wild animal prey and plant species; however, in winter consumption of domestic prey increases.

Ancient divergence of Indian and Tibetan wolves revealed by recombination-aware phylogenomics.

The grey wolf (Canis lupus) expanded its range across Holarctic regions during the late Pleistocene. Consequently, most grey wolves share recent (<100,000 years ago) maternal origins corresponding to a widespread Holarctic clade. However, two deeply divergent (200,000-700,000 years ago) mitochondrial clades are restricted, respectively, to the Indian subcontinent and the Tibetan Plateau, where remaining wolves are endangered. No genome-wide analysis had previously included wolves corresponding to the mitochondrial Indian clade or attempted to parse gene flow and phylogeny. We sequenced four Indian and two Tibetan wolves and included 31 additional canid genomes to resolve the phylogenomic history of grey wolves. Genomic analyses revealed Indian and Tibetan wolves to be distinct from each other and from broadly distributed wolf populations corresponding to the mitochondrial Holarctic clade. Despite gene flow, which was reflected disproportionately in high-recombination regions of the genome, analyses revealed Indian and Tibetan wolves to be basal to Holarctic grey wolves, in agreement with the mitochondrial phylogeny. In contrast to mitochondrial DNA, however, genomic findings suggest the possibility that the Indian wolf could be basal to the Tibetan wolf, a discordance potentially reflecting selection on the mitochondrial genome. Together, these findings imply that southern regions of Asia have been important centers for grey wolf evolution and that Indian and Tibetan wolves represent evolutionary significant units (ESUs). Further study is needed to assess whether these ESUs warrant recognition as distinct species. This question is especially urgent regarding the Indian wolf, which represents one of the world's most endangered wolf populations.

Pleistocene origins, western ghost lineages, and the emerging phylogeographic history of the red wolf and coyote.

The red wolf (Canis rufus) of the eastern US was driven to near-extinction by colonial-era persecution and habitat conversion, which facilitated coyote (C. latrans) range expansion and widespread hybridization with red wolves. The observation of some grey wolf (C. lupus) ancestry within red wolves sparked controversy over whether it was historically a subspecies of grey wolf with its predominant "coyote-like" ancestry obtained from post-colonial coyote hybridization (2-species hypothesis) versus a distinct species closely related to the coyote that hybridized with grey wolf (3-species hypothesis). We analysed mitogenomes sourced from before the 20th century bottleneck and coyote invasion, along with hundreds of modern amplicons, which led us to reject the 2-species model and to investigate a broader phylogeographic 3-species model suggested by the fossil record. Our findings broadly support this model, in which red wolves ranged the width of the American continent prior to arrival of the grey wolf to the mid-continent 60-80 ka; red wolves subsequently disappeared from the mid-continent, relegated to California and the eastern forests, which ushered in emergence of the coyote in their place (50-30 ka); by the early Holocene (12-10 ka), coyotes had expanded into California, where they admixed with and phenotypically replaced western red wolves in a process analogous to the 20th century coyote invasion of the eastern forests. Findings indicate that the red wolf pre-dated not only European colonization but human, and possibly coyote, presence in North America. These findings highlight the urgency of expanding conservation efforts for the red wolf.

The relative importance of phylogeny and habitat in determining the presence and prominence of a granula iridica in hooved mammals.

PURPOSE: To investigate the relationship between phylogeny and amount of shade in a species' habitat regarding the presence or absence of an iridal granula iridica (GI) in a large sample of Artiodactyl and Perissodactyl clades and using online resources. METHODS: The Comparative Ocular Pathology Laboratory of Wisconsin (COPLOW) archives were searched for glass slide material from Artiodactyl (even-toed) and Perissodactyl (odd-toed) ungulates. The slides were examined, and the presence or absence of the GI was noted. The phylogenetic tree of the ungulate species was inferred using TimeTree (http://www.timetree.org), and the habitat data are derived from Animal Diversity Web (https://animaldiversity.org/). We assessed the probability of the presence of GI occurring given the amount of shade in a species' environment using phylogenetic logistic regression. RESULTS: Forty-eight artiodactyl species were able to be evaluated and tabulated. Nine perissodactyl species were able to be evaluated. The phylogenetic logistic regression showed that the probability of GI presence was lower in artiodactyl species that inhabited shaded environments (ßshaded  = -1.774). Arctiodacyl species inhabiting a nonshaded environment were slightly more probable to have the GI present (ßnonshaded  = 0.023), with species inhabitating ambiguously shaded environments having a high probability of GI presence (ßambiguous  = 2.214). CONCLUSIONS: Our results suggest that the GI may be a common morphological feature to shade the pupil in nonshaded environments, and, in its absence, increase the amount of light reaching the retina to improve vision in shaded environments for hooved mammals. Further research on the functional optics of the GI and studies that include additional ungulate species would further elucidate phylogenetic and ecological factors influencing the occurrence of GI in hooved mammals.

Characterising the harmonic vocal repertoire of the Indian wolf (Canis lupus pallipes).

Vocal communication in social animals plays a crucial role in mate choice, maintaining social structure, and foraging strategy. The Indian grey wolf, among the least studied subspecies, is a social carnivore that lives in groups called packs and has many types of vocal communication. In this study, we characterise harmonic vocalisation types of the Indian wolf using howl survey responses and opportunistic recordings from captive and nine packs (each pack contains 2-9 individuals) of free-ranging Indian wolves. Using principal component analysis, hierarchical clustering, and discriminant function analysis, we found four distinct vocalisations using 270 recorded vocalisations (Average Silhouette width Si = 0.598) which include howls and howl-barks (N = 238), whimper (N = 2), social squeak (N = 28), and whine (N = 2). Although having a smaller body size compared to other wolf subspecies, Indian wolf howls have an average mean fundamental frequency of 422 Hz (±126), which is similar to other wolf subspecies. The whimper showed the highest frequency modulation (37.296±4.601) and the highest mean fundamental frequency (1708±524 Hz) compared to other call types. Less information is available on the third vocalisation type, i.e. 'Social squeak' or 'talking' (Mean fundamental frequency = 461±83 Hz), which is highly variable (coefficient of frequency variation = 18.778±3.587). Lastly, we identified the whine, which had a mean fundamental frequency of 906Hz (±242) and is similar to the Italian wolf (979±109 Hz). Our study's characterisation of the Indian wolf's harmonic vocal repertoire provides a first step in understanding the function and contextual use of vocalisations in this social mammal.

Howl variation across Himalayan, North African, Indian, and Holarctic wolf clades: tracing divergence in the world's oldest wolf lineages using acoustics.

Vocal divergence within species often corresponds to morphological, environmental, and genetic differences between populations. Wolf howls are long-range signals that encode individual, group, and subspecies differences, yet the factors that may drive this variation are poorly understood. Furthermore, the taxonomic division within the Canis genus remains contended and additional data are required to clarify the position of the Himalayan, North African, and Indian wolves within Canis lupus. We recorded 451 howls from the 3 most basal wolf lineages-Himalayan C. lupus chanco-Himalayan haplotype, North African C. lupus lupaster, and Indian C. lupus pallipes wolves-and present a howl acoustic description within each clade. With an additional 619 howls from 7 Holarctic subspecies, we used a random forest classifier and principal component analysis on 9 acoustic parameters to assess whether Himalayan, North African, and Indian wolf howls exhibit acoustic differences compared to each other and Holarctic wolf howls. Generally, both the North African and Indian wolf howls exhibited high mean fundamental frequency (F0) and short duration compared to the Holarctic clade. In contrast, the Himalayan wolf howls typically had lower mean F0, unmodulated frequencies, and short howls compared to Holarctic wolf howls. The Himalayan and North African wolves had the most acoustically distinct howls and differed significantly from each other and to the Holarctic wolves. Along with the influence of body size and environmental differences, these results suggest that genetic divergence and/or geographic distance may play an important role in understanding howl variation across subspecies.