Conducting a microbiome study.

Human microbiome research is an actively developing area of inquiry, with ramifications for our lifestyles, our interactions with microbes, and how we treat disease. Advances depend on carefully executed, controlled, and reproducible studies. Here, we provide a Primer for researchers from diverse disciplines interested in conducting microbiome research. We discuss factors to be considered in the design, execution, and data analysis of microbiome studies. These recommendations should help researchers to enter and contribute to this rapidly developing field.

Advancing the microbiome research community.

The human microbiome has become a recognized factor in promoting and maintaining health. We outline opportunities in interdisciplinary research, analytical rigor, standardization, and policy development for this relatively new and rapidly developing field. Advances in these aspects of the research community may in turn advance our understanding of human microbiome biology.

Facilitating bioinformatics reproducibility with QIIME 2 Provenance Replay.

Study reproducibility is essential to corroborate, build on, and learn from the results of scientific research but is notoriously challenging in bioinformatics, which often involves large data sets and complex analytic workflows involving many different tools. Additionally, many biologists are not trained in how to effectively record their bioinformatics analysis steps to ensure reproducibility, so critical information is often missing. Software tools used in bioinformatics can automate provenance tracking of the results they generate, removing most barriers to bioinformatics reproducibility. Here we present an implementation of that idea, Provenance Replay, a tool for generating new executable code from results generated with the QIIME 2 bioinformatics platform, and discuss considerations for bioinformatics developers who wish to implement similar functionality in their software.

Multi-omics data integration reveals metabolome as the top predictor of the cervicovaginal microenvironment.

Emerging evidence suggests that host-microbe interaction in the cervicovaginal microenvironment contributes to cervical carcinogenesis, yet dissecting these complex interactions is challenging. Herein, we performed an integrated analysis of multiple "omics" datasets to develop predictive models of the cervicovaginal microenvironment and identify characteristic features of vaginal microbiome, genital inflammation and disease status. Microbiomes, vaginal pH, immunoproteomes and metabolomes were measured in cervicovaginal specimens collected from a cohort (n = 72) of Arizonan women with or without cervical neoplasm. Multi-omics integration methods, including neural networks (mmvec) and Random Forest supervised learning, were utilized to explore potential interactions and develop predictive models. Our integrated analyses revealed that immune and cancer biomarker concentrations were reliably predicted by Random Forest regressors trained on microbial and metabolic features, suggesting close correspondence between the vaginal microbiome, metabolome, and genital inflammation involved in cervical carcinogenesis. Furthermore, we show that features of the microbiome and host microenvironment, including metabolites, microbial taxa, and immune biomarkers are predictive of genital inflammation status, but only weakly to moderately predictive of cervical neoplastic disease status. Different feature classes were important for prediction of different phenotypes. Lipids (e.g. sphingolipids and long-chain unsaturated fatty acids) were strong predictors of genital inflammation, whereas predictions of vaginal microbiota and vaginal pH relied mostly on alterations in amino acid metabolism. Finally, we identified key immune biomarkers associated with the vaginal microbiota composition and vaginal pH (MIF), as well as genital inflammation (IL-6, IL-10, MIP-1α).

Altered Faecal Microbiota Composition and Structure of Ghanaian Children with Acute Gastroenteritis.

Acute gastroenteritis (AGE) is a disease of global public health importance. Recent studies show that children with AGE have an altered gut microbiota relative to non-AGE controls. Yet, how the gut microbiota differs in Ghanaian children with and without AGE remains unclear. Here, we explore the 16S rRNA gene-based faecal microbiota profiles of Ghanaian children five years of age and younger, comprising 57 AGE cases and 50 healthy controls. We found that AGE cases were associated with lower microbial diversity and altered microbial sequence profiles relative to the controls. The faecal microbiota of AGE cases was enriched for disease-associated bacterial genera, including Enterococcus, Streptococcus, and Staphylococcus. In contrast, the faecal microbiota of controls was enriched for potentially beneficial genera, including Faecalibacterium, Prevotella, Ruminococcus, and Bacteroides. Lastly, distinct microbial correlation network characteristics were observed between AGE cases and controls, thereby supporting broad differences in faecal microbiota structure. Altogether, we show that the faecal microbiota of Ghanaian children with AGE differ from controls and are enriched for bacterial genera increasingly associated with diseases.

Experiences and lessons learned from two virtual, hands-on microbiome bioinformatics workshops.

