Modern microbiology: Embracing complexity through integration across scales.

Microbes were the only form of life on Earth for most of its history, and they still account for the vast majority of life's diversity. They convert rocks to soil, produce much of the oxygen we breathe, remediate our sewage, and sustain agriculture. Microbes are vital to planetary health as they maintain biogeochemical cycles that produce and consume major greenhouse gases and support large food webs. Modern microbiologists analyze nucleic acids, proteins, and metabolites; leverage sophisticated genetic tools, software, and bioinformatic algorithms; and process and integrate complex and heterogeneous datasets so that microbial systems may be harnessed to address contemporary challenges in health, the environment, and basic science. Here, we consider an inevitably incomplete list of emergent themes in our discipline and highlight those that we recognize as the archetypes of its modern era that aim to address the most pressing problems of the 21st century.

A cryptic plasmid is among the most numerous genetic elements in the human gut.

Plasmids are extrachromosomal genetic elements that often encode fitness-enhancing features. However, many bacteria carry "cryptic" plasmids that do not confer clear beneficial functions. We identified one such cryptic plasmid, pBI143, which is ubiquitous across industrialized gut microbiomes and is 14 times as numerous as crAssphage, currently established as the most abundant extrachromosomal genetic element in the human gut. The majority of mutations in pBI143 accumulate in specific positions across thousands of metagenomes, indicating strong purifying selection. pBI143 is monoclonal in most individuals, likely due to the priority effect of the version first acquired, often from one's mother. pBI143 can transfer between Bacteroidales, and although it does not appear to impact bacterial host fitness in vivo, it can transiently acquire additional genetic content. We identified important practical applications of pBI143, including its use in identifying human fecal contamination and its potential as an alternative approach to track human colonic inflammatory states.

tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes.

tRNA is the most extensively modified RNA in cells. On average, a bacterial tRNA contains 8 modifications per molecule and a eukaryotic tRNA contains 13 modifications per molecule. Recent studies reveal that tRNA modifications are highly dynamic and respond extensively to environmental conditions. Functions of tRNA modification dynamics include enhanced, on-demand decoding of specific codons in response genes and regulation of tRNA fragment biogenesis. This review summarizes recent advances in the studies of tRNA modification dynamics in biological processes, tRNA modification erasers, and human-associated bacteria. Furthermore, we use the term "metaepitranscriptomics" to describe the potential and approach of tRNA modification studies in natural biological communities such as microbiomes. tRNA is highly modified in cells, and tRNA modifications respond extensively to environmental conditions to enhance translation of specific genes and produce tRNA fragments on demand. We review recent advances in tRNA sequencing methods, tRNA modification dynamics in biological processes, and tRNA modification studies in natural communities such as the microbiomes.

Microbial-enrichment method enables high-throughput metagenomic characterization from host-rich samples.

Host-microbe interactions have been linked to health and disease states through the use of microbial taxonomic profiling, mostly via 16S ribosomal RNA gene sequencing. However, many mechanistic insights remain elusive, in part because studying the genomes of microbes associated with mammalian tissue is difficult due to the high ratio of host to microbial DNA in such samples. Here we describe a microbial-enrichment method (MEM), which we demonstrate on a wide range of sample types, including saliva, stool, intestinal scrapings, and intestinal mucosal biopsies. MEM enabled high-throughput characterization of microbial metagenomes from human intestinal biopsies by reducing host DNA more than 1,000-fold with minimal microbial community changes (roughly 90% of taxa had no significant differences between MEM-treated and untreated control groups). Shotgun sequencing of MEM-treated human intestinal biopsies enabled characterization of both high- and low-abundance microbial taxa, pathways and genes longitudinally along the gastrointestinal tract. We report the construction of metagenome-assembled genomes directly from human intestinal biopsies for bacteria and archaea at relative abundances as low as 1%. Analysis of metagenome-assembled genomes reveals distinct subpopulation structures between the small and large intestine for some taxa. MEM opens a path for the microbiome field to acquire deeper insights into host-microbe interactions by enabling in-depth characterization of host-tissue-associated microbial communities.

