Different coexistence patterns between apex carnivores and mesocarnivores based on temporal, spatial, and dietary niche partitioning analysis in Qilian Mountain National Park, China

  1. Wei Cong
  2. Jia Li
  3. Charlotte Hacker
  4. Ye Li
  5. Yu Zhang
  6. Lixiao Jin
  7. Yi Zhang
  8. Diqiang Li
  9. Yadong Xue  Is a corresponding author
  10. Yuguang Zhang  Is a corresponding author
  1. Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, China
  2. Institute of Ecological Conservation and Restoration, Chinese Academy of Forestry, China
  3. Department of Biological Sciences, Duquesne University, United States

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.sa0

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.

Table 1
Summary of occupancy rate and detection probability of different species for the optimal models (ΔAIC≤2).
SpeciesModelsNumber of parametersAICΔAICAIC Wtψp
WolfPsi (rix+ndvi+prey); P (ele+prey)7282.140.0000.0770.1800.153
Psi (rix+ndvi+prey); P (ele+hdis)7282.450.3020.0660.1830.145
Psi (rix+ndvi+prey); P (ele)6283.671.5250.0360.2170.111
Psi (rix+ndvi+prey); P (ele+hdis+prey)8284.101.9590.0290.1790.154
Psi (ele+rix+ndvi+prey); P (ele+prey)8284.121.9720.0290.1680.165
Psi (rix+ndvi+disrd+prey); P (ele+prey)8284.142.0000.0280.1780.156
Model average0.1840.147
Snow leopardPsi (disrd+ndvi); P (hdis+prey)6970.690.0000.0350.4360.419
Psi (ele+disrd); P (hdis+prey)6970.810.1130.0330.4350.421
Psi (.); P (hdis+prey)4970.990.2920.0300.4370.420
Psi (disrd+ndvi); P (ele+hdis+prey)7971.120.4290.0280.4400.413
Psi (ele+disrd+ndvi); P (hdis+prey)7971.130.4390.0280.4330.420
Psi (.); P (ele+hdis+prey)5971.210.5110.0270.4420.413
Psi (disrd); P (ele+hdis+prey)6971.220.5250.0270.4440.410
Psi (disrd); P (hdis+prey)5971.520.8250.0230.4370.420
Psi (ele); P (hdis+prey)5971.620.9280.0220.4350.421
Psi (ndvi); P (hdis+prey)5971.640.9470.0220.4370.420
Psi (ele+disrd); P (ele+hdis+prey)7971.921.2220.0190.4360.416
Psi (ndvi); P (ele+hdis+prey)6972.421.7260.0150.4400.414
Psi (ele+disrd+ndvi); P (ele+hdis+prey)8972.451.7570.0140.4350.416
Psi (disrd+ndvi+rix); P (hdis+prey)7972.671.9740.0130.4360.419
Model average0.4370.417
Eurasian lynxPsi (ndvi); P (hdis)4371.320.0000.0250.1630.321
Psi (rix+ndvi); P (hdis)5371.400.0760.0240.1610.322
Psi (ndvi); P (hdis+prey)5371.640.3250.0210.1700.296
Psi (rix+ndvi); P (hdis+prey)6371.860.5390.0190.1670.300
Psi (ndvi+hdis); P (prey)5371.940.6230.0190.1630.315
Psi (ndvi+disrd); P (hdis)5372.020.7020.0180.1630.321
Psi (rix+ndvi+hdis); P (prey)6372.230.9070.0160.1610.318
Psi (ndvi+hdis); P (.)4372.240.9160.0160.1560.343
Psi (rix+ndvi+hdis); P (.)5372.361.0440.0150.1540.345
Psi (rix+ndvi+disrd); P (hdis)6372.391.0720.0150.1610.321
Psi (ndvi+disrd); P (hdis+prey)6372.641.3180.0130.1680.300
Psi (ndvi+prey); P (hdis)5372.671.3520.0130.1630.320
Psi (rix+ndvi+prey); P (hdis)6372.991.6720.0110.1600.323
Psi (ndvi+disrd+hdis); P (.)5373.061.7440.0110.1560.344
Psi (ndvi+hdis); P (hdis)5373.091.7750.0100.1560.325
Psi (ndvi+disrd+hdis); P (prey)6373.111.7920.0100.1620.319
Psi (rix+ndvi+disrd); P (hdis+prey)7373.121.7960.0100.1650.303
Psi (rix+ndvi+hdis); P (hdis)6373.231.9130.0100.1540.326
Psi (ele+ndvi); P (hdis)5373.321.9990.0090.1640.319
Model average0.1610.320
Red foxPsi (disrd+hdis+prey); P (.)5894.440.0000.0390.3690.391
Psi (disrd+hdis+prey); P (ele)6894.590.1520.0370.3700.387
Psi (disrd+hdis+ndvi+prey); P (.)6895.300.8640.0260.3700.391
Psi (disrd+hdis+ndvi+prey); P (ele)7895.370.9340.0250.3700.387
Psi (disrd+hdis+prey+rix); P (.)6895.991.5440.0180.3690.392
Psi (disrd+hdis+prey+rix); P (ele)7896.161.7170.0170.3700.388
Psi (disrd+hdis+prey); P (prey)6896.291.8500.0160.3680.388
Psi (disrd+hdis+prey); P (hdis)6896.371.9320.0150.3680.394
Model average0.3690.390
Tibetan foxPsi (ele+disrd+prey+rix); P (ele+prey)8326.430.0000.0660.1470.236
Psi (disrd+prey+rix); P (ele+prey+hdis)8327.140.7140.0460.2770.135
Psi (ele+disrd+hdis+prey+rix); P (ele+prey)9327.200.7750.0450.1510.230
Psi (ele+disrd+rix); P (ele+prey)7327.290.8660.0430.1580.219
Psi (ele+disrd+prey+rix); P (ele+prey+hdis)9327.531.1000.0380.1530.222
Psi (ele+disrd+prey); P (ele+prey)7327.691.2630.0350.1490.234
Psi (disrd+hdis+prey+rix); P (ele+prey)8328.011.5850.0300.2600.142
Psi (ele+disrd); P (ele+prey)6328.421.9860.0240.1590.219
Model average0.1820.205
  1. disrd – distance to roads, ele – elevation, ndvi – normalized difference vegetation index, rix – roughness index, hdis – human disturbance.

