Extensive strain-level copy-number variation across human gut microbiome species

A predictive index for health status using species-level gut microbiome profiling

Providing insight into one’s health status from a gut microbiome sample is an important clinical goal in current human microbiome research. Herein, we introduce the Gut Microbiome Health Index (GMHI), a biologically-interpretable mathematical formula for predicting the likelihood of disease independent of the clinical diagnosis. GMHI is formulated upon 50 microbial species associated with healthy gut ecosystems. These species are identified through a multi-study, integrative analysis on 4347 human stool metagenomes from 34 published studies across healthy and 12 different nonhealthy conditions, i.e., disease or abnormal bodyweight. When demonstrated on our population-scale meta-dataset, GMHI is the most robust and consistent predictor of disease presence (or absence) compared to α-diversity indices. Validation on 679 samples from 9 additional studies results in a balanced accuracy of 73.7% in distinguishing healthy from non-healthy groups. Our findings suggest that gut taxonomic signatures can predict health status, and highlight how data sharing efforts can provide broadly applicable discoveries.

A biologically-interpretable and robust metric that provides insight into one’s health status from a gut microbiome sample is an important clinical goal in current human microbiome research. Herein, the authors introduce a species-level index that predicts the likelihood of having a disease.

Gut microbiota in human metabolic health and disease

Observational findings achieved during the past two decades suggest that the intestinal microbiota may contribute to the metabolic health of the human host and, when aberrant, to the pathogenesis of various common metabolic disorders including obesity, type 2 diabetes, non-alcoholic liver disease, cardio-metabolic diseases and malnutrition. However, to gain a mechanistic understanding of how the gut microbiota affects host metabolism, research is moving from descriptive microbiota census analyses to cause-and-effect studies. Joint analyses of high-throughput human multi-omics data, including metagenomics and metabolomics data, together with measures of host physiology and mechanistic experiments in humans, animals and cells hold potential as initial steps in the identification of potential molecular mechanisms behind reported associations. In this Review, we discuss the current knowledge on how gut microbiota and derived microbial compounds may link to metabolism of the healthy host or to the pathogenesis of common metabolic diseases. We highlight examples of microbiota-targeted interventions aiming to optimize metabolic health, and we provide perspectives for future basic and translational investigations within the nascent and promising research field.

Diversity within species: interpreting strains in microbiomes

Studying within-species variation has traditionally been limited to culturable bacterial isolates and low-resolution microbial community fingerprinting. Metagenomic sequencing and technical advances have enabled culture-free, high-resolution strain and subspecies analyses at high throughput and in complex environments. This holds great scientific promise but has also led to an overwhelming number of methods and terms to describe infraspecific variation. This Review aims to clarify these advances by focusing on the diversity within bacterial and archaeal species in the context of microbiomics. We cover foundational microevolutionary concepts relevant to population genetics and summarize how within-species variation can be studied and stratified directly within microbial communities with a focus on metagenomics. Finally, we describe how common applications of within-species variation can be achieved using metagenomic data. We aim to guide the selection of appropriate terms and analytical approaches to facilitate researchers in benefiting from the increasing availability of large, high-resolution microbiome genetic sequencing data.

Screening of autochthonous vaginal beneficial lactobacilli strains by their growth at high temperatures for technological applications

The pharmaceutical industry shows an emerging interest in formulas that contain live and beneficial microorganisms, also known as probiotics or pharmabiotics, which in many cases, are host-specific. The resistance to higher temperature is an essential feature of these microorganisms when working on the design of products for vaginal formula. In order to obtain a high number of viable cells and a prolonged shelf life in the designed product, it is required to apply technological procedures using high temperatures or abrupt changes of them, which result in conditions that are different from the optimal growth temperature and can affect the metabolic capabilities of the bacteria when administered to the host in order to reestablish the ecological mucosa. The aim of this work was to evaluate the behavior of 30 different species and strains of autochthonous beneficial vaginal lactobacilli (BVL) when exposed to high temperatures, determine their survival capabilities and analyze their pre-adaptation to those temperatures, in order that they still maintain their viability after technological processes and further conservation. BVL were exhibited to temperatures higher than optimal, with the purpose of evaluating their growth kinetics and parameters. Later, they were exposed to higher temperatures, and then, returned to their optimal, to determine if they were able to grow again. The strains that showed higher resistance were selected, and their viability and beneficial properties studied further. The growth kinetics of strains exposed to higher temperatures showed different patterns, which provided evidence that the thermal adaptation is strain-dependent and is not related to any particular species and/or metabolic group in which the strains were taxonomically classified. The pre-adaptive step allowed the growth of some of the strains, preserving their viability and probiotic properties after the high temperatures were applied. The results shows that BVL can be exposed to high temperatures used in different technological processes that are applied for pharmabiotic formulations, such as spray dried or vacuum rotary evaporation, and/or during the conservation period. The results obtained indicate that some specific BVL strains resist high temperatures and grow afterwards at optimal conditions.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

Vast Differences in Strain-Level Diversity in the Gut Microbiota of Two Closely Related Honey Bee Species

Most bacterial species encompass strains with vastly different gene content. Strain diversity in microbial communities is therefore considered to be of functional importance. Yet little is known about the extent to which related microbial communities differ in diversity at this level and which underlying mechanisms may constrain and maintain strain-level diversity. Here, we used shotgun metagenomics to characterize and compare the gut microbiota of two honey bee species, Apis mellifera and Apis cerana, which diverged about 6 mya. Although the host species are colonized largely by the same bacterial 16S rRNA phylotypes, we find that their communities are host specific when analyzed with genomic resolution. Moreover, despite their similar ecology, A. mellifera displayed a much higher diversity of strains and functional gene content in the microbiota compared to A. cerana, both per colony and per individual bee. In particular, the gene repertoire for polysaccharide degradation was massively expanded in the microbiota of A. mellifera relative to A. cerana. Bee management practices, divergent ecological adaptation, or habitat size may have contributed to the observed differences in microbiota genomic diversity of these key pollinator species. Our results illustrate that the gut microbiota of closely related animal hosts can differ vastly in genomic diversity while displaying similar levels of diversity based on the 16S rRNA gene. Such differences are likely to have consequences for gut microbiota functioning and host-symbiont interactions, highlighting the need for metagenomic studies to understand the ecology and evolution of microbial communities.

What Is Metagenomics Teaching Us, and What Is Missed?

Shotgun metagenomic sequencing has revolutionized our ability to detect and characterize the diversity and function of complex microbial communities. In this review, we highlight the benefits of using metagenomics as well as the breadth of conclusions that can be made using currently available analytical tools, such as greater resolution of species and strains across phyla and functional content, while highlighting challenges of metagenomic data analysis. Major challenges remain in annotating function, given the dearth of functional databases for environmental bacteria compared to model organisms, and the technical difficulties of metagenome assembly and phasing in heterogeneous environmental samples. In the future, improvements and innovation in technology and methodology will lead to lowered costs. Data integration using multiple technological platforms will lead to a better understanding of how to harness metagenomes. Subsequently, we will be able not only to characterize complex microbiomes but also to manipulate communities to achieve prosperous outcomes for health, agriculture, and environmental sustainability.

