The pollination services provided by bees are essential for supporting natural and agricultural ecosystems. However, bee population declines have been documented across the world. Many of the factors known to undermine bee health (e.g., poor nutrition) can decrease immunocompetence and, thereby, increase bees' susceptibility to diseases. Given the myriad of stressors that can exacerbate disease in wild bee populations, assessments of the relative impact of landscape habitat conditions on bee pathogen prevalence are needed to effectively conserve pollinator populations. Herein, we assess how landscape-level conditions, including various metrics of floral/nesting resources, insecticides, weather, and honey bee (Apis mellifera) abundance, drive variation in wild bumble bee (Bombus impatiens) pathogen loads. Specifically, we screened 890 bumble bee workers from varied habitats in Pennsylvania, USA for three pathogens (deformed wing virus, black queen cell virus, and Vairimorpha (= Nosema) bombi), Defensin expression, and body size. Bumble bees collected within low-quality landscapes exhibited the highest pathogen loads, with spring floral resources and nesting habitat availability serving as the main drivers. We also found higher loads of pathogens where honey bee apiaries are more abundant, a positive relationship between Vairimorpha loads and rainfall, and differences in pathogens by geographic region. Collectively, our results highlight the need to support high-quality landscapes (i.e., those with abundant floral/nesting resources) to maintain healthy wild bee populations.

Wild bees, like many other taxa, are threatened by land-use and climate change, which, in turn, jeopardizes pollination of crops and wild plants. Understanding how land-use and climate factors interact is critical to predicting and managing pollinator populations and ensuring adequate pollination services, but most studies have evaluated either land-use or climate effects, not both. Furthermore, bee species are incredibly variable, spanning an array of behavioral, physiological, and life-history traits that can increase or decrease resilience to land-use or climate change. Thus, there are likely bee species that benefit, while others suffer, from changing climate and land use, but few studies have documented taxon-specific trends. To address these critical knowledge gaps, we analyzed a long-term dataset of wild bee occurrences from Maryland, Delaware, and Washington DC, USA, examining how different bee genera and functional groups respond to landscape composition, quality, and climate factors. Despite a large body of literature documenting land-use effects on wild bees, in this study, climate factors emerged as the main drivers of wild-bee abundance and richness. For wild-bee communities in spring and summer/fall, temperature and precipitation were more important predictors than landscape composition, landscape quality, or topography. However, relationships varied substantially between wild-bee genera and functional groups. In the Northeast USA, past trends and future predictions show a changing climate with warmer winters, more intense precipitation in winter and spring, and longer growing seasons with higher maximum temperatures. In almost all of our analyses, these conditions were associated with lower abundance of wild bees. Wild-bee richness results were more mixed, including neutral and positive relationships with predicted temperature and precipitation patterns. Thus, in this region and undoubtedly more broadly, changing climate poses a significant threat to wild-bee communities.

Regulation of worker behavior by dominant queens or workers is a hallmark of insect societies, but the underlying molecular mechanisms and their evolutionary conservation are not well understood. Honey bee and bumble bee colonies consist of a single reproductive queen and facultatively sterile workers. The queens' influences on the workers are mediated largely via inhibition of juvenile hormone titers, which affect division of labor in honey bees and worker reproduction in bumble bees. Studies in honey bees identified a transcription factor, Krüppel-homolog 1 (Kr-h1), whose expression in worker brains is significantly downregulated in the presence of a queen or queen pheromone and higher in forager bees, making this gene an ideal candidate for examining the evolutionary conservation of socially regulated pathways in Hymenoptera.

