Next Article in Journal
Taxonomic Revision of the Nearctic Genus Drepanaphis Del Guercio (Hemiptera, Aphididae: Drepanosiphinae)
Previous Article in Journal
Can the Necrophagous Blow Fly Calliphora vicina (Diptera: Calliphoridae) Be Reared on Plant-Based Meal?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 15
Issue 7
10.3390/insects15070552
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Patterns and Drivers of Bumblebee Diversity in Gansu

by
Muhammad Naeem
1,2,†,
Huanhuan Chen
1,2,3,†,
Wenbo Li
1,
Alice C. Hughes
4,
Paul H. Williams
5,
Nawaz Haider Bashir
1,
Zhengying Miao
6,
Jiaxing Huang
2,* and
Jiandong An
2,*
1
College of Biological Resource and Food Engineering, Qujing Normal University, Qujing 655011, China
2
State Key Laboratory of Resource Insects, Key Laboratory of Insect-Pollinator Biology of Ministry of Agriculture and Rural Affairs, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100193, China
3
Key Laboratory of Yunnan Provincial Department of Education of the Deep-Time Evolution on Biodiversity from the Origin of the Pearl River, Qujing Normal University, Qujing 655011, China
4
School of Biological Sciences, University of Hong Kong, Hong Kong SAR 999077, China
5
The Natural History Museum, London SW75BD, UK
6
Gansu Institute of Apiculture, Tianshui 741022, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Insects 2024, 15(7), 552; https://doi.org/10.3390/insects15070552
Submission received: 15 May 2024 / Revised: 6 July 2024 / Accepted: 17 July 2024 / Published: 21 July 2024
(This article belongs to the Section Insect Societies and Sociality)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Bumblebee species are crucial pollinators for both crops and natural ecosystems. Managing diversity and securing the services they provide requires sufficient data to manage these systems effectively. In this study, seventeen environmental factors were selected to understand their contribution to shaping species assemblages and to provide baseline information for understanding drivers of community assembly. The study focused on Gansu Province because it has some of the highest richness of bumblebees within China, and China represents a major data gap for bee distribution data. Net primary productivity followed by water vapor pressure were top-ranked drivers of species distributions. The results enhance our understanding of the importance of environmental factors in shaping bumblebee community assemblages.

Abstract

Understanding the influence of factors responsible for shaping community assemblage is crucial for biodiversity management and conservation. Gansu is one of the richest regions for bumblebee species in the world. We explored the distribution data of 52 bumblebee species collected in Gansu and its surroundings between 2002 and 2022, predicting habitat suitability based on 17 environmental variables using MaxEnt. The factors influencing community assemblage were assessed using canonical correspondence analysis. Net primary productivity, water vapor pressure, temperature seasonality, annual precipitation, and precipitation seasonality were some of the most influential drivers of species distributions. Based on Ward’s agglomerative cluster analysis, four biogeographic zones are described: the Southern humid zone, the Western Qilian snow mountain zone, the Eastern Loess plateau zone, and the Western dry mountain zone. In the clusters of grid cells based on beta diversity values, the Southern humid zone comprised 42.5% of the grid cells, followed by the Eastern Loess plateau zone (32.5%), the Western dry mountain zone (20%), and the Western Qilian snow mountain zone (5%). Almost all the environmental factors showed a significant contribution to the assemblages of bumblebees of different groups. Our findings highlight the need for better data to understand species biogeography and diversity patterns, and they provide key baseline data for refining conservation strategies.

1. Introduction

Species assemblages may be shaped by a combination of biogeographic, climatic, and biotic factors, but the relative influence of each is likely to vary considerably [1,2,3,4,5,6]. Climatic factors are known to shape the community assemblage of dung beetles, but few studies exist for other invertebrate taxa [7,8]. Invertebrates provide important services, yet the drivers of their distribution require further research. Bumblebees are key wild pollinators, yet various species are showing population declines [9,10]. Bumblebee community assemblages are impacted by landcover, yet the effect of other factors is still little known [11,12,13].
Gansu Province provides an ideal place to explore drivers of distribution patterns in bumblebees, as 45% of the 125 species recorded in China have been recorded there, one of the highest provincial counts globally [14,15]. For example, the Qinghai-Tibetan plateau has 57 species, whereas South America has 44, the USA has 49, and Canada has 40 bumblebee species [16,17,18,19]. Thus, these data fill a vital gap as China is a hotspot of bee diversity [20], yet, the majority of prior research in Gansu has been limited to basic lists of specimens and species of bumblebees, and higher-resolution analysis is clearly needed [21]. Before 2011, reports on bumblebees of Gansu were reliant on small collections from a few specific parts of the province [22,23,24,25]. Later, a comprehensive report was compiled based on systematic field surveys of bumblebee species collections, which resulted in an increased number of records for Gansu province [26].
Here, we explore bumblebee faunal composition based on bioregionalization at species level across Gansu [25,27]. We analyze (1) patterns of richness of bumblebees in Gansu and (2) map biogeographic regions of bumblebees in Gansu (3) to estimate the contribution of different environmental factors in the assembly of bumblebee communities within each biogeographic zone.

2. Materials and Methods

2.1. Study Site and Bumblebee Species

