Climate-land-use interactions shape tropical mountain biodiversity and ecosystem functions.

Agriculture and the exploitation of natural resources have transformed tropical mountain ecosystems across the world, and the consequences of these transformations for biodiversity and ecosystem functioning are largely unknown1-3. Conclusions that are derived from studies in non-mountainous areas are not suitable for predicting the effects of land-use changes on tropical mountains because the climatic environment rapidly changes with elevation, which may mitigate or amplify the effects of land use4,5. It is of key importance to understand how the interplay of climate and land use constrains biodiversity and ecosystem functions to determine the consequences of global change for mountain ecosystems. Here we show that the interacting effects of climate and land use reshape elevational trends in biodiversity and ecosystem functions on Africa's largest mountain, Mount Kilimanjaro (Tanzania). We find that increasing land-use intensity causes larger losses of plant and animal species richness in the arid lowlands than in humid submontane and montane zones. Increases in land-use intensity are associated with significant changes in the composition of plant, animal and microorganism communities; stronger modifications of plant and animal communities occur in arid and humid ecosystems, respectively. Temperature, precipitation and land use jointly modulate soil properties, nutrient turnover, greenhouse gas emissions, plant biomass and productivity, as well as animal interactions. Our data suggest that the response of ecosystem functions to land-use intensity depends strongly on climate; more-severe changes in ecosystem functioning occur in the arid lowlands and the cold montane zone. Interactions between climate and land use explained-on average-54% of the variation in species richness, species composition and ecosystem functions, whereas only 30% of variation was related to single drivers. Our study reveals that climate can modulate the effects of land use on biodiversity and ecosystem functioning, and points to a lowered resistance of ecosystems in climatically challenging environments to ongoing land-use changes in tropical mountainous regions.

Hill-Chao numbers allow decomposing gamma multifunctionality into alpha and beta components.

Biodiversity-ecosystem functioning (BEF) research has provided strong evidence and mechanistic underpinnings to support positive effects of biodiversity on ecosystem functioning, from single to multiple functions. This research has provided knowledge gained mainly at the local alpha scale (i.e. within ecosystems), but the increasing homogenization of landscapes in the Anthropocene has raised the potential that declining biodiversity at the beta (across ecosystems) and gamma scales is likely to also impact ecosystem functioning. Drawing on biodiversity theory, we propose a new statistical framework based on Hill-Chao numbers. The framework allows decomposition of multifunctionality at gamma scales into alpha and beta components, a critical but hitherto missing tool in BEF research; it also allows weighting of individual ecosystem functions. Through the proposed decomposition, new BEF results for beta and gamma scales are discovered. Our novel approach is applicable across ecosystems and connects local- and landscape-scale BEF assessments from experiments to natural settings.

Ant-following behavior is correlated with plumage traits in African understory birds.

Ant-following behavior is a common phenomenon in birds of Neotropical and Afrotropical rainforests but yet little is known from Central Africa. We here report on the phenomenon in lowland rainforest in Cameroon, quantifying the strength of the interaction of different ant-following bird species with driver ants and test the hypothesis that higher levels of specialization in ant-following behavior are associated with dominance or aggression-dependent plumage and other morphological traits. Flock size varied between 1 and 11 individuals with a mean size of 5.34 ± 2.68 (mean ± SD) individuals occurring at the same time. The maximum number of species present during one raid observed was ten, whereas the minimum number was four with an overall species richness of 6.89 ± 2.1 species. The 21 attending bird species strongly varied in the degree of ant-following behavior. In an interspecific comparison, plumage traits such as the presence of a colored crown, eyespots, and bare skin around the eye, in combination with metatarsus length and weight, were significantly correlated with ant-following behavior. These results suggest that-in size and identity of species-ant-following bird assemblages in Central Africa are similar to those reported from East Africa. They also suggest that ant following favors the selection of traits that signal dominance in interactions between individuals struggling for valuable food resources in the forest understory.

Enhancing the structural diversity between forest patches-A concept and real-world experiment to study biodiversity, multifunctionality and forest resilience across spatial scales.

