Strain Structure and Dynamics Revealed by Targeted Deep Sequencing of the Honey Bee Gut Microbiome

Temporospatial dynamics and host specificity of honeybee gut bacteria

Honeybees are important pollinators worldwide, with their gut microbiota playing a crucial role in maintaining their health. The gut bacteria of honeybees consist of primarily five core lineages that are spread through social interactions. Previous studies have provided a basic understanding of the composition and function of the honeybee gut microbiota, with recent advancements focusing on analyzing diversity at the strain level and changes in bacterial functional genes. Research on honeybee gut microbiota across different regions globally has provided insights into microbial ecology. Additionally, recent findings have shed light on the mechanisms of host specificity of honeybee gut bacteria. This review explores the temporospatial dynamics in honeybee gut microbiota, discussing the reasons and mechanisms behind these fluctuations. This synopsis provides insights into host-microbe interactions and is invaluable for honeybee health.

Shotgun metagenomics reveals the gut microbial diversity and functions in Vespa mandarinia (Hymenoptera: Vespidae) at multiple life stages

Wasps play important roles as predators and pollinators in the ecosystem. The Jingpo minority residing in Yunnan Province, China, has a traditional practice of using wine infused with mature wasps as a customary remedy for managing rheumatoid arthritis. The larva of the wasp is also a tasteful folk dish that has created a tremendous market. There is a paucity of survival knowledge, which has greatly restricted their potential applications in food and healthcare. Recent research has highlighted the importance of gut microbiota in insect growth. Nevertheless, there is still a lack of understanding regarding the composition, changes, and functions of the gut microbiota in Vespa mandarinia during development. In this research, the gut microbiota were investigated across three growth stages of Vespa mandarinia using a metagenomic technology. The result revealed that there are significant variations in the proportion of main gut microbes during the metamorphosis of Vespa mandarinia. Tenericutes were found to dominate during the larval stage, while Proteobacteria emerged as the dominant group post-pupation. Through a comprehensive analysis of the gut microbiota metagenome, this study revealed functional differences in the wasp gut microbiota at various growth stages. During the larval stage, the gut microbiota plays a central role in promoting metabolism. Following pupation, the gut microbiota exhibited diversified functions, likely due to the complex environments and diverse food sources encountered after metamorphosis. These functions included amino acid metabolism, compound degradation, and defense mechanisms. This research provides an extensive dataset on the gut microbiota during the metamorphosis of Vespa mandarinia, contributing to a deeper understanding of the influence of gut microbiota on wasp growth. Furthermore, this study uncovers a unique microbial treasure within insect guts, which is important for advancing the application of wasps in the fields of food and medicine.

The honeybee microbiota and its impact on health and disease

Honeybees (Apis mellifera) are key pollinators that support global agriculture and are long-established models for developmental and behavioural research. Recently, they have emerged as models for studying gut microbial communities. Earlier research established that hindguts of adult worker bees harbour a conserved set of host-restricted bacterial species, each showing extensive strain variation. These bacteria can be cultured axenically and introduced to gnotobiotic hosts, and some have basic genetic tools available. In this Review, we explore the most recent research showing how the microbiota establishes itself in the gut and impacts bee biology and health. Microbiota members occupy specific niches within the gut where they interact with each other and the host. They engage in cross-feeding and antagonistic interactions, which likely contribute to the stability of the community and prevent pathogen invasion. An intact gut microbiota provides protection against diverse pathogens and parasites and contributes to the processing of refractory components of the pollen coat and dietary toxins. Absence or disruption of the microbiota results in altered expression of genes that underlie immunity, metabolism, behaviour and development. In the field, such disruption by agrochemicals may negatively impact bees. These findings demonstrate a key developmental and protective role of the microbiota, with broad implications for bee health.

Pollen Diet Diversity does not Affect Gut Bacterial Communities or Melanization in a Social and Solitary Bee Species