Table 1—source data 1

Data captured by camera traps on carnivore species.

https://cdn.elifesciences.org/articles/90559/elife-90559-table1-data1-v1.xlsx
Table 2
Covariates influencing carnivore occupancy rate and detection probability based on the optimal models (ΔAIC≤2).
SpeciesModel componentCovariatesEstimate (β)SEZp
WolfOccupancyIntercept–2.0310.4974.083<0.001***
ndvi–0.8430.2872.9350.003**
prey1.3180.5342.4700.014*
rix–1.1660.5702.0460.041*
DetectionIntercept–1.9200.4574.198<0.001***
ele0.6180.2662.3200.020*
Snow leopardOccupancyIntercept–0.2560.1431.7830.075
disrd–0.1630.1720.9530.341
ndvi–0.0950.1450.6530.514
DetectionIntercept–0.4510.1243.628<0.001***
hdis–1.8300.5203.519<0.001***
ele0.6610.2472.6780.007**
Eurasian lynxOccupancyIntercept–1.8790.3165.941<0.001***
ndvi–0.5830.2012.8990.004**
DetectionIntercept–0.9990.3752.6650.008**
hdis–1.7851.6901.0560.291
Red foxOccupancyIntercept–0.4410.1772.4880.013*
disrd0.3310.1781.8660.062
hdis1.3330.6112.1840.029*
prey0.6650.3192.0830.037*
DetectionIntercept–0.4510.1213.729<0.001***
ele0.0570.0930.6060.545
Tibetan foxOccupancyIntercept–2.0870.6433.247<0.001***
disrd0.9550.3212.9730.003**
ele0.7680.5471.4030.161
DetectionIntercept–1.6790.6852.4530.014*
ele0.8820.4172.1120.035*
prey0.3690.1672.2160.027*
  1. The different superscript letters represent significance, ***p<0.001, **0.001<p<0.01, *0.01<p<0.05.

The relationship between carnivore species occupancy probability and covariates by the optimal models (∆AIC≤2).

The solid line represents the fitted polynomial regression and the gray area represents 95% confidence intervals.

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.