In October of 2020, in response to the Coronavirus Disease 2019 (COVID-19) pandemic, our team hosted our first fully online workshop teaching the QIIME 2 microbiome bioinformatics platform. We had 75 enrolled participants who joined from at least 25 different countries on 6 continents, and we had 22 instructors on 4 continents. In the 5-day workshop, participants worked hands-on with a cloud-based shared compute cluster that we deployed for this course. The event was well received, and participants provided feedback and suggestions in a postworkshop questionnaire. In January of 2021, we followed this workshop with a second fully online workshop, incorporating lessons from the first. Here, we present details on the technology and protocols that we used to run these workshops, focusing on the first workshop and then introducing changes made for the second workshop. We discuss what worked well, what didn't work well, and what we plan to do differently in future workshops.

Species-level functional profiling of metagenomes and metatranscriptomes.

Functional profiles of microbial communities are typically generated using comprehensive metagenomic or metatranscriptomic sequence read searches, which are time-consuming, prone to spurious mapping, and often limited to community-level quantification. We developed HUMAnN2, a tiered search strategy that enables fast, accurate, and species-resolved functional profiling of host-associated and environmental communities. HUMAnN2 identifies a community's known species, aligns reads to their pangenomes, performs translated search on unclassified reads, and finally quantifies gene families and pathways. Relative to pure translated search, HUMAnN2 is faster and produces more accurate gene family profiles. We applied HUMAnN2 to study clinal variation in marine metabolism, ecological contribution patterns among human microbiome pathways, variation in species' genomic versus transcriptional contributions, and strain profiling. Further, we introduce 'contributional diversity' to explain patterns of ecological assembly across different microbial community types.

Qiita: rapid, web-enabled microbiome meta-analysis.

Multi-omic insights into microbiome function and composition typically advance one study at a time. However, in order for relationships across studies to be fully understood, data must be aggregated into meta-analyses. This makes it possible to generate new hypotheses by finding features that are reproducible across biospecimens and data layers. Qiita dramatically accelerates such integration tasks in a web-based microbiome-comparison platform, which we demonstrate with Human Microbiome Project and Integrative Human Microbiome Project (iHMP) data.

Leptospira in river and soil in a highly endemic area of Ecuador.

BACKGROUND: Leptospira are shed into the environment via urine of infected animals. Rivers are thought to be an important risk factor for transmission to humans, though much is unknown about the types of environment or characteristics that favor survival. To address this, we screened for Leptospira DNA in two rivers in rural Ecuador where Leptospirosis is endemic. RESULTS: We collected 112 longitudinal samples and recorded pH, temperature, river depth, precipitation, and dissolved oxygen. We also performed a series of three experiments designed to provide insight into Leptospira presence in the soil. In the first soil experiment, we characterized prevalence and co-occurrence of Leptospira with other bacterial taxa in the soil at dispersed sites along the rivers (n = 64). In the second soil experiment, we collected 24 river samples and 48 soil samples at three points along eight transects to compare the likelihood of finding Leptospira in the river and on the shore at different distances from the river. In a third experiment, we tested whether Leptospira presence is associated with soil moisture by collecting 25 soil samples from two different sites. In our river experiment, we found pathogenic Leptospira in only 4 (3.7%) of samples. In contrast, pathogenic Leptospira species were found in 22% of shore soil at dispersed sites, 16.7% of soil samples (compared to 4.2% of river samples) in the transects, and 40% of soil samples to test for associations with soil moisture. CONCLUSIONS: Our data are limited to two sites in a highly endemic area, but the scarcity of Leptospira DNA in the river is not consistent with the widespread contention of the importance of river water for leptospirosis transmission. While Leptospira may be shed directly into the river, onto the shores, or washed into the river from more remote sites, massive dilution and limited persistence in rivers may reduce the environmental load and therefore, the epidemiological significance of such sources. It is also possible that transmission may occur more frequently on shores where people are liable to be barefoot. Molecular studies that further explore the role of rivers and water bodies in the epidemiology of leptospirosis are needed.

The international sinonasal microbiome study: A multicentre, multinational characterization of sinonasal bacterial ecology.

