Using Bayesian stable isotope mixing models to estimate wolf diet in a multi-prey ecosystem

Land cover and space use influence coyote carnivory: evidence from stable-isotope analysis

For many species, the relationship between space use and diet composition is complex, with individuals adopting varying space use strategies such as territoriality to facilitate resource acquisition. Coyotes (Canis latrans) exhibit two disparate types of space use; defending mutually exclusive territories (residents) or moving nomadically across landscapes (transients). Resident coyotes have increased access to familiar food resources, thus improved foraging opportunities to compensate for the energetic costs of defending territories. Conversely, transients do not defend territories and are able to redirect energetic costs of territorial defense towards extensive movements in search of mates and breeding opportunities. These differences in space use attributed to different behavioral strategies likely influence foraging and ultimately diet composition, but these relationships have not been well studied. We investigated diet composition of resident and transient coyotes in the southeastern United States by pairing individual space use patterns with analysis of stable carbon (δ13C) and nitrogen (δ15N) isotope values to assess diet. During 2016-2017, we monitored 41 coyotes (26 residents, 15 transients) with GPS radio-collars along the Savannah River area in the southeastern United States. We observed a canopy effect on δ13C values and little anthropogenic food in coyote diets, suggesting 13C enrichment is likely more influenced by reduced canopy cover than consumption of human foods. We also observed other land cover effects, such as agricultural cover and road density, on δ15N values as well as reduced space used by coyotes, suggesting that cover types and localized, resident-like space use can influence the degree of carnivory in coyotes. Finally, diets and niche space did not differ between resident and transient coyotes despite differences observed in the proportional contribution of potential food sources to their diets. Although our stable isotope mixing models detected differences between the diets of resident and transient coyotes, both relied mostly on mammalian prey (52.8%, SD = 15.9 for residents, 42.0%, SD = 15.6 for transients). Resident coyotes consumed more game birds (21.3%, SD = 11.6 vs 13.7%, SD = 8.8) and less fruit (10.5%, SD = 6.9 vs 21.3%, SD = 10.7) and insects (7.2%, SD = 4.7 vs 14.3%, SD = 8.5) than did transients. Our findings indicate that coyote populations fall on a feeding continuum of omnivory to carnivory in which variability in feeding strategies is influenced by land cover characteristics and space use behaviors.

The mercury flow through a terrestrial songbird food chain in subtropical pine forest: Elucidated by Bayesian isotope mixing model and stable mercury isotopes

In order to comprehend the transfer of inorganic mercury (IHg) and methylmercury (MeHg) within food chains in terrestrial pine forests, we collected samples of Great Tit nestlings, common invertebrates, plants, and soil in a subtropical pine forest and used Bayesian isotope mixing model analysis, Hg daily intake, and stable Hg isotopes to elucidate the flow of MeHg and IHg in these food chains. Results indicate that caterpillars and cockroaches are the predominant prey items for nestlings, accounting for a combined contribution of 81.5%. Furthermore, caterpillars, cockroaches, and spiders were found to contribute the most (∼80%) of both IHg and MeHg that dietary accumulated in nestlings. The provisoned invertebrates tend to supply more IHg and diluting the proportion of MeHg as total Hg (MeHg%). Notably, nestling feathers displayed the highest Δ199Hg values but a relatively lower MeHg%, suggesting an imbalanced incorporation of Hg from maternal transfer and dietary accumulation during the nestling stage. This study highlights the efficacy of nestlings as indicators for identifying Hg sources and transfers in avian species and food chains. However, caution must be exercised when using Hg isotope compositions in growing feathers, and the contribution of maternally transferred Hg should not be ignored.

Temporal variation in translocated Isle Royale wolf diet