Table 3
Spatial overlap (Sørensen’s index), diel activity overlap (Δ), and dietary overlap (Pianka’s index), as well as confidence intervals for carnivore species.
Sørensen’s indexΔPianka’s index
Wolf – Snow leopard0.2770.676 (0.562–0.756)0.892 (0.804–0.982)
Wolf – Eurasian lynx0.2720.661 (0.541–0.759)0.585 (0.141–0.881)
Snow leopard – Eurasian lynx0.3050.900 (0.854–0.992)0.458 (0.160–0.886)
Wolf – Pallas’s cat0.658 (0.053–0.950)
Wolf – Red fox0.3650.497 (0.359–0.563)0.811 (0.497–0.962)
Wolf – Tibetan fox0.3500.855 (0.777–0.937)0.689 (0.456–0.967)
Snow leopard – Pallas’s cat0.354 (0.092–0.827)
Snow leopard – Red fox0.4770.814 (0.756–0.863)0.586 (0.568–0.941)
Snow leopard – Tibetan fox0.2920.711 (0.629–0.744)0.390 (0.299–0.905)
Eurasian lynx – Pallas’s cat0.376 (0–0.921)
Eurasian lynx – Red fox0.2050.800 (0.728–0.878)0.536 (0.442–0.910)
Eurasian lynx – Tibetan fox0.1980.695 (0.601–0.756)0.385 (0.078–0.919)
Pallas’s cat – Red fox0.945 (0.369–0.988)
Pallas’s cat – Tibetan fox0.997 (0–1)
Red fox – Tibetan fox0.3490.576 (0.467–0.580)0.949 (0.279–0.988)

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).

Daily activity patterns of carnivore species.

Shades of lavender indicate temporal areas of overlap. The p-values are derived from Wald tests. A p-value less than 0.05 indicates that the parameter significantly contributes to the model.

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.

Figure 3 with 2 supplements see all
The food web of carnivore species (SL – snow leopard, EL – Eurasian lynx, PC – Pallas’s cat, RF – red fox, TF – Tibetan fox).

The heights of the left bars represent the frequency of occurrence of the taxonomic order of prey species, the middle bars represent the number of samples for each carnivore, and the heights of the right bars represent the frequency of occurrence of prey species in scats. The colors of prey match the taxonomic orders. The connecting line widths represent the prey frequency of occurrence in the diet of each carnivore species.

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 protocol

The 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).

Locations of camera trap stations and scat collection sites in this study.

Camera trap monitoring and noninvasive sampling

Request a detailed protocol

The 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 protocol

To 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 protocol

Estimates 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 protocol

Host 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

  1. Book
    1. Boitani L
    (2016)
    Camera Trapping for Wildlife Research
    Pelagic Publishing Ltd.
    1. Dormann CF
    (2011)
    How to be a specialist? Quantifying specialisation in pollination networks
    Network Biology 1:1–20.
  2. Book
    1. Gause GF
    (1934)
    The Struggle for Existence
    Baltimore: Williams and Wilkins.
    1. Gong M
    2. Hu J
    (2003)
    The summer microhabitat selection of Tibetan Fox in the Northwest Plateau of Sichuan
    Acta Theriologica Sinica 23:266–269.
  3. Book
    1. MacKenzie DI
    (2018)
    Occupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species Occurrence
    Academic Press, an imprint of Elsevier.
  4. Thesis
    1. Ross S
    (2009)
    Providing an ecological basis for the conservation of the Pallas’s cat (Otocolobus manul) PhD thesis
    University of Bristol.
  5. Book
    1. Ross S
    2. Barashkova A
    3. Kirilyuk V
    4. Naidenko S
    (2019)
    The Behaviour and Ecology of the Manul
    Cat Specialist Group.
  6. Book
    1. Sorensen TJ
    (1948)
    A 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 Commons
    Royal Academy press.
  7. Book
    1. Zheng J
    (2011)
    The Study of Qinghai Nature Reserve
    Qinghai People’s Publishing House.

Article and author information

