skip to main content


Title: Strain Structure and Dynamics Revealed by Targeted Deep Sequencing of the Honey Bee Gut Microbiome
ABSTRACT Host-associated microbiomes can be critical for the health and proper development of animals and plants. The answers to many fundamental questions regarding the modes of acquisition and microevolution of microbiome communities remain to be established. Deciphering strain-level dynamics is essential to fully understand how microbial communities evolve, but the forces shaping the strain-level dynamics of microbial communities remain largely unexplored, mostly because of methodological issues and cost. Here, we used targeted strain-level deep sequencing to uncover the strain dynamics within a host-associated microbial community using the honey bee gut microbiome as a model system. Our results revealed that amplicon sequencing of conserved protein-coding gene regions using species-specific primers is a cost-effective and accurate method for exploring strain-level diversity. In fact, using this method we were able to confirm strain-level results that have been obtained from whole-genome shotgun sequencing of the honey bee gut microbiome but with a much higher resolution. Importantly, our deep sequencing approach allowed us to explore the impact of low-frequency strains (i.e., cryptic strains) on microbiome dynamics. Results show that cryptic strain diversity is not responsible for the observed variations in microbiome composition across bees. Altogether, the findings revealed new fundamental insights regarding strain dynamics of host-associated microbiomes. IMPORTANCE The factors driving fine-scale composition and dynamics of gut microbial communities are poorly understood. In this study, we used metagenomic amplicon deep sequencing to decipher the strain dynamics of two key members of the honey bee gut microbiome. Using this high-throughput and cost-effective approach, we were able to confirm results from previous large-scale whole-genome shotgun (WGS) metagenomic sequencing studies while also gaining additional insights into the community dynamics of two core members of the honey bee gut microbiome. Moreover, we were able to show that cryptic strains are not responsible for the observed variations in microbiome composition across bees.  more » « less
Award ID(s):
1930776
NSF-PAR ID:
10250797
Author(s) / Creator(s):
; ;
Editor(s):
Campbell, Barbara J.
Date Published:
Journal Name:
mSphere
Volume:
5
Issue:
4
ISSN:
2379-5042
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    One of the best indicators of colony health for the European honey bee ( Apis mellifera ) is its performance in the production of honey. Recent research into the microbial communities naturally populating the bee gut raise the question as to whether there is a correlation between microbial community structure and colony productivity. In this work, we used 16S rRNA amplicon sequencing to explore the microbial composition associated with forager bees from honey bee colonies producing large amounts of surplus honey (productive) and compared them to colonies producing less (unproductive). As supported by previous work, the honey bee microbiome was found to be dominated by three major phyla: the Proteobacteria, Bacilli and Actinobacteria, within which we found a total of 23 different bacterial genera, including known “core” honey bee microbiome members. Using discriminant function analysis and correlation-based network analysis, we identified highly abundant members (such as Frischella and Gilliamella ) as important in shaping the bacterial community; libraries from colonies with high quantities of these Orbaceae members were also likely to contain fewer Bifidobacteria and Lactobacillus species (such as Firm-4). However, co-culture assays, using isolates from these major clades, were unable to confirm any antagonistic interaction between Gilliamella and honey bee gut bacteria. Our results suggest that honey bee colony productivity is associated with increased bacterial diversity, although this mechanism behind this correlation has yet to be determined. Our results also suggest researchers should not base inferences of bacterial interactions solely on correlations found using sequencing. Instead, we suggest that depth of sequencing and library size can dramatically influence statistically significant results from sequence analysis of amplicons and should be cautiously interpreted. 
    more » « less
  2. ABSTRACT Microbial communities are shaped by interactions among their constituent members. Some Gram-negative bacteria employ type VI secretion systems (T6SSs) to inject protein toxins into neighboring cells. These interactions have been theorized to affect the composition of host-associated microbiomes, but the role of T6SSs in the evolution of gut communities is not well understood. We report the discovery of two T6SSs and numerous T6SS-associated Rhs toxins within the gut bacteria of honey bees and bumble bees. We sequenced the genomes of 28 strains of Snodgrassella alvi , a characteristic bee gut microbe, and found tremendous variability in their Rhs toxin complements: altogether, these strains appear to encode hundreds of unique toxins. Some toxins are shared with Gilliamella apicola , a coresident gut symbiont, implicating horizontal gene transfer as a source of toxin diversity in the bee gut. We use data from a transposon mutagenesis screen to identify toxins with antibacterial function in the bee gut and validate the function and specificity of a subset of these toxin and immunity genes in Escherichia coli . Using transcriptome sequencing, we demonstrate that S. alvi T6SSs and associated toxins are upregulated in the gut environment. We find that S. alvi Rhs loci have a conserved architecture, consistent with the C-terminal displacement model of toxin diversification, with Rhs toxins, toxin fragments, and cognate immunity genes that are expressed and confer strong fitness effects in vivo . Our findings of T6SS activity and Rhs toxin diversity suggest that T6SS-mediated competition may be an important driver of coevolution within the bee gut microbiota. IMPORTANCE The structure and composition of host-associated bacterial communities are of broad interest, because these communities affect host health. Bees have a simple, conserved gut microbiota, which provides an opportunity to explore interactions between species that have coevolved within their host over millions of years. This study examined the role of type VI secretion systems (T6SSs)—protein complexes used to deliver toxic proteins into bacterial competitors—within the bee gut microbiota. We identified two T6SSs and diverse T6SS-associated toxins in bacterial strains from bees. Expression of these genes is increased in bacteria in the bee gut, and toxin and immunity genes demonstrate antibacterial and protective functions, respectively, when expressed in Escherichia coli . Our results suggest that coevolution among bacterial species in the bee gut has favored toxin diversification and maintenance of T6SS machinery, and demonstrate the importance of antagonistic interactions within host-associated microbial communities. 
    more » « less
  3. Pathogens and parasites of solitary bees have been studied for decades, but the microbiome as a whole is poorly understood for most taxa. Comparative analyses of microbiome features such as composition, abundance, and specificity, can shed light on bee ecology and the evolution of host–microbe interactions. Here we study microbiomes of ground-nesting cellophane bees (Colletidae: Diphaglossinae). From a microbial point of view, the diphaglossine genus Ptiloglossa is particularly remarkable: their larval provisions are liquid and smell consistently of fermentation. We sampled larval provisions and various life stages from wild nests of Ptiloglossa arizonensis and two species of closely related genera: Caupolicana yarrowi and Crawfordapis luctuosa . We also sampled nectar collected by P. arizonensis . Using 16S rRNA gene sequencing, we find that larval provisions of all three bee species are near-monocultures of lactobacilli. Nectar communities are more diverse, suggesting ecological filtering. Shotgun metagenomic and phylogenetic data indicate that Ptiloglossa culture multiple species and strains of Apilactobacillus , which circulate among bees and flowers. Larval lactobacilli disappear before pupation, and hence are likely not vertically transmitted, but rather reacquired from flowers as adults. Thus, brood cell microbiomes are qualitatively similar between diphaglossine bees and other solitary bees: lactobacilli-dominated, environmentally acquired, and non-species-specific. However, shotgun metagenomes provide evidence of a shift in bacterial abundance. As compared with several other bee species, Ptiloglossa have much higher ratios of bacterial to plant biomass in larval provisions, matching the unusually fermentative smell of their brood cells. Overall, Ptiloglossa illustrate a path by which hosts can evolve quantitatively novel symbioses: not by acquiring or domesticating novel symbionts, but by altering the microenvironment to favor growth of already widespread and generalist microbes. 
    more » « less
  4. We introduce Operational Genomic Unit (OGU), a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent from taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldomly applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in one synthetic and two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome datasets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project dataset, and more accurate prediction of human age by the gut microbiomes in the Finnish population. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate OGU adoption in future metagenomics studies. Importance Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. However, current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution compared to 16S rRNA amplicon sequence variant analysis. To solve these challenges, we introduce Operational Genomic Units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition while (ii) permitting use of phylogeny-aware tools. Our analysis of real-world datasets shows several advantages over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGU as standard practice in metagenomic studies. 
    more » « less
  5. Gambino, Michela (Ed.)
    ABSTRACT 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. 
    more » « less