Wolves (Canis lupus) can exert top-down pressure and shape ecological communities through the predation of ungulates and beavers (Castor spp.). Therefore, understanding wolf foraging is critical to estimating their ecosystem-level effects. Specifically, if wolves are consumers that optimize tradeoffs between the cost and benefits of prey acquisition, changes in these factors may lead to prey-switching or negative-density dependent selection with potential consequences for community stability. For wolves, factors affecting cost and benefits include prey vulnerability, risk, reward, and availability, which can vary temporally. We described the wolf diet by the frequency of occurrence and percent biomass and characterized the diet using prey remains found in wolf scats on Isle Royale National Park, Michigan, USA, during May-October 2019 and 2020. We used logistic regression to estimate prey consumption over time. We predicted prey with temporal variation in cost (availability and/or vulnerability) such as adult moose (Alces alces), calf moose, and beaver (Castor canadensis) to vary in wolf diets. We analyzed 206 scats and identified 62% of remains as beaver, 26% as moose, and 12% as other species (birds, smaller mammals, and wolves). Adult moose were more likely to occur in wolf scats in May when moose are in poor condition following winter. The occurrence of moose calves peaked during June-mid-July following birth but before calf vulnerability declined as they matured. By contrast, beaver occurrence in wolf scat did not change over time, reflecting the importance of low-handling cost prey items for recently introduced lone or paired wolves. Our results demonstrate that the wolf diet is responsive to temporal changes in prey costs. Temporal fluctuation in diet may influence wolves' ecological role if prey respond to increased predation risk by altering foraging or breeding behavior.

African wild dogs (Lycaon pictus) show differences in diet composition across landscape types in Kruger National Park, South Africa

The Kruger National Park (KNP) is home to the last genetically viable, minimally managed population of African wild dogs (Lycaon pictus, wild dogs) in South Africa. Until 2004, this population remained stable, but since has been declining. In this study, we aimed to improve our understanding of the ecology of KNP wild dogs by estimating the relative contribution of different prey types to their diet across landscape types. Based on a Bayesian mixing model, we assessed wild dog diet and foraging preferences using stable isotope analysis. We sampled 73 individuals from 40 packs found in six different landscape types. In thickets, packs predominantly prey on small browsing and mixed-feeding species (accounting for ~73% of their diet), but occasionally hunt large grazers (~24%) and large browsers (~3%). In open landscape types where lions (Panthera leo) are more or less absent, such as in the Lowveld sour bushveld, wild dogs prey on large browsers and large grazers (~67%). Our results demonstrate that KNP wild dogs occupy a broader ecological niche than previously thought, with small browsers forming an integral part of their diet. We also present the first data describing differences in wild dog diet–tissue discrimination factors for tail hair and whiskers compared to respective stable nitrogen (δ15N) and carbon (δ13C) values obtained from feces of captive wild dogs, as well as from those of South Africa’s broader managed metapopulation. While these data should be considered preliminary, we suggest that until wild dog diet–tissue discrimination factors are calculated through a controlled feeding study, the discrimination factors calculated for the gray wolf (Canis lupus) should be used for wild dog-related isotope studies, rather than the often cited values for red foxes (Vulpes vulpes).

Habitat selection by wolves and mountain lions during summer in western Montana

In the Northern Rockies of the United States, predators like wolves (Canis lupus) and mountain lions (Puma concolor) have been implicated in fluctuations or declines in populations of game species like elk (Cervus canadensis) and mule deer (Odocoileus hemionus). In particular, local distributions of these predators may affect ungulate behavior, use of space, and dynamics. Our goal was to develop generalizable predictions of habitat selection by wolves and mountain lions across western Montana. We hypothesized both predator species would select habitat that maximized their chances of encountering and killing ungulates and that minimized their chances of encountering humans. We assessed habitat selection by these predators during summer using within-home range (3^rd order) resource selection functions (RSFs) in multiple study areas throughout western Montana, and tested how generalizable RSF predictions were by applying them to out-of-sample telemetry data from separate study areas. Selection for vegetation cover-types varied substantially among wolves in different study areas. Nonetheless, our predictions of 3^rd order selection by wolves were highly generalizable across different study areas. Wolves consistently selected simple topography where ungulate prey may be more susceptible to their cursorial hunting mode. Topographic features may serve as better proxies of predation risk by wolves than vegetation cover-types. Predictions of mountain lion distribution were less generalizable. Use of rugged terrain by mountain lions varied across ecosystem-types, likely because mountain lions targeted the habitats of different prey species in each study area. Our findings suggest that features that facilitate the hunting mode of a predator (i.e. simple topography for cursorial predators and hiding cover for stalking predators) may be more generalizable predictors of their habitat selection than features associated with local prey densities.

