Vikodak - A Modular Framework for Inferring Functional Potential of Microbial Communities from 16S Metagenomic Datasets

Background

The overall metabolic/functional potential of any given environmental niche is a function of the sum total of genes/proteins/enzymes that are encoded and expressed by various interacting microbes residing in that niche. Consequently, prior (collated) information pertaining to genes, enzymes encoded by the resident microbes can aid in indirectly (re)constructing/ inferring the metabolic/ functional potential of a given microbial community (given its taxonomic abundance profile). In this study, we present Vikodak—a multi-modular package that is based on the above assumption and automates inferring and/ or comparing the functional characteristics of an environment using taxonomic abundance generated from one or more environmental sample datasets. With the underlying assumptions of co-metabolism and independent contributions of different microbes in a community, a concerted effort has been made to accommodate microbial co-existence patterns in various modules incorporated in Vikodak.

Results

Validation experiments on over 1400 metagenomic samples have confirmed the utility of Vikodak in (a) deciphering enzyme abundance profiles of any KEGG metabolic pathway, (b) functional resolution of distinct metagenomic environments, (c) inferring patterns of functional interaction between resident microbes, and (d) automating statistical comparison of functional features of studied microbiomes. Novel features incorporated in Vikodak also facilitate automatic removal of false positives and spurious functional predictions.

Conclusions

With novel provisions for comprehensive functional analysis, inclusion of microbial co-existence pattern based algorithms, automated inter-environment comparisons; in-depth analysis of individual metabolic pathways and greater flexibilities at the user end, Vikodak is expected to be an important value addition to the family of existing tools for 16S based function prediction.

Availability and Implementation

A web implementation of Vikodak can be publicly accessed at: http://metagenomics.atc.tcs.com/vikodak. This web service is freely available for all categories of users (academic as well as commercial).     

Binpairs: Utilization of Illumina Paired-End Information for Improving Efficiency of Taxonomic Binning of Metagenomic Sequences

Motivation

Paired-end sequencing protocols, offered by next generation sequencing (NGS) platforms like Illumia, generate a pair of reads for every DNA fragment in a sample. Although this protocol has been utilized for several metagenomics studies, most taxonomic binning approaches classify each of the reads (forming a pair), independently. The present work explores some simple but effective strategies of utilizing pairing-information of Illumina short reads for improving the accuracy of taxonomic binning of metagenomic datasets. The strategies proposed can be used in conjunction with all genres of existing binning methods.

Results

Validation results suggest that employment of these “Binpairs” strategies can provide significant improvements in the binning outcome. The quality of the taxonomic assignments thus obtained are often comparable to those that can only be achieved with relatively longer reads obtained using other NGS platforms (such as Roche).

Availability

An implementation of the proposed strategies of utilizing pairing information is freely available for academic users at https://metagenomics.atc.tcs.com/binning/binpairs.     

New challenges in the study of the evolution of wild animals and their gut microbiome

In this viewpoint, by reviewing the recent findings on wild animals and their gut microbiomes, we found some potential new insights and challenges in the study of the evolution of wild animals and their gut microbiome. We suggested that wild animal gut microbiomes may come from microbiomes in the animals'' living habitats along with animals'' special behavior, and that the study of long‐term changes in gut microbiomes should consider both habitat and special behaviors. Also, host behavior would facilitate the gut microbiome transmission between individuals. We suggested that research should integrate the evolutionary history and physiological systems of wild animals to understand the evolution of animals and their gut microbiomes. Finally, we proposed the Noncultured‐Cultured‐Fermentation‐Model Animal pipeline to determine the function (diet digestion, physiology, and behavior) of these target strains in the wild animal gut.     

Insights into horizontal acquisition patterns of dormancy and reactivation regulon genes in mycobacterial species using a partitioning-based framework

