Transmission mode and dispersal traits correlate with host specificity in mammalian gut microbes.

Different host species associate with distinct gut microbes in mammals, a pattern sometimes referred to as phylosymbiosis. However, the processes shaping this host specificity are not well understood. One model proposes that barriers to microbial transmission promote specificity by limiting microbial dispersal between hosts. This model predicts that specificity levels measured across microbes is correlated to transmission mode (vertical vs. horizontal) and individual dispersal traits. Here, we leverage two large publicly available gut microbiota data sets (1490 samples from 195 host species) to test this prediction. We found that host specificity varies widely across bacteria (i.e., there are generalist and specialist bacteria) and depends on transmission mode and dispersal ability. Horizontally-like transmitted bacteria equipped with traits that facilitate switches between host (e.g., tolerance to oxygen) were found to be less specific (more generalist) than microbes without those traits, for example, vertically-like inherited bacteria that are intolerant to oxygen. Altogether, our findings are compatible with a model in which limited microbial dispersal abilities foster host specificity.

Response to "Vast (but avoidable) underestimation of global biodiversity".

This Formal Comment provides clarifications on the authors' recent estimates of global bacterial diversity and the current status of the field, and responds to a Formal Comment from John Wiens regarding their prior work.

Kelp and sea urchin settlement mediated by biotic interactions with benthic coralline algal species.

Species interactions can influence key ecological processes that support community assembly and composition. For example, coralline algae encompass extensive diversity and may play a major role in regime shifts from kelp forests to urchin-dominated barrens through their role in inducing invertebrate larval metamorphosis and influencing kelp spore settlement. In a series of laboratory experiments, we tested the hypothesis that different coralline communities facilitate the maintenance of either ecosystem state by either promoting or inhibiting early recruitment of kelps or urchins. Coralline algae significantly increased red urchin metamorphosis compared with a control, while they had varying effects on kelp settlement. Urchin metamorphosis and density of juvenile canopy kelps did not differ significantly across coralline species abundant in both kelp forests and urchin barrens, suggesting that recruitment of urchin and canopy kelps does not depend on specific corallines. Non-calcified fleshy red algal crusts promoted the highest mean urchin metamorphosis percentage and showed some of the lowest canopy kelp settlement. In contrast, settlement of one subcanopy kelp species was reduced on crustose corallines, but elevated on articulated corallines, suggesting that articulated corallines, typically absent in urchin barrens, may need to recover before this subcanopy kelp could return. Coralline species differed in surface bacterial microbiome composition; however, urchin metamorphosis was not significantly different when microbiomes were removed with antibiotics. Our results clarify the role played by coralline algal species in kelp forest community assembly and could have important implications for kelp forest recovery.

Hierarchical phylogenetic community assembly of soil protists in a temperate agricultural field.

Protists are abundant, diverse and perform essential functions in soils. Protistan community structure and its change across time or space are traditionally studied at the species level but the relative importance of the processes shaping these patterns depends on the taxon phylogenetic resolution. Using 18S rDNA amplicon data of the Cercozoa, a group of dominant soil protists, from an agricultural field in western Germany, we observed a turnover of relatively closely related taxa (from sequence variants to genus-level clades) across soil depth; while across soil habitats (rhizosphere, bulk soil, drilosphere), we observed turnover of relatively distantly related taxa, confirming Paracercomonadidae as a rhizosphere-associated clade. We extended our approach to show that closely related Cercozoa encounter divergent arbuscular mycorrhizal (AM) fungi across soil depth and that distantly related Cercozoa encounter closely related AM fungi across soil compartments. This study suggests that soil Cercozoa community assembly at the field scale is driven by niche-based processes shaped by evolutionary legacy of adaptation to conditions primarily related to the soil compartment, followed by the soil layer, giving a deeper understanding on the selection pressures that shaped their evolution.

Soil protist function varies with elevation in the Swiss Alps.

Protists are abundant and play key trophic functions in soil. Documenting how their trophic contributions vary across large environmental gradients is essential to understand and predict how biogeochemical cycles will be impacted by global changes. Here, using amplicon sequencing of environmental DNA in open habitat soil from 161 locations spanning 2600 m of elevation in the Swiss Alps (from 400 to 3000 m), we found that, over the whole study area, soils are dominated by consumers, followed by parasites and phototrophs. In contrast, the proportion of these groups in local communities shows large variations in relation to elevation. While there is, on average, three times more consumers than parasites at low elevation (400-1000 m), this ratio increases to 12 at high elevation (2000-3000 m). This suggests that the decrease in protist host biomass and diversity toward mountains tops impact protist functional composition. Furthermore, the taxonomic composition of protists that infect animals was related to elevation while that of protists that infect plants or of protist consumers was related to soil pH. This study provides a first step to document and understand how soil protist functions vary along the elevational gradient.