Spatial and temporal variability in summer diet of gray wolves (Canis lupus) in the Greater Yellowstone Ecosystem

The role of predation by large carnivores in suppressing prey populations and structuring ecosystems is highly debated, calling for a detailed understanding of carnivore diets. Wolves (Canis lupus) roam across three continents and persist throughout widely different ecosystems. Their diet is flexible and may vary spatially as well as seasonally, which requires analysis of diet on different spatial and temporal scales. Few studies have investigated the summer diet of wolves, which is more variable, consists of smaller prey, and requires different methods than studying their winter diet. To better understand the summer diet of wolves, we combined three independently collected wolf scat data sets from three distinctly different portions of the Greater Yellowstone Ecosystem: Yellowstone National Park (2009), Grand Teton National Park (2003 – 2009), and the Absaroka-Beartooth Wilderness (2009 – 2010). These areas represent different ecological conditions and management regimes, which may impact wolf diet. We estimated relative biomass and compared occurrence of different prey species among packs, years, as well as the three regions. In total, we analyzed 1,906 wolf scats and found that neonate cervids, adult elk, and adult deer were the most important prey species in the summer diet of the wolves. We found dietary variation among packs residing in the same area, as well as across years. The occurrence of neonate cervids displayed the most variation, and low occurrence of this prey type often was associated with a more diverse diet. Wolf packs within the national parks had a higher occurrence of medium-sized prey (~ 50 – 70 kg) and lower occurrence of small-sized prey (≤ 20 kg) compared to wolves in the Absaroka-Beartooth Wilderness. These results demonstrate flexibility in summer diet across packs, years, and between regions within the Greater Yellowstone Ecosystem.

Experimental removals reveal dietary niche partitioning facilitates coexistence between native and introduced species

Niche overlap between native species and ecologically similar invaders can lead to competitive exclusion of threatened native species, but if two such species also co-occur naturally elsewhere, interactions between native and introduced populations may mirror coevolved niche partitioning that reduces competition and promotes coexistence.A single, insular population of Fremont's squirrel (Tamiasciurus fremonti) the Mount Graham red squirrel (MGRS; T. f. grahamensis) in the Pinaleño Mountains, Arizona, USA, is critically endangered and resource competition with introduced Abert's squirrels (Sciurus aberti) may threaten its long-term persistence. The species are naturally synoptic in other mountain sites, and both consume diets comprised primarily of conifer seeds and fungi.We conducted experimental removals of introduced Abert's squirrels and used stable isotope analysis of diets before and after removals, and of diets in naturally syntopic populations to test the hypothesis that dietary niche partitioning can facilitate coexistence between native and introduced species. We also developed a novel approach to determine the influence of fluctuating food availability on carbon enrichment in consumers.Mount Graham red squirrels and introduced Abert's squirrels partitioned the dietary niche similarly to naturally syntopic populations. Removals had no apparent effect. Diet of MGRS was more closely linked to availability of resources than to presence of Abert's squirrels.Flexible dietary niche of introduced Abert's squirrels may have allowed them to exploit a resource opportunity in syntopy with MGRS. Variable food production of MGRS habitat may intensify competition in poor years, and territorial defense against non-native Abert's squirrels likely imposes fitness costs on individual MGRS. Similarity in our model species' diets may make MGRS more vulnerable to competition if climate change eliminates the advantages of larder-hoarding. Where introduced populations of ecologically similar species are better adapted to changing conditions, they may ultimately replace southern peripheral populations of native species.

The problems with pooling poop: confronting sampling method biases in wolf (Canis lupus) diet studies

