Variation in gut microbial contribution of essential amino acids to host protein metabolism in a wild small mammal community.

Herbivory is a dominant feeding strategy among animals, yet herbivores are often protein limited. The gut microbiome is hypothesized to help maintain host protein balance by provisioning essential macromolecules, but this has never been tested in wild consumers. Using amino acid carbon (δ13 C) and nitrogen (δ15 N) isotope analysis, we estimated the proportional contributions of essential amino acids (AAESS ) synthesized by gut microbes to five co-occurring desert rodents representing herbivorous, omnivorous and insectivorous functional groups. We found that herbivorous rodents occupying lower trophic positions (Dipodomys spp.) routed a substantial proportion (~40%-50%) of their AAESS from gut microbes, while higher trophic level omnivores (Peromyscus spp.) and insectivores (Onychomys arenicola) obtained most of their AAESS (~58%) from plant-based energy channels but still received ~20% of their AAESS from gut microbes. These findings empirically demonstrate that gut microbes play a key functional role in host protein metabolism in wild animals.

Limits to the three domains of life: lessons from community assembly along an Antarctic salinity gradient.

Extremophiles exist among all three domains of life; however, physiological mechanisms for surviving harsh environmental conditions differ among Bacteria, Archaea and Eukarya. Consequently, we expect that domain-specific variation of diversity and community assembly patterns exist along environmental gradients in extreme environments. We investigated inter-domain community compositional differences along a high-elevation salinity gradient in the McMurdo Dry Valleys, Antarctica. Conductivity for 24 soil samples collected along the gradient ranged widely from 50 to 8355 µS cm-1. Taxonomic richness varied among domains, with a total of 359 bacterial, 2 archaeal, 56 fungal, and 69 non-fungal eukaryotic operational taxonomic units (OTUs). Richness for bacteria, archaea, fungi, and non-fungal eukaryotes declined with increasing conductivity (all P < 0.05). Principal coordinate ordination analysis (PCoA) revealed significant (ANOSIM R = 0.97) groupings of low/high salinity bacterial OTUs, while OTUs from other domains were not significantly clustered. Bacterial beta diversity was unimodally distributed along the gradient and had a nested structure driven by species losses, whereas in fungi and non-fungal eukaryotes beta diversity declined monotonically without strong evidence of nestedness. Thus, while increased salinity acts as a stressor in all domains, the mechanisms driving community assembly along the gradient differ substantially between the domains.

Isotopic and genetic methods reveal the role of the gut microbiome in mammalian host essential amino acid metabolism.

Intestinal microbiota perform many functions for their host, but among the most important is their role in metabolism, especially the conversion of recalcitrant biomass that the host is unable to digest into bioavailable compounds. Most studies have focused on the assistance gut microbiota provide in the metabolism of carbohydrates, however, their role in host amino acid metabolism is poorly understood. We conducted an experiment on Mus musculus using 16S rRNA gene sequencing and carbon isotope analysis of essential amino acids (AAESS) to quantify the community composition of gut microbiota and the contribution of carbohydrate carbon used by the gut microbiome to synthesize AAESS that are assimilated by mice to build skeletal muscle tissue. The relative abundances of Firmicutes and Bacteroidetes inversely varied as a function of dietary macromolecular content, with Firmicutes dominating when mice were fed low-protein diets that contained the highest proportions of simple carbohydrates (sucrose). Mixing models estimated that the microbial contribution of AAESS to mouse muscle varied from less than 5% (threonine, lysine, and phenylalanine) to approximately 60% (valine) across diet treatments, with the Firmicute-dominated microbiome associated with the greatest contribution. Our results show that intestinal microbes can provide a significant source of the AAESS their host uses to synthesize structural tissues. The role that gut microbiota play in the amino acid metabolism of animals that consume protein-deficient diets is likely a significant but under-recognized aspect of foraging ecology and physiology.

Patterns of Symbiodinium (Dinophyceae) diversity and assemblages among diverse hosts and the coral reef environment of Lizard Island, Australia.

