Sean M. Gibbons

Subsampled open-reference clustering creates consistent, comprehensive OTU definitions and scales to billions of sequences

We present a performance-optimized algorithm, subsampled open-reference OTU picking, for assigning marker gene (e.g., 16S rRNA) sequences generated on next-generation sequencing platforms to operational taxonomic units (OTUs) for microbial community analysis. This algorithm provides benefits over de novo OTU picking (clustering can be performed largely in parallel, reducing runtime) and closed-reference OTU picking (all reads are clustered, not only those that match a reference database sequence with high similarity). Because more of our algorithm can be run in parallel relative to “classic” open-reference OTU picking, it makes open-reference OTU picking tractable on massive amplicon sequence data sets (though on smaller data sets, “classic” open-reference OTU clustering is often faster). We illustrate that here by applying it to the first 15,000 samples sequenced for the Earth Microbiome Project (1.3 billion V4 16S rRNA amplicons). To the best of our knowledge, this is the largest OTU picking run ever performed, and we estimate that our new algorithm runs in less than 1/5 the time than would be required of “classic” open reference OTU picking. We show that subsampled open-reference OTU picking yields results that are highly correlated with those generated by “classic” open-reference OTU picking through comparisons on three well-studied datasets. An implementation of this algorithm is provided in the popular QIIME software package, which uses uclust for read clustering. All analyses were performed using QIIME’s uclust wrappers, though we provide details (aided by the open-source code in our GitHub repository) that will allow implementation of subsampled open-reference OTU picking independently of QIIME (e.g., in a compiled programming language, where runtimes should be further reduced). Our analyses should generalize to other implementations of these OTU picking algorithms. Finally, we present a comparison of parameter settings in QIIME’s OTU picking workflows and make recommendations on settings for these free parameters to optimize runtime without reducing the quality of the results. These optimized parameters can vastly decrease the runtime of uclust-based OTU picking in QIIME.

Effects of diurnal variation of gut microbes and high-fat feeding on host circadian clock function and metabolism.

Circadian clocks and metabolism are inextricably intertwined, where central and hepatic circadian clocks coordinate metabolic events in response to light-dark and sleep-wake cycles. We reveal an additional key element involved in maintaining host circadian rhythms, the gut microbiome. Despite persistence of light-dark signals, germ-free mice fed low or high-fat diets exhibit markedly impaired central and hepatic circadian clock gene expression and do not gain weight compared to conventionally raised counterparts. Examination of gut microbiota in conventionally raised mice showed differential diurnal variation in microbial structure and function dependent upon dietary composition. Additionally, specific microbial metabolites induced under low- or high-fat feeding, particularly short-chain fatty acids, but not hydrogen sulfide, directly modulate circadian clock gene expression within hepatocytes. These results underscore the ability of microbially derived metabolites to regulate or modify central and hepatic circadian rhythm and host metabolic function, the latter following intake of a Westernized diet.

Evidence for a persistent microbial seed bank throughout the global ocean

Do bacterial taxa demonstrate clear endemism, like macroorganisms, or can one site’s bacterial community recapture the total phylogenetic diversity of the world’s oceans? Here we compare a deep bacterial community characterization from one site in the English Channel (L4-DeepSeq) with 356 datasets from the International Census of Marine Microbes (ICoMM) taken from around the globe (ranging from marine pelagic and sediment samples to sponge-associated environments). At the L4-DeepSeq site, increasing sequencing depth uncovers greater phylogenetic overlap with the global ICoMM data. This site contained 31.7–66.2% of operational taxonomic units identified in a given ICoMM biome. Extrapolation of this overlap suggests that 1.93 × 1011 sequences from the L4 site would capture all ICoMM bacterial phylogenetic diversity. Current technology trends suggest this limit may be attainable within 3 y. These results strongly suggest the marine biosphere maintains a previously undetected, persistent microbial seed bank.

Severe plant invasions can increase mycorrhizal fungal abundance and diversity