Vast Differences in Strain-Level Diversity in the Gut Microbiota of Two Closely Related Honey Bee Species

Most bacterial species encompass strains with vastly different gene content. Strain diversity in microbial communities is therefore considered to be of functional importance. Yet little is known about the extent to which related microbial communities differ in diversity at this level and which underlying mechanisms may constrain and maintain strain-level diversity. Here, we used shotgun metagenomics to characterize and compare the gut microbiota of two honey bee species, Apis mellifera and Apis cerana, which diverged about 6 mya. Although the host species are colonized largely by the same bacterial 16S rRNA phylotypes, we find that their communities are host specific when analyzed with genomic resolution. Moreover, despite their similar ecology, A. mellifera displayed a much higher diversity of strains and functional gene content in the microbiota compared to A. cerana, both per colony and per individual bee. In particular, the gene repertoire for polysaccharide degradation was massively expanded in the microbiota of A. mellifera relative to A. cerana. Bee management practices, divergent ecological adaptation, or habitat size may have contributed to the observed differences in microbiota genomic diversity of these key pollinator species. Our results illustrate that the gut microbiota of closely related animal hosts can differ vastly in genomic diversity while displaying similar levels of diversity based on the 16S rRNA gene. Such differences are likely to have consequences for gut microbiota functioning and host-symbiont interactions, highlighting the need for metagenomic studies to understand the ecology and evolution of microbial communities.

•

Metagenomics reveals differences in gut microbiota diversity beyond the 16S rRNA gene

•

Apis cerana and Apis mellifera harbor distinct species and strains in their gut

•

Diversity is much higher in A. mellifera per individual bee and within colonies

•

Major differences in functions are related to polysaccharide degradation

A Genomic Toolkit for the Mechanistic Dissection of Intractable Human Gut Bacteria

Despite the remarkable microbial diversity found within humans, our ability to link genes to phenotypes is based upon a handful of model microorganisms. We report a comparative genomics platform for Eggerthella lenta and other Coriobacteriia, a neglected taxon broadly relevant to human health and disease. We uncover extensive genetic and metabolic diversity and validate a tool for mapping phenotypes to genes and sequence variants. We also present a tool for the quantification of strains from metagenomic sequencing data, enabling the identification of genes that predict bacterial fitness. Competitive growth is reproducible under laboratory conditions and attributable to intrinsic growth rates and resource utilization. Unique signatures of in vivo competition in gnotobiotic mice include an adhesin enriched in poor colonizers. Together, these computational and experimental resources represent a strong foundation for the continued mechanistic dissection of the Coriobacteriia and a template that can be applied to study other genetically intractable taxa.

The gut microbiome and metabolic syndrome

The metabolic syndrome (MetS) is a constellation of risk factors that, if left untreated, will often progress to greater metabolic defects such as type 2 diabetes and nonalcoholic fatty liver disease. While these risk factors have been established for over 40 years, the definition of MetS warrants reconsideration in light of the substantial data that have emerged from studies of the gut microbiome. In this Review we present the existing recent literature that supports the gut microbiome's potential influence on the various risk factors of MetS. The interplay of the intestinal microbiota with host metabolism has been shown to be mediated by a myriad of factors, including a defective gut barrier, bile acid metabolism, antibiotic use, and the pleiotropic effects of microbially produced metabolites. These data show that events that start in the gut, often in response to external cues such as diet and circadian disruption, have far-reaching effects beyond the gut.

MSPminer: abundance-based reconstitution of microbial pan-genomes from shotgun metagenomic data

Motivation

Analysis toolkits for shotgun metagenomic data achieve strain-level characterization of complex microbial communities by capturing intra-species gene content variation. Yet, these tools are hampered by the extent of reference genomes that are far from covering all microbial variability, as many species are still not sequenced or have only few strains available. Binning co-abundant genes obtained from de novo assembly is a powerful reference-free technique to discover and reconstitute gene repertoire of microbial species. While current methods accurately identify species core parts, they miss many accessory genes or split them into small gene groups that remain unassociated to core clusters.

Results

We introduce MSPminer, a computationally efficient software tool that reconstitutes Metagenomic Species Pan-genomes (MSPs) by binning co-abundant genes across metagenomic samples. MSPminer relies on a new robust measure of proportionality coupled with an empirical classifier to group and distinguish not only species core genes but accessory genes also. Applied to a large scale metagenomic dataset, MSPminer successfully delineates in a few hours the gene repertoires of 1661 microbial species with similar specificity and higher sensitivity than existing tools. The taxonomic annotation of MSPs reveals microorganisms hitherto unknown and brings coherence in the nomenclature of the species of the human gut microbiota. The provided MSPs can be readily used for taxonomic profiling and biomarkers discovery in human gut metagenomic samples. In addition, MSPminer can be applied on gene count tables from other ecosystems to perform similar analyses.

Availability and implementation

The binary is freely available for non-commercial users at www.enterome.com/downloads.

Supplementary information

Supplementary data are available at Bioinformatics online.

Apis cerana gut microbiota contribute to host health though stimulating host immune system and strengthening host resistance to Nosema ceranae

Gut microbial communities play vital roles in the modulation of many insects' immunity, including Apis mellifera. However, little is known about the interaction of Apis cerana gut bacteria and A. cerana immune system. Here in this study, we conducted a comparison between germ-free gut microbiota deficient (GD) workers and conventional gut community (CV) workers, to reveal the possible impact of gut microbiota on the expression of A. cerana antimicrobial peptides and immune regulate pathways. We also test whether A. cerana gut microbiota can strengthen host resistance to Nosema ceranae. We find that the expression of apidaecin, abaecin and hymenoptaecin were significantly upregulated with the presence of gut bacteria, and JNK pathway was activated; in the meanwhile, the existence of gut bacteria inhibited the proliferation of Nosema ceranae. These demonstrated the essential role of A. cerana gut microbiota to host health and provided critical insight into the honeybee host-microbiome interaction.

Adaptation of the Gut Microbiota to Modern Dietary Sugars and Sweeteners

The consumption of sugar has become central to the Western diet. Cost and health concerns associated with sucrose spurred the development and consumption of other sugars and sweeteners, with the average American consuming 10 times more sugar than 100 y ago. In this review, we discuss how gut microbes are affected by changes in the consumption of sugars and other sweeteners through transcriptional, abundance, and genetic adaptations. We propose that these adaptations result in microbes taking on different metabolic, ecological, and genetic profiles along the intestinal tract. We suggest novel approaches to assess the consequences of these changes on host-microbe interactions to determine the safety of novel sugars and sweeteners.