Intensification of land use by humans has led to a homogenization of landscapes and decreasing resilience of ecosystems globally due to a loss of biodiversity, including the majority of forests. Biodiversity-ecosystem functioning (BEF) research has provided compelling evidence for a positive effect of biodiversity on ecosystem functions and services at the local (α-diversity) scale, but we largely lack empirical evidence on how the loss of between-patch ß-diversity affects biodiversity and multifunctionality at the landscape scale (γ-diversity). Here, we present a novel concept and experimental framework for elucidating BEF patterns at α-, ß-, and γ-scales in real landscapes at a forest management-relevant scale. We examine this framework using 22 temperate broadleaf production forests, dominated by Fagus sylvatica. In 11 of these forests, we manipulated the structure between forest patches by increasing variation in canopy cover and deadwood. We hypothesized that an increase in landscape heterogeneity would enhance the ß-diversity of different trophic levels, as well as the ß-functionality of various ecosystem functions. We will develop a new statistical framework for BEF studies extending across scales and incorporating biodiversity measures from taxonomic to functional to phylogenetic diversity using Hill numbers. We will further expand the Hill number concept to multifunctionality allowing the decomposition of γ-multifunctionality into α- and ß-components. Combining this analytic framework with our experimental data will allow us to test how an increase in between patch heterogeneity affects biodiversity and multifunctionality across spatial scales and trophic levels to help inform and improve forest resilience under climate change. Such an integrative concept for biodiversity and functionality, including spatial scales and multiple aspects of diversity and multifunctionality as well as physical and environmental structure in forests, will go far beyond the current widely applied approach in forestry to increase resilience of future forests through the manipulation of tree species composition.

The role of diversity, body size and climate in dung removal: A correlative and experimental approach.

The mechanisms by which climatic changes influence ecosystem functions, that is, by a direct climatic control of ecosystem processes or by modifying richness and trait compositions of species communities, remain unresolved. This study is a contribution to this discourse by elucidating the linkages between climate, land use, biodiversity, body size and ecosystem functions. We disentangled direct climatic from biodiversity-mediated effects by using dung removal by dung beetles as a model system and by combining correlative field data and exclosure experiments along an extensive elevational gradient on Mt. Kilimanjaro, Tanzania. Dung removal declined with increasing elevation, being associated with a strong reduction in the richness and body size traits of dung beetle communities. Climate influenced dung removal rates by modifying biodiversity rather than by direct effects. The biodiversity-ecosystem effect was driven by a change in the mean body size of dung beetles. Dung removal rates were strongly reduced when large dung beetles were experimentally excluded. This study underscores that climate influences ecosystem functions mainly by modifying biodiversity and underpins the important role of body size for dung removal.

Increasing the phylogenetic coverage for understanding broad-scale diversity gradients.

Despite decades of scientific effort, there is still no consensus on the determinants of broad-scale gradients of animal diversity. We argue that general drivers of diversity are unlikely to be found among the narrowly defined taxa which are typically analyzed in studies of broad-scale diversity gradients because ecological niches evolve largely conservatively. This causes constraints in the use of available niche space leading to systematic differences in diversity gradients among taxa. We instead advocate studies of phylogenetically diverse animal communities along broad environmental gradients. Such multi-taxa communities are less constrained in resource use and diversification and may be better targets for testing major classical hypotheses on diversity gradients. Besides increasing the spatial scale in analyses, expanding the phylogenetic coverage may be a second way to achieve higher levels of generality in studies of broad-scale diversity gradients.

Primary productivity and habitat protection predict elevational species richness and community biomass of large mammals on Mt. Kilimanjaro.