Gansu is situated at the northeastern edge of the Qinghai-Tibetan plateau (Figure 1), in a transition zone between the Neimenggu plateau, the Qinghai-Tibetan plateau, and the Loess plateau in China. No bumblebee records were found in northwestern Gansu, so we expanded our study area to include surrounding regions in order to obtain bumblebee records from northwestern regions, increasing the number of species included from 57 within Gansu to 62 (Figure 1).
More than 8000 specimens of 62 total bumblebee species were recorded from 2295 accessible sampling sites within Gansu and its surrounding regions (Figure 1 and Figure S1 and Table 1). The bumblebee species data used in this study are based on the collection of Institute of Apicultural Research, Chinese Academy of Agricultural Sciences (IAR-CAAS), obtained across Gansu between 2002 and 2022. For the collection of bumblebee species, a range of different environmental conditions was surveyed, including most protected areas and nature reserves across Gansu, to ensure high-quality habitats were included. Collection was conducted from June to September each year, with each site visited by 3-4 non-specialists for 1–2 h on sunny days to maximize bumblebee activity. All specimens were pinned and labeled with individual identifier numbers. Specimens were identified using both a morphological key [24] and COI barcode sequences. All collected specimens were stored in entomological boxes and deposited in the collection of the Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing, China. The biogeographic zones assessment was based on habitat suitability ranges of modeled bumblebee species within Gansu.

2.2. Statistical Analysis

The distribution ranges of bumblebee species calculated the pixels of 0.5° longitude × 0.5° latitude within the boundary of Gansu, and there was total of 240 grid cells throughout Gansu (Figure S1) [28]. The values of habitat suitability per grid cell were assessed for all 240 grid cells, and a species vs. sites (grid cells) matrix was developed, where rows represent the grid cells or sites and columns represent the species. Principal component analysis was performed on all 17 variables to better cluster and assess distributional patterns (Figure 2).
For the community assemblage of bumblebee species, Simpson’s dissimilarity was determined for 240 grid cells based on predicted habitat suitability values per grid cell within the Gansu of each of the bumblebee species using the “vegan” and “betapart” packages in R version 4.3.3 [29,30]. This calculation was made based on the following formula:
β s i m = min ⁡ ( b ,   c ) a + min ⁡ ( b ,   c )
Here, the number of species shared by the two grid cells is represented by “a”, the unique number of species to one grid cell is represented by “b”, and the unique number of species to another grid cell is represented by “c”. Based on the dissimilarity values, Ward’s agglomerative clustering analysis was used to classify the grid cells into different groups [31] (Figure 3a,b). To explore the impact of environmental heterogeneity, we used Canonical Correspondence Analysis (CCA) in the “vegan” package of R 3.3.3 [29,30] (Figure 4). CCA was used here for the analysis of environmental heterogeneity because it can assess that multiple variables are influenced simultaneously on the species assemblage. The permutation test is used to assess the statistical significance of the CCA axes. Furthermore, the strength and direction of the relationship between environmental factors and community assemblage can be calculated [32,33].
Whilst our study provides an exploration of distribution and zones of bumblebees in Gansu, it was not without limitations. Given persistent data gaps [34], modeling can provide a key mechanism to overcome spatial sampling gaps, when used with caution. The sampling of bumblebees is difficult in most of Gansu because in more arid environments species are patchily distributed and much more difficult to collect [19,20,21].

2.3. Model Development and Assessment

The suitability ranges were assessed using MaxEnt software v. 3.4.4 [35]. MaxEnt is widely used to assess species distribution. It works on the major assumption of presence-only data, meaning it does not require absence data, and requires environmental variables that are most relevant to the species [35]. To model the distribution of bumblebee species, we selected environmental variables that are most likely to drive species distributions and their reproductive success of species, and these are cross-referenced with the distribution of each species to assess the habitat requirements of each species and then to assess where each species may occur [19,20] (see Table S1 for variables and sources).
For autocorrelation, SDMtoolbox in ArcMap v. 10.3 was used to calculate Pearson’s correlation coefficients (r < 0.9) [36], and a set of 17 variables was selected. These included bioclimatic factors, tree density, canopy height, soil pH, water vapor, and solar radiation factors, using similar factors as previous studies [37]. We did not include elevation because it is used to downscale bioclimatic factors [38]. The model evaluation was assessed on the basis of the area under the curve (AUC) of the receiver operating characteristic (ROC) curves [39].

3. Results

3.1. The Contribution of Environmental Factors in the Habitat Suitability of Bumblebee Species