Association of HLA-dependent islet autoimmunity with systemic antibody responses to intestinal commensal bacteria in children

Microbiome sequence analyses have suggested that changes in gut bacterial composition are associated with autoimmune disease in humans and animal models. However, little is known of the mechanisms through which the gut microbiota influences autoimmune responses to distant tissues. Here, we evaluated systemic antibody responses against cultured human gut bacterial strains to determine whether observed patterns of anticommensal antibody (ACAb) responses are associated with type 1 diabetes (T1D) in two cohorts of pediatric study participants. In the first cohort, ACAb responses in sera collected from participants within 6 months of T1D diagnosis were compared with age-matched healthy controls and also with patients with recent onset Crohn's disease. ACAb responses against multiple bacterial species discriminated among these three groups. In the second cohort, we asked whether ACAb responses present before diagnosis were associated with later T1D development and with HLA genotype in participants who were discordant for subsequent progression to diabetes. Serum IgG2 antibodies against Roseburia faecis and against a bacterial consortium were associated with future T1D diagnosis in an HLA DR3/DR4 haplotype-dependent manner. These analyses reveal associations between antibody responses to intestinal microbes and HLA-DR genotype and islet autoantibody specificity and with a future diagnosis of T1D. Further, we present a platform to investigate antibacterial antibodies in biological fluids that is applicable to studies of autoimmune diseases and responses to therapeutic interventions.

Spaceship Earth Revisited: The Co-Benefits of Overcoming Biological Extinction of Experience at the Level of Person, Place and Planet

Extensive research underscores that we interpret the world through metaphors; moreover, common metaphors are a useful means to enhance the pursuit of personal and collective goals. In the context of planetary health—defined as the interdependent vitality of all natural and anthropogenic ecosystems (social, political and otherwise)—one enduring metaphor can be found in the concept of “Spaceship Earth”. Although not without criticism, the term “Spaceship Earth” has been useful to highlight both resource limitations and the beauty and fragility of delicate ecosystems that sustain life. Rene Dubos, who helped popularize the term, underscored the need for an exposome perspective, one that examines the total accumulated environmental exposures (both detrimental and beneficial) that predict the biological responses of the “total organism to the total environment” over time. In other words, how large-scale environmental changes affect us all personally, albeit in individualized ways. This commentary focuses the ways in which microbes, as an essential part of all ecosystems, provide a vital link between personal and planetary systems, and mediate the biopsychosocial aspects of our individualized experience—and thus health—over our life course journey. A more fine-grained understanding of these dynamics and our power to change them, personally and collectively, lies at the core of restoring “ecosystems balance” for person, place and planet. In particular, restoring human connectedness to the natural world, sense of community and shared purpose must occur in tandem with technological solutions, and will enhance individual empowerment for personal well-being, as well as our collective potential to overcome our grand challenges. Such knowledge can help shape the use of metaphor and re-imagine solutions and novel ways for restoration or rewilding of ecosystems, and the values, behaviors and attitudes to light the path toward exiting the Anthropocene.

Host-Specific Evolutionary and Transmission Dynamics Shape the Functional Diversification of Staphylococcus epidermidis in Human Skin

Metagenomic inferences of bacterial strain diversity and infectious disease transmission studies largely assume a dominant, within-individual haplotype. We hypothesize that within-individual bacterial population diversity is critical for homeostasis of a healthy microbiome and infection risk. We characterized the evolutionary trajectory and functional distribution of Staphylococcus epidermidis-a keystone skin microbe and opportunistic pathogen. Analyzing 1,482 S. epidermidis genomes from 5 healthy individuals, we found that skin S. epidermidis isolates coalesce into multiple founder lineages rather than a single colonizer. Transmission events, natural selection, and pervasive horizontal gene transfer result in population admixture within skin sites and dissemination of antibiotic resistance genes within-individual. We provide experimental evidence for how admixture can modulate virulence and metabolism. Leveraging data on the contextual microbiome, we assess how interspecies interactions can shape genetic diversity and mobile gene elements. Our study provides insights into how within-individual evolution of human skin microbes shapes their functional diversification.

High-resolution sweep metagenomics using fast probabilistic inference

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

The antibody/microbiota interface in health and disease

The human intestine is densely colonized with commensal microbes that stimulate the immune system. While secretory Immunoglobulin (Ig) A is known to play a crucial role in gut microbiota compartmentalization, secretory IgM, and systemic IgG have recently been highlighted in host-microbiota interactions as well. In this review, we discuss important aspects of secretory IgA biology, but rather than focusing on mechanistic aspects of IgA impact on microbiota, we stress the current knowledge of systemic antibody responses to whole gut microbiota, in particular their generation, specificities, and function. We also provide a comprehensive picture of secretory IgM biology. Finally, therapeutic and diagnostic implications of these novel findings for the treatment of various diseases are outlined.

Comprehensive analysis of chromosomal mobile genetic elements in the gut microbiome reveals phylum-level niche-adaptive gene pools

Mobile genetic elements (MGEs) drive extensive horizontal transfer in the gut microbiome. This transfer could benefit human health by conferring new metabolic capabilities to commensal microbes, or it could threaten human health by spreading antibiotic resistance genes to pathogens. Despite their biological importance and medical relevance, MGEs from the gut microbiome have not been systematically characterized. Here, we present a comprehensive analysis of chromosomal MGEs in the gut microbiome using a method that enables the identification of the mobilizable unit of MGEs. We curated a database of 5,219 putative MGEs encompassing seven MGE classes called ImmeDB. We observed that many MGEs carry genes that could confer an adaptive advantage to the gut environment including gene families involved in antibiotic resistance, bile salt detoxification, mucus degradation, capsular polysaccharide biosynthesis, polysaccharide utilization, and sporulation. We find that antibiotic resistance genes are more likely to be spread by conjugation via integrative conjugative elements or integrative mobilizable elements than transduction via prophages. Horizontal transfer of MGEs is extensive within phyla but rare across phyla, supporting phylum level niche-adaptive gene pools in the gut microbiome. ImmeDB will be a valuable resource for future studies on the gut microbiome and MGE communities.

Population Genetics in the Human Microbiome

While the human microbiome's structure and function have been extensively studied, its within-species genetic diversity is less well understood. However, genetic mutations in the microbiome can confer biomedically relevant traits, such as the ability to extract nutrients from food, metabolize drugs, evade antibiotics, and communicate with the host immune system. The population genetic processes by which these traits evolve are complex, in part due to interacting ecological and evolutionary forces in the microbiome. Advances in metagenomic sequencing, coupled with bioinformatics tools and population genetic models, facilitate quantification of microbiome genetic variation and inferences about how this diversity arises, evolves, and correlates with traits of both microbes and hosts. In this review, we explore the population genetic forces (mutation, recombination, drift, and selection) that shape microbiome genetic diversity within and between hosts, as well as efforts towards predictive models that leverage microbiome genetics.

