A pile of pipelines: An overview of the bioinformatics software for metabarcoding data analyses.

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high-throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon-specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty-two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects.

Microbial diversity and antimicrobial resistance in faecal samples from acute medical patients assessed through metagenomic sequencing.

Antimicrobial resistance (AMR) is a threat to global public health. However, unsatisfactory approaches to directly measuring the AMR burden carried by individuals has hampered efforts to assess interventions aimed at reducing selection for AMR. Metagenomics can provide accurate detection and quantification of AMR genes within an individual person's faecal flora (their gut "resistome"). Using this approach, we aimed to test the hypothesis that differences in antimicrobial use across different hospitals in the United Kingdom will result in observable differences in the resistome of individual patients. Three National Health Service acute Hospital Trusts with markedly different antibiotic use and Clostridioides difficile infection rates collected faecal samples from anonymous patients which were discarded after C. difficile testing over a period of 9 to 15 months. Metagenomic DNA was extracted from these samples and sequenced using an Illumina NovaSeq 6000 platform. The resulting sequencing reads were analysed for taxonomic composition and for the presence of AMR genes. Among 683 faecal metagenomes we found huge variation between individuals in terms of taxonomic diversity (Shannon Index range 0.10-3.99) and carriage of AMR genes (Median 1.50 genes/cell/sample overall). We found no statistically significant differences in diversity (median Shannon index 2.16 (IQR 1.71-2.56), 2.15 (IQR 1.62-2.50) and 2.26 (IQR 1.55-2.51)) or carriage of AMR genes (median 1.37 genes/cell/sample (IQR 0.70-3.24), 1.70 (IQR 0.70-4.52) and 1.43 (IQR 0.55-3.71)) at the three trusts respectively. This was also the case across the sample collection period within the trusts. While we have not demonstrated differences over place or time using metagenomic sequencing of faecal discards, other sampling frameworks may be more suitable to determine whether organisational level differences in antibiotic use are associated with individual-level differences in burden of AMR carriage.

The interactions and hierarchical effects of long-term agricultural stressors on soil bacterial communities.

Soils are subjected to multiple anthropogenic modifications, but the synergistic impacts of simultaneous environmental stressors on below-ground communities are poorly understood. We used a large-scale (1152 plots), long-term (26 years), multi-factorial grassland experiment to assess the impact of five common agricultural practises (pesticides, herbicide, liming, fertilizers and grazing exclusion) and their interactive effects on the composition and activity of soil microbial communities. We confirmed that pH strongly impacts belowground communities, but further demonstrate that pH strongly mediates the impacts of other management factors. Notably, there was a significant interaction between liming and the effect of pesticide application, with only half of the taxa responding to pesticide being shared in both limed and unlimed treatments. Likewise, nutrient amendments significantly altered bacterial community structure in acidic soils. Not only do these results highlight an hierarchy of effect of commonly used agricultural practices but also the widespread interactions between treatments: many taxa were significantly affected by interactions between treatments, even in the absence of significant main effects. Furthermore, the results demonstrated that chemical amendments may not percolate deeply into physically unperturbed soils with effects concentrated between 0 and 30 cm, despite 20+ years of treatment. The research shows that future changes to agricultural practices will need to consider interactions among multiple factors.

Accumulation of nylon microplastics and polybrominated diphenyl ethers and effects on gut microbial community of Chironomus sancticaroli.

Microplastics (MP) are emerging contaminants with the capacity to bind and transport hydrophobic organic compounds of environmental concern, such as polybrominated diphenyl ethers (PBDEs). The aim of this study was to investigate the ingestion of nylon (polyamide) MP alone and when associated with PBDEs and their effects on Chironomus sancticaroli larvae survival and microbiome structure. Survival, PBDE uptake and microbial community composition were measured in fourth instar larvae exposed for 96 h to BDEs- 47, 99, 100 and 153 in the presence and absence of 1% w/w MP in sediment. Microbiome community structures were determined through high throughput sequencing of 16S small subunit ribosomal RNA gene (16S rRNA). Initial experiments showed that larvae ingested MP faster at 0.5% w/w MP, while depuration was more efficient at 1% w/w MP, although retention of MP was seen even after 168 h depuration. No mortality was observed as a result of PBDEs and MP exposure. MP had a negative effect on PBDE concentration within larvae (η2 = 0.94) and a negative effect on sediment concentrations (η2 = 0.48). In all samples, microbial communities were dominated by Alphaproteobacteria, Betaproteobacteria, Actinobacteria and Gammaproteobacteria. Bacterial alpha diversity was not significantly affected by PBDEs or MP exposure. However, the abundance of discrete bacterial taxa was more sensitive to MP (X2 = 45.81, p = 0.02), and PBDE exposure. Our results highlight that C. sancticaroli showed no acute response to MPs and PBDEs, but that MPs influenced bacterial microbiome structure even after only short-term acute exposure.

