Evolutionary Origins and Patterns of Diversification in Animal Brood Parasitism.

AbstractBrood parasitism involves the exploitation of host parental care rather than the extraction of resources directly from hosts. We identify defining characteristics of this strategy and consider its position along continua with adjacent behaviors but focus on canonical brood parasites, where parasitism is obligate and hosts are noneusocial (thereby distinguishing from social parasitism). A systematic literature survey revealed 59 independently derived brood parasitic lineages with most origins (49) in insects, particularly among bees and wasps, and other origins in birds (seven) and fish (three). Insects account for more than 98% of brood parasitic species, with much of that diversity reflecting ancient (≥100-million-year-old) brood parasitic lineages. Brood parasites usually, but not always, evolve from forms that show parental care. In insects, brood parasitism often first evolves through exploitation of a closely related species, following Emery's rule, but this is less typical in birds, which we discuss. We conducted lineage-level comparisons between brood parasitic clades and their sister groups, finding mixed results but an overall neutral to negative effect of brood parasitism on species richness and diversification. Our review of brood parasites reveals many unanswered questions requiring new research, including further modeling of the coevolutionary dynamics of brood parasites and their hosts.

Phylogenetic relationships and the evolution of host preferences in the largest clade of brood parasitic bees (Apidae: Nomadinae).

Brood parasites (also known as cleptoparasites) represent a substantial fraction of global bee diversity. Rather than constructing their own nests, these species instead invade those of host bees to lay their eggs. Larvae then hatch and consume the food provisions intended for the host's offspring. While this life history strategy has evolved numerous times across the phylogeny of bees, the oldest and most speciose parasitic clade is the subfamily Nomadinae (Apidae). However, the phylogenetic relationships among brood parasitic apids both within and outside the Nomadinae have not been fully resolved. Here, we present new findings on the phylogeny of this diverse group of brood parasites based on ultraconserved element (UCE) sequence data and extensive taxon sampling with 114 nomadine species representing all tribes. We suggest a broader definition of the subfamily Nomadinae to describe a clade that includes almost all parasitic members of the family Apidae. The tribe Melectini forms the sister group to all other Nomadinae, while the remainder of the subfamily is composed of two sister clades: a "nomadine line" representing the former Nomadinae sensu stricto, and an "ericrocidine line" that unites several mostly Neotropical lineages. We find the tribe Osirini Handlirsch to be polyphyletic, and divide it into three lineages, including the newly described Parepeolini trib. nov. In addition to our taxonomic findings, we use our phylogeny to explore the evolution of different modes of parasitism, detecting two independent transitions from closed-cell to open-cell parasitism. Finally, we examine how nomadine host-parasite associations have evolved over time. In support of Emery's rule, which suggests close relationships between hosts and parasites, we confirm that the earliest nomadines were parasites of their close free-living relatives within the family Apidae, but that over time their host range broadened to include more distantly related hosts spanning the diversity of bees. This expanded breadth of host taxa may also be associated with the transition to open-cell parasitism.

Gene Tree Estimation Error with Ultraconserved Elements: An Empirical Study on Pseudapis Bees.

Summarizing individual gene trees to species phylogenies using two-step coalescent methods is now a standard strategy in the field of phylogenomics. However, practical implementations of summary methods suffer from gene tree estimation error, which is caused by various biological and analytical factors. Greatly understudied is the choice of gene tree inference method and downstream effects on species tree estimation for empirical data sets. To better understand the impact of this method choice on gene and species tree accuracy, we compare gene trees estimated through four widely used programs under different model-selection criteria: PhyloBayes, MrBayes, IQ-Tree, and RAxML. We study their performance in the phylogenomic framework of $>$800 ultraconserved elements from the bee subfamily Nomiinae (Halictidae). Our taxon sampling focuses on the genus Pseudapis, a distinct lineage with diverse morphological features, but contentious morphology-based taxonomic classifications and no molecular phylogenetic guidance. We approximate topological accuracy of gene trees by assessing their ability to recover two uncontroversial, monophyletic groups, and compare branch lengths of individual trees using the stemminess metric (the relative length of internal branches). We further examine different strategies of removing uninformative loci and the collapsing of weakly supported nodes into polytomies. We then summarize gene trees with ASTRAL and compare resulting species phylogenies, including comparisons to concatenation-based estimates. Gene trees obtained with the reversible jump model search in MrBayes were most concordant on average and all Bayesian methods yielded gene trees with better stemminess values. The only gene tree estimation approach whose ASTRAL summary trees consistently produced the most likely correct topology, however, was IQ-Tree with automated model designation (ModelFinder program). We discuss these findings and provide practical advice on gene tree estimation for summary methods. Lastly, we establish the first phylogeny-informed classification for Pseudapis s. l. and map the distribution of distinct morphological features of the group. [ASTRAL; Bees; concordance; gene tree estimation error; IQ-Tree; MrBayes, Nomiinae; PhyloBayes; RAxML; phylogenomics; stemminess].