Exposure to carbon dioxide (CO2) has pleiotropic effects in many insect species, ranging from eliciting rapid behavioral responses such as attraction, to dramatic physiological changes, including ovary activation. In bumble bees, CO2 narcosis causes queens to bypass diapause and initiate egg laying, but its mode of action is not well-understood. Here, we evaluated the effects of CO2 narcosis on the behavior, physiology and immune function of virgin bumble bee queens (Bombus impatiens). We tested the hypothesis that CO2 induces these changes by stimulating oxidative stress response pathways. We found that CO2 stimulates ovarian activation and egg production and suppresses lipid (but not glycogen) accumulation in virgin queens. Additionally, CO2 treated queens were more active (particularly in terms of flight) and performed, but did not receive, more aggressive behaviors compared to controls. Moreover, CO2 positively affected immune function in queens, reduced transcript levels of 5/6 antioxidant enzyme genes and had no effect on longevity. Thus, although CO2 treatment stimulated reproduction, we did not observe any evidence of a trade-off in queen health parameters, aside from a reduction in lipids. Overall CO2 narcosis does not appear to stimulate a typical stress response in virgin bumble bee queens. On the contrary, CO2 narcosis appears to stimulate changes that prepare queens to cope with the nutritional, metabolic and behavioral challenges associated with reproduction and colony-founding.

Populations of pollinators are in decline worldwide. These declines are best documented in honey bees and are due to a combination of stressors. In particular, pesticides have been linked to decreased longevity and performance in honey bees; however, the molecular and physiological pathways mediating sensitivity and resistance to pesticides are not well characterized. We explored the impact of coumaphos and fluvalinate, the two most abundant and frequently detected pesticides in the hive, on genome-wide gene expression patterns of honey bee workers. We found significant changes in 1118 transcripts, including genes involved in detoxification, behavioral maturation, immunity, and nutrition. Since behavioral maturation is regulated by juvenile hormone III (JH), we examined effects of these miticides on hormone titers; while JH titers were unaffected, titers of methyl farnesoate (MF), the precursor to JH, were decreased. We further explored the association between nutrition- and pesticide-regulated gene expression patterns and demonstrated that bees fed a pollen-based diet exhibit reduced sensitivity to a third pesticide, chlorpyrifos. Finally, we demonstrated that expression levels of several of the putative pesticide detoxification genes identified in our study and previous studies are also upregulated in response to pollen feeding, suggesting that these pesticides and components in pollen modulate similar molecular response pathways. Our results demonstrate that pesticide exposure can substantially impact expression of genes involved in several core physiological pathways in honey bee workers. Additionally, there is substantial overlap in responses to pesticides and pollen-containing diets at the transcriptional level, and subsequent analyses demonstrated that pollen-based diets reduce workers' pesticide sensitivity. Thus, providing honey bees and other pollinators with high quality nutrition may improve resistance to pesticides.

Variation in individual behavior within social groups can affect the fitness of the group as well as the individual, and can be caused by a combination of genetic and environmental factors. However, the molecular factors associated with individual variation in social behavior remain relatively unexplored. We used honey bees (Apis mellifera) as a model to examine differences in socially-regulated behavior among individual workers, and used transcriptional profiling to determine if specific gene expression patterns are associated with these individual differences. In honey bees, the reproductive queen produces a pheromonal signal that regulates many aspects of worker behavior and physiology and maintains colony organization.

Each year, millions of kilograms of insecticides are applied to crops in the US. While insecticide use supports food, fuel, and fiber production, it can also threaten non-target organisms, a concern underscored by mounting evidence of widespread decline of pollinator populations. Here, we integrate several public datasets to generate county-level annual estimates of total 'bee toxic load' (honey bee lethal doses) for insecticides applied in the US between 1997-2012, calculated separately for oral and contact toxicity. To explore the underlying components of the observed changes, we divide bee toxic load into extent (area treated) and intensity (application rate x potency). We show that while contact-based bee toxic load remained relatively steady, oral-based bee toxic load increased roughly 9-fold, with reductions in application rate outweighed by disproportionate increases in potency (toxicity/kg) and extent. This pattern varied markedly by region, with the greatest increase seen in Heartland (121-fold increase), likely driven by use of neonicotinoid seed treatments in corn and soybean. In this "potency paradox", farmland in the central US has become more hazardous to bees despite lower volumes of insecticides applied, raising concerns about insect conservation and highlighting the importance of integrative approaches to pesticide use monitoring.