Invasions by non-native plants can alter ecosystem functions and reduce native plant diversity, but relatively little is known about their effect on belowground microbial communities. We show that invasions by knapweed (Centaurea stoebe) and leafy spurge (Euphorbia esula, hereafter spurge)—but not cheatgrass (Bromus tectorum)—support a higher abundance and diversity of symbiotic arbuscular mycorrhizal fungi (AMF) than multi-species native plant communities. The higher AMF richness associated with knapweed and spurge is unlikely due to a co-invasion by AMF, because a separate sampling showed that individual native forbs hosted a similar AMF abundance and richness as exotic forbs. Native grasses associated with fewer AMF taxa, which could explain the reduced AMF richness in native, grass-dominated communities. The three invasive plant species harbored distinct AMF communities, and analyses of co-occurring native and invasive plants indicate that differences were partly driven by the invasive plants and were not the result of pre-invasion conditions. Our results suggest that invasions by mycotrophic plants that replace poorer hosts can increase AMF abundance and richness. The high AMF richness in monodominant plant invasions also indicates that the proposed positive relationship between above and belowground diversity is not always strong. Finally, the disparate responses among exotic plants and consistent results between grasses and forbs suggest that AMF respond more to plant functional group than plant provenance.

Whole-grain wheat consumption reduces inflammation in a randomized controlled trial on overweight and obese subjects with unhealthy dietary and lifestyle behaviors: role of polyphenols bound to cereal dietary fiber

BACKGROUND Epidemiology associates whole-grain (WG) consumption with several health benefits. Mounting evidence suggests that WG wheat polyphenols play a role in mechanisms underlying health benefits. OBJECTIVE The objective was to assess circulating concentration, excretion, and the physiologic role of WG wheat polyphenols in subjects with suboptimal dietary and lifestyle behaviors. DESIGN A placebo-controlled, parallel-group randomized trial with 80 healthy overweight/obese subjects with low intake of fruit and vegetables and sedentary lifestyle was performed. Participants replaced precise portions of refined wheat (RW) with a fixed amount of selected WG wheat or RW products for 8 wk. At baseline and every 4 wk, blood, urine, feces, and anthropometric and body composition measures were collected. Profiles of phenolic acids in biological samples, plasma markers of metabolic disease and inflammation, and fecal microbiota composition were assessed. RESULTS WG consumption for 4-8 wk determined a 4-fold increase in serum dihydroferulic acid (DHFA) and a 2-fold increase in fecal ferulic acid (FA) compared with RW consumption (no changes). Similarly, urinary FA at 8 wk doubled the baseline concentration only in WG subjects. Concomitant reduction in plasma tumor necrosis factor-α (TNF-α) after 8 wk and increased interleukin (IL)-10 only after 4 wk with WG compared with RW (P = 0.04) were observed. No significant change in plasma metabolic disease markers over the study period was observed, but a trend toward lower plasma plasminogen activator inhibitor 1 with higher excretion of FA and DHFA in the WG group was found. Fecal FA was associated with baseline low Bifidobacteriales and Bacteroidetes abundances, whereas after WG consumption, it correlated with increased Bacteroidetes and Firmicutes but reduced Clostridium. TNF-α reduction correlated with increased Bacteroides and Lactobacillus. No effect of dietary interventions on anthropometric measurements and body composition was found. CONCLUSIONS WG wheat consumption significantly increased excreted FA and circulating DHFA. Bacterial communities influenced fecal FA and were modified by WG wheat consumption. This trial was registered at clinicaltrials.gov as NCT01293175.

Investigating the Impact of Storage Conditions on Microbial Community Composition in Soil Samples

Recent advances in DNA sequencing technologies have allowed scientists to probe increasingly complex biological systems, including the diversity of bacteria in the environment. However, despite a multitude of recent studies incorporating these methods, many questions regarding how environmental samples should be collected and stored still persist. Here, we assess the impact of different soil storage conditions on microbial community composition using Illumina-based 16S rRNA V4 amplicon sequencing. Both storage time and temperature affected bacterial community composition and structure. Frozen samples maintained the highest alpha diversity and differed least in beta diversity, suggesting the utility of cold storage for maintaining consistent communities. Samples stored for intermediate times (three and seven days) had both the highest alpha diversity and the largest differences in overall beta diversity, showing the degree of community change after sample collection. These divergences notwithstanding, differences in neither storage time nor storage temperature substantially altered overall communities relative to more than 500 previously examined soil samples. These results systematically support previous studies and stress the importance of methodological consistency for accurate characterization and comparison of soil microbiological assemblages.

Human and environmental impacts on river sediment microbial communities.