We modeled 52 out of 62 bumblebee species, excluding 10 species with less than three records of distribution. Species varied from 0 to 23 (Figure 5; Table 1). The area under the curve values of training data vary from 0.94 to 0.99, whereas the range of values of test data varies from 0.526 to 1. The 10th-percentile training presence threshold values showed that the predicted suitable areas vary from 0.003 to 0.188 (Table S2). However, the factor with the highest contribution varies for each species. For example, for 87% of bumblebees, the most important factor was “Net primary productivity”, and in 65% of species, this factor showed top contribution in the distribution of species. Water vapor pressure was the top contributing factor for 10% of bumblebee species. Additionally, “Vegetation height”, “Solar radiation”, and “Tree density” contributed 2% each in the distribution of bumblebee species. Various bioclimatic factors were also important; “TS” (temperature seasonality), “AP” (annual precipitation), and “PS” (precipitation seasonality) accounted for 5% each in the distribution of species. Similarly, “MaxTWM” (maximum temperature of warmest month) exhibited a 4% contribution to the distribution of bumblebee species, whereas “MinTCM” (minimum temperature of coldest month) and “PWM” (precipitation of wettest month) contributed 2% each to the overall species distribution pattern. Further details are present in the Supplementary Materials (Table S2).
Using principal component analysis (PCA), the first three axes included 79.24% of total variation (Figure 2). PC1 gave a substantial proportion variation, almost 48.41% among 240 grid cells of Gansu. PC2 accounts for 22.82% environmental heterogeneity, and PC3 accounts for 8.01% environmental heterogeneity. Distinct clustering patterns emerged in 3D scatter plot representations of PCA (Figure 2). The spatial distribution of the grid cells of different groups in scatter plot representations show similar kinds of clustering, as displayed in CCA analysis (Figure 4).

3.2. Biogeographic Zones of Bumblebee Community Assemblage

Based on predicted values of habitat suitability modeling, Ward’s agglomerative cluster analysis segregated the species into four distinct groups (Figure 3a). In the clusters of grid cells based on beta diversity values, the Southern humid zone (green) comprised 42.5% of the grid cells, followed by the Eastern Loess plateau zone (pink) (32.5%), the Western dry mountain zone (blue) (20%) and the Western Qilian snow mountain zone (red) (5% of the grid cells). A phenon line was also drawn on the clusters to visualize the four major bumblebee biogeographic regions based on the cophenetic correlation coefficient value (r = 0.65) (Figure 3a). The relative positions of the four bumblebee biogeographic regions of bumblebees in Gansu are displayed in Figure 3b. The bumblebee species found only in the Southern humid zone are Bombus atripes, B. bicoloratus, B. trifasciatus, B. flavescens, B. breviceps, B. hengduanensis, B. remotus, B. grahami, B. festivus, etc., which is the relative moist group of South China. B. lantschouensis, B. ganjsuensis, B. ignitus, B. pyrosoma, B. melanurus, B. sibiricus, B. opulentus, B. longipes, B. hedini, etc., are mainly distributed in the Eastern Loess plateau zone, which belongs to the temperate group of North China. The species limited to the Western Qilian snow mountain zone are B. validus, B. wangae, B. waltoni, B. tibetanus, B. prshewalskyi, B. supremus, B. personatus, B. qilianensis, B. kashmirensis, B. minshanicola, and B. convexus, most of which are Qinghai-Tibetan plateau species with high elevation. Compared to the above three groups, few species, such as B. asiaticus and B. difficillimus, are restricted to the Western dry mountain zone of the Central Asian group.

3.3. Bumblebee Community Assemblage Explained by Environmental Heterogeneity

In our set of 17 environmental variables (please see Table S1), simple linear regression analysis indicates that bee presence in the Southern humid zone is significantly (p < 0.05) positively related to NPP, isothermality, AP, PWM, PDM, and PS. In the Eastern Loess plateau zone, the presence of bees is significantly related to NPP, AP, PWM, PDM, MinTCM, tree density, and AMT. In the Western dry mountain zone, the presence of bumblebee species is positively and significantly related to NPP, AP, PWM, MinTCM, AMT, and vapr. Finally, the presence of species in the Western Qilian snow mountain zone is positively associated with NPP, AP, and PWM.
There are no bumblebees found in the northwest part of Gansu. The main factor is that this area is a very dry desert. The relationship between bumblebee community assemblage and environmental variables is also shown in Figure 4 based on Canonical Correspondence Analysis (CCA) constrained space. The eigenvalues indicate that axis 1, representing 38.35% of the variance, plays a substantial role in explaining the observed patterns of community assemblage. Axis 2, with 18.83% of the variance, contributes further to the understanding of species–environment relationships. Together, these two axes explain a significant portion (57.18%) of the total variance present between the species assemblage of different sampling sites. In our results, the permutation “p” value for both axis 1 and axis 2 is 0.001. This low p-value suggests that the observed relationship between bumblebee species assemblage and environmental variables on these axes is statistically significant (Figure 4).
The most important variables found to be highly related in the regionalization of the Southern humid zone and the Western Qilian snow mountain zone are soil pH, temperature seasonality, and solar radiation. Regarding the regionalization of the Eastern Loess plateau zone, the most important variables are average mean temperature (AMT), tree density, maximum temperature of warmest month (MaxTWM), and Net primary productivity (NPP). Similarly, for the region of Western dry mountain, the most important variables found were precipitation seasonality (PS) and isothermality. Axis 1 of the CCA plot showed a strong positive relationship between the species and environmental factors such as wind (with the values of 0.87), mean diurnal range (0.74), solar radiation (0.70), temperature seasonality (0.48), isothermality (0.36), precipitation seasonality (0.33), and soil pH (0.26). However, the variables average precipitation (0.86), average precipitation of wettest month (0.85), Net primary productivity (0.76), water vapor pressure (0.47), tree density (0.45), minimum temperature of coldest month (0.40), canopy height (0.38), precipitation of driest month (0.37), isothermality (0.30), and precipitation seasonality (0.08) showed a strong positive relationship on axis 2 of the CCA plot (Figure 4).