Horizontal Gene Transfer (HGT) events, initially thought to be rare in Mycobacterium tuberculosis, have recently been shown to be involved in the acquisition of virulence operons in M. tuberculosis. We have developed a new partitioning framework based HGT prediction algorithm, called Grid3M, and applied the same for the prediction of HGTs in Mycobacteria. Validation and testing using simulated and real microbial genomes indicated better performance of Grid3M as compared with other widely used HGT prediction methods. Specific analysis of the genes belonging to dormancy/reactivation regulons across 14 mycobacterial genomes indicated that horizontal acquisition is specifically restricted to important accessory proteins. The results also revealed Burkholderia species to be a probable source of HGT genes belonging to these regulons. The current study provides a basis for similar analyses investigating the functional/evolutionary aspects of HGT genes in other pathogens. A database of Grid3M predicted HGTs in completely sequenced genomes is available at https://metagenomics.atc.tcs.com/Grid3M/ .     

Mammalian gut metabolomes mirror microbiome composition and host phylogeny

In the past decade, studies on the mammalian gut microbiome have revealed that different animal species have distinct gut microbial compositions. The functional ramifications of this variation in microbial composition remain unclear: do these taxonomic differences indicate microbial adaptations to host-specific functionality, or are these diverse microbial communities essentially functionally redundant, as has been indicated by previous metagenomics studies? Here, we examine the metabolic content of mammalian gut microbiomes as a direct window into ecosystem function, using an untargeted metabolomics platform to analyze 101 fecal samples from a range of 25 exotic mammalian species in collaboration with a zoological center. We find that mammalian metabolomes are chemically diverse and strongly linked to microbiome composition, and that metabolome composition is further correlated to the phylogeny of the mammalian host. Specific metabolites enriched in different animal species included modified and degraded host and dietary compounds such as bile acids and triterpenoids, as well as fermentation products such as lactate and short-chain fatty acids. Our results suggest that differences in microbial taxonomic composition are indeed translated to host-specific metabolism, indicating that taxonomically distant microbiomes are more functionally diverse than redundant.Subject terms: Metabolomics, Microbiome, Microbial ecology     

Analysis of Stomach and Gut Microbiomes of the Eastern Oyster (Crassostrea virginica) from Coastal Louisiana,USA

We used high throughput pyrosequencing to characterize stomach and gut content microbiomes of Crassostrea virginica, the Easter oyster, obtained from two sites, one in Barataria Bay (Hackberry Bay) and the other in Terrebonne Bay (Lake Caillou), Louisiana, USA. Stomach microbiomes in oysters from Hackberry Bay were overwhelmingly dominated by Mollicutes most closely related to Mycoplasma; a more rich community dominated by Planctomyctes occurred in Lake Caillou oyster stomachs. Gut communities for oysters from both sites differed from stomach communities, and harbored a relatively diverse assemblage of phylotypes. Phylotypes most closely related to Shewanella and a Chloroflexi strain dominated the Lake Caillou and Hackberry Bay gut microbiota, respectively. While many members of the stomach and gut microbiomes appeared to be transients or opportunists, a putative core microbiome was identified based on phylotypes that occurred in all stomach or gut samples only. The putative core stomach microbiome comprised 5 OTUs in 3 phyla, while the putative core gut microbiome contained 44 OTUs in 12 phyla. These results collectively revealed novel microbial communities within the oyster digestive system, the functions of the oyster microbiome are largely unknown. A comparison of microbiomes from Louisiana oysters with bacterial communities reported for other marine invertebrates and fish indicated that molluscan microbiomes were more similar to each other than to microbiomes of polychaetes, decapods and fish.     

Exploring Vertical Transmission of Bifidobacteria from Mother to Child

Passage through the birth canal and consequent exposure to the mother''s microbiota is considered to represent the initiating event for microbial colonization of the gastrointestinal tract of the newborn. However, a precise evaluation of such suspected vertical microbiota transmission has yet to be performed. Here, we evaluated the microbiomes of four sample sets, each consisting of a mother''s fecal and milk samples and the corresponding infant''s fecal sample, by means of amplicon-based profiling supported by shotgun metagenomics data for two key samples. Notably, targeted genome reconstruction from microbiome data revealed vertical transmission of a Bifidobacterium breve strain and a Bifidobacterium longum subsp. longum strain from mother to infant, a notion confirmed by strain isolation and genome sequencing. Furthermore, PCR analyses targeting unique genes from these two strains highlighted their persistence in the infant gut at 6 months. Thus, this study demonstrates the existence of specific bifidobacterial strains that are common to mother and child and thus indicative of vertical transmission and that are maintained in the infant for at least relatively short time spans.     