Metatranscriptomic response of deep ocean microbial populations to infusions of oil and/or synthetic chemical dispersant.

Oil spills are a frequent perturbation to the marine environment that has rapid and significant impacts on the local microbiome. Previous studies have shown that exposure to synthetic dispersant alone did not enhance heterotrophic microbial activity or oxidation rates of specific hydrocarbon components but increased the abundance of some taxa (e.g., Colwellia). In contrast, exposure to oil, but not dispersants, increased the abundance of other taxa (e.g., Marinobacter) and stimulated hydrocarbon oxidation rates. Here, we advance these findings by interpreting metatranscriptomic data from this experiment to explore how and why specific components of the microbial community responded to distinct organic carbon exposure regimes. Dispersant alone was selected for a unique community and for dominant organisms that reflected treatment- and time-dependent responses. Dispersant amendment also led to diverging functional profiles among the different treatments. Similarly, oil alone was selected for a community that was distinct from treatments amended with dispersants. The presence of oil and dispersants with added nutrients led to substantial differences in microbial responses, likely suggesting increased fitness driven by the presence of additional inorganic nutrients. The oil-only additions led to a marked increase in the expression of phages, prophages, transposable elements, and plasmids (PPTEPs), suggesting that aspects of microbial community response to oil are driven by the "mobilome," potentially through viral-associated regulation of metabolic pathways in ciliates and flagellates that would otherwise throttle the microbial community through grazing.IMPORTANCEMicrocosm experiments simulated the April 2010 Deepwater Horizon oil spill by applying oil and synthetic dispersants (Corexit EC9500A and EC9527A) to deep ocean water samples. The exposure regime revealed severe negative alterations in the treatments' heterotrophic microbial activity and hydrocarbon oxidation rates. We expanded these findings by exploring metatranscriptomic signatures of the microbial communities during the chemical amendments in the microcosm experiments. Here we report how dominant organisms were uniquely associated with treatment- and time-dependent trajectories during the exposure regimes; nutrient availability was a significant factor in driving changes in metatranscriptomic responses. Remarkable signals associated with PPTEPs showed the potential role of mobilome and viral-associated survival responses. These insights underscore the time-dependent environmental perturbations of fragile marine environments under oil and anthropogenic stress.

Microbial signals drive pre-leukaemic myeloproliferation in a Tet2-deficient host.

Somatic mutations in tet methylcytosine dioxygenase 2 (TET2), which encodes an epigenetic modifier enzyme, drive the development of haematopoietic malignancies1-7. In both humans and mice, TET2 deficiency leads to increased self-renewal of haematopoietic stem cells with a net developmental bias towards the myeloid lineage1,4,8,9. However, pre-leukaemic myeloproliferation (PMP) occurs in only a fraction of Tet2-/- mice8,9 and humans with TET2 mutations1,3,5-7, suggesting that extrinsic non-cell-autonomous factors are required for disease onset. Here we show that bacterial translocation and increased interleukin-6 production, resulting from dysfunction of the small-intestinal barrier, are critical for the development of PMP in mice that lack Tet2 expression in haematopoietic cells. Furthermore, in symptom-free Tet2-/- mice, PMP can be induced by disrupting intestinal barrier integrity, or in response to systemic bacterial stimuli such as the toll-like receptor 2 agonist. PMP was reversed by antibiotic treatment and failed to develop in germ-free Tet2-/- mice, which illustrates the importance of microbial signals in the development of this condition. Our findings demonstrate the requirement for microbial-dependent inflammation in the development of PMP and provide a mechanistic basis for the variation in PMP penetrance observed in Tet2-/- mice. This study will prompt new lines of investigation that may profoundly affect the prevention and management of haematopoietic malignancies.

Accurate and complete genomes from metagenomes.

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

Early-Life Microbial Restitution Reduces Colitis Risk Promoted by Antibiotic-Induced Gut Dysbiosis in Interleukin 10-/- Mice.