Large-scale environmental disturbances may impact both partners in coral host-Symbiodinium systems. Elucidation of the assembly patterns in such complex and interdependent communities may enable better prediction of environmental impacts across coral reef ecosystems. In this study, we investigated how the community composition and diversity of dinoflagellate symbionts in the genus Symbiodinium were distributed among 12 host species from six taxonomic orders (Actinaria, Alcyonacea, Miliolida, Porifera, Rhizostoma, Scleractinia) and in the reef water and sediments at Lizard Island, Great Barrier Reef before the 3rd Global Coral Bleaching Event. 454 pyrosequencing of the ITS2 region of Symbiodinium yielded 83 operational taxonomic units (OTUs) at a 97% similarity cut-off. Approximately half of the Symbiodinium OTUs from reef water or sediments were also present in symbio. OTUs belonged to six clades (A-D, F-G), but community structure was uneven. The two most abundant OTUs (100% matches to types C1 and A3) comprised 91% of reads and OTU C1 was shared by all species. However, sequence-based analysis of these dominant OTUs revealed host species specificity, suggesting that genetic similarity cut-offs of Symbiodinium ITS2 data sets need careful evaluation. Of the less abundant OTUs, roughly half occurred at only one site or in one species and the background Symbiodinium communities were distinct between individual samples. We conclude that sampling multiple host taxa with differing life history traits will be critical to fully understand the symbiont diversity of a given system and to predict coral ecosystem responses to environmental change and disturbance considering the differential stress response of the taxa within.

Decoupled responses of soil bacteria and their invertebrate consumer to warming, but not freeze-thaw cycles, in the Antarctic Dry Valleys.

Altered temperature profiles resulting in increased warming and freeze-thaw cycle (FTC) frequency pose great ecological challenges to organisms in alpine and polar ecosystems. We performed a laboratory microcosm experiment to investigate how temperature variability affects soil bacterial cell numbers, and abundance and traits of soil microfauna (the microbivorous nematode Scottnema lindsayae) from McMurdo Dry Valleys, Antarctica. FTCs and constant freezing shifted nematode body size distribution towards large individuals, driven by higher mortality among smaller individuals. FTCs reduced both bacterial and nematode abundance, but bacterial cell numbers also declined under warming, demonstrating decoupled consumer-prey responses. We predict that higher occurrence of FTCs in cold ecosystems will select for large body size within soil microinvertebrates and overall reduce their abundance. In contrast, warm temperatures without FTCs could lead to divergent responses in soil bacteria and their microinvertebrate consumers, potentially affecting energy and nutrient transfer rates in soil food webs of cold ecosystems.

Niche specialization of bacteria in permanently ice-covered lakes of the McMurdo Dry Valleys, Antarctica.

Perennially ice-covered lakes in the McMurdo Dry Valleys, Antarctica, are chemically stratified with depth and have distinct biological gradients. Despite long-term research on these unique environments, data on the structure of the microbial communities in the water columns of these lakes are scarce. Here, we examined bacterial diversity in five ice-covered Antarctic lakes by 16S rRNA gene-based pyrosequencing. Distinct communities were present in each lake, reflecting the unique biogeochemical characteristics of these environments. Further, certain bacterial lineages were confined exclusively to specific depths within each lake. For example, candidate division WM88 occurred solely at a depth of 15 m in Lake Fryxell, whereas unknown lineages of Chlorobi were found only at a depth of 18 m in Lake Miers, and two distinct classes of Firmicutes inhabited East and West Lobe Bonney at depths of 30 m. Redundancy analysis revealed that community variation of bacterioplankton could be explained by the distinct conditions of each lake and depth; in particular, assemblages from layers beneath the chemocline had biogeochemical associations that differed from those in the upper layers. These patterns of community composition may represent bacterial adaptations to the extreme and unique biogeochemical gradients of ice-covered lakes in the McMurdo Dry Valleys.

Microbial community analysis of pH 4 thermal springs in Yellowstone National Park.

The pH of the majority of thermal springs in Yellowstone National Park (YNP) is from 1 to 3 and 6 to 10; relatively few springs (~5%) have a pH range of 4-5. We used 16S rRNA gene pyrosequencing to investigate microbial communities sampled from four pH 4 thermal springs collected from four regions of YNP that differed in their fluid temperature and geochemistry. Our results revealed that the composition of bacterial communities varied among the sites, despite sharing similar pH values. The taxonomic composition and metabolic functional potential of the site with the lowest temperature (55 °C), a thermal spring from the Seven Mile Hole (SMH) area, were further investigated using shotgun metagenome sequencing. The taxonomic classification, based on 372 Mbp of unassembled metagenomic reads, indicated that this community included a high proportion of Chloroflexi, Bacteroidetes, Proteobacteria, and Firmicutes. Functional comparison with other YNP thermal spring metagenomes indicated that the SMH metagenome was enriched in genes related to energy production and conversion, transcription, and carbohydrate transport. Analysis of genes involved in nitrogen metabolism revealed assimilatory and dissimilatory nitrate reduction pathways, whereas the majority of genes involved in sulfur metabolism were related to the reduction of sulfate to adenylylsulfate, sulfite, and H2S. Given that pH 4 thermal springs are relatively less common in YNP and thermal areas worldwide, they may harbor novel microbiota and the communities that inhabit them deserve further investigation.