The sinonasal microbiome remains poorly defined, with our current knowledge based on a few cohort studies whose findings are inconsistent. Furthermore, the variability of the sinus microbiome across geographical divides remains unexplored. We characterize the sinonasal microbiome and its geographical variations in both health and disease using 16S rRNA gene sequencing of 410 individuals from across the world. Although the sinus microbial ecology is highly variable between individuals, we identify a core microbiome comprised of Corynebacterium, Staphylococcus, Streptococcus, Haemophilus and Moraxella species in both healthy and chronic rhinosinusitis (CRS) cohorts. Corynebacterium (mean relative abundance = 44.02%) and Staphylococcus (mean relative abundance = 27.34%) appear particularly dominant in the majority of patients sampled. Amongst patients suffering from CRS with nasal polyps, a statistically significant reduction in relative abundance of Corynebacterium (40.29% vs 50.43%; P = .02) was identified. Despite some measured differences in microbiome composition and diversity between some of the participating centres in our cohort, these differences would not alter the general pattern of core organisms described. Nevertheless, atypical or unusual organisms reported in short-read amplicon sequencing studies and that are not part of the core microbiome should be interpreted with caution. The delineation of the sinonasal microbiome and standardized methodology described within our study will enable further characterization and translational application of the sinus microbiota.

Domestic canines do not display evidence of gut microbial dysbiosis in the presence of Clostridioides (Clostridium) difficile, despite cellular susceptibility to its toxins.

Clostridioides difficile infection (CDI) is an emerging public health threat and C. difficile is the most common cause of antimicrobial-associated diarrhea worldwide and the leading cause of hospital-associated infections in the US, yet the burden of community-acquired infections (CAI) is poorly understood. Characterizing C. difficile isolated from canines is important for understanding the role that canines may play in CAI. In addition, several studies have suggested that canines carry toxigenic C. difficile asymptomatically, which may imply that there are mechanisms responsible for resistance to CDI in canines that could be exploited to help combat human CDI. To assess the virulence potential of canine-derived C. difficile, we tested whether toxins TcdA and TcdB (hereafter toxins) derived from a canine isolate were capable of causing tight junction disruptions to colonic epithelial cells. Additionally, we addressed whether major differences exist between human and canine cells regarding C. difficile pathogenicity by exposing them to identical toxins. We then examined the canine gut microbiome associated with C. difficile carriage using 16S rRNA gene sequencing and searched for deviations from homeostasis as an indicator of CDI. Finally, we queried 16S rRNA gene sequences for bacterial taxa that may be associated with resistance to CDI in canines. Clostridioides difficile isolated from a canine produced toxins that reduced tight junction integrity in both human and canine cells in vitro. However, canine guts were not dysbiotic in the presence of C. difficile. These findings support asymptomatic carriage in canines and, furthermore, suggest that there are features of the gut microbiome and/or a canine-specific immune response that may protect canines against CDI. We identified two biologically relevant bacteria that may aid in CDI resistance in canines: 1) Clostridium hiranonis, which synthesizes secondary bile acids that have been shown to provide resistance to CDI in mice; and 2) Sphingobacterium faecium, which produces sphingophospholipids that may be associated with regulating homeostasis in the canine gut. Our findings suggest that canines may be cryptic reservoirs for C. difficile and, furthermore, that mechanisms of CDI resistance in the canine gut could provide insights into targeted therapeutics for human CDI.

Human gut microbiome viewed across age and geography.

Gut microbial communities represent one source of human genetic and metabolic diversity. To examine how gut microbiomes differ among human populations, here we characterize bacterial species in fecal samples from 531 individuals, plus the gene content of 110 of them. The cohort encompassed healthy children and adults from the Amazonas of Venezuela, rural Malawi and US metropolitan areas and included mono- and dizygotic twins. Shared features of the functional maturation of the gut microbiome were identified during the first three years of life in all three populations, including age-associated changes in the genes involved in vitamin biosynthesis and metabolism. Pronounced differences in bacterial assemblages and functional gene repertoires were noted between US residents and those in the other two countries. These distinctive features are evident in early infancy as well as adulthood. Our findings underscore the need to consider the microbiome when evaluating human development, nutritional needs, physiological variations and the impact of westernization.

Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing.

High-throughput sequencing has revolutionized microbial ecology, but read quality remains a considerable barrier to accurate taxonomy assignment and α-diversity assessment for microbial communities. We demonstrate that high-quality read length and abundance are the primary factors differentiating correct from erroneous reads produced by Illumina GAIIx, HiSeq and MiSeq instruments. We present guidelines for user-defined quality-filtering strategies, enabling efficient extraction of high-quality data and facilitating interpretation of Illumina sequencing results.

From Rare to Dominant: a Fine-Tuned Soil Bacterial Bloom during Petroleum Hydrocarbon Bioremediation.