Temporal variation selects for diet–microbe co-metabolic traits in the gut of Gorilla spp

Although the critical role that our gastrointestinal microbes play in host physiology is now well established, we know little about the factors that influenced the evolution of primate gut microbiomes. To further understand current gut microbiome configurations and diet–microbe co-metabolic fingerprints in primates, from an evolutionary perspective, we characterized fecal bacterial communities and metabolomic profiles in 228 fecal samples of lowland and mountain gorillas (G. g. gorilla and G. b. beringei, respectively), our closest evolutionary relatives after chimpanzees. Our results demonstrate that the gut microbiomes and metabolomes of these two species exhibit significantly different patterns. This is supported by increased abundance of metabolites and bacterial taxa associated with fiber metabolism in mountain gorillas, and enrichment of markers associated with simple sugar, lipid and sterol turnover in the lowland species. However, longitudinal sampling shows that both species'' microbiomes and metabolomes converge when hosts face similar dietary constraints, associated with low fruit availability in their habitats. By showing differences and convergence of diet–microbe co-metabolic fingerprints in two geographically isolated primate species, under specific dietary stimuli, we suggest that dietary constraints triggered during their adaptive radiation were potential factors behind the species-specific microbiome patterns observed in primates today.     

The effect of captivity on the primate gut microbiome varies with host dietary niche

Despite careful attention to animal nutrition and wellbeing, gastrointestinal distress remains relatively common in captive non‐human primates (NHPs), particularly dietary specialists such as folivores. These patterns may be a result of marked dietary differences between captive and wild settings and associated impacts on the gut microbiome. However, given that most existing studies target NHP dietary specialists, it is unclear if captive environments have distinct impacts on the gut microbiome of NHPs with different dietary niches. To begin to examine this question, we used 16S ribosomal RNA gene amplicon sequences to compare the gut microbiomes of five NHP genera categorized either as folivores (Alouatta, Colobus) or non‐folivores (Cercopithecus, Gorilla, Pan) sampled both in captivity and in the wild. Though captivity affected the gut microbiomes of all NHPs in this study, the effects were largest in folivorous NHPs. Shifts in gut microbial diversity and in the relative abundances of fiber‐degrading microbial taxa suggest that these findings are driven by marked dietary shifts for folivorous NHPs in captive settings. We propose that zoos and other captive care institutions consider including more natural browse in folivorous NHP diets and regularly bank fecal samples to further explore the relationship between NHP diet, the gut microbiome, and health outcomes.     

Primate microbiomes over time: Longitudinal answers to standing questions in microbiome research

To date, most insights into the processes shaping vertebrate gut microbiomes have emerged from studies with cross‐sectional designs. While this approach has been valuable, emerging time series analyses on vertebrate gut microbiomes show that gut microbial composition can change rapidly from 1 day to the next, with consequences for host physical functioning, health, and fitness. Hence, the next frontier of microbiome research will require longitudinal perspectives. Here we argue that primatologists, with their traditional focus on tracking the lives of individual animals and familiarity with longitudinal fecal sampling, are well positioned to conduct research at the forefront of gut microbiome dynamics. We begin by reviewing some of the most important ecological processes governing microbiome change over time, and briefly summarizing statistical challenges and approaches to microbiome time series analysis. We then introduce five questions of general interest to microbiome science where we think field‐based primate studies are especially well positioned to fill major gaps: (a) Do early life events shape gut microbiome composition in adulthood? (b) Do shifting social landscapes cause gut microbial change? (c) Are gut microbiome phenotypes heritable across variable environments? (d) Does the gut microbiome show signs of host aging? And (e) do gut microbiome composition and dynamics predict host health and fitness? For all of these questions, we highlight areas where primatologists are uniquely positioned to make substantial contributions. We review preliminary evidence, discuss possible study designs, and suggest future directions.     