BACKGROUND & AIMS: Perturbations in the early-life gut microbiome are associated with increased risk for complex immune disorders like inflammatory bowel diseases. We previously showed that maternal antibiotic-induced gut dysbiosis vertically transmitted to offspring increases experimental colitis risk in interleukin (IL) 10 gene deficient (IL10-/-) mice, a finding that may result from the loss/lack of essential microbes needed for appropriate immunologic education early in life. Here, we aimed to identify key microbes required for proper development of the early-life gut microbiome that decrease colitis risk in genetically susceptible animals. METHODS: Metagenomic sequencing followed by reconstruction of metagenome-assembled genomes was performed on fecal samples of IL10-/- mice with and without antibiotic-induced dysbiosis to identify potential missing microbial members needed for immunologic education. One high-value target strain was then engrafted early and/or late into the gut microbiomes of IL10-/- mice with antibiotic-induced dysbiosis. RESULTS: Early-, but not late-, life engraftment of a single dominant Bacteroides strain of non-antibiotic-treated IL10-/- mice was sufficient to restore the development of the gut microbiome, promote immune tolerance, and prevent colitis in IL10-/- mice that had antibiotic-induced dysbiosis. CONCLUSIONS: Restitution of a keystone microbial strain missing in the early-life antibiotic-induced gut dysbiosis results in recovery of the microbiome, proper development of immune tolerance, and reduced risk for colitis in genetically prone hosts.

Biochemical diversity of glycosphingolipid biosynthesis as a driver of Coccolithovirus competitive ecology.

Coccolithoviruses (EhVs) are large, double-stranded DNA-containing viruses that infect the single-celled, marine coccolithophore Emiliania huxleyi. Given the cosmopolitan nature and global importance of E. huxleyi as a bloom-forming, calcifying, photoautotroph, E. huxleyi-EhV interactions play a key role in oceanic carbon biogeochemistry. Virally-encoded glycosphingolipids (vGSLs) are virulence factors that are produced by the activity of virus-encoded serine palmitoyltransferase (SPT). Here, we characterize the dynamics, diversity and catalytic production of vGSLs in an array of EhV strains in relation to their SPT sequence composition and explore the hypothesis that they are a determinant of infectivity and host demise. vGSL production and diversity was positively correlated with increased virulence, virus replication rate and lytic infection dynamics in laboratory experiments, but they do not explain the success of less-virulent EhVs in natural EhV communities. The majority of EhV-derived SPT amplicon sequences associated with infected cells in the North Atlantic derived from slower infecting, less virulent EhVs. Our lab-, field- and mathematical model-based data and simulations support ecological scenarios whereby slow-infecting, less-virulent EhVs successfully compete in North Atlantic populations of E. huxleyi, through either the preferential removal of fast-infecting, virulent EhVs during active infection or by having access to a broader host range.

Oligotyping analysis of the human oral microbiome.

The Human Microbiome Project provided a census of bacterial populations in healthy individuals, but an understanding of the biomedical significance of this census has been hindered by limited taxonomic resolution. A high-resolution method termed oligotyping overcomes this limitation by evaluating individual nucleotide positions using Shannon entropy to identify the most information-rich nucleotide positions, which then define oligotypes. We have applied this method to comprehensively analyze the oral microbiome. Using Human Microbiome Project 16S rRNA gene sequence data for the nine sites in the oral cavity, we identified 493 oligotypes from the V1-V3 data and 360 oligotypes from the V3-V5 data. We associated these oligotypes with species-level taxon names by comparison with the Human Oral Microbiome Database. We discovered closely related oligotypes, differing sometimes by as little as a single nucleotide, that showed dramatically different distributions among oral sites and among individuals. We also detected potentially pathogenic taxa in high abundance in individual samples. Numerous oligotypes were preferentially located in plaque, others in keratinized gingiva or buccal mucosa, and some oligotypes were characteristic of habitat groupings such as throat, tonsils, tongue dorsum, hard palate, and saliva. The differing habitat distributions of closely related oligotypes suggest a level of ecological and functional biodiversity not previously recognized. We conclude that the Shannon entropy approach of oligotyping has the capacity to analyze entire microbiomes, discriminate between closely related but distinct taxa and, in combination with habitat analysis, provide deep insight into the microbial communities in health and disease.