Despite their diversity and their large functional and cultural importance, the patterns and predictors of large mammal diversity along elevational gradients on tropical mountains remain poorly understood. Today, large mammals are threatened by human disturbances such as habitat destruction and hunting and may increasingly depend on the conservation of protected areas. Here, we use field data on the diversity of large mammals along a 3.6 km elevational gradient on Mt. Kilimanjaro to evaluate the importance of climate, net primary productivity and human impact for the distribution, species richness and community biomass of wild mammals. Mammal species richness was explored with camera traps on 66 study plots along an elevational gradient from 870 to 4,550 m a.s.l.. We applied path analysis and variance partitioning analysis to unravel the direct and indirect effects of temperature, precipitation, primary productivity, land use, land area, the protection of habitats and the occurrence of domestic mammals on the diversity of wild mammals. Both species richness and community biomass of wild mammals showed a unimodal distribution with elevation, peaking in the montane zone of Mt. Kilimanjaro. However, the peak shifted significantly to lower elevations when only protected habitats were considered. Wild mammal diversity increased with net primary productivity, protection of habitats and temperature. Our study underscores the importance of energy resources for the establishment of diversity gradients in large mammals. While temperature has been revealed as a direct predictor of diversity in most ectothermic taxa, in endothermic organisms temperature has stronger indirect effects, via a modulation of net primary productivity. Moreover, our study reveals how patterns of diversity on tropical mountains are influenced by human impact, pointing to the pivotal role of protected areas for the long-term conservation of mountain biodiversity.

Leaf traits mediate changes in invertebrate herbivory along broad environmental gradients on Mt. Kilimanjaro, Tanzania.

Temperature, primary productivity, plant functional traits, and herbivore abundances are considered key predictors of leaf herbivory but their direct and indirect contributions to community-level herbivory are not well understood along broad climatic gradients. Here, we determined elevational herbivory patterns and used a path analytical approach to disentangle the direct and indirect effects of climate, land use, net primary productivity (NPP), herbivore abundance, and plant functional traits on community-level invertebrate herbivory along the extensive elevational and land use gradients at Mt. Kilimanjaro, Tanzania. We recorded standing leaf herbivory caused by leaf chewers, leaf miners and leaf gallers on 55 study sites distributed in natural and anthropogenic habitats along a 3,060 m elevation gradient. We related the total community-level herbivory to climate (temperature and precipitation), NPP, plant functional traits (specific leaf area, leaf carbon-to-nitrogen [CN] ratio and leaf nitrogen-to-phosphorus [NP] ratio) and herbivore abundances. Leaf herbivory ranged from 5% to 11% along the elevation gradient. Total leaf herbivory showed unimodal pattern in natural habitats but a strongly contrasting bimodal pattern in anthropogenic habitats. We also detected some variation in the patterns of leaf herbivory along environmental gradients across feeding guilds with leaf chewers being responsible for a disproportionally large part of herbivory. Path analyses indicated that the variation in leaf herbivory was mainly driven by changes in leaf CN and NP ratios which were closely linked to changes in NPP in natural habitats. Similarly, patterns of leaf herbivory in anthropogenic habitats were best explained by variation in leaf CN ratios and a negative effect of land use. Our study elucidates the strong role of leaf nutrient stoichiometry and its linkages to climate and NPP for explaining the variation in leaf herbivory along broad climatic gradients. Furthermore, the study suggests that climatic changes and nutrient inputs in the course of land use change may alter leaf herbivory and consequently energy and nutrient fluxes in terrestrial habitats.

Small-scale topography modulates elevational α-, ß- and γ-diversity of Andean leaf beetles.

Elevational diversity gradients are typically studied without considering the complex small-scale topography of large mountains, which generates habitats of strongly different environmental conditions within the same elevational zones. Here we analyzed the importance of small-scale topography for elevational diversity patterns of hyperdiverse tropical leaf beetles (Coleoptera: Chrysomelidae). We compared patterns of elevational diversity and species composition of beetles in two types of forests (on mountain ridges and in valleys) and analyzed whether differences in the rate of species turnover among forest habitats lead to shifts in patterns of elevational diversity when scaling up from the local study site to the elevational belt level. We sampled beetle assemblages at 36 sites in the Podocarpus National Park, Ecuador, which were equally distributed over two forest habitats and three elevational levels. DNA barcoding and Poisson tree processes modelling were used to delimitate putative species. On average, local leaf beetle diversity showed a clear hump-shaped pattern. However, only diversity in forests on mountain ridges peaked at mid-elevation, while beetle diversity in valleys was similarly high at low- and mid-elevation and only declined at highest elevations. A higher turnover of species assemblages at lower than at mid-elevations caused a shift from a hump-shaped diversity pattern found at the local level to a low-elevation plateau pattern (with similar species numbers at low and mid-elevation) at the elevational belt level. Our study reveals an important role of small-scale topography and spatial scale for the inference on gradients of elevational species diversity.