BEEM-Static: Accurate inference of ecological interactions from cross-sectional microbiome data

The structure and function of diverse microbial communities is underpinned by ecological interactions that remain uncharacterized. With rapid adoption of next-generation sequencing for studying microbiomes, data-driven inference of microbial interactions based on abundance correlations is widely used, but with the drawback that ecological interpretations may not be possible. Leveraging cross-sectional microbiome datasets for unravelling ecological structure in a scalable manner thus remains an open problem. We present an expectation-maximization algorithm (BEEM-Static) that can be applied to cross-sectional datasets to infer interaction networks based on an ecological model (generalized Lotka-Volterra). The method exhibits robustness to violations in model assumptions by using statistical filters to identify and remove corresponding samples. Benchmarking against 10 state-of-the-art correlation based methods showed that BEEM-Static can infer presence and directionality of ecological interactions even with relative abundance data (AUC-ROC>0.85), a task that other methods struggle with (AUC-ROC<0.63). In addition, BEEM-Static can tolerate a high fraction of samples (up to 40%) being not at steady state or coming from an alternate model. Applying BEEM-Static to a large public dataset of human gut microbiomes (n = 4,617) identified multiple stable equilibria that better reflect ecological enterotypes with distinct carrying capacities and interactions for key species.ConclusionBEEM-Static provides new opportunities for mining ecologically interpretable interactions and systems insights from the growing corpus of microbiome data.     

Abdominal microbial communities in ants depend on colony membership rather than caste and are linked to colony productivity

Gut bacteria aid their host in digestion and pathogen defense, and bacterial communities that differ in diversity or composition may vary in their ability to do so. Typically, the gut microbiomes of animals living in social groups converge as members share a nest environment and frequently interact. Social insect colonies, however, consist of individuals that differ in age, physiology, and behavior, traits that could affect gut communities or that expose the host to different bacteria, potentially leading to variation in the gut microbiome within colonies. Here we asked whether bacterial communities in the abdomen of Temnothorax nylanderi ants, composed largely of the gut microbiome, differ between different reproductive and behavioral castes. We compared microbiomes of queens, newly eclosed workers, brood carers, and foragers by high‐throughput 16S rRNA sequencing. Additionally, we sampled individuals from the same colonies twice, in the field and after 2 months of laboratory housing. To disentangle the effects of laboratory environment and season on microbial communities, additional colonies were collected at the same location after 2 months. There were no large differences between ant castes, although queens harbored more diverse microbial communities than workers. Instead, we found effects of colony, environment, and season on the abdominal microbiome. Interestingly, colonies with more diverse communities had produced more brood. Moreover, the queens' microbiome composition was linked to egg production. Although long‐term coevolution between social insects and gut bacteria has been repeatedly evidenced, our study is the first to find associations between abdominal microbiome characteristics and colony productivity in social insects.     

Human reproductive system microbiomes exhibited significantly different heterogeneity scaling with gut microbiome,but the intra-system scaling is invariant

Maintaining sexual reproduction in a highly competitive world is still one of the major mysteries of biology given the apparently high efficiency of asexual reproduction. Co-evolutionary theories such as the Red Queen hypothesis would suggest that the microbiomes in human reproductive systems, specifically the microbiomes contained in semen and vaginal fluids, should reach some level of homogeneity thanks to arguably the most conspicuous microbiome transmission between two sexes. The long-term sexual coevolution should favor the dynamic homogeneity or stability, which should also be beneficial for sexual reproduction such as sperm survival or fertilization on physiological/ecological time scale. We present a piece of quantitative evidence in the form of microbial community spatial heterogeneity to support the stability notion by analyzing three big datasets of the human vaginal, semen and gut microbiome. Methodologically, we applied a recent community-level extension to the classic Taylor's power law, which reached the rare status of ecological law and has found applications beyond biology. Both ecological and evolutionary theories, such as hologenome/holobiont and Red Queen, as well as consideration of first principles, would predict that microbiome transmissions between two sexes should have homogenizing effects on the composition and stability of the microbiomes in human reproductive systems, and therefore have similar variance structures. This is supported by the finding that the power law analysis revealed human vaginal and semen microbiomes exhibited the same scaling parameter size in their community spatial (inter-individual) heterogeneities, while the both exhibited significantly different heterogeneity scaling parameter size with the human gut microbiome.     