Wolf (Canis lupus L., 1758) diet is commonly estimated via scat analysis. Several researchers have concluded that scat collection method can bias diet estimates, but none of these studies properly accounted for interpack, age class, and temporal variability, all of which could bias diet estimates. We tested whether different scat collection methods yielded different wolf diet estimates after accounting for these other potential biases. We collected scats (n = 2406) monthly from four packs via three scat collection methods (at home sites, at clusters of GPS locations, and opportunistically) in and adjacent to Voyageurs National Park, Minnesota, USA, during April–October 2015. Diet estimates were not affected by scat collection method but did vary temporally, among packs, and by age class. To more accurately estimate wolf population diets, researchers should collect 10–20 adult scats/pack per month from home sites and (or) opportunistically from packs that are representative of the population of interest. Doing so will minimize the potential biases associated with temporal, interpack, and age-class variability.

Dietary niche partitioning among black bears, grizzly bears, and wolves in a multiprey ecosystem

Identifying mechanisms that promote coexistence of sympatric species is important for predicting ecological effects of anthropogenic change. Many caribou (Rangifer tarandus (L., 1758)) populations are declining, and it is unclear to what extent sympatric predators consume caribou or how alternative prey affect caribou–predator relationships. We used stable isotope mixing models to estimate diets of black bear (Ursus americanus Pallas, 1780), grizzly bear (Ursus arctos L., 1758), and grey wolves (Canis lupus L., 1758) during early, middle, and late summer of 2009–2010 in northwestern British Columbia, Canada. Although we expected wolf diet to be primarily composed of moose (Alces alces (L., 1758)) — as they exist at twice the density of caribou — wolf diet consisted principally of caribou, and to a lesser extent moose and beaver (Castor canadensis Kuhl, 1820), with little change occurring throughout summer. Black bear diet consisted mainly of vegetation and moose, shifting from moose to vegetation through summer. Grizzly bear diet consisted primarily of vegetation and moose, and did not change throughout summer. Our results demonstrate the role of dietary niche partitioning in bear and wolf coexistence, and that caribou may be primary prey for wolves in an ecosystem with relatively high moose abundance and low human development.

Alternative prey use affects helminth parasite infections in grey wolves

Predators affect prey populations not only through direct predation, but also by acting as definitive hosts for their parasites and completing parasite life cycles. Understanding the affects of parasitism on prey population dynamics requires knowing how their predators' parasite community is affected by diet and prey availability. Ungulates, such as moose (Alces americanus) and white-tailed deer (Odocoileus virginianus), are often important prey for wolves (Canis lupus), but wolves also consume a variety of alternative prey, including beaver (Castor canadensis) and snowshoe hare (Lepus americanus). The use of alternative prey, which may host different or fewer parasites than ungulates, could potentially reduce overall abundance of ungulate parasites within the ecosystem, benefiting both wolves and ungulate hosts. We examined parasites in wolf carcasses from eastern Manitoba and estimated wolf diet using stable isotope analysis. Taeniidae cestodes were present in most wolves (75%), reflecting a diet primarily comprised of ungulates, but nematodes were unexpectedly rare. Cestode abundance was negatively related to the wolf's δ(13) C value, indicating diet affects parasite abundance. Wolves that consumed a higher proportion of beaver and caribou (Rangifer tarandus), estimated using Bayesian mixing models, had lower cestode abundance, suggesting the use of these alternative prey can reduce parasite loads. Long-term consumption of beavers may lower the abundance of adult parasites in wolves, eventually lowering parasite density in the region and ultimately benefiting ungulates that serve as intermediate hosts. Thus, alternative prey can affect both predator-prey and host-parasite interactions and potentially affect food web dynamics.

Experimentally derived δ¹³C and δ¹⁵N discrimination factors for gray wolves and the impact of prior information in Bayesian mixing models

