Gut Microbiome Activity Contributes to Prediction of Individual Variation in Glycemic Response in Adults.

Limiting postprandial glycemic response (PPGR) is an important intervention in reducing the risk of chronic metabolic diseases and has been shown to impart significant health benefits in people with elevated levels of blood sugar. In this study, we collected gut microbiome activity data by assessing the metatranscriptome, and we measured the glycemic responses of 550 adults who consumed more than 30,000 meals, collectively, from omnivore or vegetarian/gluten-free diets. We demonstrate that gut microbiome activity, anthropometric factors, and food macronutrients modulate individual variation in glycemic response. We employ two predictive models, including a mixed-effects linear regression model (R = 0.77) and a gradient boosting machine model (Rtrain = 0.80/R2train = 0.64; Rtest = 0.64/R2test = 0.40), which demonstrate variation in PPGR between individuals when ingesting the same foods. All features in the final mixed-effects linear regression model were significant (p < 0.05) except for two features which were retained as suggestive: glutamine production pathways (p = 0.08) and the interaction between tyrosine metabolizers and carbs (p = 0.06). We introduce molecular functions as features in these two models, aggregated from microbial activity data, and show their statistically significant contributions to glycemic control. In summary, we demonstrate for the first time that metatranscriptomic activity of the gut microbiome is correlated with PPGR among adults.

Blood sugar dysregulation is caused by various underlying conditions, including type 2 diabetes, and this may lead to extended periods of hypoglycemia or hyperglycemia, which can be harmful or deadly. Clinically, glycemic control is a primary therapeutic target for dysglycemia, and food and nutrition are frequent interventions used to reduce postprandial blood glucose excursions. Primary determinants of postprandial glycemic response (PPGR) include dietary carbohydrates, individual phenotypes, and individual molecular characteristics which include the gut microbiome. Typical investigations of gut microbiomes depend on analysis methods which have poor taxonomic resolution, cannot identify certain microorganisms, and are prone to errors. In this study, each RNA molecule was identified and counted, allowing quantitative strain-level taxonomic classification and molecular pathway analysis. The primary goal of the study was to assess the impact of microbial functional activity on PPGR. The study was conducted in the USA and involved a multiethnic population of healthy adults with HbA1c levels below 6.5. All participants received 14-day omnivore diets or vegetarian/gluten-free diets, depending on nutritional requirements (omnivore diets include meat while vegetarian/gluten-free diets exclude both gluten and meat). Over this timeframe, blood glucose levels were measured in 15-min intervals, 24 h per day, capturing postprandial responses for more than 27,000 meals, including more than 18,000 provided meals which spanned a wide range of foods and macronutrient characteristics. Computational modeling demonstrated the statistical significance of all features and identified new features which may be relevant to glycemic control. These results show, for the first time, that a person's glycemic response depends on individual traits, including both their anthropometrics and their gut metatranscriptome, representing the activity of gut microbiomes.

A Robust Metatranscriptomic Technology for Population-Scale Studies of Diet, Gut Microbiome, and Human Health.

A functional readout of the gut microbiome is necessary to enable precise control of the gut microbiome's functions, which support human health and prevent or minimize a wide range of chronic diseases. Stool metatranscriptomic analysis offers a comprehensive functional view of the gut microbiome, but despite its usefulness, it has rarely been used in clinical studies due to its complexity, cost, and bioinformatic challenges. This method has also received criticism due to potential intrasample variability, rapid changes, and RNA degradation. Here, we describe a robust and automated stool metatranscriptomic method, called Viomega, which was specifically developed for population-scale studies. Viomega includes sample collection, ambient temperature sample preservation, total RNA extraction, physical removal of ribosomal RNAs (rRNAs), preparation of directional Illumina libraries, Illumina sequencing, taxonomic classification based on a database of >110,000 microbial genomes, and quantitative microbial gene expression analysis using a database of ~100 million microbial genes. We applied this method to 10,000 human stool samples and performed several small-scale studies to demonstrate sample stability and consistency. In summary, Viomega is an inexpensive, high-throughput, automated, and accurate sample-to-result stool metatranscriptomic technology platform for large-scale studies and a wide range of applications.

Gut Microbiome-Based Metagenomic Signature for Non-invasive Detection of Advanced Fibrosis in Human Nonalcoholic Fatty Liver Disease.