Niche and metabolic principles explain patterns of diversity and distribution: theory and a case study with soil bacterial communities.

The causes of biodiversity patterns are controversial and elusive due to complex environmental variation, covarying changes in communities, and lack of baseline and null theories to differentiate straightforward causes from more complex mechanisms. To address these limitations, we developed general diversity theory integrating metabolic principles with niche-based community assembly. We evaluated this theory by investigating patterns in the diversity and distribution of soil bacteria taxa across four orders of magnitude variation in spatial scale on an Antarctic mountainside in low complexity, highly oligotrophic soils. Our theory predicts that lower temperatures should reduce taxon niche widths along environmental gradients due to decreasing growth rates, and the changing niche widths should lead to contrasting α- and ß-diversity patterns. In accord with the predictions, α-diversity, niche widths and occupancies decreased while ß-diversity increased with increasing elevation and decreasing temperature. The theory also successfully predicts a hump-shaped relationship between α-diversity and pH and a negative relationship between α-diversity and salinity. Thus, a few simple principles explained systematic microbial diversity variation along multiple gradients. Such general theory can be used to disentangle baseline effects from more complex effects of temperature and other variables on biodiversity patterns in a variety of ecosystems and organisms.

Detection and analysis of elusive members of a novel and diverse archaeal community within a thermal spring streamer consortium.

Recent metagenomic analyses of Yellowstone National Park (YNP) thermal spring communities suggested the presence of minor archaeal populations that simultaneous PCR-based assays using traditional 'universal' 16S rRNA gene primers failed to detect. Here we use metagenomics to identify PCR primers effective at detecting elusive members of the Archaea, assess their efficacy, and describe the diverse and novel archaeal community from a circum-neutral thermal spring from the Bechler region of YNP. We determined that a less commonly used PCR primer, Arch349F, captured more diversity in this spring than the widely used A21F primer. A search of the PCR primers against the RDP 16S rRNA gene database indicated that Arch349F also captured the largest percentage of Archaea, including 41 % more than A21F. Pyrosequencing using the Arch349F primer recovered all of the phylotypes present in the clone-based portion of the study and the metagenome of this spring in addition to several other populations of Archaea, some of which are phylogenetically novel. In contrast to the lack of amplification with traditional 16S rRNA gene primers, our comprehensive analyses suggested a diverse archaeal community in the Bechler spring, with implications for recently discovered groups such as the Geoarchaeota and other undescribed archaeal groups.

Soil microbial responses to increased moisture and organic resources along a salinity gradient in a polar desert.

Microbial communities in extreme environments often have low diversity and specialized physiologies suggesting a limited resistance to change. The McMurdo Dry Valleys (MDV) are a microbially dominated, extreme ecosystem currently undergoing climate change-induced disturbances, including the melting of massive buried ice, cutting through of permafrost by streams, and warming events. These processes are increasing moisture across the landscape, altering conditions for soil communities by mobilizing nutrients and salts and stimulating autotrophic carbon inputs to soils. The goal of this study was to determine the effects of resource addition (water/organic matter) on the composition and function of microbial communities in the MDV along a natural salinity gradient representing an additional gradient of stress in an already extreme environment. Soil respiration and the activity of carbon-acquiring extracellular enzymes increased significantly (P < 0.05) with the addition of resources at the low- and moderate-salinity sites but not the high-salinity site. The bacterial community composition was altered, with an increase in Proteobacteria and Firmicutes with water and organic matter additions at the low- and moderate-salinity sites and a near dominance of Firmicutes at the high-salinity site. Principal coordinate analyses of all samples using a phylogenetically informed distance matrix (UniFrac) demonstrated discrete clustering among sites (analysis of similarity [ANOSIM], P < 0.05 and R > 0.40) and among most treatments within sites. The results from this experimental work suggest that microbial communities in this environment will undergo rapid change in response to the altered resources resulting from climate change impacts occurring in this region.

Bacterial community structure along moisture gradients in the parafluvial sediments of two ephemeral desert streams.