Wild insect diversity increases inter-annual stability in global crop pollinator communities.

While an increasing number of studies indicate that the range, diversity and abundance of many wild pollinators has declined, the global area of pollinator-dependent crops has significantly increased over the last few decades. Crop pollination studies to date have mainly focused on either identifying different guilds pollinating various crops, or on factors driving spatial changes and turnover observed in these communities. The mechanisms driving temporal stability for ecosystem functioning and services, however, remain poorly understood. Our study quantifies temporal variability observed in crop pollinators in 21 different crops across multiple years at a global scale. Using data from 43 studies from six continents, we show that (i) higher pollinator diversity confers greater inter-annual stability in pollinator communities, (ii) temporal variation observed in pollinator abundance is primarily driven by the three-most dominant species, and (iii) crops in tropical regions demonstrate higher inter-annual variability in pollinator species richness than crops in temperate regions. We highlight the importance of recognizing wild pollinator diversity in agricultural landscapes to stabilize pollinator persistence across years to protect both biodiversity and crop pollination services. Short-term agricultural management practices aimed at dominant species for stabilizing pollination services need to be considered alongside longer term conservation goals focussed on maintaining and facilitating biodiversity to confer ecological stability.

Omnivory in Bees: Elevated Trophic Positions among All Major Bee Families.

As pollen and nectar foragers, bees have long been considered strictly herbivorous. Their pollen provisions, however, are host to abundant microbial communities, which feed on the pollen before and/or while it is consumed by bee larvae. In the process, microbes convert pollen into a complex of plant and microbial components. Since microbes are analogous to metazoan consumers within trophic hierarchies, the pollen-eating microbes are, functionally, herbivores. When bee larvae consume a microbe-rich pollen complex, they ingest proteins from plant and microbial sources and thus should register as omnivores on the trophic "ladder." We tested this hypothesis by examining the isotopic compositions of amino acids extracted from native bees collected in North America over multiple years. We measured bee trophic position across the six major bee families. Our findings indicate that bee trophic identity was consistently and significantly higher than that of strict herbivores, providing the first evidence that omnivory is ubiquitous among bee fauna. Such omnivory suggests that pollen-borne microbes represent an important protein source for larval bees, which introduces new questions as to the link between floral fungicide residues and bee development.

Viral transmission in honey bees and native bees, supported by a global black queen cell virus phylogeny.

In recent decades, we have realized that honey bee viruses are not, in fact, exclusive to honey bees. The potential impact of Apis-affiliated viruses on native pollinators is prompting concern. Our research addresses the issue of virus crossover between honey bees and native bees foraging in the same localities. We measured the presence of black queen cell virus (BQCV), deformed wing virus (DWV) and sacbrood virus (SBV) in managed Apis mellifera (honey bees) and native Andrena spp. (subgenus Melandrena) bee populations in five commercial orchards. We identified viral presence across sites and bees and related these data to measures of bee community diversity. All viruses were found in both managed and native bees, and BQCV was the most common virus in each. To establish evidence for viral crossover between taxa, we undertook an additional examination of BQCV where 74 samples were sequenced and placed in a global phylogenic framework of hundreds of BQCV strains. We demonstrate pathogen sharing across managed honey bees and distantly related wild bees. This phylogenetic analysis contributes to growing evidence for host switching and places local incidence patterns in a worldwide context, revealing multispecies viral transmission.