Diet-microbiome-disease: Investigating diet's influence on infectious disease resistance through alteration of the gut microbiome

Abiotic and biotic factors can affect host resistance to parasites. Host diet and host gut microbiomes are two increasingly recognized factors influencing disease resistance. In particular, recent studies demonstrate that (1) particular diets can reduce parasitism; (2) diets can alter the gut microbiome; and (3) the gut microbiome can decrease parasitism. These three separate relationships suggest the existence of indirect links through which diets reduce parasitism through an alteration of the gut microbiome. However, such links are rarely considered and even more rarely experimentally validated. This is surprising because there is increasing discussion of the therapeutic potential of diets and gut microbiomes to control infectious disease. To elucidate these potential indirect links, we review and examine studies on a wide range of animal systems commonly used in diet, microbiome, and disease research. We also examine the relative benefits and disadvantages of particular systems for the study of these indirect links and conclude that mice and insects are currently the best animal systems to test for the effect of diet-altered protective gut microbiomes on infectious disease. Focusing on these systems, we provide experimental guidelines and highlight challenges that must be overcome. Although previous studies have recommended these systems for microbiome research, here we specifically recommend these systems because of their proven relationships between diet and parasitism, between diet and the microbiome, and between the microbiome and parasite resistance. Thus, they provide a sound foundation to explore the three-way interaction between diet, the microbiome, and infectious disease.

A library of human gut bacterial isolates paired with longitudinal multiomics data enables mechanistic microbiome research

Our understanding of how the gut microbiome interacts with its human host has been restrained by limited access to longitudinal datasets to examine stability and dynamics, and by having only a few isolates to test mechanistic hypotheses. Here, we present the Broad Institute-OpenBiome Microbiome Library (BIO-ML), a comprehensive collection of 7,758 gut bacterial isolates paired with 3,632 genome sequences and longitudinal multi-omics data. We show that microbial species maintain stable population sizes within and across humans and that commonly used 'omics' survey methods are more reliable when using averages over multiple days of sampling. Variation of gut metabolites within people over time is associated with amino acid levels, and differences across people are associated with differences in bile acids. Finally, we show that genomic diversification can be used to infer eco-evolutionary dynamics and in vivo selection pressures for strains within individuals. The BIO-ML is a unique resource designed to enable hypothesis-driven microbiome research.

Pharmacomicrobiomics informs clinical pharmacogenomics

Aim: Microbiota-host-xenobiotics interactions in humans become of prime interest when clinical pharmacogenomics is to be implemented. Despite the advent of technology, information still needs to be translated into knowledge for optimum patient stratification and disease management. Material & methods: Herein, we mined metagenomic, pharmacometagenomic and pharmacomicrobiomic datasets to map microbiota-host-drugs networks. Results: Datasets were multifaceted and voluminous. Interoperability, data sparsity and scarcity remain a challenge. Mapping microbiota-host-drugs networks allowed the prediction of drug response/toxicity and modulation of the microbiota-host-drugs interplay. Conclusion: Our approach triangulated microbiota, host and drug networks revealing the need for contextual data and open science via microattribution to accelerate knowledge growth. Our findings may serve as a data storehouse for a user-friendly query system, coupled with databanks and databases.

Polar bear evolution is marked by rapid changes in gene copy number in response to dietary shift

Polar bear (Ursus maritimus) and brown bear (Ursus arctos) are recently diverged species that inhabit vastly differing habitats. Thus, analysis of the polar bear and brown bear genomes represents a unique opportunity to investigate the evolutionary mechanisms and genetic underpinnings of rapid ecological adaptation in mammals. Copy number (CN) differences in genomic regions between closely related species can underlie adaptive phenotypes and this form of genetic variation has not been explored in the context of polar bear evolution. Here, we analyzed the CN profiles of 17 polar bears, 9 brown bears, and 2 black bears (Ursus americanus). We identified an average of 318 genes per individual that showed evidence of CN variation (CNV). Nearly 200 genes displayed species-specific CN differences between polar bear and brown bear species. Principal component analysis of gene CN provides strong evidence that CNV evolved rapidly in the polar bear lineage and mainly resulted in CN loss. Olfactory receptors composed 47% of CN differentiated genes, with the majority of these genes being at lower CN in the polar bear. Additionally, we found significantly fewer copies of several genes involved in fatty acid metabolism as well as AMY1B, the salivary amylase-encoding gene in the polar bear. These results suggest that natural selection shaped patterns of CNV in response to the transition from an omnivorous to primarily carnivorous diet during polar bear evolution. Our analyses of CNV shed light on the genomic underpinnings of ecological adaptation during polar bear evolution.

An evolving view of copy number variants

Copy number variants (CNVs) are regions of the genome that vary in integer copy number. CNVs, which comprise both amplifications and deletions of DNA sequence, have been identified across all domains of life, from bacteria and archaea to plants and animals. CNVs are an important source of genetic diversity, and can drive rapid adaptive evolution and progression of heritable and somatic human diseases, such as cancer. However, despite their evolutionary importance and clinical relevance, CNVs remain understudied compared to single-nucleotide variants (SNVs). This is a consequence of the inherent difficulties in detecting CNVs at low-to-intermediate frequencies in heterogeneous populations of cells. Here, we discuss molecular methods used to detect CNVs, the limitations associated with using these techniques, and the application of new and emerging technologies that present solutions to these challenges. The goal of this short review and perspective is to highlight aspects of CNV biology that are understudied and define avenues for further research that address specific gaps in our knowledge of these complex alleles. We describe our recently developed method for CNV detection in which a fluorescent gene functions as a single-cell CNV reporter and present key findings from our evolution experiments in Saccharomyces cerevisiae. Using a CNV reporter, we found that CNVs are generated at a high rate and undergo selection with predictable dynamics across independently evolving replicate populations. Many CNVs appear to be generated through DNA replication-based processes that are mediated by the presence of short, interrupted, inverted-repeat sequences. Our results have important implications for the role of CNVs in evolutionary processes and the molecular mechanisms that underlie CNV formation. We discuss the possible extension of our method to other applications, including tracking the dynamics of CNVs in models of human tumors.

Structural variation in the gut microbiome associates with host health

Differences in the presence of even a few genes between otherwise identical bacterial strains may result in critical phenotypic differences. Here we systematically identify microbial genomic structural variants (SVs) and find them to be prevalent in the human gut microbiome across phyla and to replicate in different cohorts. SVs are enriched for CRISPR-associated and antibiotic-producing functions and depleted from housekeeping genes, suggesting that they have a role in microbial adaptation. We find multiple associations between SVs and host disease risk factors, many of which replicate in an independent cohort. Exploring genes that are clustered in the same SV, we uncover several possible mechanistic links between the microbiome and its host, including a region in Anaerostipes hadrus that encodes a composite inositol catabolism-butyrate biosynthesis pathway, the presence of which is associated with lower host metabolic disease risk. Overall, our results uncover a nascent layer of variability in the microbiome that is associated with microbial adaptation and host health.