Sediment microbial communities are responsible for a majority of the metabolic activity in river and stream ecosystems. Understanding the dynamics in community structure and function across freshwater environments will help us to predict how these ecosystems will change in response to human land-use practices. Here we present a spatiotemporal study of sediments in the Tongue River (Montana, USA), comprising six sites along 134 km of river sampled in both spring and fall for two years. Sequencing of 16S rRNA amplicons and shotgun metagenomes revealed that these sediments are the richest (∼65,000 microbial ‘species’ identified) and most novel (93% of OTUs do not match known microbial diversity) ecosystems analyzed by the Earth Microbiome Project to date, and display more functional diversity than was detected in a recent review of global soil metagenomes. Community structure and functional potential have been significantly altered by anthropogenic drivers, including increased pathogenicity and antibiotic metabolism markers near towns and metabolic signatures of coal and coalbed methane extraction byproducts. The core (OTUs shared across all samples) and the overall microbial community exhibited highly similar structure, and phylogeny was weakly coupled with functional potential. Together, these results suggest that microbial community structure is shaped by environmental drivers and niche filtering, though stochastic assembly processes likely play a role as well. These results indicate that sediment microbial communities are highly complex and sensitive to changes in land use practices.

Forensic analysis of the microbiome of phones and shoes

BackgroundMicrobial interaction between human-associated objects and the environments we inhabit may have forensic implications, and the extent to which microbes are shared between individuals inhabiting the same space may be relevant to human health and disease transmission. In this study, two participants sampled the front and back of their cell phones, four different locations on the soles of their shoes, and the floor beneath them every waking hour over a 2-day period. A further 89 participants took individual samples of their shoes and phones at three different scientific conferences.ResultsSamples taken from different surface types maintained significantly different microbial community structures. The impact of the floor microbial community on that of the shoe environments was strong and immediate, as evidenced by Procrustes analysis of shoe replicates and significant correlation between shoe and floor samples taken at the same time point. Supervised learning was highly effective at determining which participant had taken a given shoe or phone sample, and a Bayesian method was able to determine which participant had taken each shoe sample based entirely on its similarity to the floor samples. Both shoe and phone samples taken by conference participants clustered into distinct groups based on location, though much more so when an unweighted distance metric was used, suggesting sharing of low-abundance microbial taxa between individuals inhabiting the same space.ConclusionsCorrelations between microbial community sources and sinks allow for inference of the interactions between humans and their environment.

Ecological succession and viability of human-associated microbiota on restroom surfaces

ABSTRACT Human-associated bacteria dominate the built environment (BE). Following decontamination of floors, toilet seats, and soap dispensers in four public restrooms, in situ bacterial communities were characterized hourly, daily, and weekly to determine their successional ecology. The viability of cultivable bacteria, following the removal of dispersal agents (humans), was also assessed hourly. A late-successional community developed within 5 to 8 h on restroom floors and showed remarkable stability over weeks to months. Despite late-successional dominance by skin- and outdoor-associated bacteria, the most ubiquitous organisms were predominantly gut-associated taxa, which persisted following exclusion of humans. Staphylococcus represented the majority of the cultivable community, even after several hours of human exclusion. Methicillin-resistant Staphylococcus aureus (MRSA)-associated virulence genes were found on floors but were not present in assembled Staphylococcus pan-genomes. Viral abundances, which were predominantly enterophages, human papilloma virus, and herpesviruses, were significantly correlated with bacterial abundances and showed an unexpectedly low virus-to-bacterium ratio in surface-associated samples, suggesting that bacterial hosts are mostly dormant on BE surfaces.

Spatial scale drives patterns in soil bacterial diversity

Soil microbial communities are essential for ecosystem function, but linking community composition to biogeochemical processes is challenging because of high microbial diversity and large spatial variability of most soil characteristics. We investigated soil bacterial community structure in a switchgrass stand planted on soil with a history of grassland vegetation at high spatial resolution to determine whether biogeographic trends occurred at the centimeter scale. Moreover, we tested whether such heterogeneity, if present, influenced community structure within or among ecosystems. Pronounced heterogeneity was observed at centimeter scales, with abrupt changes in relative abundance of phyla from sample to sample. At the ecosystem scale (> 10 m), however, bacterial community composition and structure were subtly, but significantly, altered by fertilization, with higher alpha diversity in fertilized plots. Moreover, by comparing these data with data from 1772 soils from the Earth Microbiome Project, it was found that 20% of bacterial taxa were shared between their site and diverse globally sourced soil samples, while grassland soils shared approximately 40% of their operational taxonomic units with the current study. By spanning several orders of magnitude, the analysis suggested that extreme patchiness characterized community structure at smaller scales but that coherent patterns emerged at larger length scales.