A census-based estimate of Earth's bacterial and archaeal diversity.

The global diversity of Bacteria and Archaea, the most ancient and most widespread forms of life on Earth, is a subject of intense controversy. This controversy stems largely from the fact that existing estimates are entirely based on theoretical models or extrapolations from small and biased data sets. Here, in an attempt to census the bulk of Earth's bacterial and archaeal ("prokaryotic") clades and to estimate their overall global richness, we analyzed over 1.7 billion 16S ribosomal RNA amplicon sequences in the V4 hypervariable region obtained from 492 studies worldwide, covering a multitude of environments and using multiple alternative primers. From this data set, we recovered 739,880 prokaryotic operational taxonomic units (OTUs, 16S-V4 gene clusters at 97% similarity), a commonly used measure of microbial richness. Using several statistical approaches, we estimate that there exist globally about 0.8-1.6 million prokaryotic OTUs, of which we recovered somewhere between 47%-96%, representing >99.98% of prokaryotic cells. Consistent with this conclusion, our data set independently "recaptured" 91%-93% of 16S sequences from multiple previous global surveys, including PCR-independent metagenomic surveys. The distribution of relative OTU abundances is consistent with a log-normal model commonly observed in larger organisms; the total number of OTUs predicted by this model is also consistent with our global richness estimates. By combining our estimates with the ratio of full-length versus partial-length (V4) sequence diversity in the SILVA sequence database, we further estimate that there exist about 2.2-4.3 million full-length OTUs worldwide. When restricting our analysis to the Americas, while controlling for the number of studies, we obtain similar richness estimates as for the global data set, suggesting that most OTUs are globally distributed. Qualitatively similar results are also obtained for other 16S similarity thresholds (90%, 95%, and 99%). Our estimates constrain the extent of a poorly quantified rare microbial biosphere and refute recent predictions that there exist trillions of prokaryotic OTUs.

An integrated, modular approach to data science education in microbiology.

We live in an increasingly data-driven world, where high-throughput sequencing and mass spectrometry platforms are transforming biology into an information science. This has shifted major challenges in biological research from data generation and processing to interpretation and knowledge translation. However, postsecondary training in bioinformatics, or more generally data science for life scientists, lags behind current demand. In particular, development of accessible, undergraduate data science curricula has the potential to improve research and learning outcomes as well as better prepare students in the life sciences to thrive in public and private sector careers. Here, we describe the Experiential Data science for Undergraduate Cross-Disciplinary Education (EDUCE) initiative, which aims to progressively build data science competency across several years of integrated practice. Through EDUCE, students complete data science modules integrated into required and elective courses augmented with coordinated cocurricular activities. The EDUCE initiative draws on a community of practice consisting of teaching assistants (TAs), postdocs, instructors, and research faculty from multiple disciplines to overcome several reported barriers to data science for life scientists, including instructor capacity, student prior knowledge, and relevance to discipline-specific problems. Preliminary survey results indicate that even a single module improves student self-reported interest and/or experience in bioinformatics and computer science. Thus, EDUCE provides a flexible and extensible active learning framework for integration of data science curriculum into undergraduate courses and programs across the life sciences.

The microbiota of intertidal macroalgae Fucus distichus is site-specific and resistant to change following transplant.

It is unclear how host-associated microbial communities will be affected by future environmental change. Characterizing how microbiota differ across sites with varying environmental conditions and assessing the stability of the microbiota in response to abiotic variation are critical steps towards predicting outcomes of environmental change. Intertidal organisms are valuable study systems because they experience extreme variation in environmental conditions on tractable timescales such as tide cycles and across small spatial gradients in the intertidal zone. Here we show a widespread intertidal macroalgae, Fucus distichus, hosts site-specific microbiota over small (meters to kilometres) spatial scales. We demonstrate stability of site-specific microbial associations by manipulating the host environment and microbial species pool with common garden and reciprocal transplant experiments. We hypothesized that F. distichus microbiota would readily shift to reflect the contemporary environment due to selective filtering by abiotic conditions and/or colonization by microbes from the new environment or nearby hosts. Instead, F. distichus microbiota was stable for days after transplantation in both the laboratory and field. Our findings expand the current understanding of microbiota dynamics on an intertidal foundation species. These results may also point to adaptations for withstanding short-term environmental variation, in hosts and/or microbes, facilitating stable host-microbial associations.