IgA Responses to Microbiota

Various immune mechanisms are deployed in the mucosa to confront the immense diversity of resident bacteria. A substantial fraction of the commensal microbiota is coated with immunoglobulin A (IgA) antibodies, and recent findings have established the identities of these bacteria under homeostatic and disease conditions. Here we review the current understanding of IgA biology, and present a framework wherein two distinct types of humoral immunity coexist in the gastrointestinal mucosa. Homeostatic IgA responses employ a polyreactive repertoire to bind a broad but taxonomically distinct subset of microbiota. In contrast, mucosal pathogens and vaccines elicit high-affinity, T cell-dependent antibody responses. This model raises fundamental questions including how polyreactive IgA specificities are generated, how these antibodies exert effector functions, and how they exist together with other immune responses during homeostasis and disease.

Within-Species Genomic Variation and Variable Patterns of Recombination in the Tetracycline Producer Streptomyces rimosus

Streptomyces rimosus is best known as the primary source of the tetracycline class of antibiotics, most notably oxytetracycline, which have been widely used against many gram-positive and gram-negative pathogens and protozoan parasites. However, despite the medical and agricultural importance of S. rimosus, little is known of its evolutionary history and genome dynamics. In this study, we aim to elucidate the pan-genome characteristics and phylogenetic relationships of 32 S. rimosus genomes. The S. rimosus pan-genome contains more than 22,000 orthologous gene clusters, and approximately 8.8% of these genes constitutes the core genome. A large part of the accessory genome is composed of 9,646 strain-specific genes. S. rimosus exhibits an open pan-genome (decay parameter α = 0.83) and high gene diversity between strains (genomic fluidity φ = 0.12). We also observed strain-level variation in the distribution and abundance of biosynthetic gene clusters (BGCs) and that each individual S. rimosus genome has a unique repertoire of BGCs. Lastly, we observed variation in recombination, with some strains donating or receiving DNA more often than others, strains that tend to frequently recombine with specific partners, genes that often experience recombination more than others, and variable sizes of recombined DNA sequences. We conclude that the high levels of inter-strain genomic variation in S. rimosus is partly explained by differences in recombination among strains. These results have important implications on current efforts for natural drug discovery, the ecological role of strain-level variation in microbial populations, and addressing the fundamental question of why microbes have pan-genomes.

Employing toxin-antitoxin genome markers for identification of Bifidobacterium and Lactobacillus strains in human metagenomes

Recent research has indicated that in addition to the unique genotype each individual may also have a unique microbiota composition. Difference in microbiota composition may emerge from both its species and strain constituents. It is important to know the precise composition especially for the gut microbiota (GM), since it can contribute to the health assessment, personalized treatment, and disease prevention for individuals and groups (cohorts). The existing methods for species and strain composition in microbiota are not always precise and usually not so easy to use. Probiotic bacteria of the genus Bifidobacterium and Lactobacillus make an essential component of human GM. Previously we have shown that in certain Bifidobacterium and Lactobacillus species the RelBE and MazEF superfamily of toxin-antitoxin (TA) systems may be used as functional biomarkers to differentiate these groups of bacteria at the species and strain levels. We have composed a database of TA genes of these superfamily specific for all lactobacilli and bifidobacteria species with complete genome sequence and confirmed that in all Lactobacillus and Bifidobacterium species TA gene composition is species and strain specific. To analyze composition of species and strains of two bacteria genera, Bifidobacterium and Lactobacillus, in human GM we developed TAGMA (toxin antitoxin genes for metagenomes analyses) software based on polymorphism in TA genes. TAGMA was tested on gut metagenomic samples. The results of our analysis have shown that TAGMA can be used to characterize species and strains of Lactobacillus and Bifidobacterium in metagenomes.

Chronic Inflammatory Diseases: Are We Ready for Microbiota-based Dietary Intervention?

The last 15 years have witnessed the emergence of a new field of research that focuses on the roles played by the intestinal microbiota in health and disease. This research field has produced accumulating evidence indicating that dysregulation of host-microbiota interactions contributes to a range of chronic inflammatory diseases, including inflammatory bowel diseases, colorectal cancer, and metabolic syndrome. Although dysregulation of the microbiota can take complex forms, in some cases, specific bacterial species that can drive specific clinical outcomes have been identified. Among the numerous factors influencing the intestinal microbiota composition, diet is a central actor, wherein numerous dietary factors can beneficially or detrimentally impact the host/microbiota relationship. This review will highlight recent literature that has advanced understanding of microbiota-diet-disease interplay, with a central focus on the following question: Are we ready to use intestinal microbiota composition-based personalized dietary interventions to treat chronic inflammatory diseases?

1,520 reference genomes from cultivated human gut bacteria enable functional microbiome analyses

Reference genomes are essential for metagenomic analyses and functional characterization of the human gut microbiota. We present the Culturable Genome Reference (CGR), a collection of 1,520 nonredundant, high-quality draft genomes generated from >6,000 bacteria cultivated from fecal samples of healthy humans. Of the 1,520 genomes, which were chosen to cover all major bacterial phyla and genera in the human gut, 264 are not represented in existing reference genome catalogs. We show that this increase in the number of reference bacterial genomes improves the rate of mapping metagenomic sequencing reads from 50% to >70%, enabling higher-resolution descriptions of the human gut microbiome. We use the CGR genomes to annotate functions of 338 bacterial species, showing the utility of this resource for functional studies. We also carry out a pan-genome analysis of 38 important human gut species, which reveals the diversity and specificity of functional enrichment between their core and dispensable genomes.

A resource of >1,500 bacterial reference genomes sheds light on the human gut microbiome.

Evolutionary dynamics of bacteria in the gut microbiome within and across hosts

Gut microbiota are shaped by a combination of ecological and evolutionary forces. While the ecological dynamics have been extensively studied, much less is known about how species of gut bacteria evolve over time. Here, we introduce a model-based framework for quantifying evolutionary dynamics within and across hosts using a panel of metagenomic samples. We use this approach to study evolution in approximately 40 prevalent species in the human gut. Although the patterns of between-host diversity are consistent with quasi-sexual evolution and purifying selection on long timescales, we identify new genealogical signatures that challenge standard population genetic models of these processes. Within hosts, we find that genetic differences that accumulate over 6-month timescales are only rarely attributable to replacement by distantly related strains. Instead, the resident strains more commonly acquire a smaller number of putative evolutionary changes, in which nucleotide variants or gene gains or losses rapidly sweep to high frequency. By comparing these mutations with the typical between-host differences, we find evidence that some sweeps may be seeded by recombination, in addition to new mutations. However, comparisons of adult twins suggest that replacement eventually overwhelms evolution over multi-decade timescales, hinting at fundamental limits to the extent of local adaptation. Together, our results suggest that gut bacteria can evolve on human-relevant timescales, and they highlight the connections between these short-term evolutionary dynamics and longer-term evolution across hosts.

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome in the Era of the Human Microbiome: Persistent Pathogens Drive Chronic Symptoms by Interfering With Host Metabolism, Gene Expression, and Immunity