Microorganisms inhabiting stream sediments mediate biogeochemical processes of importance to both aquatic and terrestrial ecosystems. In deserts, the lateral margins of ephemeral stream channels (parafluvial sediments) are dried and rewetted, creating periodically wet conditions that typically enhance microbial activity. However, the influence of water content on microbial community composition and diversity in desert stream sediments is unclear. We sampled stream margins along gradients of wet to dry sediments, measuring geochemistry and bacterial 16S rRNA gene composition, at streams in both a cold (McMurdo Dry Valleys, Antarctica) and hot (Chihuahuan Desert, New Mexico, USA) desert. Across the gradients, sediment water content spanned a wide range (1.6-37.9% w/w), and conductivity was highly variable (12.3-1,380 µS cm(-2)). Bacterial diversity (at 97% sequence similarity) was high and variable, but did not differ significantly between the hot and cold desert and was not correlated with sediment water content. Instead, conductivity was most strongly related to diversity. Water content was strongly related to bacterial 16S rRNA gene community composition, though samples were distributed in wet and dry clusters rather than as assemblages shifting along a gradient. Phylogenetic analyses showed that many taxa from wet sediments at the hot and cold desert site were related to, respectively, halotolerant Gammaproteobacteria, and one family within the Sphingobacteriales (Bacteroidetes), while dry sediments at both sites contained a high proportion of taxa related to the Acidobacteria. These results suggest that bacterial diversity and composition in desert stream sediments is more strongly affected by hydrology and conductivity than temperature.

Stoichiometry of soil enzyme activity at global scale.

Extracellular enzymes are the proximate agents of organic matter decomposition and measures of these activities can be used as indicators of microbial nutrient demand. We conducted a global-scale meta-analysis of the seven-most widely measured soil enzyme activities, using data from 40 ecosystems. The activities of beta-1,4-glucosidase, cellobiohydrolase, beta-1,4-N-acetylglucosaminidase and phosphatase g(-1) soil increased with organic matter concentration; leucine aminopeptidase, phenol oxidase and peroxidase activities showed no relationship. All activities were significantly related to soil pH. Specific activities, i.e. activity g(-1) soil organic matter, also varied in relation to soil pH for all enzymes. Relationships with mean annual temperature (MAT) and precipitation (MAP) were generally weak. For hydrolases, ratios of specific C, N and P acquisition activities converged on 1 : 1 : 1 but across ecosystems, the ratio of C : P acquisition was inversely related to MAP and MAT while the ratio of C : N acquisition increased with MAP. Oxidative activities were more variable than hydrolytic activities and increased with soil pH. Our analyses indicate that the enzymatic potential for hydrolyzing the labile components of soil organic matter is tied to substrate availability, soil pH and the stoichiometry of microbial nutrient demand. The enzymatic potential for oxidizing the recalcitrant fractions of soil organic material, which is a proximate control on soil organic matter accumulation, is most strongly related to soil pH. These trends provide insight into the biogeochemical processes that create global patterns in ecological stoichiometry and organic matter storage.

Volcanic calderas delineate biogeographic provinces among Yellowstone thermophiles.

It has been suggested that the distribution of microorganisms should be cosmopolitan because of their enormous capacity for dispersal. However, recent studies have revealed that geographically isolated microbial populations do exist. Geographic distance as a barrier to dispersal is most often invoked to explain these distributions. Here we show that unique and diverse sequences of the bacterial genus Sulfurihydrogenibium exist in Yellowstone thermal springs, indicating that these sites are geographically isolated. Although there was no correlation with geographic distance or the associated geochemistry of the springs, there was a strong historical signal. We found that the Yellowstone calderas, remnants of prehistoric volcanic eruptions, delineate biogeographical provinces for the Sulfurihydrogenibium within Yellowstone (chi(2): 9.7, P = 0.002). The pattern of distribution that we have detected suggests that major geological events in the past 2 million years explain more of the variation in sequence diversity in this system than do contemporary factors such as habitat or geographic distance. These findings highlight the importance of historical legacies in determining contemporary microbial distributions and suggest that the same factors that determine the biogeography of macroorganisms are also evident among bacteria.

Molecular characterization of the diversity and distribution of a thermal spring microbial community by using rRNA and metabolic genes.