Pollinators face many stressors, including reduced floral diversity. A low-diversity diet can impair organisms' ability to cope with additional stressors, such as pathogens, by altering the gut microbiome and/or immune function, but these effects are understudied for most pollinators. We investigated the impact of pollen diet diversity on two ecologically and economically important generalist pollinators, the social bumble bee (Bombus impatiens) and the solitary alfalfa leafcutter bee (Megachile rotundata). We experimentally tested the effect of one-, two-, or three-species pollen diets on gut bacterial communities in both species, and the melanization immune response in B. impatiens. Pollen diets included dandelion (Taraxacum officinale), staghorn sumac (Rhus typhina), and hawthorn (Crataegus sp.) alone, each pair-wise combination, or a mix of all three species. We fed bees their diet for 7 days and then dissected out guts and sequenced 16S rRNA gene amplicons to characterize gut bacterial communities. To assess melanization in B. impatiens, we inserted microfilament implants into the bee abdomen and measured melanin deposition on the implant. We found that pollen diet did not influence gut bacterial communities in M. rotundata. In B. impatiens, pollen diet composition, but not diversity, affected gut bacterial richness in older, but not newly-emerged bees. Pollen diet did not affect the melanization response in B. impatiens. Our results suggest that even a monofloral, low-quality pollen diet such as dandelion can support diverse gut bacterial communities in captive-reared adults of these bee species. These findings shed light on the effects of reduced diet diversity on bee health.

Social transmission of bacterial symbionts homogenizes the microbiome within and across generations of group-living spiders

Disentangling modes and fidelity of symbiont transmission are key for understanding host-symbiont associations in wild populations. In group-living animals, social transmission may evolve to ensure high-fidelity transmission of symbionts, since non-reproducing helpers constitute a dead-end for vertical transmission. We investigated symbiont transmission in the social spider Stegodyphus dumicola, which lives in family groups where the majority of females are non-reproducing helpers, females feed offspring by regurgitation, and individuals feed communally on insect prey. Group members share temporally stable microbiomes across generations, while distinct variation in microbiome composition exists between groups. We hypothesized that horizontal transmission of symbionts is enhanced by social interactions, and investigated transmission routes within (horizontal) and across (vertical) generations using bacterial 16S rRNA gene amplicon sequencing in three experiments: (i) individuals were sampled at all life stages to assess at which life stage the microbiome is acquired. (ii) a cross-fostering design was employed to test whether offspring carry the microbiome from their natal nest, or acquire the microbiome of the foster nest via social transmission. (iii) adult spiders with different microbiome compositions were mixed to assess whether social transmission homogenizes microbiome composition among group members. We demonstrate that offspring hatch symbiont-free, and bacterial symbionts are transmitted vertically across generations by social interactions with the onset of regurgitation feeding by (foster)mothers in an early life stage. Social transmission governs horizontal inter-individual mixing and homogenization of microbiome composition among nest mates. We conclude that temporally stable host-symbiont associations in social species can be facilitated and maintained by high-fidelity social transmission.

A phylogenomic and comparative genomic analysis of Commensalibacter, a versatile insect symbiont

BACKGROUND

To understand mechanisms of adaptation and plasticity of pollinators and other insects a better understanding of diversity and function of their key symbionts is required. Commensalibacter is a genus of acetic acid bacterial symbionts in the gut of honey bees and other insect species, yet little information is available on the diversity and function of Commensalibacter bacteria. In the present study, whole-genome sequences of 12 Commensalibacter isolates from bumble bees, butterflies, Asian hornets and rowan berries were determined, and publicly available genome assemblies of 14 Commensalibacter strains were used in a phylogenomic and comparative genomic analysis.

RESULTS

The phylogenomic analysis revealed that the 26 Commensalibacter isolates represented four species, i.e. Commensalibacter intestini and three novel species for which we propose the names Commensalibacter melissae sp. nov., Commensalibacter communis sp. nov. and Commensalibacter papalotli sp. nov. Comparative genomic analysis revealed that the four Commensalibacter species had similar genetic pathways for central metabolism characterized by a complete tricarboxylic acid cycle and pentose phosphate pathway, but their genomes differed in size, G + C content, amino acid metabolism and carbohydrate-utilizing enzymes. The reduced genome size, the large number of species-specific gene clusters, and the small number of gene clusters shared between C. melissae and other Commensalibacter species suggested a unique evolutionary process in C. melissae, the Western honey bee symbiont.

CONCLUSION