The illness ME/CFS has been repeatedly tied to infectious agents such as Epstein Barr Virus. Expanding research on the human microbiome now allows ME/CFS-associated pathogens to be studied as interacting members of human microbiome communities. Humans harbor these vast ecosystems of bacteria, viruses and fungi in nearly all tissue and blood. Most well-studied inflammatory conditions are tied to dysbiosis or imbalance of the human microbiome. While gut microbiome dysbiosis has been identified in ME/CFS, microbes and viruses outside the gut can also contribute to the illness. Pathobionts, and their associated proteins/metabolites, often control human metabolism and gene expression in a manner that pushes the body toward a state of illness. Intracellular pathogens, including many associated with ME/CFS, drive microbiome dysbiosis by directly interfering with human transcription, translation, and DNA repair processes. Molecular mimicry between host and pathogen proteins/metabolites further complicates this interference. Other human pathogens disable mitochondria or dysregulate host nervous system signaling. Antibodies and/or clonal T cells identified in patients with ME/CFS are likely activated in response to these persistent microbiome pathogens. Different human pathogens have evolved similar survival mechanisms to disable the host immune response and host metabolic pathways. The metabolic dysfunction driven by these organisms can result in similar clusters of inflammatory symptoms. ME/CFS may be driven by this pathogen-induced dysfunction, with the nature of dysbiosis and symptom presentation varying based on a patient's unique infectious and environmental history. Under such conditions, patients would benefit from treatments that support the human immune system in an effort to reverse the infectious disease process.

Single-cell copy number variant detection reveals the dynamics and diversity of adaptation

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1's proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102-104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18-21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution.

The Human Microbiota and Obesity: A Literature Systematic Review of In Vivo Models and Technical Approaches

Obesity is a noncommunicable disease that affects a considerable part of humanity. Recently, it has been recognized that gut microbiota constitutes a fundamental factor in the triggering and development of a large number of pathologies, among which obesity is one of the most related to the processes of dysbiosis. In this review, different animal model approaches, methodologies, and genome scale metabolic databases were revisited to study the gut microbiota and its relationship with metabolic disease. As a data source, PubMed for English-language published material from 1 January 2013, to 22 August 2018, were screened. Some previous studies were included if they were considered classics or highly relevant. Studies that included innovative technical approaches or different in vivo or in vitro models for the study of the relationship between gut microbiota and obesity were selected after a 16-different-keyword exhaustive search. A clear panorama of the current available options for the study of microbiota’s influence on obesity, both for animal model election and technical approaches, is presented to the researcher. All the knowledge generated from the study of the microbiota opens the possibility of considering fecal transplantation as a relevant therapeutic alternative for obesity and other metabolic disease treatment.

Prebiotics for Lactose Intolerance: Variability in Galacto-Oligosaccharide Utilization by Intestinal Lactobacillus rhamnosus

Lactose intolerance, characterized by a decrease in host lactase expression, affects approximately 75% of the world population. Galacto-oligosaccharides (GOS) are prebiotics that have been shown to alleviate symptoms of lactose intolerance and to modulate the intestinal microbiota, promoting the growth of beneficial microorganisms. We hypothesized that mechanisms of GOS utilization by intestinal bacteria are variable, impacting efficacy and response, with differences occurring at the strain level. This study aimed to determine the mechanisms by which human-derived Lactobacillus rhamnosus strains metabolize GOS. Genomic comparisons between strains revealed differences in carbohydrate utilization components, including transporters, enzymes for degradation, and transcriptional regulation, despite a high overall sequence identity (>95%) between strains. Physiological and transcriptomics analyses showed distinct differences in carbohydrate metabolism profiles and GOS utilization between strains. A putative operon responsible for GOS utilization was identified and characterized by genetic disruption of the 6-phospho-β-galactosidase, which had a critical role in GOS utilization. Our findings highlight the importance of strain-specific bacterial metabolism in the selection of probiotics and synbiotics to alleviate symptoms of gastrointestinal disorders including lactose intolerance.

AluScanCNV2: An R package for copy number variation calling and cancer risk prediction with next-generation sequencing data

The usage of next-generation sequencing (NGS) to detect copy number variation (CNV) is widely accepted in cancer research. Based on an AluScanCNV software developed by us previously, an AluScanCNV2 software has been developed in the present study as an R package that performs CNV detection from NGS data obtained through AluScan, whole-genome sequencing or other targeted NGS platforms. Its applications would include the expedited usage of somatic CNVs for cancer subtyping, and usage of recurrent germline CNVs to perform machine learning-assisted prediction of a test subject's susceptibility to cancer.

GePMI: A statistical model for personal intestinal microbiome identification

Human gut microbiomes consist of a large number of microbial genomes, which vary by diet and health conditions and from individual to individual. In the present work, we asked whether such variation or similarity could be measured and, if so, whether the results could be used for personal microbiome identification (PMI). To address this question, we herein propose a method to estimate the significance of similarity among human gut metagenomic samples based on reference-free, long k-mer features. Using these features, we find that pairwise similarities between the metagenomes of any two individuals obey a beta distribution and that a p value derived accordingly well characterizes whether two samples are from the same individual or not. We develop a computational framework called GePMI (Generating inter-individual similarity distribution for Personal Microbiome Identification) and apply it to several human gut metagenomic datasets (>300 individuals and >600 samples in total). From the results of GePMI, most of the human gut microbiomes can be identified (auROC = 0.9470, auPRC = 0.8702). Even after antibiotic treatment or fecal microbiota transplantation, the individual k-mer signature still maintains a certain specificity.

The model squid-vibrio symbiosis provides a window into the impact of strain- and species-level differences during the initial stages of symbiont engagement

Among horizontally acquired symbioses, the mechanisms underlying microbial strain- and species-level specificity remain poorly understood. Here, confocal-microscopy analyses and genetic manipulation of the squid-vibrio association revealed quantitative differences in a symbiont's capacity to interact with the host during initial engagement. Specifically, dominant strains of Vibrio fischeri, 'D-type', previously named for their dominant, single-strain colonization of the squid's bioluminescent organ, were compared with 'S-type', or 'sharing', strains, which can co-colonize the organ. These D-type strains typically: (i) formed aggregations of 100s-1000s of cells on the light-organ surface, up to 3 orders of magnitude larger than those of S-type strains; (ii) showed dominance in co-aggregation experiments, independent of inoculum size or strain proportion; (iii) perturbed larger areas of the organ's ciliated surface; and, (iv) appeared at the pore of the organ approximately 4×s more quickly than S-type strains. At least in part, genes responsible for biofilm synthesis control the hyperaggregation phenotype of a D-type strain. Other marine vibrios produced relatively small aggregations, while an array of marine Gram-positive and -negative species outside of the Vibrionaceae did not attach to the organ's surface. These studies provide insight into the impact of strain variation on early events leading to establishment of an environmentally acquired symbiosis.