Combining transcriptomes and ultraconserved elements to illuminate the phylogeny of Apidae.

Two increasingly popular approaches to reconstruct the Tree of Life involve whole transcriptome sequencing and the target capture of ultraconserved elements (UCEs). Both methods can be used to generate large, multigene datasets for analysis of phylogenetic relationships in non-model organisms. While targeted exon sequencing across divergent lineages is now a standard method, it is still not clear if UCE data can be readily combined with published transcriptomes. In this study, we evaluate the combination of UCEs and transcriptomes in a single analysis using genome-, transcriptome-, and UCE data for 79 bees in the largest and most biologically diverse bee family, Apidae. Using existing tools, we first developed a workflow to assemble phylogenomic data from different sources and produced two large nucleotide matrices of combined data. We then reconstructed the phylogeny of the Apidae using concatenation- and coalescent-based methods, and critically evaluated the resulting phylogenies in the context of previously published genetic, genomic, and morphological data sets. Our estimated phylogenetic trees are robustly supported and largely congruent with previous molecular hypotheses, from deep nodes to shallow species-level phylogenies. Moreover, the combined approach allows us to resolve controversial nodes of the apid Tree of Life, by clarifying the relationships among the genera of orchid bees (Euglossini) and the monophyly of the Centridini. Additionally, we present novel phylogenetic evidence supporting the monophyly of the diverse clade of cleptoparasitic Apidae and the placement of two enigmatic, oil-collecting genera (Ctenoplectra and Tetrapedia). Lastly, we propose a revised classification of the family Apidae that reflects our improved understanding of apid higher-level relationships.

Pollinivory and the diversification dynamics of bees.

Pollinivory-the consumption of pollen rather than arthropod prey-is a defining feature of bees (Anthophila; the flower lovers). In virtually all bee species, larvae consume a diet composed of pollen mixed with nectar or floral oils. Bees arose from within a group of solitary, carnivorous, apoid wasps in the Early to Mid-Cretaceous, coincident with the rapid rise of flowering plants. It is assumed that the switch from carnivory to pollen-feeding was a key innovation that led to the rapid diversification of bees, but this has never been examined empirically. Here, we explore the hypothesis that pollinivory led to the increased diversification of bees. In contrast to common perception, we find that the switch to pollen-feeding per se does not explain their extensive diversification. Rather, our results indicate that pollinivory was a necessary but not sufficient condition for diversification, and that other complementary innovations, such as a broadening of host-plant diet, allowed the diversification of the major bee lineages. Our results have broad implications for understanding tempo and mode of bee diversification dynamics in light of their floral resources.

Long-legged bees make adaptive leaps: linking adaptation to coevolution in a plant-pollinator network.

Adaptation is evolution in response to natural selection. Hence, an adaptation is expected to originate simultaneously with the acquisition of a particular selective environment. Here we test whether long legs evolve in oil-collecting Rediviva bees when they come under selection by long-spurred, oil-secreting flowers. To quantify the selective environment, we drew a large network of the interactions between Rediviva species and oil-secreting plant species. The selective environment of each bee species was summarized as the average spur length of the interacting plant species weighted by interaction frequency. Using phylogenetically independent contrasts, we calculated divergence in selective environment and evolutionary divergence in leg length between sister species (and sister clades) of Rediviva We found that change in the selective environment explained 80% of evolutionary change in leg length, with change in body size contributing an additional 6% of uniquely explained variance. The result is one of four proposed steps in testing for plant-pollinator coevolution.

The impact of GC bias on phylogenetic accuracy using targeted enrichment phylogenomic data.