The presence of advanced fibrosis in nonalcoholic fatty liver disease (NAFLD) is the most important predictor of liver mortality. There are limited data on the diagnostic accuracy of gut microbiota-derived signature for predicting the presence of advanced fibrosis. In this prospective study, we characterized the gut microbiome compositions using whole-genome shotgun sequencing of DNA extracted from stool samples. This study included 86 uniquely well-characterized patients with biopsy-proven NAFLD, of which 72 had mild/moderate (stage 0-2 fibrosis) NAFLD, and 14 had advanced fibrosis (stage 3 or 4 fibrosis). We identified a set of 40 features (p < 0.006), which included 37 bacterial species that were used to construct a Random Forest classifier model to distinguish mild/moderate NAFLD from advanced fibrosis. The model had a robust diagnostic accuracy (AUC 0.936) for detecting advanced fibrosis. This study provides preliminary evidence for a fecal-microbiome-derived metagenomic signature to detect advanced fibrosis in NAFLD.

Utilization of defined microbial communities enables effective evaluation of meta-genomic assemblies.

BACKGROUND: Metagenomics is the study of the microbial genomes isolated from communities found on our bodies or in our environment. By correctly determining the relation between human health and the human associated microbial communities, novel mechanisms of health and disease can be found, thus enabling the development of novel diagnostics and therapeutics. Due to the diversity of the microbial communities, strategies developed for aligning human genomes cannot be utilized, and genomes of the microbial species in the community must be assembled de novo. However, in order to obtain the best metagenomic assemblies, it is important to choose the proper assembler. Due to the rapidly evolving nature of metagenomics, new assemblers are constantly created, and the field has not yet agreed on a standardized process. Furthermore, the truth sets used to compare these methods are either too simple (computationally derived diverse communities) or complex (microbial communities of unknown composition), yielding results that are hard to interpret. In this analysis, we interrogate the strengths and weaknesses of five popular assemblers through the use of defined biological samples of known genomic composition and abundance. We assessed the performance of each assembler on their ability to reassemble genomes, call taxonomic abundances, and recreate open reading frames (ORFs). RESULTS: We tested five metagenomic assemblers: Omega, metaSPAdes, IDBA-UD, metaVelvet and MEGAHIT on known and synthetic metagenomic data sets. MetaSPAdes excelled in diverse sets, IDBA-UD performed well all around, metaVelvet had high accuracy in high abundance organisms, and MEGAHIT was able to accurately differentiate similar organisms within a community. At the ORF level, metaSPAdes and MEGAHIT had the least number of missing ORFs within diverse and similar communities respectively. CONCLUSIONS: Depending on the metagenomics question asked, the correct assembler for the task at hand will differ. It is important to choose the appropriate assembler, and thus clearly define the biological problem of an experiment, as different assemblers will give different answers to the same question.

A robust ambient temperature collection and stabilization strategy: Enabling worldwide functional studies of the human microbiome.

As reports on possible associations between microbes and the host increase in number, more meaningful interpretations of this information require an ability to compare data sets across studies. This is dependent upon standardization of workflows to ensure comparability both within and between studies. Here we propose the standard use of an alternate collection and stabilization method that would facilitate such comparisons. The DNA Genotek OMNIgene∙Gut Stool Microbiome Kit was compared to the currently accepted community standard of freezing to store human stool samples prior to whole genome sequencing (WGS) for microbiome studies. This stabilization and collection device allows for ambient temperature storage, automation, and ease of shipping/transfer of samples. The device permitted the same data reproducibility as with frozen samples, and yielded higher recovery of nucleic acids. Collection and stabilization of stool microbiome samples with the DNA Genotek collection device, combined with our extraction and WGS, provides a robust, reproducible workflow that enables standardized global collection, storage, and analysis of stool for microbiome studies.

Microbial diversity in individuals and their household contacts following typical antibiotic courses.

BACKGROUND: Antibiotics are a mainstay of treatment for bacterial infections worldwide, yet the effects of typical antibiotic prescriptions on human indigenous microbiota have not been thoroughly evaluated. We examined the effects of the two most commonly prescribed antibiotics (amoxicillin and azithromycin) in the USA to discern whether short-term antibiotic courses may have prolonged effects on human microbiota. RESULTS: We sampled the feces, saliva, and skin specimens from a cohort of unrelated, cohabitating individuals over 6 months. An individual in each household was given an antibiotic, and the other a placebo to discern antibiotic impacts on microbiota, as well as determine whether antibiotic use might reshape the microbiota of each household. We observed household-specific patterns of microbiota on each body surface, which persevered despite antibiotic perturbations. While the gut microbiota within an individual became more dissimilar over time, there was no evidence that the use of antibiotics accelerated this process when compared to household members. There was a significant change in microbiota diversity in the gut and mouth in response to antibiotics, but analogous patterns were not observed on the skin. Those who received 7 days of amoxicillin generally had greater reductions in diversity compared to those who received 3 days, in contrast to those who received azithromycin. CONCLUSIONS: As few as 3 days of treatment with the most commonly prescribed antibiotics can result in sustained reductions in microbiota diversity, which could have implications for the maintenance of human health and resilience to disease.