Stunted childhood growth is associated with decompartmentalization of the gastrointestinal tract and overgrowth of oropharyngeal taxa.

Linear growth delay (stunting) affects roughly 155 million children under the age of 5 years worldwide. Treatment has been limited by a lack of understanding of the underlying pathophysiological mechanisms. Stunting is most likely associated with changes in the microbial community of the small intestine, a compartment vital for digestion and nutrient absorption. Efforts to better understand the pathophysiology have been hampered by difficulty of access to small intestinal fluids. Here, we describe the microbial community found in the upper gastrointestinal tract of stunted children aged 2-5 y living in sub-Saharan Africa. We studied 46 duodenal and 57 gastric samples from stunted children, as well as 404 fecal samples from stunted and nonstunted children living in Bangui, Central African Republic, and in Antananarivo, Madagascar, using 16S Illumina Amplicon sequencing and semiquantitative culture methods. The vast majority of the stunted children showed small intestinal bacterial overgrowth dominated by bacteria that normally reside in the oropharyngeal cavity. There was an overrepresentation of oral bacteria in fecal samples of stunted children, opening the way for developing noninvasive diagnostic markers. In addition, Escherichia coli/Shigella sp. and Campylobacter sp. were found to be more prevalent in stunted children, while Clostridia, well-known butyrate producers, were reduced. Our data suggest that stunting is associated with a microbiome "decompartmentalization" of the gastrointestinal tract characterized by an increased presence of oropharyngeal bacteria from the stomach to the colon, hence challenging the current view of stunting arising solely as a consequence of small intestine overstimulation through recurrent infections by enteric pathogens.

The relative importance of ecological drivers of arbuscular mycorrhizal fungal distribution varies with taxon phylogenetic resolution.

The phylogenetic depth at which arbuscular mycorrhizal (AM) fungi harbor a coherent ecological niche is unknown, which has consequences for operational taxonomic unit (OTU) delineation from sequence data and the study of their biogeography. We tested how changes in AM fungi community composition across habitats (beta diversity) vary with OTU phylogenetic resolution. We inferred exact sequence variants (ESVs) to resolve phylotypes at resolutions finer than provided by traditional sequence clustering and analyzed beta diversity profiles up to order-level sequence clusters. At the ESV level, we detected the environmental predictors revealed with traditional OTUs or at higher genetic distances. However, the correlation between environmental predictors and community turnover steeply increased at a genetic distance of c. 0.03 substitutions per site. Furthermore, we observed a turnover of either closely or distantly related taxa (respectively at or above 0.03 substitutions per site) along different environmental gradients. This study suggests that different axes of AM fungal ecological niche are conserved at different phylogenetic depths. Delineating AM fungal phylotypes using DNA sequences should screen different phylogenetic resolutions to better elucidate the factors that shape communities and predict the fate of AM symbioses in a changing environment.

Living the high life: Could gut microbiota matter for adaptation to high altitude?

At high altitude, the reduced availability of thermal energy and oxygen poses major challenges to organisms. Different species or populations have evolved similar solutions to these challenges, such as blood flow regulation in animals (Bouverot, 1985). Previous studies investigating such convergent adaptations have primarily looked at changes in host genomes (e.g., see Scheinfeldt & Tishkoff, 2010), but have rarely considered the potential role of the gut microbiome in mediating host adaptation. As gut microbes can indirectly regulate host blood pressure (Pluznick, 2014) and energy intake efficiency, it has been hypothesized that they could help maintain normal energy production and/or optimize nutritional assimilation in high-altitude hypoxic environments (e.g., Li & Zhao, 2015). However, it has been hard to (a) show that there is a direct effect of altitude on the gut microbiota, because of the many potential confounding effects of altitude (e.g., diet is correlated to altitude, as well as to the microbiome) and to (b) understand the mechanisms by which the microbiota could mediate host hypoxic and thermoregulatory stresses. In this issue of Molecular Ecology, Suzuki, Martins, and Nachman (2018) show that, independently of diet, taxonomic composition and functions of mouse gut microbiota converge in independent high-altitude environments and propose the intriguing hypothesis that some of these functional convergences might be beneficial to their host.