The genus Commensalibacter is a widely distributed insect symbiont that consists of multiple species, each contributing in a species specific manner to the physiology of the holobiont host.

Fidelity varies in the symbiosis between a gutless marine worm and its microbial consortium

BACKGROUND

Many animals live in intimate associations with a species-rich microbiome. A key factor in maintaining these beneficial associations is fidelity, defined as the stability of associations between hosts and their microbiota over multiple host generations. Fidelity has been well studied in terrestrial hosts, particularly insects, over longer macroevolutionary time. In contrast, little is known about fidelity in marine animals with species-rich microbiomes at short microevolutionary time scales, that is at the level of a single host population. Given that natural selection acts most directly on local populations, studies of microevolutionary partner fidelity are important for revealing the ecological and evolutionary processes that drive intimate beneficial associations within animal species.

RESULTS

In this study on the obligate symbiosis between the gutless marine annelid Olavius algarvensis and its consortium of seven co-occurring bacterial symbionts, we show that partner fidelity varies across symbiont species from strict to absent over short microevolutionary time. Using a low-coverage sequencing approach that has not yet been applied to microbial community analyses, we analysed the metagenomes of 80 O. algarvensis individuals from the Mediterranean and compared host mitochondrial and symbiont phylogenies based on single-nucleotide polymorphisms across genomes. Fidelity was highest for the two chemoautotrophic, sulphur-oxidizing symbionts that dominated the microbial consortium of all O. algarvensis individuals. In contrast, fidelity was only intermediate to absent in the sulphate-reducing and spirochaetal symbionts with lower abundance. These differences in fidelity are likely driven by both selective and stochastic forces acting on the consistency with which symbionts are vertically transmitted.

CONCLUSIONS

We hypothesize that variable degrees of fidelity are advantageous for O. algarvensis by allowing the faithful transmission of their nutritionally most important symbionts and flexibility in the acquisition of other symbionts that promote ecological plasticity in the acquisition of environmental resources. Video Abstract.

Strain-level resolution and pneumococcal carriage dynamics by single-molecule real-time (SMRT) sequencing of the plyNCR marker: a longitudinal study in Swiss infants

BACKGROUND

Pneumococcal carriage has often been studied from a serotype perspective; however, little is known about the strain-specific carriage and inter-strain interactions. Here, we examined the strain-level carriage and co-colonization dynamics of Streptococcus pneumoniae in a Swiss birth cohort by PacBio single-molecule real-time (SMRT) sequencing of the plyNCR marker.

METHODS

A total of 872 nasal swab (NS) samples were included from 47 healthy infants during the first year of life. Pneumococcal carriage was determined based on the quantitative real-time polymerase chain reaction (qPCR) targeting the lytA gene. The plyNCR marker was amplified from 214 samples having lytA-based carriage for pneumococcal strain resolution. Amplicons were sequenced using SMRT technology, and sequences were analyzed with the DADA2 pipeline. In addition, pneumococcal serotypes were determined using conventional, multiplex PCR (cPCR).

RESULTS

PCR-based plyNCR amplification demonstrated a 94.2% sensitivity and 100% specificity for Streptococcus pneumoniae if compared to lytA qPCR. The overall carriage prevalence was 63.8%, and pneumococcal co-colonization (≥ 2 plyNCR amplicon sequence variants (ASVs)) was detected in 38/213 (17.8%) sequenced samples with the relative proportion of the least abundant strain(s) ranging from 1.1 to 48.8% (median, 17.2%; IQR, 5.8-33.4%). The median age to first acquisition was 147 days, and having ≥ 2 siblings increased the risk of acquisition.

CONCLUSION

The plyNCR amplicon sequencing is species-specific and enables pneumococcal strain resolution. We therefore recommend its application for longitudinal strain-level carriage studies of Streptococcus pneumoniae. Video Abstract.

Significant compositional and functional variation reveals the patterns of gut microbiota evolution among the widespread Asian honeybee populations