Library preparation methodology can influence genomic and functional predictions in human microbiome research.

Observations from human microbiome studies are often conflicting or inconclusive. Many factors likely contribute to these issues including small cohort sizes, sample collection, and handling and processing differences. The field of microbiome research is moving from 16S rDNA gene sequencing to a more comprehensive genomic and functional representation through whole-genome sequencing (WGS) of complete communities. Here we performed quantitative and qualitative analyses comparing WGS metagenomic data from human stool specimens using the Illumina Nextera XT and Illumina TruSeq DNA PCR-free kits, and the KAPA Biosystems Hyper Prep PCR and PCR-free systems. Significant differences in taxonomy are observed among the four different next-generation sequencing library preparations using a DNA mock community and a cell control of known concentration. We also revealed biases in error profiles, duplication rates, and loss of reads representing organisms that have a high %G+C content that can significantly impact results. As with all methods, the use of benchmarking controls has revealed critical differences among methods that impact sequencing results and later would impact study interpretation. We recommend that the community adopt PCR-free-based approaches to reduce PCR bias that affects calculations of abundance and to improve assemblies for accurate taxonomic assignment. Furthermore, the inclusion of a known-input cell spike-in control provides accurate quantitation of organisms in clinical samples.

Isolation of a significant fraction of non-phototroph diversity from a desert Biological Soil Crust.

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

A rapid molecular approach for chromosomal phasing.

Determining the chromosomal phase of pairs of sequence variants - the arrangement of specific alleles as haplotypes - is a routine challenge in molecular genetics. Here we describe Drop-Phase, a molecular method for quickly ascertaining the phase of pairs of DNA sequence variants (separated by 1-200 kb) without cloning or manual single-molecule dilution. In each Drop-Phase reaction, genomic DNA segments are isolated in tens of thousands of nanoliter-sized droplets together with allele-specific fluorescence probes, in a single reaction well. Physically linked alleles partition into the same droplets, revealing their chromosomal phase in the co-distribution of fluorophores across droplets. We demonstrated the accuracy of this method by phasing members of trios (revealing 100% concordance with inheritance information), and demonstrate a common clinical application by phasing CFTR alleles at genomic distances of 11-116 kb in the genomes of cystic fibrosis patients. Drop-Phase is rapid (requiring less than 4 hours), scalable (to hundreds of samples), and effective at long genomic distances (200 kb).

Multiplexed target detection using DNA-binding dye chemistry in droplet digital PCR.

Two years ago, we described the first droplet digital PCR (ddPCR) system aimed at empowering all researchers with a tool that removes the substantial uncertainties associated with using the analogue standard, quantitative real-time PCR (qPCR). This system enabled TaqMan hydrolysis probe-based assays for the absolute quantification of nucleic acids. Due to significant advancements in droplet chemistry and buoyed by the multiple benefits associated with dye-based target detection, we have created a "second generation" ddPCR system compatible with both TaqMan-probe and DNA-binding dye detection chemistries. Herein, we describe the operating characteristics of DNA-binding dye based ddPCR and offer a side-by-side comparison to TaqMan probe detection. By partitioning each sample prior to thermal cycling, we demonstrate that it is now possible to use a DNA-binding dye for the quantification of multiple target species from a single reaction. The increased resolution associated with partitioning also made it possible to visualize and account for signals arising from nonspecific amplification products. We expect that the ability to combine the precision of ddPCR with both DNA-binding dye and TaqMan probe detection chemistries will further enable the research community to answer complex and diverse genetic questions.

The ability of flux balance analysis to predict evolution of central metabolism scales with the initial distance to the optimum.