Drivers of Microbiome Biodiversity: A Review of General Rules, Feces, and Ignorance

The alpha diversity of ecologic communities is affected by many biotic and abiotic drivers and, in turn, affects ecosystem functioning. Yet, patterns of alpha diversity in host-associated microbial communities (microbiomes) are poorly studied and the appropriateness of general theory is untested. Expanding diversity theory to include microbiomes is essential as diversity is a frequently cited metric of their status. Here, we review and newly analyze reports of alpha diversity for animal gut microbiomes. We demonstrate that both diet and body size affect diversity in the gut but that gut physiology (fermenter versus simple) is the most important driver. We also assess the advantages of various diversity metrics. The importance of diversity in microbiomes is often assumed but has not been tested outright. Therefore, we close by discussing how to integrate microbiomes into the field of biodiversity-ecosystem functioning to more clearly understand when and why a host supports diverse microbial communities.

Profiling microbial strains in urban environments using metagenomic sequencing data

Background

The microbial communities populating human and natural environments have been extensively characterized with shotgun metagenomics, which provides an in-depth representation of the microbial diversity within a sample. Microbes thriving in urban environments may be crucially important for human health, but have received less attention than those of other environments. Ongoing efforts started to target urban microbiomes at a large scale, but the most recent computational methods to profile these metagenomes have never been applied in this context. It is thus currently unclear whether such methods, that have proven successful at distinguishing even closely related strains in human microbiomes, are also effective in urban settings for tasks such as cultivation-free pathogen detection and microbial surveillance. Here, we aimed at a) testing the currently available metagenomic profiling tools on urban metagenomics; b) characterizing the organisms in urban environment at the resolution of single strain and c) discussing the biological insights that can be inferred from such methods.

Results

We applied three complementary methods on the 1614 metagenomes of the CAMDA 2017 challenge. With MetaMLST we identified 121 known sequence-types from 15 species of clinical relevance. For instance, we identified several Acinetobacter strains that were close to the nosocomial opportunistic pathogen A. nosocomialis. With StrainPhlAn, a generalized version of the MetaMLST approach, we inferred the phylogenetic structure of Pseudomonas stutzeri strains and suggested that the strain-level heterogeneity in environmental samples is higher than in the human microbiome. Finally, we also probed the functional potential of the different strains with PanPhlAn. We further showed that SNV-based and pangenome-based profiling provide complementary information that can be combined to investigate the evolutionary trajectories of microbes and to identify specific genetic determinants of virulence and antibiotic resistances within closely related strains.

Conclusion

We show that strain-level methods developed primarily for the analysis of human microbiomes can be effective for city-associated microbiomes. In fact, (opportunistic) pathogens can be tracked and monitored across many hundreds of urban metagenomes. However, while more effort is needed to profile strains of currently uncharacterized species, this work poses the basis for high-resolution analyses of microbiomes sampled in city and mass transportation environments.

Reviewers

This article was reviewed by Alexandra Bettina Graf, Daniel Huson and Trevor Cickovski.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0211-z) contains supplementary material, which is available to authorized users.

Complementary Mechanisms for Degradation of Inulin-Type Fructans and Arabinoxylan Oligosaccharides among Bifidobacterial Strains Suggest Bacterial Cooperation

Inulin-type fructans (ITF) and arabinoxylan oligosaccharides (AXOS) are broken down to different extents by various bifidobacterial strains present in the human colon. To date, phenotypic heterogeneity in the consumption of these complex oligosaccharides at the strain level remains poorly studied. To examine mechanistic variations in ITF and AXOS constituent preferences present in one individual, ITF and AXOS consumption by bifidobacterial strains isolated from the simulator of the human intestinal microbial ecosystem (SHIME) after inoculation with feces from one healthy individual was investigated. Among the 18 strains identified, four species-independent clusters displaying different ITF and AXOS degradation mechanisms and preferences were found. Bifidobacterium bifidum B46 showed limited growth on all substrates, whereas B. longum B24 and B. longum B18 could grow better on short-chain-length fractions of fructooligosaccharides (FOS) than on fructose. B. longum B24 could cleave arabinose substituents of AXOS extracellularly, without using the AXOS-derived xylose backbones, whereas B. longum B18 was able to consume oligosaccharides (up to xylotetraose) preferentially and consumed AXOS to a limited extent. B. adolescentis B72 degraded all fractions of FOS simultaneously, partially degraded inulin, and could use xylose backbones longer than xylotetraose extracellularly. The strain-specific degradation mechanisms were suggested to be complementary and indicated resource partitioning. Specialization in the degradation of complex carbohydrates by bifidobacteria present on the individual level could have in vivo implications for the successful implementation of ITF and AXOS, aiming at bifidogenic and/or butyrogenic effects. Finally, this work shows the importance of taking microbial strain-level differences into account in gut microbiota research.IMPORTANCE It is well known that bifidobacteria degrade undigestible complex polysaccharides, such as ITF and AXOS, in the human colon. However, this process has never been studied for strains coexisting in the same individual. To examine strain-dependent mechanistic variations in ITF and AXOS constituent preferences present in one individual, ITF and AXOS consumption by bifidobacterial strains isolated from the SHIME after inoculation with feces from one healthy individual was investigated. Among the 18 bifidobacterial strains identified, four species-independent clusters displaying different ITF and AXOS degradation mechanisms and preferences were found, indicating that such strains can coexist in the human colon. Such specialization in the degradation of complex carbohydrates by bifidobacteria present on the individual level could have in vivo implications for the successful implementation of ITF and AXOS, aiming at bifidogenic and/or butyrogenic effects.

Taxa-function robustness in microbial communities

BACKGROUND

The species composition of a microbial community is rarely fixed and often experiences fluctuations of varying degrees and at varying frequencies. These perturbations to a community's taxonomic profile naturally also alter the community's functional profile-the aggregate set of genes encoded by community members-ultimately altering the community's overall functional capacities. The magnitude of such functional changes and the specific shift that will occur in each function, however, are strongly dependent on how genes are distributed across community members' genomes. This gene distribution, in turn, is determined by the taxonomic composition of the community and would markedly differ, for example, between communities composed of species with similar genomic content vs. communities composed of species whose genomes encode relatively distinct gene sets. Combined, these observations suggest that community functional robustness to taxonomic perturbations could vary widely across communities with different compositions, yet, to date, a systematic study of the inherent link between community composition and robustness is lacking.