Mapping the imprint of biotic interactions on ß-diversity.

Investigating how trophic interactions influence the ß-diversity of meta-communities is of paramount importance to understanding the processes shaping biodiversity distribution. Here, we apply a statistical method for inferring the strength of spatial dependencies between pairs of species groups. Using simulated community data generated from a multi-trophic model, we showed that this method can approximate biotic interactions in multi-trophic communities based on ß-diversity patterns across groups. When applied to soil multi-trophic communities along an elevational gradient in the French Alps, we found that fungi make a major contribution to the structuring of ß-diversity across trophic groups. We also demonstrated that there were strong spatial dependencies between groups known to interact specifically (e.g. plant-symbiotic fungi, bacteria-nematodes) and that the influence of environment was less important than previously reported in the literature. Our method paves the way for a better understanding and mapping of multi-trophic communities through space and time.

Conserving Phylogenetic Diversity Can Be a Poor Strategy for Conserving Functional Diversity.

For decades, academic biologists have advocated for making conservation decisions in light of evolutionary history. Specifically, they suggest that policy makers should prioritize conserving phylogenetically diverse assemblages. The most prominent argument is that conserving phylogenetic diversity (PD) will also conserve diversity in traits and features (functional diversity [FD]), which may be valuable for a number of reasons. The claim that PD-maximized ("maxPD") sets of taxa will also have high FD is often taken at face value and in cases where researchers have actually tested it, they have done so by measuring the phylogenetic signal in ecologically important functional traits. The rationale is that if traits closely mirror phylogeny, then saving the maxPD set of taxa will tend to maximize FD and if traits do not have phylogenetic structure, then saving the maxPD set of taxa will be no better at capturing FD than criteria that ignore PD. Here, we suggest that measuring the phylogenetic signal in traits is uninformative for evaluating the effectiveness of using PD in conservation. We evolve traits under several different models and, for the first time, directly compare the FD of a set of taxa that maximize PD to the FD of a random set of the same size. Under many common models of trait evolution and tree shapes, conserving the maxPD set of taxa will conserve more FD than conserving a random set of the same size. However, this result cannot be generalized to other classes of models. We find that under biologically plausible scenarios, using PD to select species can actually lead to less FD compared with a random set. Critically, this can occur even when there is phylogenetic signal in the traits. Predicting exactly when we expect using PD to be a good strategy for conserving FD is challenging, as it depends on complex interactions between tree shape and the assumptions of the evolutionary model. Nonetheless, if our goal is to maintain trait diversity, the fact that conserving taxa based on PD will not reliably conserve at least as much FD as choosing randomly raises serious concerns about the general utility of PD in conservation.

Improving phylogenetic regression under complex evolutionary models.

Phylogenetic Generalized Least Square (PGLS) is the tool of choice among phylogenetic comparative methods to measure the correlation between species features such as morphological and life-history traits or niche characteristics. In its usual form, it assumes that the residual variation follows a homogenous model of evolution across the branches of the phylogenetic tree. Since a homogenous model of evolution is unlikely to be realistic in nature, we explored the robustness of the phylogenetic regression when this assumption is violated. We did so by simulating a set of traits under various heterogeneous models of evolution, and evaluating the statistical performance (type I error [the percentage of tests based on samples that incorrectly rejected a true null hypothesis] and power [the percentage of tests that correctly rejected a false null hypothesis]) of classical phylogenetic regression. We found that PGLS has good power but unacceptable type I error rates. This finding is important since this method has been increasingly used in comparative analyses over the last decade. To address this issue, we propose a simple solution based on transforming the underlying variance-covariance matrix to adjust for model heterogeneity within PGLS. We suggest that heterogeneous rates of evolution might be particularly prevalent in large phylogenetic trees, while most current approaches assume a homogenous rate of evolution. Our analysis demonstrates that overlooking rate heterogeneity can result in inflated type I errors, thus misleading comparative analyses. We show that it is possible to correct for this bias even when the underlying model of evolution is not known a priori.

Mammalian phylogenetic diversity-area relationships at a continental scale.