The field of sequence based phylogenetic analyses is currently being transformed by novel hybrid-based targeted enrichment methods, such as the use of ultraconserved elements (UCEs). Rather than analyzing relationships among organisms using a small number of genes, these methods now allow us to evaluate relationships with many hundreds to thousands of individual gene loci. However, the inclusion of thousands of loci does not necessarily overcome the long-standing challenge of incongruence among phylogenetic trees derived from different genes or gene regions. One factor that impacts the level of incongruence in phylogenomic data sets is the level of GC bias. GC rich gene regions are prone to higher recombination rates than AT rich regions, driven by a process referred to as "GC biased gene conversion". As a result, high GC content can be negatively associated with phylogenetic accuracy, but the extent to which this impacts incongruence among UCEs is currently unstudied. We investigated the impact of GC content on phylogeny reconstruction using in silico captured UCE data for the corbiculate bees (Hymenoptera: Apidae). The phylogeny of this group has been the subject of extensive study, and incongruence among gene trees is thought to be a source of phylogenetic error. We conducted coalescent- and concatenation-based analyses of 810 individual gene loci from all 13 currently available bee genomes, including 8 corbiculate taxa. Both coalescent- and concatenation-based methods converged on a single topology for the corbiculate tribes. In contrast to concatenation, the coalescent-based methods revealed significant topological conflict at nodes involving the orchid bees (Euglossini) and honeybees (Apini). Partitioning the loci by GC content reveals decreasing support for the inferred topology with increasing GC bias. Based on the results of this study, we report the first evidence that GC biased gene conversion may contribute to topological incongruence in studies based on ultraconserved elements.

Playing with extremes: Origins and evolution of exaggerated female forelegs in South African Rediviva bees.

Despite close ecological interactions between plants and their pollinators, only some highly specialised pollinators adapt to a specific host plant trait by evolving a bizarre morphology. Here we investigated the evolution of extremely elongated forelegs in females of the South African bee genus Rediviva (Hymenoptera: Melittidae), in which long forelegs are hypothesised to be an adaptation for collecting oils from the extended spurs of their Diascia host flowers. We first reconstructed the phylogeny of the genus Rediviva using seven genes and inferred an origin of Rediviva at around 29MYA (95% HPD=19.2-40.5), concurrent with the origin and radiation of the Succulent Karoo flora. The common ancestor of Rediviva was inferred to be a short-legged species that did not visit Diascia. Interestingly, all our analyses strongly supported at least two independent origins of long legs within Rediviva. Leg length was not correlated with any variable we tested (ecological specialisation, Diascia visitation, geographic distribution, pilosity type) but seems to have evolved very rapidly. Overall, our results indicate that foreleg length is an evolutionary highly labile, rapidly evolving trait that might enable Rediviva bees to respond quickly to changing floral resource availability.

Crop domestication facilitated rapid geographical expansion of a specialist pollinator, the squash bee Peponapis pruinosa.

Squash was first domesticated in Mexico and is now found throughout North America (NA) along with Peponapis pruinosa, a pollen specialist bee species of the squash genus Cucurbita The origin and spread of squash cultivation is well-studied archaeologically and phylogenetically; however, no study has documented how cultivation of this or any other crop has influenced species in mutualistic interactions. We used molecular markers to reconstruct the demographic range expansion and colonization routes of P. pruinosa from its native range into temperate NA. Populations east of the Rocky Mountains expanded from the wild host plant's range in Mexico and were established by a series of founder events. Eastern North America was most likely colonized from squash bee populations in the present-day continental Midwest USA and not from routes that followed the Gulf and Atlantic coasts from Mexico. Populations of P. pruinosa west of the Rockies spread north from the warm deserts much more recently, showing two genetically differentiated populations with no admixture: one in California and the other one in eastern Great Basin. These bees have repeatedly endured severe bottlenecks as they colonized NA, following human spread of their Cucurbita pollen hosts during the Holocene.

Phylogenetic systematics and a revised generic classification of anthidiine bees (Hymenoptera: Megachilidae).