Complementary ecosystem services provided by pest predators and pollinators increase quantity and quality of coffee yields.

Wild animals substantially support crop production by providing ecosystem services, such as pollination and natural pest control. However, the strengths of synergies between ecosystem services and their dependencies on land-use management are largely unknown. Here, we took an experimental approach to test the impact of land-use intensification on both individual and combined pollination and pest control services in coffee production systems at Mount Kilimanjaro. We established a full-factorial pollinator and vertebrate exclosure experiment along a land-use gradient from traditional homegardens (agroforestry systems), shaded coffee plantations to sun coffee plantations (total sample size = 180 coffee bushes). The exclusion of vertebrates led to a reduction in fruit set of ca 9%. Pollinators did not affect fruit set, but significantly increased fruit weight of coffee by an average of 7.4%. We found no significant decline of these ecosystem services along the land-use gradient. Pest control and pollination service were thus complementary, contributing to coffee production by affecting the quantity and quality of a major tropical cash crop across different coffee production systems at Mount Kilimanjaro.

Elevated ozone and carbon dioxide affects the composition of volatile organic compounds emitted by Vicia faba (L.) and visitation by European orchard bee (Osmia cornuta).

Recent studies link increased ozone (O3) and carbon dioxide (CO2) levels to alteration of plant performance and plant-herbivore interactions, but their interactive effects on plant-pollinator interactions are little understood. Extra floral nectaries (EFNs) are essential organs used by some plants for stimulating defense against herbivory and for the attraction of insect pollinators, e.g., bees. The factors driving the interactions between bees and plants regarding the visitation of bees to EFNs are poorly understood, especially in the face of global change driven by greenhouse gases. Here, we experimentally tested whether elevated levels of O3 and CO2 individually and interactively alter the emission of Volatile Organic Compound (VOC) profiles in the field bean plant (Vicia faba, L., Fabaceae), EFN nectar production and EFN visitation by the European orchard bee (Osmia cornuta, Latreille, Megachilidae). Our results showed that O3 alone had significant negative effects on the blends of VOCs emitted while the treatment with elevated CO2 alone did not differ from the control. Furthermore, as with O3 alone, the mixture of O3 and CO2 also had a significant difference in the VOCs' profile. O3 exposure was also linked to reduced nectar volume and had a negative impact on EFN visitation by bees. Increased CO2 level, on the other hand, had a positive impact on bee visits. Our results add to the knowledge of the interactive effects of O3 and CO2 on plant volatiles emitted by Vicia faba and bee responses. As greenhouse gas levels continue to rise globally, it is important to take these findings into consideration to better prepare for changes in plant-insect interactions.

Seasonal and elevational changes of plant-pollinator interaction networks in East African mountains.

Across an elevation gradient, several biotic and abiotic factors influence community assemblages of interacting species leading to a shift in species distribution, functioning, and ultimately topologies of species interaction networks. However, empirical studies of climate-driven seasonal and elevational changes in plant-pollinator networks are rare, particularly in tropical ecosystems. Eastern Afromontane Biodiversity Hotspots in Kenya, East Africa. We recorded plant-bee interactions at 50 study sites between 515 and 2600 m asl for a full year, following all four major seasons in this region. We analysed elevational and seasonal network patterns using generalised additive models (GAMs) and quantified the influence of climate, floral resource availability, and bee diversity on network structures using a multimodel inference framework. We recorded 16,741 interactions among 186 bee and 314 plant species of which a majority involved interactions with honeybees. We found that nestedness and bee species specialisation of plant-bee interaction networks increased with elevation and that the relationships were consistent in the cold-dry and warm-wet seasons respectively. Link rewiring increased in the warm-wet season with elevation but remained indifferent in the cold-dry seasons. Conversely, network modularity and plant species were more specialised at lower elevations during both the cold-dry and warm-wet seasons, with higher values observed during the warm-wet seasons. We found flower and bee species diversity and abundance rather than direct effects of climate variables to best predict modularity, specialisation, and link rewiring in plant-bee-interaction networks. This study highlights changes in network architectures with elevation suggesting a potential sensitivity of plant-bee interactions with climate warming and changes in rainfall patterns along the elevation gradients of the Eastern Afromontane Biodiversity Hotspot.

