It is presently unclear how much individual community members contribute to the overall metabolic output of a gut microbiota. To address this question, we used the honey bee, which harbors a relatively simple and remarkably conserved gut microbiota with striking parallels to the mammalian system and importance for bee health. Using untargeted metabolomics, we profiled metabolic changes in gnotobiotic bees that were colonized with the complete microbiota reconstituted from cultured strains. We then determined the contribution of individual community members in mono-colonized bees and recapitulated our findings using in vitro cultures. Our results show that the honey bee gut microbiota utilizes a wide range of pollen-derived substrates, including flavonoids and outer pollen wall components, suggesting a key role for degradation of recalcitrant secondary plant metabolites and pollen digestion. In turn, multiple species were responsible for the accumulation of organic acids and aromatic compound degradation intermediates. Moreover, a specific gut symbiont, Bifidobacterium asteroides, stimulated the production of host hormones known to impact bee development. While we found evidence for cross-feeding interactions, approximately 80% of the identified metabolic changes were also observed in mono-colonized bees, with Lactobacilli being responsible for the largest share of the metabolic output. These results show that, despite prolonged evolutionary associations, honey bee gut bacteria can independently establish and metabolize a wide range of compounds in the gut. Our study reveals diverse bacterial functions that are likely to contribute to bee health and provide fundamental insights into how metabolic activities are partitioned within gut communities. Honey bees are important pollinators that harbor a relatively simple gut microbiota with striking parallels to the mammalian system. This makes them relevant models to study gut microbiota functions and their impact on host health. We applied untargeted metabolomics to characterize metabolic changes induced by the gut microbiota and to characterize the contributions of the major community members. We find that the gut microbiota digests recalcitrant substrates derived from the pollen-based diet of bees. Most metabolic changes could be explained by the activity of individual community members, suggesting substrate specificity and independent metabolic functions. We did identify some cross-feeding interactions between species, including for pyruvate. Our study provides novel insights into the functional understanding of the bee gut microbiota and provides a framework for applying untargeted metabolomics to disentangle metabolic functions of gut bacteria.

The structure and distribution of genomic diversity in natural microbial communities is largely unexplored. Here, we used shotgun metagenomics to assess the diversity of the honey bee gut microbiota, a community consisting of few bacterial phylotypes. Our results show that most phylotypes are composed of sequence-discrete populations, which co-exist in individual bees and show age-specific abundance profiles. In contrast, strains present within these sequence-discrete populations were found to segregate into individual bees. Consequently, despite a conserved phylotype composition, each honey bee harbors a distinct community at the functional level. While ecological differentiation seems to facilitate coexistence at higher taxonomic levels, our findings suggest that, at the level of strains, priority effects during community assembly result in individualized profiles, despite the social lifestyle of the host. Our study underscores the need to move beyond phylotype-level characterizations to understand the function of this community, and illustrates its potential for strain-level analysis.

The importance of host-specialization to speciation processes in obligate host-associated bacteria is well known, as is also the ability of recombination to generate cohesion in bacterial populations. However, whether divergent strains of highly recombining intracellular bacteria, such as Wolbachia, can maintain their genetic distinctness when infecting the same host is not known. We first developed a protocol for the genome sequencing of uncultivable endosymbionts. Using this method, we have sequenced the complete genomes of the Wolbachia strains wHa and wNo, which occur as natural double infections in Drosophila simulans populations on the Seychelles and in New Caledonia. Taxonomically, wHa belong to supergroup A and wNo to supergroup B. A comparative genomics study including additional strains supported the supergroup classification scheme and revealed 24 and 33 group-specific genes, putatively involved in host-adaptation processes. Recombination frequencies were high for strains of the same supergroup despite different host-preference patterns, leading to genomic cohesion. The inferred recombination fragments for strains of different supergroups were of short sizes, and the genomes of the co-infecting Wolbachia strains wHa and wNo were not more similar to each other and did not share more genes than other A- and B-group strains that infect different hosts. We conclude that Wolbachia strains of supergroup A and B represent genetically distinct clades, and that strains of different supergroups can co-exist in the same arthropod host without converging into the same species. This suggests that the supergroups are irreversibly separated and that barriers other than host-specialization are able to maintain distinct clades in recombining endosymbiont populations. Acquiring a good knowledge of the barriers to genetic exchange in Wolbachia will advance our understanding of how endosymbiont communities are constructed from vertically and horizontally transmitted genes. Speciation in sexual organisms is defined as the inability of two populations to get viable offspring. Speciation in asexual, obligate endosymbionts is thought to be an indirect consequence of host-specialization. An important question is if divergent endosymbionts would start blending if the host barrier isolating them were removed. Here, we have studied Wolbachia, an abundant group of bacteria in the insect world. Wolbachia is classified into supergroups based on multi-locus sequence typing. We have sequenced the genomes from the Wolbachia strains wNo and wHa. These are particularly interesting since they belong to different supergroups yet co-occur as a double-infection in natural populations of Drosophila simulans. A comparative genomics study showed that wHa and wNo contain no uniquely shared genes. Instead, each strain shares unique gene functions with members of the same supergroup that infect other hosts. This unexpected finding suggests an alternative means of ecological speciation, indicating that speciation is not restricted to host-specialization but rather that related endosymbionts can coexist as separate species in the same host. Our study sheds light on the genomic divergence between different partners inhabiting the intracellular niche of the same host organism.