The bee tribe Anthidiini (Hymenoptera: Megachilidae) is a large, cosmopolitan group of solitary bees that exhibit intriguing nesting behavior. We present the first molecular-based phylogenetic analysis of relationships within Anthidiini using model-based methods and a large, multi-locus dataset (five nuclear genes, 5081 base pairs), as well as a combined analysis using our molecular dataset in conjunction with a previously published morphological matrix. We discuss the evolution of nesting behavior in Anthidiini and the relationship between nesting material and female mandibular morphology. Following an examination of the morphological characters historically used to recognize anthidiine genera, we recommend the use of a molecular-based phylogenetic backbone to define taxonomic groups prior to the assignment of diagnostic morphological characters for these groups. Finally, our results reveal the paraphyly of numerous genera and have significant consequences for anthidiine classification. In order to promote a classification system based on stable, monophyletic clades, we hereby make the following changes to Michener's (2007) classification: The subgenera Afranthidium (Zosteranthidium) Michener and Griswold, 1994, Afranthidium (Branthidium) Pasteels, 1969 and Afranthidium (Immanthidium) Pasteels, 1969 are moved into the genus Pseudoanthidium, thus forming the new combinations Pseudoanthidium (Zosteranthidium), Pseudoanthidium (Branthidium), and Pseudoanthidium (Immanthidium). The genus Neanthidium Pasteels, 1969 is also moved into the genus Pseudoanthidium, thus forming the new combination Pseudoanthidium (Neanthidium). Based on morphological characters shared with our new definition of the genus Pseudoanthidium, the subgenus Afranthidium (Mesanthidiellum) Pasteels, 1969 and the genus Gnathanthidium Pasteels, 1969 are also moved into the genus Pseudoanthidium, thus forming the new combinations Pseudoanthidium (Mesanthidiellum) and Pseudoanthidium (Gnathanthidium).

Historical changes in northeastern US bee pollinators related to shared ecological traits.

Pollinators such as bees are essential to the functioning of terrestrial ecosystems. However, despite concerns about a global pollinator crisis, long-term data on the status of bee species are limited. We present a long-term study of relative rates of change for an entire regional bee fauna in the northeastern United States, based on >30,000 museum records representing 438 species. Over a 140-y period, aggregate native species richness weakly decreased, but richness declines were significant only for the genus Bombus. Of 187 native species analyzed individually, only three declined steeply, all of these in the genus Bombus. However, there were large shifts in community composition, as indicated by 56% of species showing significant changes in relative abundance over time. Traits associated with a declining relative abundance include small dietary and phenological breadth and large body size. In addition, species with lower latitudinal range boundaries are increasing in relative abundance, a finding that may represent a response to climate change. We show that despite marked increases in human population density and large changes in anthropogenic land use, aggregate native species richness declines were modest outside of the genus Bombus. At the same time, we find that certain ecological traits are associated with declines in relative abundance. These results should help target conservation efforts focused on maintaining native bee abundance and diversity and therefore the important ecosystems services that they provide.

Negative effects of pesticides on wild bee communities can be buffered by landscape context.

Wild bee communities provide underappreciated but critical agricultural pollination services. Given predicted global shortages in pollination services, managing agroecosystems to support thriving wild bee communities is, therefore, central to ensuring sustainable food production. Benefits of natural (including semi-natural) habitat for wild bee abundance and diversity on farms are well documented. By contrast, few studies have examined toxicity of pesticides on wild bees, let alone effects of farm-level pesticide exposure on entire bee communities. Whether beneficial natural areas could mediate effects of harmful pesticides on wild bees is also unknown. Here, we assess the effect of conventional pesticide use on the wild bee community visiting apple (Malus domestica) within a gradient of percentage natural area in the landscape. Wild bee community abundance and species richness decreased linearly with increasing pesticide use in orchards one year after application; however, pesticide effects on wild bees were buffered by increasing proportion of natural habitat in the surrounding landscape. A significant contribution of fungicides to observed pesticide effects suggests deleterious properties of a class of pesticides that was, until recently, considered benign to bees. Our results demonstrate extended benefits of natural areas for wild pollinators and highlight the importance of considering the landscape context when weighing up the costs of pest management on crop pollination services.

Climate, physiological tolerance and sex-biased dispersal shape genetic structure of Neotropical orchid bees.