Seasonal and environmental dynamics of intra-urban freshwater habitats and their influence on the abundance of Bulinus snail host of Schistosoma haematobium in the Tiko endemic focus, Mount Cameroon region.

Urogenital schistosomiasis (UGS) caused by Schistosoma haematobium is endemic in the South West Region of Cameroon. An understanding of the abundance and distribution of the Bulinus snail, intermediate host can inform strategic snail control programmes at a local scale. This study investigated seasonal dynamics and environmental factors influencing occurrence and abundance of freshwater snail intermediate hosts in Tiko, a semi-urban endemic focus in the Mount Cameroon area. A longitudinal malacological field survey was conducted between December 2019 and December 2020 in the Tiko municipality. Snails were collected for one year monthly at 12 different human water contact sites along a stretch of the Ndongo stream using a standardized sampling technique. Freshwater snails were identified using shell morphological features. In addition, water temperature, pH, electrical conductivity, total dissolved solutes, salinity, water depth, width and flow velocity were measured, and vegetation cover as well as substrate type were determined. Bayesian regression models were used to identify the main environmental factors affecting the occurrence and abundance of Bulinus intermediate host. In total, 2129 fresh water snails were collected during the study period. Physa (51.4%) was the most abundant genus followed by Melanoides (28.6%) then, Bulinus (15.5%), Lymnaea (4.2%), Indoplanorbis (0.2%) and Potadoma (0.1%). Seasonality in abundance was significant in Bulinus sp as well as other genera, with greater numbers in the dry season (peaks between December and February). Water temperature, a rocky or sandy substrate type associated positively with Bulinus sp, meanwhile a higher water flow rate and medium vegetation negatively influenced the snail intermediate host population. These findings underscore the importance of timing behavioural and snail control interventions against schistosomiasis as well as increase vigilance of other trematode diseases in the study area. The continuous spread of planorbid snail hosts is a major concern.

A dataset of occurrence of wild bees and their interaction with foraging plants along a livestock grazing gradient of northern Tanzania.

A dataset describing the occurrence of wild bees and their interaction with forage plants along livestock grazing gradient is critical in understanding bee-plant interaction networks and in developing conservation plans to ensure ecosystem services in human-modified landscapes. Despite this need, bee-plant datasets are scarce in Africa, and Tanzania is no exception. Therefore, in this article, we present a dataset of wild bee species richness, occurrence, and distribution collected across sites with different levels of livestock grazing intensity and forage resources thereby. The data presented in this paper supports a research article by Lasway et al., 2022 describing the effects of grazing intensity on East African bee assemblages. The paper presents primary data on bee species, collection method, date of collection, bee family, identifier, plant forage resource, forage plant life form, forage plant family, location (GPS coordinates), grazing intensity category, mean annual temperature (°C), and elevation (m asl). The data were collected intermittently between August 2018 and March 2020 from 24 study sites distributed along three levels of livestock grazing intensity with eight replicates for each: low, moderate, and high livestock grazing intensity. In each study site, two 50 × 50 m study plots were set from which bees and floral resources were sampled and quantified. The two plots were placed in a way to capture the overall structural heterogeneity of the respective habitat by placing the two plots in contrasting microhabitats where possible. For example, in moderately livestock-grazed habitats, plots were placed on sites with and without tree or shrub cover to ensure representativeness. This paper presents a dataset comprising 2,691 bee individuals from 183 species representing 55 genera of the five bee families: Halictidae (74), Apidae (63), Megachilidae (40), Andrenidae (5), and Colletidae (1). In addition, the dataset comprises 112 species of flowering plants that were identified as potential forage resources for bees. This paper supplements rare but critical data on bee pollinators in Northern Tanzania and advances our knowledge of the potential drivers of bee-pollinator whose populations diversity are declining globally. The dataset will also promote collaborations among researchers who would wish to combine and extend their data for further analysis to gain a broader understanding of the phenomenon on a larger spatial scale.