In the honeybee Apis mellifera, the bacterial gut community is consistently colonized by eight distinct phylotypes of bacteria. Managed bee colonies are of considerable economic interest and it is therefore important to elucidate the diversity and role of this microbiota in the honeybee. In this study, we have sequenced the genomes of eleven strains of lactobacilli and bifidobacteria isolated from the honey crop of the honeybee A. mellifera. Single gene phylogenies confirmed that the isolated strains represent the diversity of lactobacilli and bifidobacteria in the gut, as previously identified by 16S rRNA gene sequencing. Core genome phylogenies of the lactobacilli and bifidobacteria further indicated extensive divergence between strains classified as the same phylotype. Phylotype-specific protein families included unique surface proteins. Within phylotypes, we found a remarkably high level of gene content diversity. Carbohydrate metabolism and transport functions contributed up to 45% of the accessory genes, with some genomes having a higher content of genes encoding phosphotransferase systems for the uptake of carbohydrates than any previously sequenced genome. These genes were often located in highly variable genomic segments that also contained genes for enzymes involved in the degradation and modification of sugar residues. Strain-specific gene clusters for the biosynthesis of exopolysaccharides were identified in two phylotypes. The dynamics of these segments contrasted with low recombination frequencies and conserved gene order structures for the core genes. Hits for CRISPR spacers were almost exclusively found within phylotypes, suggesting that the phylotypes are associated with distinct phage populations. The honeybee gut microbiota has been described as consisting of a modest number of phylotypes; however, the genomes sequenced in the current study demonstrated a very high level of gene content diversity within all three described phylotypes of lactobacilli and bifidobacteria, particularly in terms of metabolic functions and surface structures, where many features were strain-specific. Together, these results indicate niche differentiation within phylotypes, suggesting that the honeybee gut microbiota is more complex than previously thought.

Interactions with microbes affect many aspects of animal biology, including immune system development, nutrition and health. In vertebrates, the gut microbiota is dominated by a small subset of phyla, but the species composition within these phyla is typically not conserved. Moreover, several recent studies have shown that bacterial species in the gut are composed of a multitude of strains, which frequently co-exist in their host, and may be host-specific. However, since the study of intra-species diversity is challenging, particularly in the setting of complex, host-associated microbial communities, our current understanding of the distribution, evolution and functional relevance of intra-species diversity in the gut is scarce. In order to unravel how genomic diversity translates into phenotypic diversity, community analyses going beyond 16S rRNA profiling, in combination with experimental approaches, are needed. Recently, the honeybee has emerged as a promising model for studying gut bacterial communities, particularly in terms of strain-level diversity. Unlike most other invertebrates, the honeybee gut is colonized by a remarkably consistent and specific core microbiota, which is dominated by only 8 bacterial species. As for the vertebrate gut microbiota, these species are composed of highly diverse strains suggesting that similar evolutionary forces shape gut community structures in vertebrates and social insects. In this review, we outline current knowledge on the evolution and functional relevance of strain diversity within the gut microbiota, including recent insights gained from mammals and other animals such as the honeybee. We discuss methodological approaches and propose possible future avenues for studying strain diversity in complex bacterial communities.

Most bacterial species encompass strains with vastly different gene content. Strain diversity in microbial communities is therefore considered to be of functional importance. Yet little is known about the extent to which related microbial communities differ in diversity at this level and which underlying mechanisms may constrain and maintain strain-level diversity. Here, we used shotgun metagenomics to characterize and compare the gut microbiota of two honey bee species, Apis mellifera and Apis cerana, which diverged about 6 mya. Although the host species are colonized largely by the same bacterial 16S rRNA phylotypes, we find that their communities are host specific when analyzed with genomic resolution. Moreover, despite their similar ecology, A. mellifera displayed a much higher diversity of strains and functional gene content in the microbiota compared to A. cerana, both per colony and per individual bee. In particular, the gene repertoire for polysaccharide degradation was massively expanded in the microbiota of A. mellifera relative to A. cerana. Bee management practices, divergent ecological adaptation, or habitat size may have contributed to the observed differences in microbiota genomic diversity of these key pollinator species. Our results illustrate that the gut microbiota of closely related animal hosts can differ vastly in genomic diversity while displaying similar levels of diversity based on the 16S rRNA gene. Such differences are likely to have consequences for gut microbiota functioning and host-symbiont interactions, highlighting the need for metagenomic studies to understand the ecology and evolution of microbial communities.