Understanding how vectors alter the interactions between viruses and their hosts is a fundamental question in virology and disease ecology. In honey bees, transmission of deformed wing virus (DWV) by parasitic Varroa mites has been associated with elevated disease and host mortality, and Varroa transmission has been hypothesized to lead to increased viral titres or select for more virulent variants. Here, we mimicked Varroa transmission by serially passaging a mixed population of two DWV variants, A and B, by injection through in vitro reared honey bee pupae and tracking these viral populations through five passages. The DWV-A and DWV-B variant proportions shifted dynamically through passaging, with DWV-B outcompeting DWV-A after one passage, but levels of both variants becoming equivalent by Passage 5. Sequencing analysis revealed a dominant, recombinant DWV-B strain (DWV-A derived 5' IRES region with the rest of the genome DWV-B), with low nucleotide diversity that decreased through passaging. DWV-A populations had higher nucleotide diversity compared to DWV-B, but this also decreased through passaging. Selection signatures were found across functional regions of the DWV-A and DWV-B genomes, including amino acid mutations in the putative capsid protein region. Simulated vector transmission differentially impacted two closely related viral variants which could influence viral interactions with the host, demonstrating surprising plasticity in vector-host-viral dynamics.

With the development of inexpensive, high-throughput sequencing technologies, it has become feasible to examine questions related to population genetics and molecular evolution of non-model species in their ecological contexts on a genome-wide scale. Here, we employed a newly developed suite of integrated, web-based programs to examine population dynamics and signatures of selection across the genome using several well-established tests, including F ST, pN/pS, and McDonald-Kreitman. We applied these techniques to study populations of honey bees (Apis mellifera) in East Africa. In Kenya, there are several described A. mellifera subspecies, which are thought to be localized to distinct ecological regions.

Parent-of-origin methylation arises when the methylation patterns of a particular allele are dependent on the parent it was inherited from. Previous work in honey bees has shown evidence of parent-of-origin-specific expression, yet the mechanisms regulating such pattern remain unknown in honey bees. In mammals and plants, DNA methylation is known to regulate parent-of-origin effects such as genomic imprinting. Here, we utilize genotyping of reciprocal European and Africanized honey bee crosses to study genome-wide allele-specific methylation patterns in sterile and reproductive individuals. Our data confirm the presence of allele-specific methylation in honey bees in lineage-specific contexts but also importantly, though to a lesser degree, parent-of-origin contexts. We show that the majority of allele-specific methylation occurs due to lineage rather than parent-of-origin factors, regardless of the reproductive state. Interestingly, genes affected by allele-specific DNA methylation often exhibit both lineage and parent-of-origin effects, indicating that they are particularly labile in terms of DNA methylation patterns. Additionally, we re-analyzed our previous study on parent-of-origin-specific expression in honey bees and found little association with parent-of-origin-specific methylation. These results indicate strong genetic background effects on allelic DNA methylation and suggest that although parent-of-origin effects are manifested in both DNA methylation and gene expression, they are not directly associated with each other.

In social groups, dominant individuals may socially inhibit reproduction of subordinates using aggressive interactions or, in the case of highly eusocial insects, pheromonal communication. It has been hypothesized these two modes of reproductive inhibition utilize conserved pathways. Here, we use a comparative framework to investigate the chemical and genomic underpinnings of reproductive dominance in the primitively eusocial wasp Polistes metricus. Our goals were to first characterize transcriptomic and chemical correlates of reproductive dominance and second, to test whether dominance-associated mechanisms in paper wasps overlapped with aggression or pheromone-related gene expression patterns in other species. To explore whether conserved molecular pathways relate to dominance, we compared wasp transcriptomic data to previous studies of gene expression associated with pheromonal communication and queen-worker differences in honey bees, and aggressive behavior in bees, Drosophila, and mice.

Social insects, such as honey bees, use molecular, physiological and behavioral responses to combat pathogens and parasites. The honey bee genome contains all of the canonical insect immune response pathways, and several studies have demonstrated that pathogens can activate expression of immune effectors. Honey bees also use behavioral responses, termed social immunity, to collectively defend their hives from pathogens and parasites. These responses include hygienic behavior (where workers remove diseased brood) and allo-grooming (where workers remove ectoparasites from nestmates). We have previously demonstrated that immunostimulation causes changes in the cuticular hydrocarbon profiles of workers, which results in altered worker-worker social interactions. Thus, cuticular hydrocarbons may enable workers to identify sick nestmates, and adjust their behavior in response. Here, we test the specificity of behavioral, chemical and genomic responses to immunostimulation by challenging workers with a panel of different immune stimulants (saline, Sephadex beads and Gram-negative bacteria E. coli).

RNA interference has been frequently applied to modulate gene function in organisms where the production and maintenance of mutants is challenging, as in our model of study, the honey bee, Apis mellifera. A green fluorescent protein (GFP)-derived double-stranded RNA (dsRNA-GFP) is currently commonly used as control in honey bee RNAi experiments, since its gene does not exist in the A. mellifera genome. Although dsRNA-GFP is not expected to trigger RNAi responses in treated bees, undesirable effects on gene expression, pigmentation or developmental timing are often observed. Here, we performed three independent experiments using microarrays to examine the effect of dsRNA-GFP treatment (introduced by feeding) on global gene expression patterns in developing worker bees. Our data revealed that the expression of nearly 1,400 genes was altered in response to dsRNA-GFP, representing around 10% of known honey bee genes. Expression changes appear to be the result of both direct off-target effects and indirect downstream secondary effects; indeed, there were several instances of sequence similarity between putative siRNAs generated from the dsRNA-GFP construct and genes whose expression levels were altered. In general, the affected genes are involved in important developmental and metabolic processes associated with RNA processing and transport, hormone metabolism, immunity, response to external stimulus and to stress. These results suggest that multiple dsRNA controls should be employed in RNAi studies in honey bees. Furthermore, any RNAi studies involving these genes affected by dsRNA-GFP in our studies should use a different dsRNA control.

Foraging behavior in honey bees (Apis mellifera) is a complex phenotype that is regulated by physiological state and social signals. How these factors are integrated at the molecular level to modulate foraging behavior has not been well characterized. The transition of worker bees from nursing to foraging behaviors is mediated by large-scale changes in brain gene expression, which are influenced by pheromones produced by the queen and larvae. Larval pheromones can also stimulate foragers to leave the colony to collect pollen. However, the mechanisms underpinning this rapid behavioral plasticity in foragers that specialize in collecting pollen over nectar, and how larval pheromones impact these different behavioral states, remains to be determined. Here, we investigated the patterns of gene expression related to rapid behavioral plasticity and task allocation among honey bee foragers exposed to two larval pheromones, brood pheromone (BP) and (E)-beta-ocimene (EBO). We hypothesized that both pheromones would alter expression of genes in the brain related to foraging and would differentially impact brain gene expression depending on foraging specialization.

The ectoparasitic mite Varroa destructor has emerged as the primary pest of domestic honey bees (Apis mellifera). Here we present an initial survey of the V. destructor genome carried out to advance our understanding of Varroa biology and to identify new avenues for mite control. This sequence survey provides immediate resources for molecular and population-genetic analyses of Varroa-Apis interactions and defines the challenges ahead for a comprehensive Varroa genome project.

The viral ecology of bee communities is complex, where viruses are readily shared among co-foraging bee species. Additionally, in honey bees (Apis mellifera), many viruses are transmitted - and their impacts exacerbated - by the parasitic Varroa destructor mite. Thus far, the viruses found to be shared across bee species and transmitted by V. destructor mites are positive-sense single-stranded RNA viruses. Recently, a negative-sense RNA enveloped virus, Apis rhabdovirus-1 (ARV-1), was found in A. mellifera honey bees in Africa, Europe, and islands in the Pacific. Here, we describe the identification - using a metagenomics approach - of ARV-1 in two bee species (A. mellifera and Bombus impatiens) and in V. destructor mites from populations collected in the United States and Israel. We confirmed the presence of ARV-1 in pools of A. mellifera, B. impatiens, and V. destructor from Israeli and U.S. populations by RT-PCR and found that it can reach high titers in individual honey bees and mites (107-108 viral genomic copies per individual). To estimate the prevalence of ARV-1 in honey bee populations, we screened 104 honey bee colonies across Israel, with 21 testing ARV-1-positive. Tagged-primer-mediated RT-PCR analysis detected the presence of the positive-sense ARV-1 RNA in A. mellifera and V. destructor, indicating that ARV-1 replicates in both hosts. This is the first report of the presence of ARV-1 in B. impatiens and of the replication of a rhabdovirus in A. mellifera and V. destructor. Our data suggest that Varroa mites could act as an ARV-1 vector; however, the presence of ARV-1 in B. impatiens (which are not parasitized by Varroa) suggests that it may not require the mite for transmission and ARV-1 may be shared among co-foraging bee species. Given that ARV-1 is found in non-Apis bee species, and because "ARV" is used for the Adelaide River virus, we propose that this virus should be called bee rhabdovirus 1 and abbreviated BRV-1. These results greatly expand our understanding of the diversity of viruses that can infect bee communities, though further analysis is required to determine how infection with this virus impacts these different hosts.

The molecular mechanisms underlying the post-mating behavioral and physiological transitions undergone by females have not been explored in great detail. Honey bees represent an excellent model system in which to address these questions because they exhibit a range of "mating states," with two extremes (virgins and egg-laying, mated queens) that differ dramatically in their behavior, pheromone profiles, and physiology. We used an incompletely-mated mating-state to understand the molecular processes that underlie the transition from a virgin to a mated, egg-laying queen. We used same-aged virgins, queens that mated once but did not initiate egg-laying, and queens that mated once and initiated egg-laying.

Understanding the ecological and evolutionary processes that drive host-pathogen interactions is critical for combating epidemics and conserving species. The Varroa destructor mite and deformed wing virus (DWV) are two synergistic threats to Western honeybee (Apis mellifera) populations across the globe. Distinct honeybee populations have been found to self-sustain despite Varroa infestations, including colonies within the Arnot Forest outside Ithaca, NY, USA. We hypothesized that in these bee populations, DWV has been selected to produce an avirulent infection phenotype, allowing for the persistence of both host and disease-causing agents. To investigate this, we assessed the titre of viruses in bees from the Arnot Forest and managed apiaries, and assessed genomic variation and virulence differences between DWV isolates. Across groups, we found viral abundance was similar, but DWV genotypes were distinct. We also found that infections with isolates from the Arnot Forest resulted in higher survival and lower rates of symptomatic deformed wings, compared to analogous isolates from managed colonies, providing preliminary evidence to support the hypothesis of adaptive decreased viral virulence. Overall, this multi-level investigation of virus genotype and phenotype indicates that host ecological context can be a significant driver of viral evolution and host-pathogen interactions in honeybees.

Bee viral ecology is a fascinating emerging area of research: viruses exert a range of effects on their hosts, exacerbate impacts of other environmental stressors, and, importantly, are readily shared across multiple bee species in a community. However, our understanding of bee viral communities is limited, as it is primarily derived from studies of North American and European Apis mellifera populations. Here, we examined viruses in populations of A. mellifera and 11 other bee species from 9 countries, across 4 continents and Oceania. We developed a novel pipeline to rapidly and inexpensively screen for bee viruses. This pipeline includes purification of encapsulated RNA/DNA viruses, sequence-independent amplification, high throughput sequencing, integrated assembly of contigs, and filtering to identify contigs specifically corresponding to viral sequences. We identified sequences for (+)ssRNA, (-)ssRNA, dsRNA, and ssDNA viruses. Overall, we found 127 contigs corresponding to novel viruses (i.e. previously not observed in bees), with 27 represented by >0.1% of the reads in a given sample, and 7 contained an RdRp or replicase sequence which could be used for robust phylogenetic analysis. This study provides a sequence-independent pipeline for viral metagenomics analysis, and greatly expands our understanding of the diversity of viruses found in bee communities.

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.