Author details

  1. Wei Cong

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Investigation, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0009-0005-6071-9564
  2. Jia Li

    Institute of Ecological Conservation and Restoration, Chinese Academy of Forestry, Beijing, China
    Contribution
    Funding acquisition, Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
  3. Charlotte Hacker

    Department of Biological Sciences, Duquesne University, Pittsburgh, United States
    Contribution
    Investigation, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
  4. Ye Li

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  5. Yu Zhang

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Funding acquisition, Investigation
    Competing interests
    No competing interests declared
  6. Lixiao Jin

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  7. Yi Zhang

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  8. Diqiang Li

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Funding acquisition, Methodology, Project administration
    Competing interests
    No competing interests declared
  9. Yadong Xue

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Data curation, Investigation, Project administration, Writing – review and editing
    For correspondence
    xueyadong334@163.com
    Competing interests
    No competing interests declared
  10. Yuguang Zhang

    Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Key Laboratory of Biodiversity Conservation, State Forestry and Grassland Administration, Beijing, China
    Contribution
    Data curation, Funding acquisition, Investigation, Methodology, Project administration, Writing – review and editing
    For correspondence
    yugzhang@sina.com.cn
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9801-8556

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

  1. Sent for peer review:
  2. Preprint posted:
  3. Reviewed Preprint version 1:
  4. Reviewed Preprint version 2:
  5. 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

  • 528
    views
  • 83
    downloads
  • 0
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Wei Cong
  2. Jia Li
  3. Charlotte Hacker
  4. Ye Li
  5. Yu Zhang
  6. Lixiao Jin
  7. Yi Zhang
  8. Diqiang Li
  9. Yadong Xue
  10. Yuguang Zhang
(2024)
Different coexistence patterns between apex carnivores and mesocarnivores based on temporal, spatial, and dietary niche partitioning analysis in Qilian Mountain National Park, China
eLife 13:RP90559.
https://doi.org/10.7554/eLife.90559.3

Share this article

https://doi.org/10.7554/eLife.90559

Further reading

    1. Ecology
    Ivan Pokrovsky, Teja Curk ... Martin Wikelski
    Research Article

    Advances in tracking technologies have revealed the diverse migration patterns of birds, which are critical for range mapping and population estimation. Population trends are usually estimated in breeding ranges where birds remain stationary, but for species that breed in remote areas like the Arctic, these trends are often assessed in over-wintering ranges. Assessing population trends during the wintering season is challenging due to the extensive movements of birds in these ranges, which requires a deep understanding of the movement dynamics. However, these movements remain understudied, particularly in the mid-latitudes, where many Arctic breeders overwinter, increasing uncertainty in their ranges and numbers. Here, we show that the Arctic breeding raptor Rough-legged buzzard, which overwinters in the mid-latitudes, has a specific wintering strategy. After migrating ca. 1500 km from the Arctic to mid-latitudes, the birds continue to move throughout the entire over-wintering period, traveling another 1000 km southwest and then back northeast as the snowline advances. This continuous movement makes their wintering range dynamic throughout the season. In essence, this movement represents an extension of the quick migration process, albeit at a slower pace, and we have termed this migration pattern ‘foxtrot migration’, drawing an analogy to the alternating fast and slow movements of the foxtrot dance. These results highlight the potential errors in range mapping from single mid-winter surveys and emphasize the importance of this migration pattern in assessing the conservation status of bird species. Understanding this migration pattern could help to correctly estimate bird populations in over-wintering ranges, which is especially important for species that nest in hard-to-reach regions such as the Arctic.

    1. Ecology
    2. Evolutionary Biology
    Zhixian Zhang, Jianying Li ... Songdou Zhang
    Research Article

    Seasonal polyphenism enables organisms to adapt to environmental challenges by increasing phenotypic diversity. Cacopsylla chinensis exhibits remarkable seasonal polyphenism, specifically in the form of summer-form and winter-form, which have distinct morphological phenotypes. Previous research has shown that low temperature and the temperature receptor CcTRPM regulate the transition from summer-form to winter-form in C. chinensis by impacting cuticle content and thickness. However, the underling neuroendocrine regulatory mechanism remains largely unknown. Bursicon, also known as the tanning hormone, is responsible for the hardening and darkening of the insect cuticle. In this study, we report for the first time on the novel function of Bursicon and its receptor in the transition from summer-form to winter-form in C. chinensis. Firstly, we identified CcBurs-α and CcBurs-β as two typical subunits of Bursicon in C. chinensis, which were regulated by low temperature (10 °C) and CcTRPM. Subsequently, CcBurs-α and CcBurs-β formed a heterodimer that mediated the transition from summer-form to winter-form by influencing the cuticle chitin contents and cuticle thickness. Furthermore, we demonstrated that CcBurs-R acts as the Bursicon receptor and plays a critical role in the up-stream signaling of the chitin biosynthesis pathway, regulating the transition from summer-form to winter-form. Finally, we discovered that miR-6012 directly targets CcBurs-R, contributing to the regulation of Bursicon signaling in the seasonal polyphenism of C. chinensis. In summary, these findings reveal the novel function of the neuroendocrine regulatory mechanism underlying seasonal polyphenism and provide critical insights into the insect Bursicon and its receptor.