4. Discussion

Studying the patterns and drivers of assemblage of bumblebee species on a local and fine scale is crucial in enabling local-scale management [20]. Climate is the main factor driving species distributions and community assembly in bees in Gansu (Table S2 and Figure 4) [40,41]. Within our analysis, all bumblebee species showed major responses to climate variables, except for Bombus grahami, though this may be due to undersampling, as this area represents the northern limit of its range [41].
The variable importance varied between species, for example, the most influential factor, Net primary productivity (NPP), was found to be influential in 87% of bumblebee species and the top driver for 65% of bumblebee species. However, NPP was not an important factor for B. atripes, B. difficillimus, B. festivus, B. semenovi, B. supremus, B. tanguticus, and B. tibetanus. Some of these species, such as B. atripes (a rare species in South China) and B. festivus, are southern species that inhabit humid regions with forests. Thus, the data here may be insufficient for analysis, as we are only exploring the distribution on the edge of their range [40]. However, whilst NPP is important, relationships are nuanced, especially in areas without significant tree-cover such as Gansu (and non-forest NPP is known to be a more significant driver than NPP overall) [20].
Solar radiation was selected as a variable for 72% of bumblebee species, possibly due to its influence on plant growth and bee thermoregulation [42]. Wind speed also contributes to the distribution of 77% of bumblebee species, and water vapor pressure contributes to the distribution of 79% of bumblebee species and may impact flight control of bumblebee species [43]. Tree canopy height was also incorporated in our study, because dense wooded systems have different floral provisioning and structures than open habitats [44].
The highest species richness was observed in Lintan County, in the transition zone between the western Qinghai-Tibetan Plateau and the Eastern Loess plateau. This area is characterized by varied topography and diverse habitats, encompassing extensive natural environments with mountains and grasslands to agricultural areas. This combination likely contributes to a high turnover in species presence. Areas with the highest species richness (>20) come under the climatic zone of temperate monsoon, and highly seasonal habitats can often support high diversity [45]. Further, it has been found that temperate areas support higher bumblebee richness than tropical areas [14]. Gansu is a temperate region. But part of southernmost Gansu is sub-tropical.
In our study, the Gansu was divided into four bumblebee biogeographic zones (Figure 3b). Similar studies have divided the whole of China into four bumblebee biogeographic regions, though we increased resolution to explore distribution on more regional scales [46]. Within China, the four biogeographic zones of Gansu fall under the North China and Tibetan Plateau regions of bumblebee groups. At a larger scale, the world was divided into 10 principal biogeographic regions of bumblebee species [1], and our four zones of Gansu fall under the category of the northern Oriental region. Another study categorized the bumblebee species of the Tibetan plateau into three principal groups, a distinct group of Himalayan faunas in the south, an interior Tibetan fauna, and the Qinghai and Gansu fauna of the northeast [19]. The Himalayan faunas and Qinghai and Gansu faunas showed closer association with climatic factors, whereas the interior Tibetan fauna showed weaker association [19].
Sampling bumblebees in Gansu is challenging due to the region’s complex topography and arid environment in the northwest, which limited sampling in some areas. In arid areas, species tend to be patchily distributed, making collection efforts more difficult [19]. Consequently, certain grid cells of dry mountains and deserts lacked bumblebee samples, and further research is needed for a more comprehensive understanding of regional distributions (Figure 1).
This study explains the most basic elements of distributional dynamics and the assemblage of bumblebee species at a regional scale, which are unknown for many insects [47]. The information about the faunal composition and contribution of different factors enables us to consider those factors while designing conservation strategies and management decisions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects15070552/s1, Figure S1: Habitat suitability modeling of bumblebee species of Gansu. X-axis represent the longitude and Y-axis are latitude. Here the values closer to 1 showed more suitability range of a species; Table S1: The environmental variables utilized in the investigation of the bumblebee biotic com-munity assemblage. The spatial resolution is approximately 1 km2; Table S2: Area under the curve values at training and test data, fractional predicted area and the highest contribution of environmental factors in the spatial distribution modeling of modelled bumblebee species of Gansu and its surrounding areas. Here, ‘PS’ stands for precipitation seasonality, ‘NPP’, Net primary productivity, ‘vapr’, vapor pressure, ‘TS’, temperature seasonality, ‘MaxTWM’, maximum temperature of warmest month, ‘AP’ average precipitation, ‘PWM’ precipitation of wettest month. Refs. [39,48] are mentioned in Supplementary Materials.

Author Contributions

Conceptualization, M.N., J.H. and J.A.; methodology, M.N., J.A. and A.C.H.; software, M.N., A.C.H., N.H.B. and Z.M.; validation, J.A., H.C., M.N., W.L., Z.M. and N.H.B.; formal analysis, H.C. and P.H.W.; investigation, J.A., Z.M. and J.H.; resources, J.A.; species taxonomy and identification, J.H. and P.H.W.; data curation, J.H.; writing—original draft preparation, M.N. and H.C.; writing—review and editing, M.N., J.A., H.C., A.C.H., N.H.B., J.H., P.H.W., W.L. and Z.M.; visualization, J.A.; supervision, J.A.; project administration, J.A.; funding acquisition, J.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Project of Discovering and Pollination Use of Bumblebees of Gansu (GIA-IAR-2017), National Natural Science Foundation of China (31672500), the China Agriculture Research System-Bee (CARS-44-KXJ5), the Yunnan Provincial Department of Science and Technology “Yunnan Talent Program” (202303AM140032), the China Postdoctoral Science Foundation (2022M723447), the Special Basic Cooperative Research Innovation Programs of QSTB and QNU (KJLH2023ZD05), and Innovative Research Team at Qujing Normal University.