The most powerful genome-scale framework to model metabolism, flux balance analysis (FBA), is an evolutionary optimality model. It hypothesizes selection upon a proposed optimality criterion in order to predict the set of internal fluxes that would maximize fitness. Here we present a direct test of the optimality assumption underlying FBA by comparing the central metabolic fluxes predicted by multiple criteria to changes measurable by a (13)C-labeling method for experimentally-evolved strains. We considered datasets for three Escherichia coli evolution experiments that varied in their length, consistency of environment, and initial optimality. For ten populations that were evolved for 50,000 generations in glucose minimal medium, we observed modest changes in relative fluxes that led to small, but significant decreases in optimality and increased the distance to the predicted optimal flux distribution. In contrast, seven populations evolved on the poor substrate lactate for 900 generations collectively became more optimal and had flux distributions that moved toward predictions. For three pairs of central metabolic knockouts evolved on glucose for 600-800 generations, there was a balance between cases where optimality and flux patterns moved toward or away from FBA predictions. Despite this variation in predictability of changes in central metabolism, two generalities emerged. First, improved growth largely derived from evolved increases in the rate of substrate use. Second, FBA predictions bore out well for the two experiments initiated with ancestors with relatively sub-optimal yield, whereas those begun already quite optimal tended to move somewhat away from predictions. These findings suggest that the tradeoff between rate and yield is surprisingly modest. The observed positive correlation between rate and yield when adaptation initiated further from the optimum resulted in the ability of FBA to use stoichiometric constraints to predict the evolution of metabolism despite selection for rate.

Dynamic cyanobacterial response to hydration and dehydration in a desert biological soil crust.

Biological soil crusts (BSCs) cover extensive portions of the earth's deserts. In order to survive desiccation cycles and utilize short periods of activity during infrequent precipitation, crust microorganisms must rely on the unique capabilities of vegetative cells to enter a dormant state and be poised for rapid resuscitation upon wetting. To elucidate the key events involved in the exit from dormancy, we performed a wetting experiment of a BSC and followed the response of the dominant cyanobacterium, Microcoleus vaginatus, in situ using a whole-genome transcriptional time course that included two diel cycles. Immediate, but transient, induction of DNA repair and regulatory genes signaled the hydration event. Recovery of photosynthesis occurred within 1 h, accompanied by upregulation of anabolic pathways. Onset of desiccation was characterized by the induction of genes for oxidative and photo-oxidative stress responses, osmotic stress response and the synthesis of C and N storage polymers. Early expression of genes for the production of exopolysaccharides, additional storage molecules and genes for membrane unsaturation occurred before drying and hints at preparedness for desiccation. We also observed signatures of preparation for future precipitation, notably the expression of genes for anaplerotic reactions in drying crusts, and the stable maintenance of mRNA through dormancy. These data shed light on possible synchronization between this cyanobacterium and its environment, and provides key mechanistic insights into its metabolism in situ that may be used to predict its response to climate, and or, land-use driven perturbations.

Deep sequencing of the oral microbiome reveals signatures of periodontal disease.

The oral microbiome, the complex ecosystem of microbes inhabiting the human mouth, harbors several thousands of bacterial types. The proliferation of pathogenic bacteria within the mouth gives rise to periodontitis, an inflammatory disease known to also constitute a risk factor for cardiovascular disease. While much is known about individual species associated with pathogenesis, the system-level mechanisms underlying the transition from health to disease are still poorly understood. Through the sequencing of the 16S rRNA gene and of whole community DNA we provide a glimpse at the global genetic, metabolic, and ecological changes associated with periodontitis in 15 subgingival plaque samples, four from each of two periodontitis patients, and the remaining samples from three healthy individuals. We also demonstrate the power of whole-metagenome sequencing approaches in characterizing the genomes of key players in the oral microbiome, including an unculturable TM7 organism. We reveal the disease microbiome to be enriched in virulence factors, and adapted to a parasitic lifestyle that takes advantage of the disrupted host homeostasis. Furthermore, diseased samples share a common structure that was not found in completely healthy samples, suggesting that the disease state may occupy a narrow region within the space of possible configurations of the oral microbiome. Our pilot study demonstrates the power of high-throughput sequencing as a tool for understanding the role of the oral microbiome in periodontal disease. Despite a modest level of sequencing (~2 lanes Illumina 76 bp PE) and high human DNA contamination (up to ~90%) we were able to partially reconstruct several oral microbes and to preliminarily characterize some systems-level differences between the healthy and diseased oral microbiomes.

Detection of transcriptional triggers in the dynamics of microbial growth: application to the respiratorily versatile bacterium Shewanella oneidensis.

The capacity of microorganisms to respond to variable external conditions requires a coordination of environment-sensing mechanisms and decision-making regulatory circuits. Here, we seek to understand the interplay between these two processes by combining high-throughput measurement of time-dependent mRNA profiles with a novel computational approach that searches for key genetic triggers of transcriptional changes. Our approach helped us understand the regulatory strategies of a respiratorily versatile bacterium with promising bioenergy and bioremediation applications, Shewanella oneidensis, in minimal and rich media. By comparing expression profiles across these two conditions, we unveiled components of the transcriptional program that depend mainly on the growth phase. Conversely, by integrating our time-dependent data with a previously available large compendium of static perturbation responses, we identified transcriptional changes that cannot be explained solely by internal network dynamics, but are rather triggered by specific genes acting as key mediators of an environment-dependent response. These transcriptional triggers include known and novel regulators that respond to carbon, nitrogen and oxygen limitation. Our analysis suggests a sequence of physiological responses, including a coupling between nitrogen depletion and glycogen storage, partially recapitulated through dynamic flux balance analysis, and experimentally confirmed by metabolite measurements. Our approach is broadly applicable to other systems.