The gut microbiome is a crucial element that facilitates a host’s adaptation to a changing environment. Compared to the western honeybee Apis mellifera, the Asian honeybee, Apis cerana populations across its natural range remain mostly semi-feral and are less affected by bee management, which provides a good system to investigate how gut microbiota evolve under environmental heterogeneity on large geographic scales. We compared and analyzed the gut microbiomes of 99 Asian honeybees, from genetically diverged populations covering 13 provinces across China. Bacterial composition varied significantly across populations at phylotype, sequence-discrete population (SDP), and strain levels, but with extensive overlaps, indicating that the diversity of microbial community among A. cerana populations is driven by nestedness. Pollen diets were significantly correlated with both the composition and function of the gut microbiome. Core bacteria, Gilliamella and Lactobacillus Firm-5, showed antagonistic turnovers and contributed to the enrichment in carbohydrate transport and metabolism. By feeding and inoculation bioassays, we confirmed that the variations in pollen polysaccharide composition contributed to the trade-off of these core bacteria. Progressive change, i.e., nestedness, is the foundation of gut microbiome evolution among the Asian honeybee. Such a transition during the co-diversification of gut microbiomes is affected by environmental factors, diets in general, and pollen polysaccharides in particular.

The Role of Feeding Characteristics in Shaping Gut Microbiota Composition and Function of Ensifera (Orthoptera)

Feeding habits were the primary factor affecting the gut bacterial communities in Ensifera. However, the interaction mechanism between the gut microbiota and feeding characteristics is not precisely understood. Here, the gut microbiota of Ensifera with diverse feeding habits was analyzed by shotgun metagenomic sequencing to further clarify the composition and function of the gut microbiota and its relationship with feeding characteristics. Our results indicate that under the influence of feeding habits, the gut microbial communities of Ensifera showed specific characteristics. Firstly, the gut microbial communities of the Ensifera with different feeding habits differed significantly, among which the gut microbial diversity of the herbivorous Mecopoda niponensis was the highest. Secondly, the functional genes related to feeding habits were in high abundance. Thirdly, the specific function of the gut microbial species in the omnivorous Gryllotalpa orientalis showed that the more diverse the feeding behavior of Ensifera, the worse the functional specificity related to the feeding characteristics of its gut microbiota. However, feeding habits were not the only factors affecting the gut microbiota of Ensifera. Some microorganisms' genes, whose functions were unrelated to feeding characteristics but were relevant to energy acquisition and nutrient absorption, were detected in high abundance. Our results were the first to report on the composition and function of the gut microbiota of Ensifera based on shotgun metagenomic sequencing and to explore the potential mechanism of the gut microbiota's association with diverse feeding habits.

Chlorothalonil alters the gut microbiota and reduces the survival of immature honey bees reared in vitro

BACKGROUND

Chlorothalonil is a nonsystemic fungicide, and it is one of the most widely detected pesticides in bee hives. The effect of chlorothalonil on the survival, weight, and gut microbiota of immature Apis mellifera L. reared in vitro was studied.

RESULTS

Larvae were fed 1, 2, 4, 8, and 16 μg/mL chlorothalonil and compared with larvae fed the negative control (diet without any additives), positive control (45 mg/L dimethoate), and solvent control (2% acetone). Compared with the control groups, the survival of the 2, 4, 8, and 16 μg/mL chlorothalonil treatments was significantly reduced. The no-observed-adverse-effect concentration of chlorothalonil was 1 μg/mL. Chlorothalonil had no significant effect on larval weight. The gut bacterial community composition of newly emerged bees was determined by PacBio 16S rDNA gene sequencing. linear discriminant analysis effect size (LEFSe) analysis showed that Pseudomonadales and Burkholderiales were affected by exposure to chlorothalonil.

CONCLUSION

Chlorothalonil reduced the survival of honey bee larvae and altered the gut microbiota of newly emerged bees. The risk of pesticides to honey bees is related to their toxicity and exposure dose.

Amplicon Sequencing of Single-Copy Protein-Coding Genes Reveals Accurate Diversity for Sequence-Discrete Microbiome Populations