The diversity and distribution of a bacterial community from Coffee Pots Hot Spring, a thermal spring in Yellowstone National Park with a temperature range of 39.3 to 74.1 degrees C and pH range of 5.75 to 6.91, were investigated by sequencing cloned PCR products and quantitative PCR (qPCR) of 16S rRNA and metabolic genes. The spring was inhabited by three Aquificae genera--Thermocrinis, Hydrogenobaculum, and Sulfurihydrogenibium--and members of the Alpha-, Beta-, and Gammaproteobacteria, Firmicutes, Acidobacteria, Deinococcus-Thermus, and candidate division OP5. The in situ chemical affinities were calculated for 41 potential metabolic reactions using measured environmental parameters and a range of hydrogen and oxygen concentrations. Reactions that use oxygen, ferric iron, sulfur, and nitrate as electron acceptors were predicted to be the most energetically favorable, while reactions using sulfate were expected to be less favorable. Samples were screened for genes used in ammonia oxidation (amoA, bacterial gene only), the reductive tricarboxylic acid (rTCA) cycle (aclB), the Calvin cycle (cbbM), sulfate reduction (dsrAB), nitrogen fixation (nifH), nitrite reduction (nirK), and sulfide oxidation (soxEF1) by PCR. Genes for carbon fixation by the rTCA cycle and nitrogen fixation were detected. All aclB sequences were phylogenetically related and spatially correlated to Sulfurihydrogenibium 16S rRNA gene sequences using qPCR (R(2) = 0.99). This result supports the recent finding of citrate cleavage by enzymes other than ATP citrate lyase in the rTCA cycle of the Aquificaceae family. We briefly consider potential biochemical mechanisms that may allow Sulfurihydrogenibium and Thermocrinis to codominate some hydrothermal environments.

Local and Regional Scale Heterogeneity Drive Bacterial Community Diversity and Composition in a Polar Desert.

The distribution of organisms in an environment is neither uniform nor random but is instead spatially patterned. The factors that control this patterning are complex and the underlying mechanisms are poorly understood. Soil microbes are critical to ecosystem function but exhibit highly complex distributions and community dynamics due in large part to the scale-dependent effects of environmental heterogeneity. To better understand the impact of environmental heterogeneity on the distribution of soil microbes, we sequenced the 16S rRNA gene from bacterial communities in the microbe-dominated polar desert ecosystem of the McMurdo Dry Valleys (MDV), Antarctica. Significant differences in key edaphic variables and alpha diversity were observed among the three lake basins of the Taylor Valley (Kruskal-Wallis; pH: χ2 = 68.89, P < 0.001, conductivity: χ2 = 35.03, P < 0.001, observed species: χ2 = 7.98, P = 0.019 and inverse Simpson: χ2 = 18.52, P < 0.001) and each basin supported distinctive microbial communities (ANOSIM R = 0.466, P = 0.001, random forest ratio of 14.1). However, relationships between community structure and edaphic characteristics were highly variable and contextual, ranging in magnitude and direction across regional, basin, and local scales. Correlations among edaphic factors (pH and soil conductivity) and the relative abundance of specific phyla were most pronounced along local environmental gradients in the Lake Fryxell basin where Acidobacteria, Bacteroidetes, and Proteobacteria declined while Deinococcus-Thermus and Gemmatimonadetes increased with soil conductivity (all P < 0.1). Species richness was most strongly related to the soil conductivity gradient present within this study system. We suggest that the relative importance of pH versus soil conductivity in structuring microbial communities is related to the length of edaphic gradients and the spatial scale of sampling. These results highlight the importance of conducting studies over large ranges of key environmental gradients and across multiple spatial scales to assess the influence of environmental heterogeneity on the composition and diversity of microbial communities.

Leaf endophytic fungus interacts with precipitation to alter belowground microbial communities in primary successional dunes.

Understanding interactions between above- and belowground components of ecosystems is an important next step in community ecology. These interactions may be fundamental to predicting ecological responses to global change because indirect effects occurring through altered species interactions can outweigh or interact with the direct effects of environmental drivers. In a multiyear field experiment (2010-2015), we tested how experimental addition of a mutualistic leaf endophyte (Epichloë amarillans) associated with American beachgrass (Ammophila breviligulata) interacted with an altered precipitation regime (±30%) to affect the belowground microbial community. Epichloë addition increased host root biomass at the plot scale, but reduced the length of extraradical arbuscular mycorrhizal (AM) fungal hyphae in the soil. Under ambient precipitation alone, the addition of Epichloë increased root biomass per aboveground tiller and reduced the diversity of AM fungi in A. breviligulata roots. Furthermore, with Epichloë added, the diversity of root-associated bacteria declined with higher soil moisture, whereas in its absence, bacterial diversity increased with higher soil moisture. Thus, the aboveground fungal mutualist not only altered the abundance and composition of belowground microbial communities but also affected how belowground communities responded to climate, suggesting that aboveground microbes have potential for cascading influences on community dynamics and ecosystem processes that occur belowground.