DRISEE overestimates errors in metagenomic sequencing data.

The extremely high error rates reported by Keegan et al. in 'A platform-independent method for detecting errors in metagenomic sequencing data: DRISEE' (PLoS Comput Biol 2012; 8: :e1002541) for many next-generation sequencing datasets prompted us to re-examine their results. Our analysis reveals that the presence of conserved artificial sequences, e.g. Illumina adapters, and other naturally occurring sequence motifs accounts for most of the reported errors. We conclude that DRISEE reports inflated levels of sequencing error, particularly for Illumina data. Tools offered for evaluating large datasets need scrupulous review before they are implemented.

Extensive Modulation of the Fecal Metagenome in Children With Crohn's Disease During Exclusive Enteral Nutrition.

OBJECTIVES: Exploring associations between the gut microbiota and colonic inflammation and assessing sequential changes during exclusive enteral nutrition (EEN) may offer clues into the microbial origins of Crohn's disease (CD). METHODS: Fecal samples (n=117) were collected from 23 CD and 21 healthy children. From CD children fecal samples were collected before, during EEN, and when patients returned to their habitual diets. Microbiota composition and functional capacity were characterized using sequencing of the 16S rRNA gene and shotgun metagenomics. RESULTS: Microbial diversity was lower in CD than controls before EEN (P=0.006); differences were observed in 36 genera, 141 operational taxonomic units (OTUs), and 44 oligotypes. During EEN, the microbial diversity of CD children further decreased, and the community structure became even more dissimilar than that of controls. Every 10 days on EEN, 0.6 genus diversity equivalents were lost; 34 genera decreased and one increased during EEN. Fecal calprotectin correlated with 35 OTUs, 14 of which accounted for 78% of its variation. OTUs that correlated positively or negatively with calprotectin decreased during EEN. The microbiota of CD patients had a broader functional capacity than healthy controls, but diversity decreased with EEN. Genes involved in membrane transport, sulfur reduction, and nutrient biosynthesis differed between patients and controls. The abundance of genes involved in biotin (P=0.005) and thiamine biosynthesis decreased (P=0.017), whereas those involved in spermidine/putrescine biosynthesis (P=0.031), or the shikimate pathway (P=0.058), increased during EEN. CONCLUSIONS: Disease improvement following treatment with EEN is associated with extensive modulation of the gut microbiome.

Comparison of Sewage and Animal Fecal Microbiomes by Using Oligotyping Reveals Potential Human Fecal Indicators in Multiple Taxonomic Groups.

Most DNA-based microbial source tracking (MST) approaches target host-associated organisms within the order Bacteroidales, but the gut microbiota of humans and other animals contain organisms from an array of other taxonomic groups that might provide indicators of fecal pollution sources. To discern between human and nonhuman fecal sources, we compared the V6 regions of the 16S rRNA genes detected in fecal samples from six animal hosts to those found in sewage (as a proxy for humans). We focused on 10 abundant genera and used oligotyping, which can detect subtle differences between rRNA gene sequences from ecologically distinct organisms. Our analysis showed clear patterns of differential oligotype distributions between sewage and animal samples. Over 100 oligotypes of human origin occurred preferentially in sewage samples, and 99 human oligotypes were sewage specific. Sequences represented by the sewage-specific oligotypes can be used individually for development of PCR-based assays or together with the oligotypes preferentially associated with sewage to implement a signature-based approach. Analysis of sewage from Spain and Brazil showed that the sewage-specific oligotypes identified in U.S. sewage have the potential to be used as global alternative indicators of human fecal pollution. Environmental samples with evidence of prior human fecal contamination had consistent ratios of sewage signature oligotypes that corresponded to the trends observed for sewage. Our methodology represents a promising approach to identifying new bacterial taxa for MST applications and further highlights the potential of the family Lachnospiraceae to provide human-specific markers. In addition to source tracking applications, the patterns of the fine-scale population structure within fecal taxa suggest a fundamental relationship between bacteria and their hosts.

