Different coexistence patterns between apex carnivores and mesocarnivores based on temporal, spatial, and dietary niche partitioning analysis in Qilian Mountain National Park, China
eLife assessment
This paper provides an important assessment of competition dynamics allowing coexistence of the carnivore guild within a large national park in China. Multiple surveying techniques (camera traps and DNA metabarcoding) provide convincing evidence that spatial segregation represents the main strategy of coexistence, while species have a certain degree of temporal and dietary overlap. Altogether, the manuscript provides information critical to the conservation and management agenda of the park.
https://doi.org/10.7554/eLife.90559.3.sa0Important: Findings that have theoretical or practical implications beyond a single subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Carnivores play key roles in maintaining ecosystem structure and function as well as ecological processes. Understanding how sympatric species coexist in natural ecosystems is a central research topic in community ecology and biodiversity conservation. In this study, we explored intra- and interspecific niche partitioning along spatial, temporal, and dietary niche partitioning between apex carnivores (wolf Canis lupus, snow leopard Panthera uncia, Eurasian lynx Lynx lynx) and mesocarnivores (Pallas’s cat Otocolobus manul, red fox Vulpes vulpes, Tibetan fox Vulpes ferrilata) in Qilian Mountain National Park, China, using camera trapping data and DNA metabarcoding sequencing data. Our study showed that apex carnivore species had more overlap temporally (coefficients of interspecific overlap ranging from 0.661 to 0.900) or trophically (Pianka’s index ranging from 0.458 to 0.892), mesocarnivore species had high dietary overlap with each other (Pianka’s index ranging from 0.945 to 0.997), and apex carnivore and mesocarnivore species had high temporal overlap (coefficients of interspecific overlap ranging from 0.497 to 0.855). Large dietary overlap was observed between wolf and snow leopard (Pianka’s index = 0.892) and Pallas’s cat and Tibetan fox (Pianka’s index = 0.997), suggesting the potential for increased resource competition for these species pairs. We concluded that spatial niche partitioning is likely to key driver in facilitating the coexistence of apex carnivore species, while spatial and temporal niche partitioning likely facilitate the coexistence of mesocarnivore species, and spatial and dietary niche partitioning facilitate the coexistence between apex and mesocarnivore species. Our findings consider partitioning across temporal, spatial, and dietary dimensions while examining diverse coexistence patterns of carnivore species in Qilian Mountain National Park, China. These findings will contribute substantially to current understanding of carnivore guilds and effective conservation management in fragile alpine ecosystems.
Introduction
Carnivores play key roles in maintaining ecosystem structure and function, as well as ecological processes (Ripple et al., 2014). Apex carnivores are often classically specialized hunters occupying top trophic positions that work to suppress the number of herbivores and mesocarnivores through predation, competition, and trophic cascades (Ripple et al., 2014; Newsome et al., 2017; Ritchie and Johnson, 2009). Understanding how sympatric carnivore species coexist in natural ecosystems is a central focus in community ecology and biodiversity conservation (Chu et al., 2017). Human activities and climate change are causing large-scale habitat loss and fragmentation, posing significant challenges to carnivore interspecies relationships (Li et al., 2021; Manlick and Pauli, 2020; Ripple et al., 2014). Such relationships are becoming increasingly intricate, leading to greater research efforts aimed at determining the strategies that facilitate coexistence despite intense competition (Kuijper et al., 2016; Smith et al., 2017).
The competitive exclusion principle dictates that species with similar ecological requirements are unable to successfully coexist (Hardin, 1960; Gause, 1934). Thus, carnivores within a guild occupy different ecological niches based on a combination of three niche dimensions, i.e., spatial, temporal, and trophic (Schoener, 1974). Spatially, carnivore species within the same geographical area exhibit distinct distributions that minimize overlap in resource use and competition. For example, carnivores can partition habitats based on habitat feature preferences and availability of prey (de Satgé et al., 2017; Garrote and Pérez de Ayala, 2019; Goldyn et al., 2003; Strampelli et al., 2023). Temporally, differences in seasonal or daily activity patterns among sympatric carnivores can reduce competitive interactions and facilitate coexistence. For example, carnivores can exhibit temporal segregation in their foraging behaviors, such as diurnal versus nocturnal activity, to avoid direct competition (Finnegan et al., 2021; Nasanbat et al., 2021; Searle et al., 2021). Trophically, carnivore species can diversify their diets to exploit different prey species or sizes, thereby reducing competition for food resources. For example, carnivores can exhibit dietary specialization to optimize their foraging efficiency and minimize competitive pressures (Steinmetz et al., 2021). Currently, research on the niche partitioning of sympatric carnivores primarily focuses on one or two dimensions, with little attention paid to all three (Alexander et al., 2016; Karanth et al., 2017; Li et al., 2022a; Santos et al., 2019; Shao et al., 2021; Strampelli et al., 2023; Tsunoda et al., 2020; Vilella et al., 2020). This limitation restricts the multi-scale understanding of coexistence mechanisms among carnivores.
In recent years, camera trapping and DNA metabarcoding technology has been widely used in wildlife monitoring and research. Camera trapping enables the monitoring of elusive species in remote areas (Boitani, 2016; Palencia et al., 2022), and has become a helpful tool for accumulating large amounts of time-recorded data that can provide detailed information surrounding animal behavior and activity patterns, which is widely used to investigate species-interspecies relationships (Frey et al., 2017). For example, Li et al., 2019, used camera trap data to conclude that temporal segregation is a key mechanism for promoting the coexistence of tigers (Panthera tigris) and leopards (Panthera pardus). DNA metabarcoding provides a noninvasive molecular tool that is more accurate than traditional dietary analysis methods, such as microhistology (Deiner et al., 2017; Newmaster et al., 2013). DNA metabarcoding technology has been applied to many studies and has offered profound insight into the ecology, conservation, and biological monitoring of rare and endangered species (Deagle et al., 2019; Kartzinel et al., 2015). For example, Shao et al., 2021, found that dietary niche partitioning promoted the coexistence of species in the mountains of southwestern China based on DNA metabarcoding. The combination of camera trapping and DNA metabarcoding can result in robust data for exploring the mechanisms of coexistence surrounding carnivore guilds.
The Qilian Mountains constitute a biodiversity hotspot with one of the richest carnivore assemblages in the world (Di Minin et al., 2016). To understand the competition and coexistence among different carnivore species in this area, we explored the habitat use, activity patterns, and prey item composition of sympatric carnivore species comprised of apex carnivores and mesocarnivores across Qilian Mountain National Park using camera trap data and DNA metabarcoding data. Based on theories surrounding resource partitioning and niche differentiation, as well as studies on interactions and coexistence among carnivorous species (Haswell et al., 2018; Linnell and Strand, 2000), we hypothesized that differentiation along one or more niche axes is beneficial for the coexistence of the carnivorous guild in the Qilian Mountains. We expected that spatial niche differentiation promotes the coexistence of large carnivores in the Qilian Mountain region, as they are more likely than small carnivores to spatially avoid interspecific competition (Davis et al., 2018). Mesocarnivores may coexist either spatially or temporally due to increased interspecific competition for similar prey (Di Bitetti et al., 2010; Donadio and Buskirk, 2006). Dietary niche differentiation may be a significant factor for promoting coexistence between large and mesocarnivore species due to differences in body size (Gómez-Ortiz et al., 2015; Lanszki et al., 2019).
Results
Sympatric carnivore identification
Of the 480 scat samples sequenced, those which had no sequencing data were inconclusive, consisted of non-target species, or host species with low sample sizes (one Asian badger [Meles leucurus] and four upland buzzard [Buteo hemilasius]) were removed (Figure 3—figure supplement 1). The remaining 404 scat samples were composed of three apex carnivores (49 wolf, 147 snow leopard, 19 Eurasian lynx) and three mesocarnivores (63 Pallas’s cat, 87 red fox, 39 Tibetan fox).
Spatial distribution difference and overlap
A total of 322 camera trap sites were surveyed after relocating infrared cameras that did not capture any target carnivore species. A total of three cameras were considered to have failed due to loss. We analyzed data from 319 camera sites and obtained 14,316 independent detections during a total effort of 37,192 effective camera trap days. We recorded wolf in 26 sites, snow leopard in 109 sites, Eurasian lynx in 36 sites, red fox in 92 sites, and Tibetan fox in 34 sites. However, the camera detection rates of Pallas’s cat were too low to analyze species occupancy and daily activity patterns.
According to the findings derived from single-season, single-species occupancy models, the snow leopard demonstrated a notably higher probability of occupancy overall compared to other carnivore species, estimated at 0.437 (Table 1). Conversely, the Eurasian lynx exhibited a lower occupancy probability, estimated at 0.161. Further analysis revealed that the occupancy probabilities of the wolf and Eurasian lynx declined with increasing normalized difference vegetation index (ndvi) (Table 2, Figure 1). Additionally, wolf occupancy probability displayed a negative relationship with roughness index and a positive relationship with prey availability. Snow leopard occupancy probabilities exhibited a negative relationship with distance to roads and ndvi. In contrast, both red fox and Tibetan fox demonstrated a positive relationship with distance to roads. Moreover, red fox occupancy probability increased with higher human disturbance and greater prey availability.
The detection probabilities of wolf, snow leopard, red fox, and Tibetan fox exhibited an increase with elevation (Table 2). Moreover, there was a positive relationship between the detection probability of Tibetan fox and prey availability. The detection probabilities of snow leopard and Eurasian lynx declined as human disturbance increased.
The Sørensen similarity index (S) ranged from 0.1 to 0.5 (Table 3). Compared with other combinations of apex-mesocarnivore pairs, snow leopard and red fox (S=0.477) had relatively high spatial overlap, while Eurasian lynx and Tibetan fox (S=0.198) had the lowest spatial overlap. Moreover, spatial overlap of apex versus apex carnivores and mesocarnivores versus mesocarnivores was relatively low.
Daily activity patterns and differences
A total of 1444 independent records were obtained for five carnivore species, consisting of 79 records of wolf, 458 records of snow leopard, 126 records of Eurasian lynx, 421 records of red fox, and 359 records of Tibetan fox. Among apex carnivores (Table 3, Figure 2), the daily activity was similar between snow leopard and Eurasian lynx and their diel activity overlap was close to 1 (Δ4=0.900, p=0.285), their daily activity peak was at 21:00 hr and dawn. However, the wolf had a significantly different daily activity pattern with snow leopard (Δ4=0.676, p<0.001) and Eurasian lynx (Δ4=0.661, p<0.001), and its daily activity peak happened around 9:00 and 18:00 hr. Tibetan fox and red fox had different activity patterns and peaks (p<0.001). The activity peak for red fox activity peaked at 3:00 and 21:00 hr, while the Tibetan fox had a prolonged active bout between noon and dusk. Temporal activity patterns between apex carnivores and mesocarnivores were significantly different, except for wolf and Tibetan fox (Δ4=855, p=0.118).
Dietary composition, diversity, and similarity
A total of 26 unique prey species were identified from 9 taxonomic orders (Figure 3, Supplementary file 1a). Artiodactyla and lagomorpha were the most frequently detected in the diets of apex carnivores and mesocarnivores, accounting for 32.81% and 70.18% of prey counts, respectively (Figure 3—figure supplement 2). Blue sheep made up 26.50% of prey counts in apex carnivore diet, while plateau pika made up 67.11% of prey counts in mesocarnivore diet. Livestock were present in 17.98% of apex carnivore diet counts and were present in 4.82% mesocarnivore diet counts.
The dietary niche overlap among all carnivore species can be found in Table 3. Wolf and snow leopard had the highest dietary niche overlap value among apex carnivores (Ojk = 0.892). The value of Pianka’s index was generally low between apex carnivores and mesocarnivores, except wolf and red fox (Ojk = 0.811). In contrast, observed dietary overlap was greatest among the mesocarnivores, especially Pallas’s cat and Tibetan fox, with a value of 0.997.
Red fox had the greatest richness of prey with a value of 16, while Pallas’s cat and Tibetan fox had the lowest diversity of prey with richness value of 6 (Supplementary file 1b). Dietary similarity was assessed using inversed Jaccard’s distances. Diets were most similar between wolf and snow leopard with a value of 0.588 and least similar between Pallas’s cat and Tibetan fox with a value of 0.200. All other pairs fell between the values of 0.2 and 0.5 (Supplementary file 1c).
Discussion
Our study addresses, for the first time, the coexistence patterns of carnivore species present in Qilian Mountain National Park across multiple dimensions of niche partitioning. This work substantially contributes to current understanding of carnivore guilds and offers helpful information for biodiversity conservation at a regional scale. Moreover, our study provides important insights into the potential mechanisms of niche partitioning among sympatric carnivores, particularly intra- and interspecific relations between apex carnivores and mesocarnivores. Specifically, we found that the overall trend of spatial overlap across carnivores is relatively low, that apex carnivores overlap in time and diet, that mesocarnivores showed a high degree of dietary overlap, and that there was substantial similarity in diel activity patterns between apex carnivores and mesocarnivores. These results indicate that carnivores with similar ecological traits foster co-occurrence by adjusting their daily activity patterns and using differing food resources to minimize competitive interactions.
We found dietary and temporal overlap among apex carnivores, showing that spatial partitioning is responsible for their successful coexistence in this area, consistent with our hypothesis. Wolf and snow leopard had the highest dietary overlap and prey similarity between apex carnivore pairs in this study, showing that the avoidance of space and time plays an important role in their coexistence. Recent evidence proves that habitat preference facilitates the coexistence of wolf and snow leopard (Shrotriya et al., 2022). Their hunting strategies may be impacted by their habitat selection. Solitary snow leopards are more suitable for hiding in habitats with features that favor ambush predators, while wolves hunt in packs (Shrotriya et al., 2022). It is clear that wild ungulates (e.g. blue sheep) constituted the primary proportion of wolf and snow leopard diet, followed by small mammals such as plateau pika, Himalayan marmot (Marmota himalayana), and woolly hare (Lepus oiostolus). In addition, livestock consumption also contributed to the high degree of overlap in their diets (Wang et al., 2014). This supports the optimal foraging theory, in which large predators preferentially select food resources that provide maximum benefit (Brown et al., 1999), but also showed that greater competition for resources is likely to occur between wolf and snow leopard due to their use of the same prey species in this area. This may be especially true in times of habitat stress when resources are poor.
Snow leopard and Eurasian lynx had the highest temporal overlap between apex carnivore pairs in our case, showing that spatial and dietary partitioning facilitate their coexistence. Eurasian lynx is considered an opportunistic predator, and its prey varies among different regions with its primary dietary resource being ungulates and small mammals. For example, Eurasian lynx showed a strong preference for brown hare Lepus europaeus in Turkey, edible dormice Glis glis in Slovenia and Croatia, and chamois Rupicapra rupicapra or roe deer Capreolus capreolus in Switzerland (Mengüllüoğlu et al., 2018; Krofel et al., 2011; Molinari-Jobin et al., 2007). Varied prey selection may be related to sex, age, population density, and season (Mengüllüoğlu et al., 2018; Odden et al., 2006). Our results show that woolly hare make up the majority of the Eurasian lynx diet, followed by blue sheep. Woolly hare is mainly distributed in shrubland, meadow, desert, and wetland, while blue sheep tend to choose highly sheltered areas, close to bare rocks and cliffs as habitat. Therefore, prey preferences among snow leopard and Eurasian lynx also contribute to spatial avoidance.
Mesocarnivores had substantial overlap in diet, underscoring food resources as a primary competitive factor that necessitates spatial and temporal partitioning for successful coexistence. This finding aligns with previous research indicating that mesocarnivores use temporal and spatial segregation to reduce competition and the probability of antagonistic interspecific encounters (Ferreiro-Arias et al., 2021; Li et al., 2022a). Differences in habitat preference may lead to spatial niche partitioning among mesocarnivores (Wang et al., 2022). In addition, species can adjust temporal periods of behavior to respond to environmental change, competition, and predation risk (Gallo et al., 2022; Finnegan et al., 2021; van der Vinne et al., 2019). Pallas’s cat hunts during crepuscular and diurnal periods and inhabits meadow habitat with greater prey abundance (Anile et al., 2021; Greco et al., 2022; Ross et al., 2019). In contrast, red fox is primarily nocturnal and occupies diverse habitats depending on prey abundance (Goldyn et al., 2003; Pandolfi et al., 1997; Reshamwala et al., 2022). Tibetan fox is a diurnal hunter of the Tibetan plateau, preferring shrub meadow, meadow steppe, and alpine meadow steppe (Gong and Hu, 2003). It is worth noting the substantial overlap in diet between Pallas’s cat and Tibetan fox. The dietary overlap between the two was extremely high, with a Pianka’s value close to 1. Dietary analyses showed that pika contributed to more than 85% of their collective diets, with 90% of Pallas’s cat diet being pika. Pika may be an optimal prey item in the area because of size and year-round activity (Ross, 2009). Previous studies have shown that the Pallas’s cat and Tibetan fox are specialist predators of pikas (Harris et al., 2014; Ross, 2009). However, specialization on pika is facultative in that Pallas’s cat and Tibetan fox can select other prey items when pika availability is low (Harris et al., 2014; Ross, 2009). This was observed in our study, even though dietary diversity was low.
Apex carnivores and mesocarnivores exhibited considerable overlap in time overall, showing that spatial and dietary partitioning may play a large role in facilitating their coexistence. As confirmed by previous research, kit foxes (Vulpes macrotis) successfully coexisted with dominant carnivores by a combination of spatial avoidance and use of alternative resources (Lonsinger et al., 2017). Differences in body mass may play a crucial role in minimizing dietary overlap, effectively reducing interspecific competition between apex and mesocarnivores. Of exception in our study, however, was wolf and red fox, who exhibited more dietary overlap, showing that temporal and spatial avoidance may promote their coexistence. As canid generalist-opportunist species, the wolf and red fox consumed similar prey, albeit the red fox may have obtained livestock and ungulate species via scavenging or by preying on very young individuals (Hacker et al., 2022). In addition, the wolf and red fox had different peak activity times, suggesting temporal segregation as a potentially strong driver of coexistence. Our occupancy model also revealed that mesocarnivores prefer habitats that coincide with human disturbance, possibly contributing to spatial niche differentiation with other apex carnivore species like the wolf. Our findings echo other recent research which supports the idea that red foxes can coexist with wolves by exploiting a broader ecological niche (Shrotriya et al., 2022).
Our study highlights the effectiveness of combining camera trapping with DNA metabarcoding for detecting and identifying both cryptic and rare species within a sympatric carnivore guild. This integrated approach allowed us to capture a more comprehensive view of species presence and interactions compared to traditional visual surveys, whereas it is important to acknowledge the challenges associated with this technique, including the high costs of equipment and the need for specialized training and computational resources to manage and analyze the large volumes of sequence data. Despite these challenges, the benefits of this combined method in improving biodiversity assessments and understanding species coexistence outweigh the drawbacks. However, several restrictions remain for this research. The first limitation involves differences in samples sizes. Although the scat samples of Tibetan fox were relatively low, the accuracy of DNA metabarcoding in informing species presence in diet ensures that data are informative and thus still important for species conservation management decisions (Hacker et al., 2022). Second, the methodology of foraging (e.g. predation or scavenging) and the condition of the prey item (e.g. age or size) cannot be identified in dietary studies (Hacker et al., 2022). Pika is a prime component of diet among mesocarnivores, especially in the diet of the Pallas’s cat and the Tibetan fox. We surmise that the simultaneous dependence on pika led to partial overlap in spatial and activity patterns, resulting in increased potential competitive interactions. Due to the lack of spatial and temporal analysis of Pallas’s cat in our study, further monitoring is needed to develop a comprehensive conservation plan. Despite these limitations, our study provides a foundation from which future studies interested in niche partitioning among carnivores along spatial, temporal, and dietary dimensions can be modeled.
In summary, our study has shown that the coexistence of carnivore species in the landscapes of Qilian Mountain National Park can be facilitated along three niche axes, with spatial segregation being the most pronounced of the three. Apex carnivore species tended to overlap temporally or trophically, mesocarnivore species had high dietary overlap with each other, and apex carnivore and mesocarnivore species displayed similarity in temporal use. Pika, blue sheep, and livestock were found to make up a large proportion of carnivore diet. Resource competition between wolf and snow leopard and the interspecific competition between Pallas’s cat and Tibetan fox were strong in this area. Based on the findings presented above, we recommend targeted efforts for enhanced protection and management in several key areas. First, greater efforts are needed to protect habitat, including the establishment of habitat corridors and optimizing grassland fence layouts to safeguard migration routes. Second, resource competition should be carefully monitored between snow leopards and wolves, as well as between Pallas’s cat and Tibetan foxes. More attention is needed for pika at our study site considering the role they play in the conservation of Pallas’s cat and Tibetan fox populations. Previous poisoning campaigns targeting pikas were initiated out of concern for grassland degradation (Smith and Foggin, 1999). Recent research underscores the importance of pika population health and habitat for the distribution of Pallas’s cats (Greenspan and Giordano, 2021). It’s crucial to note that poisoning campaigns targeting small mammals may be incredibly dangerous for mesocarnivores feeding on them due to secondary poisoning, necessitating caution in such conservation strategies. Further, vigilance regarding the preponderance of smaller prey in predator diets is vital, as this may indicate severe loss of larger prey, which will increase the risk of interference competition (Steinmetz et al., 2021). Third, efforts should focus on restoring vulnerable wild prey populations, strengthening grazing area management, and supporting livelihoods of herders to mitigate livestock predation and consequent human-wildlife conflicts. Regular patrols of protected areas and core habitats are essential, along with educational outreach to conservation staff and herders. Our study corroborates and complements the findings of prior studies on these species and their coexistence mechanisms, offering insights crucial for wildlife conservation in the region.
Materials and methods
Study sites
Request a detailed protocolThe Qilian Mountains laterally span Gansu and Qinghai Provinces in China, located on the northeastern edge of the Qinghai-Tibetan Plateau (Figure 4). Qilian Mountain National Park covers an area of approximately 52,000 km2, with an average elevation of over 3000 m. The area is an alpine ecosystem with a typical plateau continental climate. The average annual temperature is below –4°C and the average annual rainfall is about 400 mm, with habitats mainly consisting of deserts, grassland, meadows, and wetland (Zheng, 2011). Wildlife present include the wolf (Canis lupus), snow leopard (Panthera uncia), Eurasian lynx (Lynx lynx), red fox (Vulpes vulpes), Tibetan fox (Vulpes ferrilata), Tibetan brown bear (Ursus arctos), Chinese mountain cat (Felis bieti), wild yak (Bos mutus), blue sheep (Pseudois nayaur), alpine musk-deer (Moschus chrysogaster), Tibetan antelope (Pantholops hodgsonii), Himalayan marmot (M. himalayana), woolly hare (L. oiostolus), and plateau pika (Ochotona curzoniae), among others (Ma et al., 2021; Xue et al., 2019).
Camera trap monitoring and noninvasive sampling
Request a detailed protocolThe study area was subdivided into sample squares of 25 km2 (5 km×5 km) as a geographical reference for placing camera survey sites and collecting scat samples (Xue et al., 2019). Species occurrence was recorded using camera trap monitoring (Model Ltl-6210; Shenzhen Ltl Acorn Electronics Co. Ltd., Shenzhen, China). Two camera traps were placed in each square with a distance of at least 1 km between them. However, due to limitations of terrain, landform, road accessibility, and other factors, the number of camera trap in some squares was adjusted in accordance with field conditions. Camera traps were set in areas believed to be important to and heavily used by wildlife, such as the bottoms of cliffs, sides of boulders, valleys, and ridges along movement corridors. Taking into account the fact that mammalian communities are sensitive to seasonality, we used camera traps to monitor animals with an extensive survey effort from December 2016 to February 2022, covering the activity of animal species in different seasons, which can reflect the overall distribution of carnivores. We placed a total of 280 infrared cameras at the study site (Figure 4), set them to be active for 4–6 months, and considered possible relocation to another position based on animal detection in an effort to improve estimates of the occupancy and detection rates for both common and rare species (Kays et al., 2020). The camera trap was set to record the time and date on a 24 hr clock when triggered, and to record a 15 s video and 1 photo with an interval of 2 min between any two consecutive triggers. The sum of camera trap effective days was defined by the total amount of trapping effort during the sampling period, which was calculated from the time the camera was placed in operation to the time the last video or photograph was taken. We visited each camera two to three times a year to download photos and check batteries.
Noninvasive sampling of scats believed to be of carnivore origin were collected within camera trapping areas. A small portion of scat (approximately 1/3) was broken off and stored in a 15 ml centrifuge tubes with silica desiccant covered by clean filter paper to separate the desiccant from the scat (Janečka et al., 2008). Gloves were replaced between sampling to avoid cross-contamination. Sampling place, date, and sample number were labeled on the tube; GPS coordinates, elevation, and nearby landscape features were recorded on the sample collection sheet (Hacker et al., 2021). A total 480 scat samples were collected from April 2019 to June 2021 (Figure 4).
Data analysis
Spatial analysis
Request a detailed protocolTo investigate the spatial distribution of carnivores, as well as the influence of environmental factors on the site occupancy of species in the study area, we performed single-season, single-species occupancy models to estimate carnivores’ occupancy (ψ) and detection (p) probability (Li et al., 2022b; MacKenzie, 2018; Moreno-Sosa et al., 2022). To ensure capture independence, only photo or video records at intervals greater than 30 min for same species were included in the data analysis (Li et al., 2010). We created a matrix recording whether each carnivore species was detected (1) or not (0) across several 30-day intervals (i.e. 0–30, 31–60, 61–90, 91–120, 121–150, >150 days) for each camera location. Based on the previous studies of habitat use of carnivores (Greenspan and Giordano, 2021; Alexander et al., 2016; Gorczynski et al., 2022), we selected terrain, vegetation, biological factors, and disturbance to construct the model. Terrain is a fundamental element of wildlife habitat and closely linked to other environmental factors (Chen et al., 2024). Terrain variables include elevation (ele) and roughness index (rix). Vegetation variables include ndvi and provide information on the level of habitat concealment. Biological variables include prey abundance (the number of independent photos of their preferred prey based on dietary analysis in this study, wolf and snow leopard: artiodactyla including livestock; Eurasian lynx and Pallas’s cat: lagomorpha; red fox and Tibetan fox: lagomorpha and rodentia) and reflect habitat preference and distribution patterns of carnivores. Disturbance variables include distance to roads (disrd) and human disturbances (hdis, the number of independent photos of herdsman and livestock) and can provide insight into the habitat selection and behavior patterns of carnivores. In addition, we used elevation, human disturbance, and prey as covariates that affect detection rate. Road data were obtained from Open Street Map (OSM, https://www.openstreetmap.org). Others environmental data were obtained from the Resource and Environment Science and Data Center (https://www.resdc.cn). We fitted all possible combinations of covariates and used Akaike’s information criterion (AIC) to rank candidate models, and selected ΔAIC≤2 model as the optimal model. If more than one optimal model resulted, then covariate estimates were obtained by using the equal-weight average.
Carnivore co-occurrence was evaluated using the Sørensen similarity index (value = 0, indicating maximum segregation and value = 1, indicating maximum co-occurrence) based on binary presence-absence data within the 5 km × 5 km camera trap grid (Torretta et al., 2021; Sorensen, 1948). Spatial analyses were performed using ArcGIS 10.8 (ESRI Inc), the ‘vegan’ packages (Oksanen et al., 2019), and ‘unmarked’ package (Fiske and Chandler, 2011) for R studio (version 1.1.463).
Temporal analysis
Request a detailed protocolEstimates of the coefficient of overlap (Δ) for activity patterns were estimated using the non-parametric kernel density method and applying time data obtained by the camera traps. Because the smallest sample had more than 50 records, we used the Δ4 estimator for pairwise comparisons between carnivore activity patterns, and used a smooth bootstrap scheme to generate 1000 resamples with 95% confidence intervals to test the reliability of the overlap value (Ridout and Linkie, 2009). Activity pattern analyses were performed using the ‘overlap’ R package. Values of the Δ4 estimates were calculated relative to 1000 simulated null models of randomized overlap data using the ‘compareCkern’ function in the ‘activity’ R package to test for differences in daily activity patterns (Ridout and Linkie, 2009; Rowcliffe et al., 2014).
Species identification and dietary analysis
Request a detailed protocolHost species and diet were identified using metabarcoding. DNA was extracted using the QIAamp Fast DNA Stool Mini Kit (QIAGEN, Hilden, Germany) following standard protocols and the MT-RNR1 (12S) and COX1 (cytochrome c oxidase subunit I) gene segments amplified using 12SV5-F/R primer and COX1 primers, respectively (Hacker et al., 2021; Riaz et al., 2011). PCR conditions followed the methods described in Hacker et al., 2021. The resulting library was sequenced on an Illumina NovaSeq platform and 250 bp paired-end reads were generated (Guangdong Magigene Biotechnology Co., Ltd., Guangzhou, China).
We used CLC Genomics Workbench version 12.0 to determine the host species as well as the prey consumed by each carnivore by mapping sequence reads to reference sequences of possible prey downloaded from GenBank and BOLD (Barcode of Life Data Systems) with representative haplotypes compiled into one.fasta file. Raw reads were required to have at least 98% similarity across at least 90% of the reference sequence for mapping (Hacker et al., 2021). Species and prey identification were made based on the reference taxa with the highest number of reads mapped and the fewest mismatches. Samples in which species could not be identified were analyzed to ensure the reference file was not incomplete by using the de nova assembly tool in CLC, then blasting the resulting contig sequence with the nucleotide databases in NCBI (https://blast.ncbi.nlm.nih.gov/Blast.cgi). As an additional precaution, the geographical range of the determined host and prey species was researched using the IUCN Red List (https://www.iucnredlist.org/) to ensure that it overlapped with the study site. For complete methods on data parameters and methods used, see Hacker et al., 2021.
Dietary data were summarized by the frequency of occurrence of prey species in scats observed. The ‘bipartite’ R package was used to construct food web networks (Dormann, 2011). Dietary diversity for each carnivore host species was assessed by calculating richness and Shannon’s index (Shannon and Weaver, 1949). Interspecific dietary niche overlap was evaluated using Pianka’s index (Ojk) (value = 0, no dietary overlap and value = 1, complete dietary overlap) and 95% confidence intervals were obtained by bootstrapping with 1000 resamples via the ‘spaa’ R package (Zhang, 2016). Dietary similarity between any two given carnivore species was assessed by calculating the inversed value of Jaccard’s index based on binary presence-absence data of prey.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
The ecology of fear: optimal foraging, game theory, and trophic interactionsJournal of Mammalogy 80:385–399.https://doi.org/10.2307/1383287
-
Coexistence mechanisms of small carnivores in a near-pristine area within the mountains of Southwest ChinaGlobal Ecology and Conservation 49:e02777.https://doi.org/10.1016/j.gecco.2023.e02777
-
Advances in species coexistence theoryBiodiversity Science 25:345–354.https://doi.org/10.17520/biods.2017034
-
Diet, morphology, and interspecific killing in carnivoraThe American Naturalist 167:524–536.https://doi.org/10.1086/501033
-
How to be a specialist? Quantifying specialisation in pollination networksNetwork Biology 1:1–20.
-
unmarked: an R package for fitting hierarchical models of wildlife occurrence and abundanceJournal of Statistical Software 43:1–23.https://doi.org/10.18637/jss.v043.i10
-
Investigating animal activity patterns and temporal niche partitioning using camera‐trap data: challenges and opportunitiesRemote Sensing in Ecology and Conservation 3:123–132.https://doi.org/10.1002/rse2.60
-
Spatial segregation between Iberian lynx and other carnivoresAnimal Biodiversity and Conservation 1:347–354.https://doi.org/10.32800/abc.2019.42.0347
-
Habitat use and diet of the red foxVulpes vulpes in an agricultural landscape in PolandZeitschrift Für Jagdwissenschaft 49:191–200.https://doi.org/10.1007/BF02189737
-
The summer microhabitat selection of Tibetan Fox in the Northwest Plateau of SichuanActa Theriologica Sinica 23:266–269.
-
Human density modulates spatial associations among tropical forest terrestrial mammal speciesGlobal Change Biology 28:7205–7216.https://doi.org/10.1111/gcb.16434
-
Spatio‐temporal occurrence and sensitivity to livestock husbandry of Pallas’s cat in the Mongolian AltaiThe Journal of Wildlife Management 86:e22150.https://doi.org/10.1002/jwmg.22150
-
Regional comparison of snow leopard (Panthera uncia) Diet using DNA metabarcodingBiodiversity and Conservation 30:797–817.https://doi.org/10.1007/s10531-021-02118-6
-
The competitive exclusion principleScience 131:1292–1297.https://doi.org/10.1126/science.131.3409.1292
-
Spatio-temporal interactions facilitate large carnivore sympatry across a resource gradientProceedings. Biological Sciences 284:20161860.https://doi.org/10.1098/rspb.2016.1860
-
An empirical evaluation of camera trap study design: How many, how long and when?Methods in Ecology and Evolution 11:700–713.https://doi.org/10.1111/2041-210X.13370
-
Paws without claws? Ecological effects of large carnivores in anthropogenic landscapesProceedings. Biological Sciences 283:20161625.https://doi.org/10.1098/rspb.2016.1625
-
Projected impacts of climate change on snow leopard habitat in Qinghai Province, ChinaEcology and Evolution 11:17202–17218.https://doi.org/10.1002/ece3.8358
-
Spatial co-occurrence and temporal activity patterns of sympatric mesocarnivores guild in Qinling MountainsGlobal Ecology and Conservation 36:e02129.https://doi.org/10.1016/j.gecco.2022.e02129
-
Interference interactions, co‐existence and conservation of mammalian carnivoresDiversity and Distributions 6:169–176.https://doi.org/10.1046/j.1472-4642.2000.00069.x
-
BookOccupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species OccurrenceAcademic Press, an imprint of Elsevier.
-
Examination of two new technologies to assess the diet of woodland caribou: video recorders attached to collars and DNA barcodingCanadian Journal of Forest Research 43:897–900.https://doi.org/10.1139/cjfr-2013-0108
-
Top predators constrain mesopredator distributionsNature Communications 8:15469.https://doi.org/10.1038/ncomms15469
-
Diet of Eurasian lynx, Lynx lynx, in the boreal forest of southeastern Norway: the relative importance of livestock and hares at low roe deer densityEuropean Journal of Wildlife Research 52:237–244.https://doi.org/10.1007/s10344-006-0052-4
-
Spatial behaviour of the red fox ( Vulpes vulpes ) in a rural area of central ItalyItalian Journal of Zoology 64:351–358.https://doi.org/10.1080/11250009709356222
-
ecoPrimers: inference of new DNA barcode markers from whole genome sequence analysisNucleic Acids Research 39:e145.https://doi.org/10.1093/nar/gkr732
-
Estimating overlap of daily activity patterns from camera trap dataJournal of Agricultural, Biological, and Environmental Statistics 14:322–337.https://doi.org/10.1198/jabes.2009.08038
-
ThesisProviding an ecological basis for the conservation of the Pallas’s cat (Otocolobus manul) PhD thesisUniversity of Bristol.
-
Quantifying levels of animal activity using camera trap dataMethods in Ecology and Evolution 5:1170–1179.https://doi.org/10.1111/2041-210X.12278
-
BookThe Mathematical Theory of CommunicationPhilosophical Review.https://doi.org/10.1063/1.3067010
-
Feeding patterns of three widespread carnivores—the wolf, snow leopard, and red fox—in the trans-himalayan landscape of IndiaFrontiers in Ecology and Evolution 10:815996.https://doi.org/10.3389/fevo.2022.815996
-
Fear of the human “super predator” reduces feeding time in large carnivoresProceedings. Biological Sciences 284:20170433.https://doi.org/10.1098/rspb.2017.0433
-
BookA Method of Establishing Groups of Equal Amplitude in Plant Sociology Based on Similarity of Species Content and Its Application to Analyses of the Vegetation on Danish CommonsRoyal Academy press.
-
Maximising survival by shifting the daily timing of activityEcology Letters 22:2097–2102.https://doi.org/10.1111/ele.13404
-
Dietary overlap of snow leopard and other carnivores in the Pamirs of Northwestern ChinaChinese Science Bulletin 59:3162–3168.https://doi.org/10.1007/s11434-014-0370-y
-
Camera-trapping survey of the mammals and birds in the Qilian Mountain National Park (Qinghai area)China. ACTA Theriologica Sinica 39:466–475.https://doi.org/10.16829/j.slxb.150307
Article and author information
Author details
Funding
National Natural Science Foundation of China (32201430)
- Jia Li
National Natural Science Foundation of China (32101409)
- Yu Zhang
Welfare Project of the National Scientific Research Institution (CAFYBB2019ZE003)
- Yuguang Zhang
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Mr. Jiong He, Yayue Gao, Duifang Ma, Liji Wu, Dazhi Hu, and other colleagues of Qilianshan National Park for their generous assistance in the field surveys. We thank Dr. Charlotte Hacker for editing the English text of this manuscript.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.90559. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Cong et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 422
- views
-
- 73
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Ecology
Most of Earth’s biomass is composed of polysaccharides. During biomass decomposition, polysaccharides are degraded by heterotrophic bacteria as a nutrient and energy source and are thereby partly remineralized into CO2. As polysaccharides are heterogeneously distributed in nature, following the colonization and degradation of a polysaccharide hotspot the cells need to reach new polysaccharide hotspots. Even though many studies indicate that these degradation-dispersal cycles contribute to the carbon flow in marine systems, we know little about how cells alternate between polysaccharide degradation and motility, and which environmental factors trigger this behavioral switch. Here, we studied the growth of the marine bacterium Vibrio cyclitrophicus ZF270 on the abundant marine polysaccharide alginate, both in its soluble polymeric form as well as on its breakdown products. We used microfluidics coupled to time-lapse microscopy to analyze motility and growth of individual cells, and RNA sequencing to study associated changes in gene expression. We found that single cells grow at reduced rate on alginate until they form large groups that cooperatively break down the polymer. Exposing cell groups to digested alginate accelerates cell growth and changes the expression of genes involved in alginate degradation and catabolism, central metabolism, ribosomal biosynthesis, and transport. However, exposure to digested alginate also triggers cells to become motile and disperse from cell groups, proportionally increasing with the group size before the nutrient switch, and this is accompanied by high expression of genes involved in flagellar assembly, chemotaxis, and quorum sensing. The motile cells chemotax toward polymeric but not digested alginate, likely enabling them to find new polysaccharide hotspots. Overall, our findings reveal cellular mechanisms that might also underlie bacterial degradation-dispersal cycles, which influence the remineralization of biomass in marine environments.
-
- Ecology
Litter decomposition is expected to be positively associated with precipitation despite evidence that decomposers of varying sizes have different moisture dependencies. We hypothesized that higher tolerance of macro-decomposers to aridity may counterbalance the effect of smaller decomposers, leading to similar decomposition rates across climatic gradients. We tested this hypothesis by placing plant litter baskets of different mesh sizes in seven sites along a sharp precipitation gradient, and by characterizing the macro-decomposer assemblages using pitfall trapping. We found that decomposers responded differently to precipitation levels based on their size. Microbial decomposition increased with precipitation in the winter while macro-decomposition peaked in arid sites during the summer. This led to similar overall decomposition rates across the gradient except in hyper-arid sites. Macro-decomposer richness, abundance, and biomass peaked in arid environments. Our findings highlight the importance of macro-decomposition in arid-lands, possibly resolving the dryland decomposition conundrum, and emphasizing the need to contemplate decomposer size when investigating zoogeochemical processes.