Host-pathogen interactome mapping for HTLV-1 and -2 retroviruses.

BACKGROUND: Human T-cell leukemia virus type 1 (HTLV-1) and type 2 both target T lymphocytes, yet induce radically different phenotypic outcomes. HTLV-1 is a causative agent of Adult T-cell leukemia (ATL), whereas HTLV-2, highly similar to HTLV-1, causes no known overt disease. HTLV gene products are engaged in a dynamic struggle of activating and antagonistic interactions with host cells. Investigations focused on one or a few genes have identified several human factors interacting with HTLV viral proteins. Most of the available interaction data concern the highly investigated HTLV-1 Tax protein. Identifying shared and distinct host-pathogen protein interaction profiles for these two viruses would enlighten how they exploit distinctive or common strategies to subvert cellular pathways toward disease progression. RESULTS: We employ a scalable methodology for the systematic mapping and comparison of pathogen-host protein interactions that includes stringent yeast two-hybrid screening and systematic retest, as well as two independent validations through an additional protein interaction detection method and a functional transactivation assay. The final data set contained 166 interactions between 10 viral proteins and 122 human proteins. Among the 166 interactions identified, 87 and 79 involved HTLV-1 and HTLV-2 -encoded proteins, respectively. Targets for HTLV-1 and HTLV-2 proteins implicate a diverse set of cellular processes including the ubiquitin-proteasome system, the apoptosis, different cancer pathways and the Notch signaling pathway. CONCLUSIONS: This study constitutes a first pass, with homogeneous data, at comparative analysis of host targets for HTLV-1 and -2 retroviruses, complements currently existing data for formulation of systems biology models of retroviral induced diseases and presents new insights on biological pathways involved in retroviral infection.

Ecosystems biology of microbial metabolism.

The metabolic capabilities of many environmentally and medically important microbes can be quantitatively explored using systems biology approaches to metabolic networks. Yet, as we learn more about the complex microbe-microbe and microbe-environment interactions in microbial communities, it is important to understand whether and how system-level approaches can be extended to the ecosystem level. Here we summarize recent work that addresses these challenges at multiple scales, starting from two-species natural and synthetic ecology models, up to biosphere-level approaches. Among the many fascinating open challenges in this field is whether the integration of high throughput sequencing methods and mathematical models will help us capture emerging principles of ecosystem-level metabolic organization and evolution.

COMBREX: a project to accelerate the functional annotation of prokaryotic genomes.

COMBREX (http://combrex.bu.edu) is a project to increase the speed of the functional annotation of new bacterial and archaeal genomes. It consists of a database of functional predictions produced by computational biologists and a mechanism for experimental biochemists to bid for the validation of those predictions. Small grants are available to support successful bids.

Environments that induce synthetic microbial ecosystems.

Interactions between microbial species are sometimes mediated by the exchange of small molecules, secreted by one species and metabolized by another. Both one-way (commensal) and two-way (mutualistic) interactions may contribute to complex networks of interdependencies. Understanding these interactions constitutes an open challenge in microbial ecology, with applications ranging from the human microbiome to environmental sustainability. In parallel to natural communities, it is possible to explore interactions in artificial microbial ecosystems, e.g. pairs of genetically engineered mutualistic strains. Here we computationally generate artificial microbial ecosystems without re-engineering the microbes themselves, but rather by predicting their growth on appropriately designed media. We use genome-scale stoichiometric models of metabolism to identify media that can sustain growth for a pair of species, but fail to do so for one or both individual species, thereby inducing putative symbiotic interactions. We first tested our approach on two previously studied mutualistic pairs, and on a pair of highly curated model organisms, showing that our algorithms successfully recapitulate known interactions, robustly predict new ones, and provide novel insight on exchanged molecules. We then applied our method to all possible pairs of seven microbial species, and found that it is always possible to identify putative media that induce commensalism or mutualism. Our analysis also suggests that symbiotic interactions may arise more readily through environmental fluctuations than genetic modifications. We envision that our approach will help generate microbe-microbe interaction maps useful for understanding microbial consortia dynamics and evolution, and for exploring the full potential of natural metabolic pathways for metabolic engineering applications.