Equine grass sickness (a multiple systems neuropathy) is associated with alterations in the gastrointestinal mycobiome.

BACKGROUND: Equine grass sickness (EGS) is a multiple systems neuropathy of grazing horses of unknown aetiology. An apparently identical disease occurs in cats, dogs, rabbits, hares, sheep, alpacas and llamas. Many of the risk factors for EGS are consistent with it being a pasture mycotoxicosis. To identify potential causal fungi, the gastrointestinal mycobiota of EGS horses were evaluated using targeted amplicon sequencing, and compared with those of two control groups. Samples were collected post mortem from up to 5 sites in the gastrointestinal tracts of EGS horses (EGS group; 150 samples from 54 horses) and from control horses that were not grazing EGS pastures and that had been euthanased for reasons other than neurologic and gastrointestinal diseases (CTRL group; 67 samples from 31 horses). Faecal samples were also collected from healthy control horses that were co-grazing pastures with EGS horses at disease onset (CoG group; 48 samples from 48 horses). RESULTS: Mycobiota at all 5 gastrointestinal sites comprised large numbers of fungi exhibiting diverse taxonomy, growth morphology, trophic mode and ecological guild. FUNGuild analysis parsed most phylotypes as ingested environmental microfungi, agaricoids and yeasts, with only 1% as gastrointestinal adapted animal endosymbionts. Mycobiota richness varied throughout the gastrointestinal tract and was greater in EGS horses. There were significant inter-group and inter-site differences in mycobiota structure. A large number of phylotypes were differentially abundant among groups. Key phylotypes (n = 56) associated with EGS were identified that had high abundance and high prevalence in EGS samples, significantly increased abundance in EGS samples, and were important determinants of the inter-group differences in mycobiota structure. Many key phylotypes were extremophiles and/or were predicted to produce cytotoxic and/or neurotoxic extrolites. CONCLUSIONS: This is the first reported molecular characterisation of the gastrointestinal mycobiota of grazing horses. Key phylotypes associated with EGS were identified. Further work is required to determine whether neurotoxic extrolites from key phylotypes contribute to EGS aetiology or whether the association of key phylotypes and EGS is a consequence of disease or is non-causal.

Gut and faecal bacterial community of the terrestrial isopod Porcellionides pruinosus: potential use for monitoring exposure scenarios.

This work aimed to characterize the gut and faeces bacterial communities (BC) of Porcellionides pruinosus using high-throughput sequencing. Isopods were collected from the field and kept in laboratory conditions similar to those normally applied in ecotoxicology tests. Faeces and purged guts of isopods (n = 3 × 30) were analysed by pyrosequencing the V3-V4 region of 16 S rRNA encoding gene. Results showed that gut and faecal BCs were dominated by Proteobacteria, particularly by an OTU (Operational Taxonomic Unit) affiliated to genus Coxiella. Diversity and richness values were statistically higher for faecal BC, mainly due to the occurrence of several low-abundance phylotypes. These results may reflect faecal carriage of bacterial groups that cannot settle in the gut. BCs of P. pruinosus comprised: (1) common members of the soil microbiota, (2) bacterial symbionts, (3) bacteria related to host metabolic/ecological features, and (4) bacterial etiological agents. Comparison of BC of this isopod species with the BC from other invertebrates revealed common bacterial groups across taxa. The baseline information provided by this work will assist the design and data interpretation of future ecotoxicological or biomonitoring assays where the analysis of P. pruinosus BC should be included as an additional indicator. CAPSULE: Terrestrial isopods bacterial communities might support ecotoxicological assays and biomonitoring processes as a valuable tool.