Data Availability Statement

All the data generated or analyzed during this study are included in this published article and its Supplementary Materials.

Acknowledgments

We are very thankful to Shiwen Zhang, Aiping Han, Wenzhong Qi, Guiqian Zhang, Jingping Xi, Guowei Zhao, Shouli Liu, Shudong Luo, Zhiyong Zhou, Yazhou Zhao, Wenfeng Chen, Xiaolong Yuan, Guangshu Zhang, Ping Liu, Hong Zhang, Yanjie Liu, Jiyao Liu, Huiyue Zhao, and Shisheng Cao for their help in the field survey of bumblebees of Gansu..

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Williams, P.H. Mapping variations in the strength and breadth of biogeographic transition zones using species turnover. Proc. R. Soc. B 1996, 263, 579–588. [Google Scholar]
  2. Williams, P.H. Can biogeography help bumblebee conservation? Eur. J. Taxon. 2023, 890, 165–183. [Google Scholar] [CrossRef]
  3. Ernst, R.; Rödel, M.-O. Patterns of community composition in two tropical tree frog assemblages: Separating spatial structure and environmental effects in disturbed and undisturbed forests. J. Trop. Ecol. 2008, 24, 111–120. [Google Scholar] [CrossRef]
  4. Gardner, T.A.; Ribeiro-Júnior, M.A.; Barlow, J.O.S.; Ávila-Pires, T.C.S.; Hoogmoed, M.S.; Peres, C.A. The value of primary, secondary, and plantation forests for a neotropical herpetofauna. Conserv. Biol. 2007, 21, 775–787. [Google Scholar] [CrossRef] [PubMed]
  5. Holt, R.D. On the evolutionary ecology of species’ ranges. Evol. Ecol. Res. 2003, 5, 159–178. [Google Scholar]
  6. Yu, D.W.; Wilson, H.B.; Pierce, N.E. An empirical model of species coexistence in a spatially structured environment. Ecology 2001, 82, 1761–1771. [Google Scholar] [CrossRef]
  7. Labidi, I.; Errouissi, F.; Nouira, S. Spatial and temporal variation in species composition, diversity, and structure of mediterranean dung beetle assemblages (Coleoptera: Scarabaeidae) across a bioclimatic gradient. Environ. Entomol. 2012, 41, 785–801. [Google Scholar] [CrossRef]
  8. Brown, C.; Burslem, D.; Illian, J.B.; Bao, L.; Brockelman, W.; Cao, M.; Chang, L.W.; Dattaraja, H.S.; Davies, S.; Gunatilleke, C.V.S. Multispecies coexistence of trees in tropical forests: Spatial signals of topographic niche differentiation increase with environmental heterogeneity. Proc. R. Soc. B Biol. Sci. 2013, 280, 20130502. [Google Scholar] [CrossRef] [PubMed]
  9. Williams, P.H.; Osborne, J.L. Bumblebee vulnerability and conservation world-wide. Apidologie 2009, 40, 367–387. [Google Scholar] [CrossRef]
  10. Cameron, S.A.; Sadd, B.M. Global trends in bumble bee health. Annu. Rev. Entomol. 2020, 65, 209–232. [Google Scholar] [CrossRef]
  11. Butler, R.G.; Lage, C.; Dobrin, S.E.; Staples, J.K.; Venturini, E.; Frank, J.; Drummond, F.A. Maine’s Bumble bee (Hymenoptera: Apidae) assemblage—Part 1: Composition, seasonal and regional distribution, and resource use. Environ. Entomol. 2021, 50, 1344–1357. [Google Scholar] [CrossRef] [PubMed]
  12. Christman, M.E.; Spears, L.R.; Strange, J.P.; Pearse, W.D.; Burchfield, E.K.; Ramirez, R.A. Land cover and climate drive shifts in Bombus assemblage composition. Agric. Ecosyst. Environ. 2022, 339, 108113. [Google Scholar] [CrossRef]
  13. Tommasi, N.; Pioltelli, E.; Biella, P.; Labra, M.; Casiraghi, M.; Galimberti, A. Effect of urbanization and its environmental stressors on the intraspecific variation of flight functional traits in two bumblebee species. Oecologia 2022, 199, 289–299. [Google Scholar] [CrossRef] [PubMed]
  14. Williams, P.H. An annotated checklist of bumble bees with an analysis of patterns of description (Hymenoptera: Apidae, Bombini). Bull. Hist. Museum Entomol. Ser. 1998, 67, 79–152. [Google Scholar]
  15. Williams, P.H.; Huang, J.; An, J. Bear wasps of the middle kingdom: A decade of discovering China’s bumblebees. Antenna 2017, 41, 21–24. [Google Scholar]
  16. USDA. Forest Service Bumblebees (Bombus Spp.). Available online: https://www.fs.usda.gov/wildflowers/pollinators/pollinator-of-the-month/bumblebees.shtml (accessed on 11 April 2024).
  17. Williams, P. Habitat use by bumble bees (Bombus Spp.). Ecol. Entomol. 1988, 13, 223–237. [Google Scholar] [CrossRef]
  18. Williams, P.H.; Thorp, R.W.; Richardson, L.L.; Colla, S.R. Bumble Bees of North America: An Identification Guide; Princeton University Press: Princeton, NJ, USA, 2014; ISBN 1400851181. [Google Scholar]
  19. Williams, P.H.; Bystriakova, N.; Huang, J.; Miao, Z.; An, J. Bumblebees, climate and glaciers across the Tibetan Plateau (Apidae: Bombus Latreille). Syst. Biodivers. 2015, 13, 164–181. [Google Scholar] [CrossRef]
  20. Orr, M.C.; Hughes, A.C.; Chesters, D.; Pickering, J.; Zhu, C.D.; Ascher, J.S. Global patterns and drivers of bee distribution. Curr. Biol. 2021, 31, 451–458. [Google Scholar] [CrossRef]
  21. An, J.; Williams, P.H.; Zhou, B.; Miao, Z.; Qi, W. The bumblebees of Gansu, Northwest China (Hymenoptera, Apidae). Zootaxa 2011, 2865, 1–36. [Google Scholar] [CrossRef]
  22. Panfilov, D.V. On the geographical distribution of bumblebees (Bombus) in China. Acta Geogr. Sin. 1957, 23, 221–239. [Google Scholar]
  23. Wang, S. Hymenoptera: Apidae-Bombus. In Insects of Xizang; Science Press: Beijing, China, 1982; pp. 427–447. [Google Scholar]
  24. Williams, P.; Tang, Y.; Yao, J.; Cameron, S. The Bumblebees of Sichuan (Hymenoptera: Apidae, Bombini). Syst. Biodivers. 2009, 7, 101–189. [Google Scholar] [CrossRef]
  25. Kreft, H.; Jetz, W. A framework for delineating biogeographical regions based on species distributions. J. Biogeogr. 2010, 37, 2029–2053. [Google Scholar] [CrossRef]
  26. Huang, J.; An, J. Species diversity, pollination application and strategy for conservation of the bumblebees of China. Biodivers. Sci. 2018, 26, 486–497. [Google Scholar] [CrossRef]
  27. He, J.; Kreft, H.; Gao, E.; Wang, Z.; Jiang, H. Patterns and drivers of zoogeographical regions of terrestrial vertebrates in China. J. Biogeogr. 2017, 44, 1172–1184. [Google Scholar] [CrossRef]
  28. Miao, Y.F.; Herrman, M.; Wu, F.L.; Yan, X.L.; Yang, S.L. What controlled mid-late miocene long-term aridifi- cation in central asia? Global cooling or Tibetan Plateau uplift: A review. Earth-Sci. Rev. 2013, 112, 155–172. [Google Scholar] [CrossRef]
  29. Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 2010, 19, 134–143. [Google Scholar] [CrossRef]
  30. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’hara, R.B.; Simpson, G.L.; Solymos, P. Vegan: Community ecology package. R Package Version 2.5-6 2019. Community Ecol. 2020, 8, 732–740. [Google Scholar]
  31. Rahman, I.U.; Khan, N.; Ali, K. Classification and ordination of understory vegetation using multivariate techniques in the Pinus wallichiana Forests of Swat. Sci. Nat. 2017, 104, 24. [Google Scholar] [CrossRef] [PubMed]
  32. Borcard, D.; Gillet, F.; Legendre, P.; Borcard, D.; Gillet, F.; Legendre, P. Canonical ordination. In Numerical Ecology with R; Springer: Cham, Switzerland, 2018; pp. 203–297. [Google Scholar]
  33. Legendre, P.; Legendre, L. Numerical Ecology; Elsevier: Amsterdam, The Netherlands, 2012; ISBN 0444538690. [Google Scholar]
  34. National Research Council; Division on Earth, Life Studies; Board on Life Sciences; Committee on the Status of Pollinators in North America. Status of Pollinators in North America; National Academies Press: Cambridge, MA, USA, 2007; ISBN 0309102898. [Google Scholar]
  35. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Modell. 2006, 190, 231–259. [Google Scholar] [CrossRef]
  36. Brown, J.L. SDMtoolbox: A python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol. Evol. 2014, 5, 694–700. [Google Scholar] [CrossRef]
  37. Stein, A.; Gerstner, K.; Kreft, H. Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales. Ecol. Lett. 2014, 17, 866–880. [Google Scholar] [CrossRef] [PubMed]
  38. Elith, J.; Leathwick, J.R. Species distribution models: Ecological explanation and prediction across space and time. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 677–697. [Google Scholar] [CrossRef]
  39. FAO (Food and Agriculture Organization of the United Nations). FAO Global Information System on Water and Agriculture (Aquastat)—Evaporation. Available online: https://data.apps.fao.org/catalog/dataset/dce0415f-fd48-468f-a197-e4603dbf88dc/resource/be9143d9-034b-46c9-8b8d-b92f82110dc2?inner_span=True (accessed on 11 April 2024).
  40. Naeem, M.; Liu, M.; Huang, J.; Ding, G.; Potapov, G. Vulnerability of east asian bumblebee species to future climate and land cover changes. Agric. Ecosyst. Environ. 2019, 277, 11–20. [Google Scholar] [CrossRef]
  41. Williams, P.H.; An, J.; Dorji, P.; Huang, J.; Jaffar, S.; Japoshvili, G.; Narah, J.; Ren, Z.; Streinzer, M.; Thanoosing, C. Bumblebees with big teeth: Revising the subgenus alpigenobombus with the good, the bad and the ugly of numts (hymenoptera: Apidae). Eur. J. Taxon. 2023, 892, 1–65. [Google Scholar] [CrossRef]
  42. Willmer, P.G.; Stone, G.N. Behavioral, ecological, and physiological determinants of the activity patterns of bees. Adv. Study Behav. 2004, 34, 347–466. [Google Scholar]
  43. Ravi, S.; Kolomenskiy, D.; Engels, T.; Schneider, K.; Wang, C.; Sesterhenn, J.; Liu, H. Bumblebees minimize control challenges by combining active and passive modes in unsteady winds. Sci. Rep. 2016, 6, 35043. [Google Scholar] [CrossRef]
  44. Donkersley, P. Trees for bees. Agric. Ecosyst. Environ. 2019, 270, 79–83. [Google Scholar] [CrossRef]
  45. Zhao, P.; He, Z. Temperature change characteristics in gansu province of China. Atmosphere 2022, 13, 728. [Google Scholar] [CrossRef]
  46. Naeem, M.; Huang, J.; Zhang, S.; Luo, S.; Liu, Y.; Zhang, H.; Luo, Q.; Zhou, Z.; Ding, G.; An, J. Diagnostic indicators of wild pollinators for biodiversity monitoring in long-term conservation. Sci. Total Environ. 2020, 708, 135231. [Google Scholar] [CrossRef]
  47. Stork, N.E. How many species of insects and other terrestrial arthropods are there on earth? Annu. Rev. Entomol. 2018, 63, 31–45. [Google Scholar] [CrossRef]
  48. Crowther, T.W.; Glick, H.B.; Covey, K.R.; Bettigole, C.; Maynard, D.S.; Thomas, S.M.; Smith, J.R.; Hintler, G.; Duguid, M.C.; Amatulli, G.; et al. Global Map of Tree Density. 2016. Available online: https://figshare.com/articles/dataset/Global_map_of_tree_density/3179986 (accessed on 11 April 2024).