In analogy to the species-area relationship (SAR), one of the few laws in ecology, the phylogenetic diversity-area relationship (PDAR) describes the tendency of phylogenetic diversity (PD) to increase with area. Although investigating PDAR has the potential to unravel the underlying processes shaping assemblages across spatial scales and to predict PD loss through habitat reduction, it has been little investigated so far. Focusing on PD has noticeable advantages compared to species richness (SR), since PD also gives insights on processes such as speciation/extinction, assembly rules and ecosystem functioning. Here we investigate the universality and pervasiveness of the PDAR at continental scale using terrestrial mammals as study case. We define the relative robustness of PD (compared to SR) to habitat loss as the area between the standardized PDAR and standardized SAR (i.e., standardized by the diversity of the largest spatial window) divided by the area under the standardized SAR only. This metric quantifies the relative increase of PD robustness compared to SR robustness. We show that PD robustness is higher than SR robustness but that it varies among continents. We further use a null model approach to disentangle the relative effect of phylogenetic tree shape and nonrandom spatial distribution of evolutionary history on the PDAR. We find that, for most spatial scales and for all continents except Eurasia, PDARs are not different from expected by a model using only the observed SAR and the shape of the phylogenetic tree at continental scale. Interestingly, we detect a strong phylogenetic structure of the Eurasian PDAR that can be predicted by a model that specifically account for a finer biogeographical delineation of this continent. In conclusion, the relative robustness of PD to habitat loss compared to species richness is determined by the phylogenetic tree shape but also depends on the spatial structure of PD.

Multifaceted diversity-area relationships reveal global hotspots of mammalian species, trait and lineage diversity.

AIM: To define biome-scale hotspots of phylogenetic and functional mammalian biodiversity (PD and FD, respectively) and compare them to 'classical' hotspots based on species richness (SR) only. LOCATION: Global. METHODS: SR, PD & FD were computed for 782 terrestrial ecoregions using distribution ranges of 4616 mammalian species. We used a set of comprehensive diversity indices unified by a recent framework that incorporates the species relative coverage in each ecoregion. We build large-scale multifaceted diversity-area relationships to rank ecoregions according to their levels of biodiversity while accounting for the effect of area on each diversity facet. Finally we defined hotspots as the top-ranked ecoregions. RESULTS: While ignoring species relative coverage led to a relative good congruence between biome top ranked SR, PD and FD hotspots, ecoregions harboring a rich and abundantly represented evolutionary history and functional diversity did not match with top ranked ecoregions defined by species richness. More importantly PD and FD hotspots showed important spatial mismatches. We also found that FD and PD generally reached their maximum values faster than species richness as a function of area. MAIN CONCLUSIONS: The fact that PD/FD reach faster their maximal value than SR may suggest that the two former facets might be less vulnerable to habitat loss than the latter. While this point is expected, it is the first time that it is quantified at global scale and should have important consequences in conservation. Incorporating species relative coverage into the delineation of multifaceted hotspots of diversity lead to weak congruence between SR, PD and FD hotspots. This means that maximizing species number may fail at preserving those nodes (in the phylogenetic or functional tree) that are relatively abundant in the ecoregion. As a consequence it may be of prime importance to adopt a multifaceted biodiversity perspective to inform conservation strategies at global scale.

Contrasted host specificity of gut and endosymbiont bacterial communities in alpine grasshoppers and crickets.

Bacteria colonize the body of macroorganisms to form associations ranging from parasitic to mutualistic. Endosymbiont and gut symbiont communities are distinct microbiomes whose compositions are influenced by host ecology and evolution. Although the composition of horizontally acquired symbiont communities can correlate to host species identity (i.e. harbor host specificity) and host phylogeny (i.e. harbor phylosymbiosis), we hypothesize that the microbiota structure of vertically inherited symbionts (e.g. endosymbionts like Wolbachia) is more strongly associated with the host species identity and phylogeny than horizontally acquired symbionts (e.g. most gut symbionts). Here, using 16S metabarcoding on 336 guts from 24 orthopteran species (grasshoppers and crickets) in the Alps, we observed that microbiota correlated to host species identity, i.e. hosts from the same species had more similar microbiota than hosts from different species. This effect was ~5 times stronger for endosymbionts than for putative gut symbionts. Although elevation correlated with microbiome composition, we did not detect phylosymbiosis for endosymbionts and putative gut symbionts: closely related host species did not harbor more similar microbiota than distantly related species. Our findings indicate that gut microbiota of studied orthopteran species is more correlated to host identity and habitat than to the host phylogeny. The higher host specificity in endosymbionts corroborates the idea that-everything else being equal-vertically transmitted microbes harbor stronger host specificity signal, but the absence of phylosymbiosis suggests that host specificity changes quickly on evolutionary time scales.