Hydrocarbons are worldwide-distributed pollutants that disturb various ecosystems. The aim of this study was to characterize the short-lapse dynamics of soil microbial communities in response to hydrocarbon pollution and different bioremediation treatments. Replicate diesel-spiked soil microcosms were inoculated with either a defined bacterial consortium or a hydrocarbonoclastic bacterial enrichment and incubated for 12 weeks. The microbial community dynamics was followed weekly in microcosms using Illumina 16S rRNA gene sequencing. Both the bacterial consortium and enrichment enhanced hydrocarbon degradation in diesel-polluted soils. A pronounced and rapid bloom of a native gammaproteobacterium was observed in all diesel-polluted soils. A unique operational taxonomic unit (OTU) related to the Alkanindiges genus represented ∼ 0.1% of the sequences in the original community but surprisingly reached >60% after 6 weeks. Despite this Alkanindiges-related bloom, inoculated strains were maintained in the community and may explain the differences in hydrocarbon degradation. This study shows the detailed dynamics of a soil bacterial bloom in response to hydrocarbon pollution, resembling microbial blooms observed in marine environments. Rare community members presumably act as a reservoir of ecological functions in high-diversity environments, such as soils. This rare-to-dominant bacterial shift illustrates the potential role of a rare biosphere facing drastic environmental disturbances. Additionally, it supports the concept of "conditionally rare taxa," in which rareness is a temporary state conditioned by environmental constraints.

Accurate Estimation of Fungal Diversity and Abundance through Improved Lineage-Specific Primers Optimized for Illumina Amplicon Sequencing.

While high-throughput sequencing methods are revolutionizing fungal ecology, recovering accurate estimates of species richness and abundance has proven elusive. We sought to design internal transcribed spacer (ITS) primers and an Illumina protocol that would maximize coverage of the kingdom Fungi while minimizing nontarget eukaryotes. We inspected alignments of the 5.8S and large subunit (LSU) ribosomal genes and evaluated potential primers using PrimerProspector. We tested the resulting primers using tiered-abundance mock communities and five previously characterized soil samples. We recovered operational taxonomic units (OTUs) belonging to all 8 members in both mock communities, despite DNA abundances spanning 3 orders of magnitude. The expected and observed read counts were strongly correlated (r = 0.94 to 0.97). However, several taxa were consistently over- or underrepresented, likely due to variation in rRNA gene copy numbers. The Illumina data resulted in clustering of soil samples identical to that obtained with Sanger sequence clone library data using different primers. Furthermore, the two methods produced distance matrices with a Mantel correlation of 0.92. Nonfungal sequences comprised less than 0.5% of the soil data set, with most attributable to vascular plants. Our results suggest that high-throughput methods can produce fairly accurate estimates of fungal abundances in complex communities. Further improvements might be achieved through corrections for rRNA copy number and utilization of standardized mock communities. IMPORTANCE: Fungi play numerous important roles in the environment. Improvements in sequencing methods are providing revolutionary insights into fungal biodiversity, yet accurate estimates of the number of fungal species (i.e., richness) and their relative abundances in an environmental sample (e.g., soil, roots, water, etc.) remain difficult to obtain. We present improved methods for high-throughput Illumina sequencing of the species-diagnostic fungal ribosomal marker gene that improve the accuracy of richness and abundance estimates. The improvements include new PCR primers and library preparation, validation using a known mock community, and bioinformatic parameter tuning.

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 × 10(11) 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.

Continental-scale variation in seaweed host-associated bacterial communities is a function of host condition, not geography.

Interactions between hosts and associated microbial communities can fundamentally shape the development and ecology of 'holobionts', from humans to marine habitat-forming organisms such as seaweeds. In marine systems, planktonic microbial community structure is mainly driven by geography and related environmental factors, but the large-scale drivers of host-associated microbial communities are largely unknown. Using 16S-rRNA gene sequencing, we characterized 260 seaweed-associated bacterial and archaeal communities on the kelp Ecklonia radiata from three biogeographical provinces spanning 10° of latitude and 35° of longitude across the Australian continent. These phylogenetically and taxonomically diverse communities were more strongly and consistently associated with host condition than geographical location or environmental variables, and a 'core' microbial community characteristic of healthy kelps appears to be lost when hosts become stressed. Microbial communities on stressed individuals were more similar to each other among locations than those on healthy hosts. In contrast to biogeographical patterns of planktonic marine microbial communities, host traits emerge as critical determinants of associated microbial community structure of these holobionts, even at a continental scale.

Quantitative microbial ecology through stable isotope probing.