Development of the Honey Bee Gut Microbiome throughout the Queen-Rearing Process

The European honey bee (Apis mellifera) is used extensively to produce hive products and for crop pollination, but pervasive concerns about colony health and population decline have sparked an interest in the microbial communities that are associated with these important insects. Currently, only the microbiome of workers has been characterized, while little to nothing is known about the bacterial communities that are associated with queens, even though their health and proper function are central to colony productivity. Here, we provide a large-scale analysis of the gut microbiome of honey bee queens during their developmental trajectory and through the multiple colonies that host them as part of modern queen-rearing practices. We found that queen microbiomes underwent a dramatic shift in size and composition as they aged and encountered different worker populations and colony environments. Queen microbiomes were dominated by enteric bacteria in early life but were comprised primarily of alphaproteobacteria at maturity. Furthermore, queen gut microbiomes did not reflect those of the workers who tended them and, indeed, they lacked many of the bacteria that are considered vital to workers. While worker gut microbiotas were consistent across the unrelated colony populations sampled, the microbiotas of the related queens were highly variable. Bacterial communities in mature queen guts were similar in size to those of mature workers and were characterized by dominant and specific alphaproteobacterial strains known to be associated with worker hypopharyngeal glands. Our results suggest a model in which queen guts are colonized by bacteria from workers'' glands, in contrast to routes of maternal inoculation for other animal microbiomes.     

Stability of the gorilla microbiome despite simian immunodeficiency virus infection

Simian immunodeficiency viruses (SIVs) have been discovered in over 45 primate species; however, the pathogenic potential of most SIV strains remains unknown due to difficulties inherent in observing wild populations. Because those SIV infections that are pathogenic have been shown to induce changes in the host's gut microbiome, monitoring the microbiota present in faecal samples can provide a noninvasive means for studying the effects of SIV infection on the health of wild‐living primates. Here, we examine the effects of SIVgor, a close relative of SIVcpz of chimpanzees and HIV‐1 of humans, on the gut bacterial communities residing within wild gorillas, revealing that gorilla gut microbiomes are exceptionally robust to SIV infection. In contrast to the microbiomes of HIV‐1‐infected humans and SIVcpz‐infected chimpanzees, SIVgor‐infected gorilla microbiomes exhibit neither rises in the frequencies of opportunistic pathogens nor elevated rates of microbial turnover within individual hosts. Regardless of SIV infection status, gorilla microbiomes assort into enterotypes, one of which is compositionally analogous to those identified in humans and chimpanzees. The other gorilla enterotype appears specialized for a leaf‐based diet and is enriched in environmentally derived bacterial genera. We hypothesize that the acquisition of this gorilla‐specific enterotype was enabled by lowered immune system control over the composition of the microbiome. Our results indicate differences between the pathology of SIVgor and SIVcpz/HIV‐1 infections, demonstrating the utility of investigating host microbial ecology as a means for studying disease in wild primates of high conservation priority.     

Bacterial and eukaryote microbiomes of mosquito habitats in dengue-endemic southern Taiwan