Temperature and livestock grazing trigger transcriptome responses in bumblebees along an elevational gradient.

Climate and land-use changes cause increasing stress to pollinators but the molecular pathways underlying stress responses are poorly understood. Here, we analyzed the transcriptomic response of Bombus lucorum workers to temperature and livestock grazing. Bumblebees sampled along an elevational gradient, and from differently managed grassland sites (livestock grazing vs unmanaged) in the German Alps did not differ in the expression of genes known for thermal stress responses. Instead, metabolic energy production pathways were upregulated in bumblebees sampled in mid- or high elevations or during cool temperatures. Extensive grazing pressure led to an upregulation of genetic pathways involved in immunoregulation and DNA-repair. We conclude that widespread bumblebees are tolerant toward temperature fluctuations in temperate mountain environments. Moderate temperature increases may even release bumblebees from metabolic stress. However, transcriptome responses to even moderate management regimes highlight the completely underestimated complexity of human influence on natural pollinators.

Interactive effects of climate and land use on pollinator diversity differ among taxa and scales.

Changes in climate and land use are major threats to pollinating insects, an essential functional group. Here, we unravel the largely unknown interactive effects of both threats on seven pollinator taxa using a multiscale space-for-time approach across large climate and land-use gradients in a temperate region. Pollinator community composition, regional gamma diversity, and community dissimilarity (beta diversity) of pollinator taxa were shaped by climate-land-use interactions, while local alpha diversity was solely explained by their additive effects. Pollinator diversity increased with reduced land-use intensity (forest < grassland < arable land < urban) and high flowering-plant diversity at different spatial scales, and higher temperatures homogenized pollinator communities across regions. Our study reveals declines in pollinator diversity with land-use intensity at multiple spatial scales and regional community homogenization in warmer and drier climates. Management options at several scales are highlighted to mitigate impacts of climate change on pollinators and their ecosystem services.

Hybridization in East African swarm-raiding army ants.

BACKGROUND: Hybridization can have complex effects on evolutionary dynamics in ants because of the combination of haplodiploid sex-determination and eusociality. While hybrid non-reproductive workers have been found in a range of species, examples of gene-flow via hybrid queens and males are rare. We studied hybridization in East African army ants (Dorylus subgenus Anomma) using morphology, mitochondrial DNA sequences, and nuclear microsatellites. RESULTS: While the mitochondrial phylogeny had a strong geographic signal, different species were not recovered as monophyletic. At our main study site at Kakamega Forest, a mitochondrial haplotype was shared between a "Dorylus molestus-like" and a "Dorylus wilverthi-like" form. This pattern is best explained by introgression following hybridization between D. molestus and D. wilverthi. Microsatellite data from workers showed that the two morphological forms correspond to two distinct genetic clusters, with a significant proportion of individuals being classified as hybrids. CONCLUSIONS: We conclude that hybridization and gene-flow between the two army ant species D. molestus and D. wilverthi has occurred, and that mating between the two forms continues to regularly produce hybrid workers. Hybridization is particularly surprising in army ants because workers have control over which males are allowed to mate with a young virgin queen inside the colony.

Cryptic species and hidden ecological interactions of halictine bees along an elevational gradient.