Beyond Taxonomic Identification: Integration of Ecological Responses to a Soil Bacterial 16S rRNA Gene Database.

High-throughput sequencing 16S rRNA gene surveys have enabled new insights into the diversity of soil bacteria, and furthered understanding of the ecological drivers of abundances across landscapes. However, current analytical approaches are of limited use in formalizing syntheses of the ecological attributes of taxa discovered, because derived taxonomic units are typically unique to individual studies and sequence identification databases only characterize taxonomy. To address this, we used sequences obtained from a large nationwide soil survey (GB Countryside Survey, henceforth CS) to create a comprehensive soil specific 16S reference database, with coupled ecological information derived from survey metadata. Specifically, we modeled taxon responses to soil pH at the OTU level using hierarchical logistic regression (HOF) models, to provide information on both the shape of landscape scale pH-abundance responses, and pH optima (pH at which OTU abundance is maximal). We identify that most of the soil OTUs examined exhibited a non-flat relationship with soil pH. Further, the pH optima could not be generalized by broad taxonomy, highlighting the need for tools and databases synthesizing ecological traits at finer taxonomic resolution. We further demonstrate the utility of the database by testing against geographically dispersed query 16S datasets; evaluating efficacy by quantifying matches, and accuracy in predicting pH responses of query sequences from a separate large soil survey. We found that the CS database provided good coverage of dominant taxa; and that the taxa indicating soil pH in a query dataset corresponded with the pH classifications of top matches in the CS database. Furthermore we were able to predict query dataset community structure, using predicted abundances of dominant taxa based on query soil pH data and the HOF models of matched CS database taxa. The database with associated HOF model outputs is released as an online portal for querying single sequences of interest (https://shiny-apps.ceh.ac.uk/ID-TaxER/), and flat files are made available for use in bioinformatic pipelines. The further development of advanced informatics infrastructures incorporating modeled ecological attributes along with new functional genomic information will likely facilitate large scale exploration and prediction of soil microbial functional biodiversity under current and future environmental change scenarios.

Integration of DNA extraction, metabarcoding and an informatics pipeline to underpin a national citizen science honey monitoring scheme.

Worldwide honeybees (Apis mellifera L.) are one of the most widely kept domesticated animals, supporting domestic and commercial livelihoods through the production of honey and wax, as well as in the delivery of pollination services to crops. Quantifying which plant species are foraged upon by honeybees provides insights into their nutritional status as well as patterns of landscape scale habitat utilization. Here we outline a rapid and reproducible methodology for identifying environmental DNA (eDNA) originating principally from pollen grains suspended within honey. The process is based on a DNA extraction incorporating vacuum filtration prior to universal eukaryotic internal transcribed spacer 2 region (ITS2) amplicon generation, sequencing and identification. To provide a pre-cursor to sequence phylotyping, we outline systems for error-corrected processing amplicon sequence variant abundance tables that removes chimeras. This methodology underpins the new UK National Honey Monitoring Scheme.•We compare the efficacy and speed of centrifugation and filtration systems for removing pollen from honey samples as a precursor to plant DNA barcoding.•We introduce the 'HONEYPI' informatics pipeline, an open access resource implemented in python 2.7, to ensure long-term reproducibility during the process of amplicon sequence variant classification.

Bacterial communities in larger islands have reduced temporal turnover.

Patterns of species diversity provide fundamental insights into the underlying mechanisms and processes that regulate biodiversity. The species-time relationship (STR) has the potential to be one such pattern; in a comparable manner to its more extensively studied spatial analogue, the species-area relationship (SAR), which has been pivotal in the development of ecological models and theories. We sought to determine the mechanisms and processes that underpin STR patterns of temporal turnover by sampling bacterial communities within ten water-filled tree-holes on the same European beech tree through the course of a year. We took this natural model system to represent an archipelago of islands of varying sizes and with shared common immigration sources. We observed an inverse relationship between STR-derived turnover rates and island size. Further, turnover was related to island size and not island isolation within the study system as indicated by a low frequency of dispersal limitation and high homogenizing dispersal. Compared to SARs, STRs are understudied, as such, the findings from the current study should provide a renewed interest in STR-based patterns and processes.

Contrasting community assembly processes structure lotic bacteria metacommunities along the river continuum.

The heterogeneous nature of lotic habitats plays an important role in the complex ecological and evolutionary processes that structure the microbial communities within them. Due to such complexity, our understanding of lotic microbial ecology still lacks conceptual frameworks for the ecological processes that shape these communities. We explored how bacterial community composition and underlying ecological assembly processes differ between lotic habitats by examining community composition and inferring community assembly processes across four major habitat types (free-living, particle-associated, biofilm on benthic stones and rocks, and sediment). This was conducted at 12 river sites from headwater streams to the main river in the River Thames, UK. Our results indicate that there are distinct differences in the bacterial communities between four major habitat types, with contrasting ecological processes shaping their community assembly processes. While the mobile free-living and particle-associated communities were consistently less diverse than the fixed sediment and biofilm communities, the latter two communities displayed higher homogeneity across the sampling sites. This indicates that the relative influence of deterministic environmental filtering is elevated in sediment and biofilm communities compared with free-living and particle-associated communities, where stochastic processes play a larger role.

Genomic surveillance of Escherichia coli and Klebsiella spp. in hospital sink drains and patients.

Escherichia coli and Klebsiella spp. are important human pathogens that cause a wide spectrum of clinical disease. In healthcare settings, sinks and other wastewater sites have been shown to be reservoirs of antimicrobial-resistant E. coli and Klebsiella spp., particularly in the context of outbreaks of resistant strains amongst patients. Without focusing exclusively on resistance markers or a clinical outbreak, we demonstrate that many hospital sink drains are abundantly and persistently colonized with diverse populations of E. coli, Klebsiella pneumoniae and Klebsiella oxytoca, including both antimicrobial-resistant and susceptible strains. Using whole-genome sequencing of 439 isolates, we show that environmental bacterial populations are largely structured by ward and sink, with only a handful of lineages, such as E. coli ST635, being widely distributed, suggesting different prevailing ecologies, which may vary as a result of different inputs and selection pressures. Whole-genome sequencing of 46 contemporaneous patient isolates identified one (2 %; 95 % CI 0.05-11 %) E. coli urine infection-associated isolate with high similarity to a prior sink isolate, suggesting that sinks may contribute to up to 10 % of infections caused by these organisms in patients on the ward over the same timeframe. Using metagenomics from 20 sink-timepoints, we show that sinks also harbour many clinically relevant antimicrobial resistance genes including blaCTX-M, blaSHV and mcr, and may act as niches for the exchange and amplification of these genes. Our study reinforces the potential role of sinks in contributing to Enterobacterales infection and antimicrobial resistance in hospital patients, something that could be amenable to intervention. This article contains data hosted by Microreact.

Accumulation of polybrominated diphenyl ethers and microbiome response in the great pond snail Lymnaea stagnalis with exposure to nylon (polyamide) microplastics.

Microplastics attract widespread attention, including for their potential to transport toxic chemicals in the form of plasticisers and associated hydrophobic organic chemicals, such as polybrominated diphenyl ethers (PBDEs). The aims of this study were to investigate how nylon (polyamide) microplastics may affect PBDE accumulation in snails, and the acute effects of nylon particles and PBDEs on survival, weight change and inherent microbiome diversity and community composition of the pond snail Lymnaea stagnalis. Snails were exposed for 96 h to BDEs-47, 99, 100 and 153 in the presence and absence of 1% w/w nylon microplastics in quartz sand sediment. No mortality was observed over the exposure period. Snails not exposed to microplastics lost significantly more weight compared to those exposed to microplastics. Increasing PBDE concentration in the sediment resulted in an increased PBDE body burden in the snails, however microplastics did not significantly influence total PBDE uptake. Based on individual congeners, uptake of BDE 47 by snails was significantly reduced in the presence of microplastics. The diversity and composition of the snail microbiome was not significantly altered by the presence of PBDEs nor by the microplastics, singly or combined. Significant effects on a few individual operational taxonomic units (OTUs) occurred when comparing the highest PBDE concentration with the control treatment, but in the absence of microplastics only. Overall within these acute experiments, only subtle effects on weight loss and slight microbiome alterations occurred. These results therefore highlight that L. stagnalis are resilient to acute exposures to microplastics and PBDEs, and that microplastics are unlikely to influence HOC accumulation or the microbiome of this species over short timescales.

Author Correction: Environment and host as large-scale controls of ectomycorrhizal fungi.

Change history: In the HTML version of this Article, author 'Filipa Cox' had no affiliation in the author list, although she was correctly associated with affiliation 3 in the PDF. In addition, the blue circles for 'oak' were missing from Extended Data Fig. 1. These errors have been corrected online.

Environment and host as large-scale controls of ectomycorrhizal fungi.

Explaining the large-scale diversity of soil organisms that drive biogeochemical processes-and their responses to environmental change-is critical. However, identifying consistent drivers of belowground diversity and abundance for some soil organisms at large spatial scales remains problematic. Here we investigate a major guild, the ectomycorrhizal fungi, across European forests at a spatial scale and resolution that is-to our knowledge-unprecedented, to explore key biotic and abiotic predictors of ectomycorrhizal diversity and to identify dominant responses and thresholds for change across complex environmental gradients. We show the effect of 38 host, environment, climate and geographical variables on ectomycorrhizal diversity, and define thresholds of community change for key variables. We quantify host specificity and reveal plasticity in functional traits involved in soil foraging across gradients. We conclude that environmental and host factors explain most of the variation in ectomycorrhizal diversity, that the environmental thresholds used as major ecosystem assessment tools need adjustment and that the importance of belowground specificity and plasticity has previously been underappreciated.

Environmental Metabarcoding Reveals Contrasting Belowground and Aboveground Fungal Communities from Poplar at a Hg Phytomanagement Site.

Characterization of microbial communities in stressful conditions at a field level is rather scarce, especially when considering fungal communities from aboveground habitats. We aimed at characterizing fungal communities from different poplar habitats at a Hg-contaminated phytomanagement site by using Illumina-based sequencing, network analysis approach, and direct isolation of Hg-resistant fungal strains. The highest diversity estimated by the Shannon index was found for soil communities, which was negatively affected by soil Hg concentration. Among the significant correlations between soil operational taxonomic units (OTUs) in the co-occurrence network, 80% were negatively correlated revealing dominance of a pattern of mutual exclusion. The fungal communities associated with Populus roots mostly consisted of OTUs from the symbiotic guild, such as members of the Thelephoraceae, thus explaining the lowest diversity found for root communities. Additionally, root communities showed the highest network connectivity index, while rarely detected OTUs from the Glomeromycetes may have a central role in the root network. Unexpectedly high richness and diversity were found for aboveground habitats, compared to the root habitat. The aboveground habitats were dominated by yeasts from the Lalaria, Davidiella, and Bensingtonia genera, not detected in belowground habitats. Leaf and stem habitats were characterized by few dominant OTUs such as those from the Dothideomycete class producing mutual exclusion with other OTUs. Aureobasidium pullulans, one of the dominating OTUs, was further isolated from the leaf habitat, in addition to Nakazawaea populi species, which were found to be Hg resistant. Altogether, these findings will provide an improved point of reference for microbial research on inoculation-based programs of tailings dumps.

Helminth burden and ecological factors associated with alterations in wild host gastrointestinal microbiota.

Infection by gastrointestinal helminths of humans, livestock and wild animals is common, but the impact of such endoparasites on wild hosts and their gut microbiota represents an important overlooked component of population dynamics. Wild host gut microbiota and endoparasites occupy the same physical niche spaces with both affecting host nutrition and health. However, associations between the two are poorly understood. Here we used the commonly parasitized European shag (Phalacrocorax aristotelis) as a model wild host. Forty live adults from the same colony were sampled. Endoscopy was employed to quantify helminth infection in situ. Microbiota from the significantly distinct proventriculus (site of infection), cloacal and faecal gastrointestinal tract microbiomes were characterised using 16S rRNA gene-targeted high-throughput sequencing. We found increasingly strong associations between helminth infection and microbiota composition progressing away from the site of infection, observing a pronounced dysbiosis in microbiota when samples were partitioned into high- and low-burden groups. We posit this dysbiosis is predominately explained by helminths inducing an anti-inflammatory environment in the proventriculus, diverting host immune responses away from themselves. This study, within live wild animals, provides a vital foundation to better understand the mechanisms that underpin the three-way relationship between helminths, microbiota and hosts.

Assessment of the bimodality in the distribution of bacterial genome sizes.

Bacterial genome sizes have previously been shown to exhibit a bimodal distribution. This phenomenon has prompted discussion regarding the evolutionary forces driving genome size in bacteria and its ecological significance. We investigated the level of inherent redundancy in the public database and the effect it has on the shape of the apparent bimodal distribution. Our study reveals that there is a significant bias in the genome sequencing efforts towards a certain group of species, and that correcting the bias using species nomenclature and clustering of the 16S rRNA gene, results in a unimodal rather than the previously published bimodal distribution. The true genome size distribution and its wider ecological implications will soon emerge as we are currently witnessing rapid growth in the number of sequenced genomes from diverse environmental niches across a range of habitats at an unprecedented rate.

Soil pH effects on the interactions between dissolved zinc, non-nano- and nano-ZnO with soil bacterial communities.

Zinc oxide nanoparticles (ZnO NPs) are used in an array of products and processes, ranging from personal care products to antifouling paints, textiles, food additives, antibacterial agents and environmental remediation processes. Soils are an environment likely to be exposed to manmade nanoparticles due to the practice of applying sewage sludge as a fertiliser or as an organic soil improver. However, understanding on the interactions between soil properties, nanoparticles and the organisms that live within soil is lacking, especially with regards to soil bacterial communities. We studied the effects of nanoparticulate, non-nanoparticulate and ionic zinc (in the form of zinc chloride) on the composition of bacterial communities in soil with a modified pH range (from pH 4.5 to pH 7.2). We observed strong pH-dependent effects on the interaction between bacterial communities and all forms of zinc, with the largest changes in bacterial community composition occurring in soils with low and medium pH levels (pH 4.8 and 5.9). The high pH soil (pH 7.2) was less susceptible to the effects of zinc exposure. At the highest doses of zinc (2500 mg/kg dw soil), both nano and non-nano particulate zinc applications elicited a similar response in the soil bacterial community, and this differed significantly to the ionic zinc salt treatment. The results highlight the importance of considering soil pH in nanotoxicology studies, although further work is needed to determine the exact mechanisms controlling the toxicity and fate and interactions of nanoparticles with soil microbial communities.

Sequence and organization of the complete mitochondrial genome of the blackfly Simulium variegatum (Diptera: Simuliidae).

The complete mitochondrial genome of the European blackfly, Simulium variegatum Meigen, 1818 was sequenced using a combined Illumina and Sanger sequencing approach. Using the known sequence of Chironomus tepperi Skuse, 1889 (Chironomidae) homologous NGS reads were identified and assembled. The genome is 15,367 bp in length and includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and a control region. Gene order resembles that of the ancestral dipteran gene arrangement. The base composition of the genome is A (37.6%), T (35.3%), C (15.8%) and G (11.3%). The control region between 12S rRNA and tRNAIle is composed of 362 bp with no obvious repetitive motifs.

Rearing and foraging affects bumblebee (Bombus terrestris) gut microbiota.

Bumblebees are ecologically and economically important as pollinators of crop and wild plants, especially in temperate systems. Species, such as the buff-tailed bumblebee (Bombus terrestris), are reared commercially to pollinate high-value crops. Their highly specific gut microbiota, characterized by low diversity, may affect nutrition and immunity and are likely to be important for fitness and colony health. However, little is known about how environmental factors affect bacterial community structure. We analysed the gut microbiota from three groups of worker bumblebees (B. terrestris) from distinct colonies that varied in rearing and foraging characteristics: commercially reared with restricted foraging (RR); commercially reared with outside foraging (RF); and wild-caught workers (W). Contrary to previous studies, which indicate that bacterial communities are highly conserved across workers, we found that RF individuals had an intermediate community structure compared with RR and W types. Further, this was shaped by differences in the abundances of common operational taxonomic units (OTUs) and the diversity of rare OTUs present, which we propose results from an increase in the variety of carbohydrates obtained through foraging.