Understanding the impact of past climatic events on the demographic history of extant species is critical for predicting species' responses to future climate change. Palaeoclimatic instability is a major mechanism of lineage diversification in taxa with low dispersal and small geographical ranges in tropical ecosystems. However, the impact of these climatic events remains questionable for the diversification of species with high levels of gene flow and large geographical distributions. In this study, we investigate the impact of Pleistocene climate change on three Neotropical orchid bee species (Eulaema bombiformis, E. meriana and E. cingulata) with transcontinental distributions and different physiological tolerances. We first generated ecological niche models to identify species-specific climatically stable areas during Pleistocene climatic oscillations. Using a combination of mitochondrial and nuclear markers, we inferred calibrated phylogenies and estimated historical demographic parameters to reconstruct the phylogeographical history of each species. Our results indicate species with narrower physiological tolerance experienced less suitable habitat during glaciations and currently exhibit strong population structure in the mitochondrial genome. However, nuclear markers with low and high mutation rates show lack of association with geography. These results combined with lower migration rate estimates from the mitochondrial than the nuclear genome suggest male-biased dispersal. We conclude that despite large effective population sizes and capacity for long-distance dispersal, climatic instability is an important mechanism of maternal lineage diversification in orchid bees. Thus, these Neotropical pollinators are susceptible to disruption of genetic connectivity in the event of large-scale climatic changes.

Climate-associated phenological advances in bee pollinators and bee-pollinated plants.

The phenology of many ecological processes is modulated by temperature, making them potentially sensitive to climate change. Mutualistic interactions may be especially vulnerable because of the potential for phenological mismatching if the species involved do not respond similarly to changes in temperature. Here we present an analysis of climate-associated shifts in the phenology of wild bees, the most important pollinators worldwide, and compare these shifts to published studies of bee-pollinated plants over the same time period. We report that over the past 130 y, the phenology of 10 bee species from northeastern North America has advanced by a mean of 10.4 ± 1.3 d. Most of this advance has taken place since 1970, paralleling global temperature increases. When the best available data are used to estimate analogous rates of advance for plants, these rates are not distinguishable from those of bees, suggesting that bee emergence is keeping pace with shifts in host-plant flowering, at least among the generalist species that we investigated.

The impact of molecular data on our understanding of bee phylogeny and evolution.

Our understanding of bee phylogeny has improved over the past fifteen years as a result of new data, primarily nucleotide sequence data, and new methods, primarily model-based methods of phylogeny reconstruction. Phylogenetic studies based on single or, more commonly, multilocus data sets have helped resolve the placement of bees within the superfamily Apoidea; the relationships among the seven families of bees; and the relationships among bee subfamilies, tribes, genera, and species. In addition, molecular phylogenies have played an important role in inferring evolutionary patterns and processes in bees. Phylogenies have provided the comparative framework for understanding the evolution of host-plant associations and pollen specialization, the evolution of social behavior, and the evolution of parasitism. In this paper, we present an overview of significant discoveries in bee phylogeny based primarily on the application of molecular data. We review the phylogenetic hypotheses family-by-family and then describe how the new phylogenetic insights have altered our understanding of bee biology.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae).

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

Symbiotic bacteria and fungi proliferate in diapause and may enhance overwintering survival in a solitary bee.

Host-microbe interactions underlie the development and fitness of many macroorganisms, including bees. Whereas many social bees benefit from vertically transmitted gut bacteria, current data suggests that solitary bees, which comprise the vast majority of species diversity within bees, lack a highly specialized gut microbiome. Here, we examine the composition and abundance of bacteria and fungi throughout the complete life cycle of the ground-nesting solitary bee Anthophora bomboides standfordiana. In contrast to expectations, immature bee stages maintain a distinct core microbiome consisting of Actinobacterial genera (Streptomyces, Nocardiodes) and the fungus Moniliella spathulata. Dormant (diapausing) larval bees hosted the most abundant and distinctive bacteria and fungi, attaining 33 and 52 times their initial copy number, respectively. We tested two adaptive hypotheses regarding microbial functions for diapausing bees. First, using isolated bacteria and fungi, we found that Streptomyces from brood cells inhibited the growth of multiple pathogenic filamentous fungi, suggesting a role in pathogen protection during overwintering, when bees face high pathogen pressure. Second, sugar alcohol composition changed in tandem with major changes in fungal abundance, suggesting links with bee cold tolerance or overwintering biology. We find that A. bomboides hosts a conserved core microbiome that may provide key fitness advantages through larval development and diapause, which raises the question of how this microbiome is maintained and faithfully transmitted between generations. Our results suggest that focus on microbiomes of mature or active insect developmental stages may overlook stage-specific symbionts and microbial fitness contributions during host dormancy.