Novel misos shape distinct microbial ecologies: opportunities for flavourful sustainable food innovation.

Fermentation is resurgent around the world as people seek healthier, more sustainable, and tasty food options. This study explores the microbial ecology of miso, a traditional Japanese fermented paste, made with novel regional substrates to develop new plant-based foods. Eight novel miso varieties were developed using different protein-rich substrates: yellow peas, Gotland lentils, and fava beans (each with two treatments: standard and nixtamalisation), as well as rye bread and soybeans. The misos were produced at Noma, a restaurant in Copenhagen, Denmark. Samples were analysed with biological and technical triplicates at the beginning and end of fermentation. We also incorporated in this study six samples of novel misos produced following the same recipe at Inua, a former affiliate restaurant of Noma in Tokyo, Japan. To analyse microbial community structure and diversity, metabarcoding (16S and ITS) and shotgun metagenomic analyses were performed. The misos contain a greater range of microbes than is currently described for miso in the literature. The composition of the novel yellow pea misos was notably similar to the traditional soybean ones, suggesting they are a good alternative, which supports our culinary collaborators' sensory conclusions. For bacteria, we found that overall substrate had the strongest effect, followed by time, treatment (nixtamalisation), and geography. For fungi, there was a slightly stronger effect of geography and a mild effect of substrate, and no significant effects for treatment or time. Based on an analysis of metagenome-assembled genomes (MAGs), strains of Staphylococccus epidermidis differentiated according to substrate. Carotenoid biosynthesis genes in these MAGs appeared in strains from Japan but not from Denmark, suggesting a possible gene-level geographical effect. The benign and possibly functional presence of S. epidermidis in these misos, a species typically associated with the human skin microbiome, suggests possible adaptation to the miso niche, and the flow of microbes between bodies and foods in certain fermentation as more common than is currently recognised. This study improves our understanding of miso ecology, highlights the potential for developing novel misos using diverse local ingredients, and suggests how fermentation innovation can contribute to studies of microbial ecology and evolution.

Deep Divergence and Genomic Diversification of Gut Symbionts of Neotropical Stingless Bees.

Social bees harbor conserved gut microbiotas that may have been acquired in a common ancestor of social bees and subsequently codiversified with their hosts. However, most of this knowledge is based on studies on the gut microbiotas of honey bees and bumblebees. Much less is known about the gut microbiotas of the third and most diverse group of social bees, the stingless bees. Specifically, the absence of genomic data from their microbiotas presents an important knowledge gap in understanding the evolution and functional diversity of the social bee microbiota. Here, we combined community profiling with culturing and genome sequencing of gut bacteria from six neotropical stingless bee species from Brazil. Phylogenomic analyses show that most stingless bee gut isolates form deep-branching sister clades of core members of the honey bee and bumblebee gut microbiota with conserved functional capabilities, confirming the common ancestry and ecology of their microbiota. However, our bacterial phylogenies were not congruent with those of the host, indicating that the evolution of the social bee gut microbiota was not driven by strict codiversification but included host switches and independent symbiont gain and losses. Finally, as reported for the honey bee and bumblebee microbiotas, we found substantial genomic divergence among strains of stingless bee gut bacteria, suggesting adaptation to different host species and glycan niches. Our study offers first insights into the genomic diversity of the stingless bee microbiota and highlights the need for broader samplings to understand the evolution of the social bee gut microbiota. IMPORTANCE Stingless bees are the most diverse group of the corbiculate bees and represent important pollinator species throughout the tropics and subtropics. They harbor specialized microbial communities in their gut that are related to those found in honey bees and bumblebees and that are likely important for bee health. Few bacteria have been cultured from the gut of stingless bees, which has prevented characterization of their genomic diversity and functional potential. Here, we established cultures of major members of the gut microbiotas of six stingless bee species and sequenced their genomes. We found that most stingless bee isolates belong to novel bacterial species distantly related to those found in honey bees and bumblebees and encoding similar functional capabilities. Our study offers a new perspective on the evolution of the social bee gut microbiota and presents a basis for characterizing the symbiotic relationships between gut bacteria and stingless bees.