An in-depth understanding of microbial function and the division of ecological niches requires accurate delineation and identification of microbes at a fine taxonomic resolution. Microbial phylotypes are typically defined using a 97% small subunit (16S) rRNA threshold. However, increasing evidence has demonstrated the ubiquitous presence of taxonomic units of distinct functions within phylotypes. These so-called sequence-discrete populations (SDPs) have used to be mainly delineated by disjunct sequence similarity at the whole-genome level. However, gene markers that could accurately identify and quantify SDPs are lacking in microbial community studies. Here, we developed a pipeline to screen single-copy protein-coding genes that could accurately characterize SDP diversity via amplicon sequencing of microbial communities. Fifteen candidate marker genes were evaluated using three criteria (extent of sequence divergence, phylogenetic accuracy, and conservation of primer regions) and the selected genes were subject to test the efficiency in differentiating SDPs within Gilliamella, a core honeybee gut microbial phylotype, as a proof-of-concept. The results showed that the 16S V4 region failed to report accurate SDP diversities due to low taxonomic resolution and changing copy numbers. In contrast, the single-copy genes recommended by our pipeline were able to successfully quantify Gilliamella SDPs for both mock samples and honeybee guts, with results highly consistent with those of metagenomics. The pipeline developed in this study is expected to identify single-copy protein coding genes capable of accurately quantifying diverse bacterial communities at the SDP level. IMPORTANCE Microbial communities can be distinguished by discrete genetic and ecological characteristics. These sequence-discrete populations are foundational for investigating the composition and functional structures of microbial communities at high resolution. In this study, we screened for reliable single-copy protein-coding marker genes to identify sequence-discrete populations through our pipeline. Using marker gene amplicon sequencing, we could accurately and efficiently delineate the population diversity in microbial communities. These results suggest that single copy protein-coding genes can be an accurate, quantitative, and economical alternative for characterizing population diversity. Moreover, the feasibility of a gene as marker for any bacterial population identification can be quickly evaluated by the pipeline proposed here.

Global Composition of the Bacteriophage Community in Honey Bees

The microbial communities in animal digestive systems are critical for host development and health. They stimulate the immune system during development, synthesize important chemical compounds like hormones, aid in digestion, competitively exclude pathogens, etc. Compared to the bacterial and fungal components of the microbiome, we know little about the temporal and spatial dynamics of bacteriophage communities in animal digestive systems. Recently, the bacteriophages of the honey bee gut were characterized in two European bee populations. Most of the bacteriophages described in these two reports were novel, harbored many metabolic genes in their genomes, and had a community structure that suggests coevolution with their bacterial hosts. To describe the conservation of bacteriophages in bees and begin to understand their role in the bee microbiome, we sequenced the virome of Apis mellifera from Austin, TX, and compared bacteriophage compositions among three locations around the world. We found that most bacteriophages from Austin are novel, sharing no sequence similarity with anything in public repositories. However, many bacteriophages are shared among the three bee viromes, indicating specialization of bacteriophages in the bee gut. Our study, along with the two previous bee virome studies, shows that the bee gut bacteriophage community is simple compared to that of many animals, consisting of several hundred types of bacteriophages that primarily infect four of the dominant bacterial phylotypes in the bee gut. IMPORTANCE Viruses that infect bacteria (bacteriophages) are abundant in the microbial communities that live on and in plants and animals. However, our knowledge of the structure, dynamics, and function of these viral communities lags far behind our knowledge of their bacterial hosts. We sequenced the first bacteriophage community of honey bees from the United States and compared the U.S. honey bee bacteriophage community to those of samples from Europe. Our work is an important characterization of an economically critical insect species and shows how bacteriophage communities can contain highly conserved individuals and be highly variable in composition across a wide geographic range.

Antibiotics Alter the Expression of Genes Related to Behavioral Development in Honey Bees (Hymenoptera: Apidae)

Honey bees, as many species of social insects, display a division of labor among colony members based on behavioral specializations related to age. Adult worker honey bees perform a series of tasks in the hive when they are young (such as brood care or nursing) and at ca. 2-3 wk of age, shift to foraging for nectar and pollen outside the hive. The transition to foraging involves changes in metabolism and neuroendocrine activities. These changes are associated with a suite of developmental genes. It was recently demonstrated that antibiotics influence behavioral development by accelerating or delaying the onset of foraging depending on timing of antibiotic exposure. To understand the mechanisms of these changes, we conducted a study on the effects of antibiotics on expression of candidate genes known to regulate behavioral development. We demonstrate a delay in the typical changes in gene expression over the lifetime of the individuals that were exposed to antibiotics during immature stage and adulthood. Additionally, we show an acceleration in the typical changes in gene expression on individuals that were expose to antibiotics only during immature stage. These results show that timing of antibiotic exposure alter the typical regulation of behavioral development by metabolic and neuroendocrine processes.