Bacteria grow and transform elements at different rates, and as yet, quantifying this variation in the environment is difficult. Determining isotope enrichment with fine taxonomic resolution after exposure to isotope tracers could help, but there are few suitable techniques. We propose a modification to stable isotope probing (SIP) that enables the isotopic composition of DNA from individual bacterial taxa after exposure to isotope tracers to be determined. In our modification, after isopycnic centrifugation, DNA is collected in multiple density fractions, and each fraction is sequenced separately. Taxon-specific density curves are produced for labeled and nonlabeled treatments, from which the shift in density for each individual taxon in response to isotope labeling is calculated. Expressing each taxon's density shift relative to that taxon's density measured without isotope enrichment accounts for the influence of nucleic acid composition on density and isolates the influence of isotope tracer assimilation. The shift in density translates quantitatively to isotopic enrichment. Because this revision to SIP allows quantitative measurements of isotope enrichment, we propose to call it quantitative stable isotope probing (qSIP). We demonstrated qSIP using soil incubations, in which soil bacteria exhibited strong taxonomic variations in (18)O and (13)C composition after exposure to [(18)O]water or [(13)C]glucose. The addition of glucose increased the assimilation of (18)O into DNA from [(18)O]water. However, the increase in (18)O assimilation was greater than expected based on utilization of glucose-derived carbon alone, because the addition of glucose indirectly stimulated bacteria to utilize other substrates for growth. This example illustrates the benefit of a quantitative approach to stable isotope probing.

Cross-biome metagenomic analyses of soil microbial communities and their functional attributes.

For centuries ecologists have studied how the diversity and functional traits of plant and animal communities vary across biomes. In contrast, we have only just begun exploring similar questions for soil microbial communities despite soil microbes being the dominant engines of biogeochemical cycles and a major pool of living biomass in terrestrial ecosystems. We used metagenomic sequencing to compare the composition and functional attributes of 16 soil microbial communities collected from cold deserts, hot deserts, forests, grasslands, and tundra. Those communities found in plant-free cold desert soils typically had the lowest levels of functional diversity (diversity of protein-coding gene categories) and the lowest levels of phylogenetic and taxonomic diversity. Across all soils, functional beta diversity was strongly correlated with taxonomic and phylogenetic beta diversity; the desert microbial communities were clearly distinct from the nondesert communities regardless of the metric used. The desert communities had higher relative abundances of genes associated with osmoregulation and dormancy, but lower relative abundances of genes associated with nutrient cycling and the catabolism of plant-derived organic compounds. Antibiotic resistance genes were consistently threefold less abundant in the desert soils than in the nondesert soils, suggesting that abiotic conditions, not competitive interactions, are more important in shaping the desert microbial communities. As the most comprehensive survey of soil taxonomic, phylogenetic, and functional diversity to date, this study demonstrates that metagenomic approaches can be used to build a predictive understanding of how microbial diversity and function vary across terrestrial biomes.

Gut microbial and short-chain fatty acid profiles in adults with chronic constipation before and after treatment with lubiprostone.

Identifying specific gut microorganisms associated with chronic constipation may be useful for diagnostic and therapeutic purposes. The objective of this study was to evaluate whether or not the gut microbial community of constipated subjects had specific microbial signatures and to assess the effects of lubiprostone treatment on the gut microbial community. Stool diaries, breath H2 and CH4 levels, and stool samples were collected from ten healthy subjects and nine patients meeting the Rome III criteria for chronic functional constipation. Constipated subjects received lubiprostone for four weeks, during which stool diaries were maintained. Stool samples were evaluated for gut microbial communities using pyrosequencing and quantitative real-time PCR (qPCR) targeting 16S-rRNA gene, along with concentrations of short-chain fatty acids (SCFAs) using high-performance liquid chromatography. Prior to treatment, gut microbial profiles were similar between constipated subjects and healthy subjects, while iso-butyrate levels were significantly higher in constipated subjects compared with healthy subjects. Despite increases in stool frequency and improvements in consistency after lubiprostone treatment, gut microbial profiles and community diversity after treatment showed no significant change compared to before treatment. While we did not observe a significant difference in either breath methane or archaeal abundance between the stool samples of healthy and constipated subjects, we confirmed a strong correlation between archaeal abundance measured by qPCR and the amount of methane gas exhaled in the fasting breath. Butyrate levels, however, were significantly higher in the stool samples of constipated subjects after lubiprostone treatment, suggesting that lubiprostone treatment had an effect on the net accumulation of SCFAs in the gut. In conclusion, lubiprostone treatment improved constipation symptoms and increased levels of butyrate without substantial modification of the gut microbial structure.