VAMPS: a website for visualization and analysis of microbial population structures.

BACKGROUND: The advent of next-generation DNA sequencing platforms has revolutionized molecular microbial ecology by making the detailed analysis of complex communities over time and space a tractable research pursuit for small research groups. However, the ability to generate 105-108 reads with relative ease brings with it many downstream complications. Beyond the computational resources and skills needed to process and analyze data, it is difficult to compare datasets in an intuitive and interactive manner that leads to hypothesis generation and testing. RESULTS: We developed the free web service VAMPS (Visualization and Analysis of Microbial Population Structures, http://vamps.mbl.edu) to address these challenges and to facilitate research by individuals or collaborating groups working on projects with large-scale sequencing data. Users can upload marker gene sequences and associated metadata; reads are quality filtered and assigned to both taxonomic structures and to taxonomy-independent clusters. A simple point-and-click interface allows users to select for analysis any combination of their own or their collaborators' private data and data from public projects, filter these by their choice of taxonomic and/or abundance criteria, and then explore these data using a wide range of analytic methods and visualizations. Each result is extensively hyperlinked to other analysis and visualization options, promoting data exploration and leading to a greater understanding of data relationships. CONCLUSIONS: VAMPS allows researchers using marker gene sequence data to analyze the diversity of microbial communities and the relationships between communities, to explore these analyses in an intuitive visual context, and to download data, results, and images for publication. VAMPS obviates the need for individual research groups to make the considerable investment in computational infrastructure and bioinformatic support otherwise necessary to process, analyze, and interpret massive amounts of next-generation sequence data. Any web-capable device can be used to upload, process, explore, and extract data and results from VAMPS. VAMPS encourages researchers to share sequence and metadata, and fosters collaboration between researchers of disparate biomes who recognize common patterns in shared data.

Diverse plasmid systems and their ecology across human gut metagenomes revealed by PlasX and MobMess.

Plasmids alter microbial evolution and lifestyles by mobilizing genes that often confer fitness in changing environments across clades. Yet our ecological and evolutionary understanding of naturally occurring plasmids is far from complete. Here we developed a machine-learning model, PlasX, which identified 68,350 non-redundant plasmids across human gut metagenomes and organized them into 1,169 evolutionarily cohesive 'plasmid systems' using our sequence containment-aware network-partitioning algorithm, MobMess. Individual plasmids were often country specific, yet most plasmid systems spanned across geographically distinct human populations. Cargo genes in plasmid systems included well-known determinants of fitness, such as antibiotic resistance, but also many others including enzymes involved in the biosynthesis of essential nutrients and modification of transfer RNAs, revealing a wide repertoire of likely fitness determinants in complex environments. Our study introduces computational tools to recognize and organize plasmids, and uncovers the ecological and evolutionary patterns of diverse plasmids in naturally occurring habitats through plasmid systems.

Wolbachia populations across organs of individual Culex pipiens: highly conserved intra-individual core pangenome with inter-individual polymorphisms.

Wolbachia is a maternally inherited intracellular bacterium that infects a wide range of arthropods including mosquitoes. The endosymbiont is widely used in biocontrol strategies due to its capacity to modulate arthropod reproduction and limit pathogen transmission. Wolbachia infections in Culex spp. are generally assumed to be monoclonal but the potential presence of genetically distinct Wolbachia subpopulations within and between individual organs has not been investigated using whole genome sequencing. Here we reconstructed Wolbachia genomes from ovary and midgut metagenomes of single naturally infected Culex pipiens mosquitoes from Southern France to investigate patterns of intra- and inter-individual differences across mosquito organs. Our analyses revealed a remarkable degree of intra-individual conservancy among Wolbachia genomes from distinct organs of the same mosquito both at the level of gene presence-absence signal and single-nucleotide polymorphisms (SNPs). Yet, we identified several synonymous and non-synonymous substitutions between individuals, demonstrating the presence of some level of genomic heterogeneity among Wolbachia that infect the same C. pipiens field population. Overall, the absence of genetic heterogeneity within Wolbachia populations in a single individual confirms the presence of a dominant Wolbachia that is maintained under strong purifying forces of evolution.

Rich microbial and depolymerising diversity in Antarctic krill gut.

With almost a quadrillion individuals, the Antarctic krill processes five million tons of organic carbon every day during austral summer. This high carbon flux requires a broad range of hydrolytic enzymes to decompose the diverse food-derived biopolymers. While krill itself possesses numerous such enzymes, it is unclear, to what extent the endogenous microbiota contribute to the hydrolytic potential of the gut environment. Here we applied amplicon sequencing, shotgun metagenomics, cultivation, and physiological assays to characterize the krill gut microbiota. The broad bacterial diversity (273 families, 919 genera, and 2,309 species) also included a complex potentially anaerobic sub-community. Plate-based assays with 198 isolated pure cultures revealed widespread capacities to utilize lipids (e.g., tributyrin), followed by proteins (casein) and to a lesser extent by polysaccharides (e.g., alginate and chitin). While most isolates affiliated with the genera Pseudoalteromonas and Psychrobacter, also Rubritalea spp. (Verrucomicrobia) were observed. The krill gut microbiota growing on marine broth agar plates possess 13,012 predicted hydrolyses; 15-fold more than previously predicted from a transcriptome-proteome compendium of krill. Cultivation-independent and -dependent approaches indicated members of the families Flavobacteriaceae and Pseudoalteromonadaceae to dominate the capacities for lipid/protein hydrolysis and to provide a plethora of carbohydrate-active enzymes, sulfatases, and laminarin- or porphyrin-depolymerizing hydrolases. Notably, also the potential to hydrolyze plastics such as polyethylene terephthalate and polylactatide was observed, affiliating mostly with Moraxellaceae. Overall, this study shows extensive microbial diversity in the krill gut, and suggests that the microbiota likely play a significant role in the nutrient acquisition of the krill by enriching its hydrolytic enzyme repertoire.IMPORTANCEThe Antarctic krill (Euphausia superba) is a keystone species of the Antarctic marine food web, connecting the productivity of phyto- and zooplankton with the nutrition of the higher trophic levels. Accordingly, krill significantly contributes to biomass turnover, requiring the decomposition of seasonally varying plankton-derived biopolymers. This study highlights the likely role of the krill gut microbiota in this ecosystem function by revealing the great number of diverse hydrolases that microbes contribute to the krill gut environment. The here resolved repertoire of hydrolytic enzymes could contribute to the overall nutritional resilience of krill and to the general organic matter cycling under changing environmental conditions in the Antarctic sea water. Furthermore, the krill gut microbiome could serve as a valuable resource of cold-adapted hydrolytic enzymes for diverse biotechnological applications.

Adaptive adjustment of significance thresholds produces large gains in microbial gene annotations and metabolic insights.

Gene function annotations enable microbial ecologists to make inferences about metabolic potential from genomes and metagenomes. However, even tools that use the same database and general approach can differ markedly in the annotations they recover. We compare three popular methods for identifying KEGG Orthologs, applying them to genomes drawn from a range of bacterial families that occupy different host-associated and free-living biomes. Our results show that by adaptively tuning sequence similarity thresholds, sensitivity can be substantially improved while maintaining accuracy. We observe the largest improvements when few reference sequences exist for a given protein family, and when annotating genomes from non-model organisms (such as gut-dwelling Lachnospiraceae). Our results suggest that straightforward heuristic adjustments can broadly improve microbial metabolic predictions.