Stable isotope analysis of diet has become a common tool in conservation research. However, the multiple sources of uncertainty inherent in this analysis framework involve consequences that have not been thoroughly addressed. Uncertainty arises from the choice of trophic discrimination factors, and for Bayesian stable isotope mixing models (SIMMs), the specification of prior information; the combined effect of these aspects has not been explicitly tested. We used a captive feeding study of gray wolves (Canis lupus) to determine the first experimentally-derived trophic discrimination factors of C and N for this large carnivore of broad conservation interest. Using the estimated diet in our controlled system and data from a published study on wild wolves and their prey in Montana, USA, we then investigated the simultaneous effect of discrimination factors and prior information on diet reconstruction with Bayesian SIMMs. Discrimination factors for gray wolves and their prey were 1.97‰ for δ13C and 3.04‰ for δ15N. Specifying wolf discrimination factors, as opposed to the commonly used red fox (Vulpes vulpes) factors, made little practical difference to estimates of wolf diet, but prior information had a strong effect on bias, precision, and accuracy of posterior estimates. Without specifying prior information in our Bayesian SIMM, it was not possible to produce SIMM posteriors statistically similar to the estimated diet in our controlled study or the diet of wild wolves. Our study demonstrates the critical effect of prior information on estimates of animal diets using Bayesian SIMMs, and suggests species-specific trophic discrimination factors are of secondary importance. When using stable isotope analysis to inform conservation decisions researchers should understand the limits of their data. It may be difficult to obtain useful information from SIMMs if informative priors are omitted and species-specific discrimination factors are unavailable.

The conflict between cheetahs and humans on Namibian farmland elucidated by stable isotope diet analysis

Large areas of Namibia are covered by farmland, which is also used by game and predator species. Because it can cause conflicts with farmers when predators, such as cheetahs (Acinonyx jubatus), hunt livestock, we assessed whether livestock constitutes a significant part of the cheetah diet by analysing the stable isotope composition of blood and tissue samples of cheetahs and their potential prey species. According to isotopic similarities, we defined three isotopic categories of potential prey: members of a C4 food web with high δ15N values (gemsbok, cattle, springhare and guinea fowl) and those with low δ15N values (hartebeest, warthog), and members of a C3 food web, namely browsers (eland, kudu, springbok, steenbok and scrub hare). We quantified the trophic discrimination of heavy isotopes in cheetah muscle in 9 captive individuals and measured an enrichment for 15N (3.2‰) but not for 13C in relation to food. We captured 53 free-ranging cheetahs of which 23 were members of groups. Cheetahs of the same group were isotopically distinct from members of other groups, indicating that group members shared their prey. Solitary males (n = 21) and males in a bachelor groups (n = 11) fed mostly on hartebeest and warthogs, followed by browsers in case of solitary males, and by grazers with high δ15N values in case of bachelor groups. Female cheetahs (n = 9) predominantly fed on browsers and used also hartebeest and warthogs. Mixing models suggested that the isotopic prey category that included cattle was only important, if at all, for males living in bachelor groups. Stable isotope analysis of fur, muscle, red blood cells and blood plasma in 9 free-ranging cheetahs identified most individuals as isotopic specialists, focussing on isotopically distinct prey categories as their food.

Mercury in gray wolves (Canis lupus) in Alaska: increased exposure through consumption of marine prey

Mercury (Hg) bioaccumulates in the tissues of organisms and biomagnifies within food-webs. Gray wolves (Canis lupus) in Alaska primarily acquire Hg through diet; therefore, comparing the extent of Hg exposure in wolves, in conjunction with stable isotopes, from interior and coastal regions of Alaska offers important insight into their feeding ecology. Liver, kidney, and skeletal muscle samples from 162 gray wolves were analyzed for total mercury (THg) concentrations and stable isotopic signatures (δ(13)C, δ(15)N, and δ(34)S). Median hepatic THg concentrations were significantly higher in wolves with coastal access compared to wolves from interior Alaska. Stable isotope ratios, in conjunction with THg concentrations, provide strong evidence that coastal wolves are utilizing marine prey representing several trophic levels. The utilization of cross-ecosystem food resources by coastal wolves is clearly contributing to increased THg exposure, and may ultimately have negative health implications for these animals.