Honey bee genetics shape the strain-level structure of gut microbiota in social transmission

BACKGROUND

Honey bee gut microbiota transmitted via social interactions are beneficial to the host health. Although the microbial community is relatively stable, individual variations and high strain-level diversity have been detected across honey bees. Although the bee gut microbiota structure is influenced by environmental factors, the heritability of the gut members and the contribution of the host genetics remains elusive. Considering bees within a colony are not readily genetically identical due to the polyandry of the queen, we hypothesize that the microbiota structure can be shaped by host genetics.

RESULTS

We used shotgun metagenomics to simultaneously profile the microbiota and host genotypes of bees from hives of four different subspecies. Gut composition is more distant between genetically different bees at both phylotype- and "sequence-discrete population" levels. We then performed a successive passaging experiment within colonies of hybrid bees generated by artificial insemination, which revealed that the microbial composition dramatically shifts across batches of bees during the social transmission. Specifically, different strains from the phylotype of Snodgrassella alvi are preferentially selected by genetically varied hosts, and strains from different hosts show a remarkably biased distribution of single-nucleotide polymorphism in the Type IV pili loci. Genome-wide association analysis identified that the relative abundance of a cluster of Bifidobacterium strains is associated with the host glutamate receptor gene specifically expressed in the bee brain. Finally, mono-colonization of Bifidobacterium with a specific polysaccharide utilization locus impacts the alternative splicing of the gluR-B gene, which is associated with an increased GABA level in the brain.

CONCLUSIONS

Our results indicated that host genetics influence the bee gut composition and suggest a gut-brain connection implicated in the gut bacterial strain preference. Honey bees have been used extensively as a model organism for social behaviors, genetics, and the gut microbiome. Further identification of host genetic function as a shaping force of microbial structure will advance our understanding of the host-microbe interactions. Video abstract.

High royal jelly production does not impact the gut microbiome of honey bees

Background

Honey bees are not only essential for pollination services, but are also economically important as a source of hive products (e.g., honey, royal jelly, pollen, wax, and propolis) that are used as foods, cosmetics, and alternative medicines. Royal jelly is a popular honey bee product with multiple potential medicinal properties. To boost royal jelly production, a long-term genetic selection program of Italian honey bees (ITBs) in China has been performed, resulting in honey bee stocks (here referred to as RJBs) that produce an order of magnitude more royal jelly than ITBs. Although multiple studies have investigated the molecular basis of increased royal jelly yields, one factor that has not been considered is the role of honey bee-associated gut microbes.

Results

Based on the behavioral, morphological, physiological, and neurological differences between RJBs and ITBs, we predicted that the gut microbiome composition of RJBs bees would differ from ITBs. To test this hypothesis, we investigated the bacterial composition of RJB and ITB workers from an urban location and RJBs from a rural location in China. Based on 16S rRNA gene profiling, we did not find any evidence that RJBs possess a unique bacterial gut community when compared to ITBs. However, we observed differences between honey bees from the urban versus rural sites.

Conclusions

Our results suggest that the environmental factors rather than stock differences are more important in shaping the bacterial composition in honey bee guts. Further studies are needed to investigate if the observed differences in relative abundance of taxa between the urban and rural bees correspond to distinct functional capabilities that impact honey bee health. Because the lifestyle, diet, and other environmental variables are different in rural and urban areas, controlled studies are needed to determine which of these factors are responsible for the observed differences in gut bacterial composition between urban and rural honeybees.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00124-1.

Field-Realistic Tylosin Exposure Impacts Honey Bee Microbiota and Pathogen Susceptibility, Which Is Ameliorated by Native Gut Probiotics