RESULTS

In this study, we examined how a community's taxonomic composition influences the robustness of that community's functional profile to taxonomic perturbation (here termed taxa-function robustness) across a wide array of environments. Using a novel simulation-based computational model to quantify this taxa-function robustness in host-associated and non-host-associated communities, we find notable differences in robustness between communities inhabiting different body sites, including significantly higher robustness in gut communities compared to vaginal communities that cannot be attributed solely to differences in species richness. We additionally find between-site differences in the robustness of specific functions, some of which are potentially related to site-specific environmental conditions. These taxa-function robustness differences are most strongly associated with differences in overall functional redundancy, though other aspects of gene distribution also influence taxa-function robustness in certain body environments, and are sufficient to cluster communities by environment. Further analysis revealed a correspondence between our robustness estimates and taxonomic and functional shifts observed across human-associated communities.

CONCLUSIONS

Our analysis approach revealed intriguing taxa-function robustness variation across environments and identified features of community and gene distribution that impact robustness. This approach could be further applied for estimating taxa-function robustness in novel communities and for informing the design of synthetic communities with specific robustness requirements.

BURRITO: An Interactive Multi-Omic Tool for Visualizing Taxa-Function Relationships in Microbiome Data

The abundance of both taxonomic groups and gene categories in microbiome samples can now be easily assayed via various sequencing technologies, and visualized using a variety of software tools. However, the assemblage of taxa in the microbiome and its gene content are clearly linked, and tools for visualizing the relationship between these two facets of microbiome composition and for facilitating exploratory analysis of their co-variation are lacking. Here we introduce BURRITO, a web tool for interactive visualization of microbiome multi-omic data with paired taxonomic and functional information. BURRITO simultaneously visualizes the taxonomic and functional compositions of multiple samples and dynamically highlights relationships between taxa and functions to capture the underlying structure of these data. Users can browse for taxa and functions of interest and interactively explore the share of each function attributed to each taxon across samples. BURRITO supports multiple input formats for taxonomic and metagenomic data, allows adjustment of data granularity, and can export generated visualizations as static publication-ready formatted figures. In this paper, we describe the functionality of BURRITO, and provide illustrative examples of its utility for visualizing various trends in the relationship between the composition of taxa and functions in complex microbiomes.

A microbiome case-control study of recurrent acute otitis media identified potentially protective bacterial genera

Background

Recurrent acute otitis media (rAOM, recurrent ear infection) is a common childhood disease caused by bacteria termed otopathogens, for which current treatments have limited effectiveness. Generic probiotic therapies have shown promise, but seem to lack specificity. We hypothesised that healthy children with no history of AOM carry protective commensal bacteria that could be translated into a specific probiotic therapy to break the cycle of re-infection. We characterised the nasopharyngeal microbiome of these children (controls) in comparison to children with rAOM (cases) to identify potentially protective bacteria. As some children with rAOM do not appear to carry any of the known otopathogens, we also hypothesised that characterisation of the middle ear microbiome could identify novel otopathogens, which may also guide the development of more effective therapies.

Results

Middle ear fluids, middle ear rinses and ear canal swabs from the cases and nasopharyngeal swabs from both groups underwent 16S rRNA gene sequencing. The nasopharyngeal microbiomes of cases and controls were distinct. We observed a significantly higher abundance of Corynebacterium and Dolosigranulum in the nasopharynx of controls. Alloiococcus, Staphylococcus and Turicella were abundant in the middle ear and ear canal of cases, but were uncommon in the nasopharynx of both groups. Gemella and Neisseria were characteristic of the case nasopharynx, but were not prevalent in the middle ear.

Conclusions

Corynebacterium and Dolosigranulum are characteristic of a healthy nasopharyngeal microbiome. Alloiococcus, Staphylococcus and Turicella are possible novel otopathogens, though their rarity in the nasopharynx and prevalence in the ear canal means that their role as normal aural flora cannot be ruled out. Gemella and Neisseria are unlikely to be novel otopathogens as they do not appear to colonise the middle ear in children with rAOM.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1154-3) contains supplementary material, which is available to authorized users.

Intra-species Genomic and Physiological Variability Impact Stress Resistance in Strains of Probiotic Potential

Large-scale microbiome studies have established that most of the diversity contained in the gastrointestinal tract is represented at the strain level; however, exhaustive genomic and physiological characterization of human isolates is still lacking. With increased use of probiotics as interventions for gastrointestinal disorders, genomic and functional characterization of novel microorganisms becomes essential. In this study, we explored the impact of strain-level genomic variability on bacterial physiology of two novel human Lactobacillus rhamnosus strains (AMC143 and AMC010) of probiotic potential in relation to stress resistance. The strains showed differences with known probiotic strains (L. rhamnosus GG, Lc705, and HN001) at the genomic level, including nucleotide polymorphisms, mutations in non-coding regulatory regions, and rearrangements of genomic architecture. Transcriptomics analysis revealed that gene expression profiles differed between strains when exposed to simulated gastrointestinal stresses, suggesting the presence of unique regulatory systems in each strain. In vitro physiological assays to test resistance to conditions mimicking the gut environment (acid, alkali, and bile stress) showed that growth of L. rhamnosus AMC143 was inhibited upon exposure to alkaline pH, while AMC010 and control strain LGG were unaffected. AMC143 also showed a significant survival advantage compared to the other strains upon bile exposure. Reverse transcription qPCR targeting the bile salt hydrolase gene (bsh) revealed that AMC143 expressed bsh poorly (a consequence of a deletion in the bsh promoter and truncation of bsh gene in AMC143), while AMC010 had significantly higher expression levels than AMC143 or LGG. Insertional inactivation of the bsh gene in AMC010 suggested that bsh could be detrimental to bacterial survival during bile stress. Together, these findings show that coupling of classical microbiology with functional genomics methods for the characterization of bacterial strains is critical for the development of novel probiotics, as variability between strains can dramatically alter bacterial physiology and functionality.

Strain Tracking Reveals the Determinants of Bacterial Engraftment in the Human Gut Following Fecal Microbiota Transplantation

Fecal microbiota transplantation (FMT) from healthy donor to patient is a treatment for microbiome-associated diseases. Although the success of FMT requires donor bacteria to engraft in the patient's gut, the forces governing engraftment in humans are unknown. Here we use an ongoing clinical experiment, the treatment of recurrent Clostridium difficile infection, to uncover the rules of engraftment in humans. We built a statistical model that predicts which bacterial species will engraft in a given host, and developed Strain Finder, a method to infer strain genotypes and track them over time. We find that engraftment can be predicted largely from the abundance and phylogeny of bacteria in the donor and the pre-FMT patient. Furthermore, donor strains within a species engraft in an all-or-nothing manner and previously undetected strains frequently colonize patients receiving FMT. We validated these findings for metabolic syndrome, suggesting that the same principles of engraftment extend to other indications.