Changes of abiotic and biotic conditions along elevational gradients represent serious challenges to organisms which may promote the turnover of species, traits and biotic interaction partners. Here, we used molecular methods to study cuticular hydrocarbon (CHC) profiles, biotic interactions and phylogenetic relationships of halictid bees of the genus Lasioglossum along a 2,900 m elevational gradient at Mt. Kilimanjaro, Tanzania. We detected a strong species turnover of morphologically indistinguishable taxa with phylogenetically clustered cryptic species at high elevations, changes in CHC profiles, pollen resource diversity, and a turnover in the gut and body surface microbiome of bees. At high elevations, increased proportions of saturated compounds in CHC profiles indicate physiological adaptations to prevent desiccation. More specialized diets with higher proportions of low-quality Asteraceae pollen imply constraints in the availability of food resources. Interactive effects of climatic conditions on gut and surface microbiomes, CHC profiles, and pollen diet suggest complex feedbacks among abiotic conditions, ecological interactions, physiological adaptations, and phylogenetic constraints as drivers of halictid bee communities at Mt. Kilimanjaro.

Plant traits mediate the effects of climate on phytophagous beetle diversity on Mt. Kilimanjaro.

Patterns of insect diversity along elevational gradients are well described in ecology. However, it remains little tested how variation in the quantity, quality, and diversity of food resources influence these patterns. Here we analyzed the direct and indirect effects of climate, food quantity (estimated by net primary productivity), quality (variation in the specific leaf area index, leaf nitrogen to phosphorus and leaf carbon to nitrogen ratio), and food diversity (diversity of leaf traits) on the species richness of phytophagous beetles along the broad elevation and land use gradients of Mt. Kilimanjaro, Tanzania. We sampled beetles at 65 study sites located in both natural and anthropogenic habitats, ranging from 866 to 4,550 m asl. We used path analysis to unravel the direct and indirect effects of predictor variables on species richness. In total, 3,154 phytophagous beetles representing 19 families and 304 morphospecies were collected. We found that the species richness of phytophagous beetles was bimodally distributed along the elevation gradient with peaks at the lowest (˜866 m asl) and upper mid-elevations (˜3,200 m asl) and sharply declined at higher elevations. Path analysis revealed temperature- and climate-driven changes in primary productivity and leaf trait diversity to be the best predictors of changes in the species richness of phytophagous beetles. Species richness increased with increases in mean annual temperature, primary productivity, and with increases in the diversity of leaf traits of local ecosystems. Our study demonstrates that, apart from temperature, the quantity and diversity of food resources play a major role in shaping diversity gradients of phytophagous insects. Drivers of global change, leading to a change of leaf traits and causing reductions in plant diversity and productivity, may consequently reduce the diversity of herbivore assemblages.

A synopsis of the Bee occurrence data of northern Tanzania.

BACKGROUND: Bees (Hymenoptera: Apoidea: Anthophila) are the most important group of pollinators with about 20,507 known species worldwide. Despite the critical role of bees in providing pollination services, studies aiming at understanding which species are present across disturbance gradients are scarce. Limited taxononomic information for the existing and unidentified bee species in Tanzania make their conservation haphazard. Here, we present a dataset of bee species records obtained from a survey in nothern Tanzania i.e. Kilimanjaro, Arusha and Manyara regions. Our findings serve as baseline data necessary for understanding the diversity and distribution of bees in the northern parts of the country, which is a critical step in devising robust conservation and monitoring strategies for their populations. NEW INFORMATION: In this paper, we present information on 45 bee species belonging to 20 genera and four families sampled using a combination of sweep-netting and pan trap methods. Most species (27, ~ 60%) belong to the family Halictidae followed by 16 species (35.5%) from the family Apidae. Megachilidae and Andrenidae were the least represented, each with only one species (2.2%). Additional species of Apidae and Megachilidae sampled during this survey are not yet published on Global Biodiversity Information Facility (GBIF), once they will be available on GBIF, they will be published in a subsequent paper. From a total of 953 occurrences, highest numbers were recorded in Kilimanjaro Region (n = 511), followed by Arusha (n = 410) and Manyara (n = 32), but this pattern reflects the sampling efforts of the research project rather than real bias in the distributions of bee species in northern Tanzania.