Insects 15 00552 g001
Figure 1. A map of Gansu, with the distribution of bumblebee sampling sites within Gansu and its surroundings. Blue dots represent the collection sites of bumblebees.
Figure 1. A map of Gansu, with the distribution of bumblebee sampling sites within Gansu and its surroundings. Blue dots represent the collection sites of bumblebees.
Insects 15 00552 g001
Insects 15 00552 g002
Figure 2. The principal component analysis of 240 grid cells based on 17 environmental variables for bumblebee distribution and assemblage represents the importance of the contributions of all variables in the development of environmental heterogeneity.
Figure 2. The principal component analysis of 240 grid cells based on 17 environmental variables for bumblebee distribution and assemblage represents the importance of the contributions of all variables in the development of environmental heterogeneity.
Insects 15 00552 g002
Insects 15 00552 g003
Figure 3. The classification of 240 grid cells dominated by different bumblebee species into four groups (zones) using Ward’s agglomerative clustering analysis, and phenon lines represent the biogeographic groups. Here, I represents the Southern humid zone, II represents the Western Qilian snow mountain zone, III represents the Eastern Loess plateau zone, and IV represents the Western dry mountain zone (a), and the hotspot regions of each group’s bumblebee species are shown with colorful grid cells in Gansu of China (b).
Figure 3. The classification of 240 grid cells dominated by different bumblebee species into four groups (zones) using Ward’s agglomerative clustering analysis, and phenon lines represent the biogeographic groups. Here, I represents the Southern humid zone, II represents the Western Qilian snow mountain zone, III represents the Eastern Loess plateau zone, and IV represents the Western dry mountain zone (a), and the hotspot regions of each group’s bumblebee species are shown with colorful grid cells in Gansu of China (b).
Insects 15 00552 g003
Insects 15 00552 g004
Figure 4. Canonical Correspondence Analysis (CCA) plot illustrating the relationship between bumblebee species composition and environmental variables within Gansu. Each point represents a sampling site within Gansu, and their positions are determined by the scores on axis 1 and axis 2. The length and direction of lines depict the importance and direction of each variable’s influence on the species composition. Here, “NPP” represents the Net primary productivity, “srad” represents the solar radiation, “vapr” represents the vapor pressure, “AMP” represents annual mean temperature, “MDR” represents mean diurnal range, “TS” represents temperature seasonality, “MaxTWM” represents maximum temperature of wettest month, “MinTCM” represents minimum temperature of coldest month, “AP” represents annual precipitation, “APWM” represents annual precipitation of wettest month, and “PDM” represents precipitation of driest month. Here, green, red, pink, and blue dots represent the Southern humid zone, the Western Qilian snow mountain zone, the Eastern Loess plateau zone, and the Western dry mountain zone, respectively.
Figure 4. Canonical Correspondence Analysis (CCA) plot illustrating the relationship between bumblebee species composition and environmental variables within Gansu. Each point represents a sampling site within Gansu, and their positions are determined by the scores on axis 1 and axis 2. The length and direction of lines depict the importance and direction of each variable’s influence on the species composition. Here, “NPP” represents the Net primary productivity, “srad” represents the solar radiation, “vapr” represents the vapor pressure, “AMP” represents annual mean temperature, “MDR” represents mean diurnal range, “TS” represents temperature seasonality, “MaxTWM” represents maximum temperature of wettest month, “MinTCM” represents minimum temperature of coldest month, “AP” represents annual precipitation, “APWM” represents annual precipitation of wettest month, and “PDM” represents precipitation of driest month. Here, green, red, pink, and blue dots represent the Southern humid zone, the Western Qilian snow mountain zone, the Eastern Loess plateau zone, and the Western dry mountain zone, respectively.
Insects 15 00552 g004
Insects 15 00552 g005
Figure 5. Species richness of bumblebees within Gansu of China.
Figure 5. Species richness of bumblebees within Gansu of China.
Insects 15 00552 g005
Table 1. Bumblebee species of different subgenera collected from Gansu and its surrounding areas. The species with “*” were excluded in the habitat suitability modeling because of their records < 3.
Table 1. Bumblebee species of different subgenera collected from Gansu and its surrounding areas. The species with “*” were excluded in the habitat suitability modeling because of their records < 3.
Sr. NoSpeciesSubgenus
1B. asiaticusSibiricobombus
2B. atripesThoracobombus
3B. bellardiiPsithyrus
4B. bicoloratusMegabombus
5B. bohemicusPsithyrus
6B. brevicepsAlpigenobombus
7B. campestrisPsithyrus
8B. chinensisPsithyrus
9B. consobrinusMegabombus
10B. convexusMendacibombus
11B. coreanus *Psithyrus
12B. cornutusPsithyrus
13B. czerskii *Megabombus
14B. deuteronymusThoracobombus
15B. difficillimusSubterraneobombus
16B. expolitus *Psithyrus
17B. festivusMelanobombus
18B. filchneraeThoracobombus
19B. flavescensPyrobombus
20B. grahamiAlpigenobombus
21B. hediniThoracobombus
22B. humilisThoracobombus
23B. hengduanensisPyrobombus
24B. ignitusBombus
25B. imitator *Thoracobombus
26B. impetuosusThoracobombus
27B. infrequens *Pyrobombus
28B. kashmirensisAlpigenobombus
29B. qilianensisMelanobombus
30B. koreanusMegabombus
31B. ladakhensis *Melanobombus
32B. laesusThoracobombus
33B. lantschouensisBombus
34B. lemniscatusPyrobombus
35B. lepidusPyrobombus
36B. minshanicolaBombus
37B. longipesMelanobombus
38B. melanurusSubterraneobombus
39B. minshanensisBombus
40B. norvegicus *Psithyrus
41B. opulentusThoracobombus
42B. ganjsuensisBombus
43B. personatusSubterraneobombus
44B. picipesPyrobombus
45B. pyrosomaMelanobombus
46B. religious *Megabombus
47B. remotusThoracobombus
48B. prshewalskyiMelanobombus
49B. rupestris *Psithyrus
50B. semenoviSibiricobombus
51B. sibiricusSibiricobombus
52B. sicheliiMelanobombus
53B. skorikoviPsithyrus
54B. supremusMegabombus
55B. sushkiniMegabombus
56B. tanguticusMelanobombus
57B. tibetanusPsithyrus
58B. trifasciatusMegabombus
59B. turneri *Psithyrus
60B. validusAlpigenobombus
61B. waltoniMendacibombus
62B. wangaePyrobombus
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Naeem, M.; Chen, H.; Li, W.; Hughes, A.C.; Williams, P.H.; Bashir, N.H.; Miao, Z.; Huang, J.; An, J. Patterns and Drivers of Bumblebee Diversity in Gansu. Insects 2024, 15, 552. https://doi.org/10.3390/insects15070552

AMA Style

Naeem M, Chen H, Li W, Hughes AC, Williams PH, Bashir NH, Miao Z, Huang J, An J. Patterns and Drivers of Bumblebee Diversity in Gansu. Insects. 2024; 15(7):552. https://doi.org/10.3390/insects15070552

Chicago/Turabian Style

Naeem, Muhammad, Huanhuan Chen, Wenbo Li, Alice C. Hughes, Paul H. Williams, Nawaz Haider Bashir, Zhengying Miao, Jiaxing Huang, and Jiandong An. 2024. "Patterns and Drivers of Bumblebee Diversity in Gansu" Insects 15, no. 7: 552. https://doi.org/10.3390/insects15070552

APA Style

Naeem, M., Chen, H., Li, W., Hughes, A. C., Williams, P. H., Bashir, N. H., Miao, Z., Huang, J., & An, J. (2024). Patterns and Drivers of Bumblebee Diversity in Gansu. Insects, 15(7), 552. https://doi.org/10.3390/insects15070552

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop