Derek S. Sikes

Predicting global change effects poses significant challenges due to the intricate interplay between climate change and anthropogenic stressors in shaping ecological communities and their function, such as pest outbreak risk. Termites are ecosystem engineers, yet some pest species are causing worldwide economic losses. While habitat fragmentation seems to drive pest-dominated termite communities, its interaction with climate change effect remains unknown. We test whether climate and habitat fragmentation interactively alter interspecific competition that may limit pest termite risk. Leveraging global termite co-occurrence including 280 pest species, we found that competitively superior termite species (e.g., large bodied) increased in large and continuous habitats solely at high precipitation. While competitive species suppressed pest species globally, habitat fragmentation drove pest termite risk only in humid biomes. Unfortunately, hu- mid tropics have experienced vast forest fragmentation and rainfall reduction over the past decades. These stressors, if not stopped, may drive pest termite risk, potentially via competitive release.

BACKGROUND Kissing bugs are the vectors of Trypanosoma cruzi, the etiological agent of Chagas disease (CD). Despite their epidemiological relevance, kissing bug species are under sampled in terms of their diversity and it is unclear what biases exist in available kissing bug data. Under climate change, range maps for kissing bugs may become less accurate as species shift their ranges to track climatic tolerance. OBJECTIVES Quantify inventory completeness in available kissing bug data. Assess how well range maps are at conveying information about current distributions and potential future distributions subject to shift under climate change. Intersect forecasted changes in kissing bug distributions with contemporary sampling gaps to identify regions for future sampling of the group. Identify whether a phylogenetic signal is present in expert range knowledge as more closely related species may be similarly well or lesser understood. METHODS We used species distribution models (SDM), specifically constructed from Bayesian additive regression trees, with Bioclim variables, to forecast kissing bug distributions into 2100 and intersect these with current sampling gaps to identify priority regions for sampling. Expert range maps were assessed by the agreement between the expert map and SDM generated occurrence probability. We used classical hypothesis testing methods as well as tests of phylogenetic signal to meet our objectives. FINDINGS Expert range maps vary in their quality of depicting current kissing bug distributions. Most expert range maps decline in their ability to convey information about kissing bug occurrence over time, especially in under sampled areas. We found limited evidence for a phylogenetic signal in expert range map performance. MAIN CONCLUSIONS Expert range maps are not a perfect account of species distributions and may degrade in their ability to accurately convey distribution knowledge under future climates. We identify regions where future sampling of kissing bugs will be crucial for completing biodiversity inventories.

We analysed the wild bee community sampled from 1921 to 2018 at a nature preserve in southern Michigan, USA, to study long-term community shifts in a protected area. During an intensive survey in 1972 and 1973, Francis C. Evans detected 135 bee species. In the most recent intensive surveys conducted in 2017 and 2018, we recorded 90 species. Only 58 species were recorded in both sampling periods, indicating a significant shift in the bee community. We found that the bee community diversity, species richness and evenness were all lower in recent samples. Additionally, 64% of the more common species exhibited a more than 30% decline in relative abundance. Neural network analysis of species traits revealed that extirpation from the reserve was most likely for oligolectic ground-nesting bees and kleptoparasitic bees, whereas polylectic cavity-nesting bees were more likely to persist. Having longer phenological ranges also increased the chance of persistence in polylectic species. Further analysis suggests a climate response as bees in the contemporary sampling period had a more southerly overall distribution compared to the historic community. Results exhibit the utility of both long-term data and machine learning in disentangling complex indicators of bee population trajectories.

Sarcophagid experts have made several efforts to associate biodiversity data and comprehend where each species occurs, but comprehensive faunal inventories remain scarce. Our aim was to provide a list of distributional patterns and endemic species and allow assessment of the sampling effort conducted within Brazilian biomes. We produced a dataset of Brazilian sarcophagids and overlaid with a biome map, to investigate distributional patterns, endemism and to build species accumulation curves. Additionally, we calculated nonparametric asymptotic species richness estimators and extrapolation of species diversity (Hill numbers). Our dataset comprised 288 sarcophagid species, which 21 were identified as endemic. The biomes with the highest species richness were the Atlantic Rainforest and the Amazon Forest, and no biome exhibited a stabilized asymptotic curve. This is the first proposal of listing Sarcophagidae species by biomes and essential to understand the spatial distribution of this family in Brazil. We present maps and richness estimators that allow identifying gaps and guiding survey planning.

Growing degree‐days (GDD), an estimate of an organism's growing season length, has been shown to be an important predictor of Lepidopteran species' distributions and could be influencing Lepidopteran range shifts to climate change. Yet, one understudied simplification in this literature is that the same thermal threshold is used in the calculations of GDD for all species instead of a species‐specific threshold. By characterizing the phenological process influenced by climate, a species‐specific estimate of GDD should improve the accuracy of species distribution models (SDMs). To test this hypothesis, we used published, experimentally estimated thermal thresholds and modeled the current geographic distribution of 30 moth species native to North America. We found that the predictive performance of models based on a species‐specific estimate of GDD was indistinguishable from models based on a standard estimate of GDD. This is likely because GDD was not an important predictor of these species' distributions. Our findings suggest that experimentally estimated thermal thresholds may not always scale up to be predictive at broad scales and that more work is needed to leverage the data from lab experiments into SDMs to accurately predict species' range shifts in response to climate change.

Species are often expected to shift their distributions either poleward or upslope to evade warming climates and colonise new suitable climatic niches. However, from 18‐years of fixed transect monitoring data on 88 species of butterfly in the midwestern United States, we show that butterflies are shifting their centroids in all directions, except towards regions that are warming the fastest (southeast).Butterflies shifted their centroids at a mean rate of 4.87 km year−1. The rate of centroid shift was significantly associated with local climate change velocity (temperature by precipitation interaction), but not with mean climate change velocity throughout the species' ranges.Species tended to shift their centroids at a faster rate towards regions that are warming at slower velocities but increasing in precipitation velocity.Surprisingly, species' thermal niche breadth (range of climates butterflies experience throughout their distribution) and wingspan (often used as metric for dispersal capability) were not correlated with the rate at which species shifted their ranges.We observed high phylogenetic signal in the direction species shifted their centroids. However, we found no phylogenetic signal in the rate species shifted their centroids, suggesting less conserved processes determine the rate of range shift than the direction species shift their ranges.This research shows important signatures of multidirectional range shifts (latitudinal and longitudinal) and uniquely shows that local climate change velocities are more important in driving range shifts than the mean climate change velocity throughout a species' entire range.

Insects are the most diverse group of organisms, but their large number of species and the lack of specialists to study them have made this group particularly vulnerable to the main limitations in biological diversity, such as the Wallacean deficit. This study will contribute to the geographical knowledge of an insect trophic guild that has been widely used as an indicator group, the Scarabaeoidea and Silphidae copro-necrophagous beetles, emphasizing their geographical distribution in two Mexican national parks (Iztaccíhuatl-Popocatepetl and La Malinche) and the intermediate region, which includes sky-island ecosystems in central Mexico. Geographic records of the 32 species that have been previously recorded in the study region were compiled and used to generate potential distribution models aiming to generate potential alpha (species richness) and beta (total beta diversity, nestedness and replacement) diversity maps. The greatest predicted species richness was found between 3,000 and 3,500 m a.s.l. in the study region. Potential species richness ranged from 2 to 24 species. Total beta diversity was low in the study region (mean 0.1), while nestedness was the most important component of beta diversity (0.8). The maximum alpha and beta diversity values were predicted outside the national parks. Consequently, we consider that the studied national parks are not able to protect completely the regional alpha and beta diversities by themselves. Implications for insect conservation: Our results show that the highest alfa and beta diversity values of copro-necrophagous beetles might occur outside the national parks, and a suitable way to protect them could be the Archipelago reserve model as an alternative to protect the regional diversity.

Understanding and mitigating the effects of anthropogenic climate change on species distributions requires the ability to track range shifts over time. This is particularly true for species occupying high-latitude regions, which are experiencing more extreme climate change than the rest of the world. In North America, the geographic ranges of many mammals reach their northernmost extent in Alaska, positioning this region at the leading edge of climate-induced distribution change. Over a decade has elapsed since the publication of the last spatial assessments of terrestrial mammals in the state. We compared public occurrence records against commonly referenced range maps to evaluate potential extralimital records and develop repeatable baseline range maps. We compared occurrence records from the Global Biodiversity Information Facility for 61 terrestrial mammal species native to mainland Alaska against a variety of range estimates (International Union for Conservation of Nature, Alaska Gap Analysis Project, and the published literature). We mapped extralimital records and calculated proportions of occurrences encompassed by range extents, measured mean direction and distance to prior range margins, evaluated predictive accuracy of published species models, and highlighted observations on federal lands in Alaska. Range comparisons identified 6,848 extralimital records for 39 of 61 (63.9%) terrestrial mainland Alaskan species. On average, 95.5% of Alaska Gap Analysis Project occurrence records and ranges were deemed accurate (i.e., > 90.0% correct) for 31 of 37 species, but overestimated extents for 13 species. The International Union for Conservation of Nature range maps encompassed 68.1% of occurrence records and were > 90% accurate for 17 of 39 species. Extralimital records represent either improved sampling and digitization or actual geographic range expansions. Here we provide new data-driven range maps, update standards for the archiving of museum-quality locational records and offer recommendations for mapping range changes for monitoring and conservation.

Background The recent discovery of emerging relapsing fever group Borrelia (RFGB) species, such as Borrelia miyamotoi, poses a growing threat to public health. However, the global distribution and associated risk burden of these species remain uncertain. We aimed to map the diversity, distribution, and potential infection risk of RFGB.MethodsWe searched PubMed, Web of Science, GenBank, CNKI, and eLibrary from Jan 1, 1874, to Dec 31, 2022, for published articles without language restriction to extract distribution data for RFGB detection in vectors, animals, and humans, and clinical information about human patients. Only articles documenting RFGB infection events were included in this study, and data for RFGB detection in vectors, animals, or humans were composed into a dataset. We used three machine learning algorithms (boosted regression trees, random forest, and least absolute shrinkage and selection operator logistic regression) to assess the environmental, ecoclimatic, biological, and socioeconomic factors associated with the occurrence of four major RFGB species: Borrelia miyamotoi, Borrelia lonestari, Borrelia crocidurae, and Borrelia hermsii; and mapped their worldwide risk level.FindingsWe retrieved 13 959 unique studies, among which 697 met the selection criteria and were used for data extraction. 29 RFGB species have been recorded worldwide, of which 27 have been identified from 63 tick species, 12 from 61 wild animals, and ten from domestic animals. 16 RFGB species caused human infection, with a cumulative count of 26 583 cases reported from Jan 1, 1874, to Dec 31, 2022. Borrelia recurrentis (17 084 cases) and Borrelia persica (2045 cases) accounted for the highest proportion of human infection. B miyamotoi showed the widest distribution among all RFGB, with a predicted environmentally suitable area of 6·92 million km2, followed by B lonestari (1·69 million km2), B crocidurae (1·67 million km2), and B hermsii (1·48 million km2). The habitat suitability index of vector ticks and climatic factors, such as the annual mean temperature, have the most significant effect among all predictive models for the geographical distribution of the four major RFGB species.InterpretationThe predicted high-risk regions are considerably larger than in previous reports. Identification, surveillance, and diagnosis of RFGB infections should be prioritised in high-risk areas, especially within low-income regions.FundingNational Key Research and Development Program of China.

This paper presents an exploratory study on the taxonomic diversity of pre-Hispanic archaeofaunas in the South-Central Andes (SCA; western South America) from the Pleistocene-Holocene boundary to the Late-Holocene. The SCA is a complex of diverse environments and has undergone distinct climate events for the last 13,000 years, such as the occurrence of warmer and drier conditions in the Middle-Holocene. The South-Central Andean area was part of the larger Andes interaction area, which was a primary center for animal and plant domestication and the emergence of agro-pastoralist economies. Since subsistence was key to these processes, the SCA provides a relevant case study on the interactions among environment, foodways and sociocultural evolution. Taxonomic diversity was used here as a proxy for diet breadth. A total of 268 archaeofaunal assemblages were sampled from the zooarchaeological literature. Reviewed variables included the cultural chronology and spatial coordinates of the assemblages, as well as the presence and abundance of taxa at the family rank. Taxonomic diversity covered two dimensions: composition (families present in each assemblage) and structure (quantitative relationships among taxa), which was measured through richness (NTAXA), ubiquity and relative abundance (NISP based rank-order). Despite the uneven distribution of samples, the analyses revealed the following trends: (1) a moderate relationship between NTAXA and distance from coastline for most of the Holocene; (2) a potential decrease in assemblage richness for coastal ecoregions during the Late-Holocene; and (3) a generalized increase in the relative abundance of Camelidae.