Mosquitoes interact with the microbiome of the waters where they oviposit in several ways. Past work suggests adult mosquitoes can detect certain microbes that stimulate oviposition. The presence or absence of certain microbes in water containers thus can attract or repel mosquito species to different containers. I hypothesized that these relationships could be detected via metagenomics. I focused on two container breeders that coexist in Southern Taiwan: the dengue vector Aedes aegypti and the less competent vector Ae. albopictus. In addition to culturing, I performed 16S and 18S rDNA metagenomics assays, the latter of which had never been applied to mosquito waters before, to identify the microbial diversity of artificial containers with and without mosquito larvae. I found no correlation between mosquito presence to any features of the containers or to their microbiomes, which instead correlated strongly with location. Microbial diversity across containers was highly variable, even within the same location, with multiple taxa only found in single containers. This variability is reasonable, because mosquito gut microbiomes are also extremely variable. The possibility remains that microbes in natural containers differ significantly from those in artificial containers, and that these differences drive Aedes preferences for human-associated containers. Broad, single-microbe experimental work is recommended to identify possible attractant or repellent microbial taxa.     

Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes

The gastrointestinal (GI) microbiome contributes significantly to host nutrition and health. However, relationships involving GI microbes, their hosts and host macrohabitats remain to be established. Here, we define clear patterns of variation in the GI microbiomes of six groups of Mexican black howler monkeys (Alouatta pigra) occupying a gradation of habitats including a continuous evergreen rainforest, an evergreen rainforest fragment, a continuous semi-deciduous forest and captivity. High throughput microbial 16S ribosomal RNA gene sequencing indicated that diversity, richness and composition of howler GI microbiomes varied with host habitat in relation to diet. Howlers occupying suboptimal habitats consumed less diverse diets and correspondingly had less diverse gut microbiomes. Quantitative real-time PCR also revealed a reduction in the number of genes related to butyrate production and hydrogen metabolism in the microbiomes of howlers occupying suboptimal habitats, which may impact host health.     

Meta-network: optimized species-species network analysis for microbial communities

Background

The explosive growth of microbiome data provides ample opportunities to gain a better understanding of the microbes and their interactions in microbial communities. Given these massive data, optimized data mining methods become important and necessary to perform deep and comprehensive analysis. Among the various priorities for microbiome data mining, the examination of species-species co-occurrence patterns becomes one of the key themes in urgent need.

Results

Hence, in this work, we propose the Meta-Network framework to lucubrate the microbial communities. Rooted in loose definitions of network (two species co-exist in a certain samples rather than all samples) as well as association rule mining (mining more complex forms of correlations like indirect correlation and mutual information), this framework outperforms other methods in restoring the microbial communities, based on two cohorts of microbial communities: (a) the loose definition strategy is capable to generate more reasonable relationships among species in the species-species co-occurrence network; (b) important species-species co-occurrence patterns could not be identified by other existing approaches, but could successfully generated by association rule mining.

Conclusions

Results have shown that the species-species co-occurrence network we generated are much more informative than those based on traditional methods. Meta-Network has consistently constructed more meaningful networks with biologically important clusters, hubs, and provides a general approach towards deciphering the species-species co-occurrence networks.

     

High-Specificity Targeted Functional Profiling in Microbial Communities with ShortBRED

Profiling microbial community function from metagenomic sequencing data remains a computationally challenging problem. Mapping millions of DNA reads from such samples to reference protein databases requires long run-times, and short read lengths can result in spurious hits to unrelated proteins (loss of specificity). We developed ShortBRED (Short, Better Representative Extract Dataset) to address these challenges, facilitating fast, accurate functional profiling of metagenomic samples. ShortBRED consists of two components: (i) a method that reduces reference proteins of interest to short, highly representative amino acid sequences (“markers”) and (ii) a search step that maps reads to these markers to quantify the relative abundance of their associated proteins. After evaluating ShortBRED on synthetic data, we applied it to profile antibiotic resistance protein families in the gut microbiomes of individuals from the United States, China, Malawi, and Venezuela. Our results support antibiotic resistance as a core function in the human gut microbiome, with tetracycline-resistant ribosomal protection proteins and Class A beta-lactamases being the most widely distributed resistance mechanisms worldwide. ShortBRED markers are applicable to other homology-based search tasks, which we demonstrate here by identifying phylogenetic signatures of antibiotic resistance across more than 3,000 microbial isolate genomes. ShortBRED can be applied to profile a wide variety of protein families of interest; the software, source code, and documentation are available for download at http://huttenhower.sph.harvard.edu/shortbred