Antibiotics have been applied to honey bee (Apis mellifera) hives for decades to treat Paenibacillus larvae, which causes American foulbrood disease and kills honey bee larvae. One of the few antibiotics approved in apiculture is tylosin tartrate. This study examined how a realistic hive treatment regimen of tylosin affected the gut microbiota of bees and susceptibility to a bacterial pathogen. Tylosin treatment reduced bacterial species richness and phylogenetic diversity and reduced the absolute abundances and strain diversity of the beneficial core gut bacteria Snodgrassella alvi and Bifidobacterium spp. Bees from hives treated with tylosin died more quickly after being fed a bacterial pathogen (Serratia marcescens) in the laboratory. We then tested whether a probiotic cocktail of core bee gut species could bolster pathogen resistance. Probiotic exposure increased survival of bees from both control and tylosin-treated hives. Finally, we measured tylosin tolerance of core bee gut bacteria by plating cultured isolates on media with different tylosin concentrations. We observed highly variable responses, including large differences among strains of both S. alvi and Gilliamella spp. Thus, probiotic treatments using cultured bee gut bacteria may ameliorate harmful perturbations of the gut microbiota caused by antibiotics or other factors. IMPORTANCE The antibiotic tylosin tartrate is used to treat honey bee hives to control Paenibacillus larvae, the bacterium that causes American foulbrood. We found that bees from tylosin-treated hives had gut microbiomes with depleted overall diversity as well as reduced absolute abundances and strain diversity of the beneficial bee gut bacteria Snodgrassella alvi and Bifidobacterium spp. Furthermore, bees from treated hives suffered higher mortality when challenged with an opportunistic pathogen. Bees receiving a probiotic treatment, consisting of a cocktail of cultured isolates of native bee gut bacteria, had increased survival following pathogen challenge. Thus, probiotic treatment with native gut bacteria may ameliorate negative effects of antibiotic exposure.

Evolution of Interbacterial Antagonism in Bee Gut Microbiota Reflects Host and Symbiont Diversification

Antagonistic interactions between bacteria affect diversity and dynamics of host-associated communities, including gut communities that are linked to host health. In many bacterial communities, including human and honey bee gut microbiotas, antagonism is mediated by type VI secretion systems (T6SSs) that deliver lethal toxins to competing strains.

Gram-negative bacteria frequently possess type VI secretion systems (T6SSs), protein complexes that are able to inject toxic proteins into nearby cells. Many aspects of T6SS structure and function have been characterized for model species, but less is known about the evolutionary processes that shape T6SS and effector (toxin) diversity in host-associated microbial communities. The bee gut microbiota is a simple community that has codiversified with bees for >80 million years. This study investigated how complements of T6SSs and effectors within the bee microbiota changed as bacteria and their hosts diversified into isolated species. We used protein homology to survey 198 isolate genomes of 9 Gram-negative species for genes encoding T6SS structural components; Rhs toxins, which are common T6SS effectors; and VgrG proteins, which are structural components associated with specific toxins. T6SS loci were present in 5 species clusters found only in bees, namely Apibacter spp., Gilliamella spp., Frischella perrara, “Candidatus Schmidhempelia bombi,” and Snodgrassella alvi. The distribution of T6SS loci suggests that at least 3 were present in the microbiota of the common ancestor of social bees and that loss of these genes in some bacterial lineages was linked to both host and bacterial speciation. Isolates differed enormously in repertoires of Rhs and VgrG proteins. We found that bacterial species employ different mechanisms for toxin acquisition and diversification and that species and strains sometimes lose the T6SS entirely, likely causing shifts in competitive dynamics within these communities.

IMPORTANCE Antagonistic interactions between bacteria affect diversity and dynamics of host-associated communities, including gut communities that are linked to host health. In many bacterial communities, including human and honey bee gut microbiotas, antagonism is mediated by type VI secretion systems (T6SSs) that deliver lethal toxins to competing strains. In this study, we explored how T6SSs and associated toxins have evolved in the simple, host-specific gut microbiota of honey bees and bumble bees. Using comparative genomics, we explored the conservation, recombination, horizontal transfer, and loss of T6SSs and effectors during 80 million years of evolution of this bee-associated community. We find that that patterns of T6SS loss and retention are linked to differences in biology across host species, while trends in effector diversification are mostly specific to bacterial lineages.