Primary productivity as a control over soil microbial diversity along environmental gradients in a polar desert ecosystem.

Primary production is the fundamental source of energy to foodwebs and ecosystems, and is thus an important constraint on soil communities. This coupling is particularly evident in polar terrestrial ecosystems where biological diversity and activity is tightly constrained by edaphic gradients of productivity (e.g., soil moisture, organic carbon availability) and geochemical severity (e.g., pH, electrical conductivity). In the McMurdo Dry Valleys of Antarctica, environmental gradients determine numerous properties of soil communities and yet relatively few estimates of gross or net primary productivity (GPP, NPP) exist for this region. Here we describe a survey utilizing pulse amplitude modulation (PAM) fluorometry to estimate rates of GPP across a broad environmental gradient along with belowground microbial diversity and decomposition. PAM estimates of GPP ranged from an average of 0.27 µmol O2/m2/s in the most arid soils to an average of 6.97 µmol O2/m2/s in the most productive soils, the latter equivalent to 217 g C/m2/y in annual NPP assuming a 60 day growing season. A diversity index of four carbon-acquiring enzyme activities also increased with soil productivity, suggesting that the diversity of organic substrates in mesic environments may be an additional driver of microbial diversity. Overall, soil productivity was a stronger predictor of microbial diversity and enzymatic activity than any estimate of geochemical severity. These results highlight the fundamental role of environmental gradients to control community diversity and the dynamics of ecosystem-scale carbon pools in arid systems.

Decadal ecosystem response to an anomalous melt season in a polar desert in Antarctica.

Amplified climate change in polar regions is significantly altering regional ecosystems, yet there are few long-term records documenting these responses. The McMurdo Dry Valleys (MDV) cold desert ecosystem is the largest ice-free area of Antarctica, comprising soils, glaciers, meltwater streams and permanently ice-covered lakes. Multi-decadal records indicate that the MDV exhibited a distinct ecosystem response to an uncharacteristic austral summer and ensuing climatic shift. A decadal summer cooling phase ended in 2002 with intense glacial melt ('flood year')-a step-change in water availability triggering distinct changes in the ecosystem. Before 2002, the ecosystem exhibited synchronous behaviour: declining stream flow, decreasing lake levels, thickening lake ice cover, decreasing primary production in lakes and streams, and diminishing soil secondary production. Since 2002, summer air temperatures and solar flux have been relatively consistent, leading to lake level rise, lake ice thinning and elevated stream flow. Biological responses varied; one stream cyanobacterial mat type immediately increased production, but another stream mat type, soil invertebrates and lake primary productivity responded asynchronously a few years after 2002. This ecosystem response to a climatic anomaly demonstrates differential biological community responses to substantial perturbations, and the mediation of biological responses to climate change by changes in physical ecosystem properties.

Distinct Bacterial Microbiomes in Sexual and Asexual Potamopyrgus antipodarum, a New Zealand Freshwater Snail.

Different reproductive strategies and the transition to asexuality can be associated with microbial symbionts. Whether such a link exists within mollusks has never been evaluated. We took the first steps towards addressing this possibility by performing pyrosequencing of bacterial 16S rRNA genes associated with Potamopyrgus antipodarum, a New Zealand freshwater snail. A diverse set of 60 tissue collections from P. antipodarum that were genetically and geographically distinct and either obligately sexual or asexual were included, which allowed us to evaluate whether reproductive mode was associated with a particular bacterial community. 2624 unique operational taxonomic units (OTU, 97% DNA similarity) were detected, which were distributed across ~30 phyla. While alpha diversity metrics varied little among individual samples, significant differences in bacterial community composition and structure were detected between sexual and asexual snails, as well as among snails from different lakes and genetic backgrounds. The mean dissimilarity of the bacterial communities between the sexual and asexual P. antipodarum was 90%, largely driven by the presence of Rickettsiales in sexual snails and Rhodobacter in asexual snails. Our study suggests that there might be a link between reproductive mode and the bacterial microbiome of P. antipodarum, though a causal connection requires additional study.