Temporal variability is a personalized feature of the human microbiome

Background

It is now apparent that the complex microbial communities found on and in the human body vary across individuals. What has largely been missing from previous studies is an understanding of how these communities vary over time within individuals. To the extent to which it has been considered, it is often assumed that temporal variability is negligible for healthy adults. Here we address this gap in understanding by profiling the forehead, gut (fecal), palm, and tongue microbial communities in 85 adults, weekly over 3 months.

Results

We found that skin (forehead and palm) varied most in the number of taxa present, whereas gut and tongue communities varied more in the relative abundances of taxa. Within each body habitat, there was a wide range of temporal variability across the study population, with some individuals harboring more variable communities than others. The best predictor of these differences in variability across individuals was microbial diversity; individuals with more diverse gut or tongue communities were more stable in composition than individuals with less diverse communities.

Conclusions

Longitudinal sampling of a relatively large number of individuals allowed us to observe high levels of temporal variability in both diversity and community structure in all body habitats studied. These findings suggest that temporal dynamics may need to be considered when attempting to link changes in microbiome structure to changes in health status. Furthermore, our findings show that, not only is the composition of an individual’s microbiome highly personalized, but their degree of temporal variability is also a personalized feature.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0531-y) contains supplementary material, which is available to authorized users.     

Comparative genomics groups phages of Negativicutes and classical Firmicutes despite different Gram-staining properties

Negativicutes are gram-negative bacteria characterized by two cell membranes, but they are phylogenetically a side-branch of gram-positive Firmicutes that contain only a single membrane. We asked whether viruses (phages) infecting Negativicutes were horizontally acquired from gram-negative Proteobacteria, given the shared outer cell structure of their bacterial hosts, or if Negativicute phages co-evolved vertically with their hosts and thus resemble gram-positive Firmicute prophages. We predicted and characterized 485 prophages (mostly Caudovirales) from gram-negative Firmicute genomes plus 2977 prophages from other bacterial clades, and we used virome sequence data from 183 human stool samples to support our predictions. The majority of identified Negativicute prophages were lambdoids closer related to prophages from other Firmicutes than Proteobacteria by sequence relationship and genome organization (position of the lysis module). Only a single Mu-like candidate prophage and no clear P2-like prophages were identified in Negativicutes, both common in Proteobacteria. Given this collective evidence, it is unlikely that Negativicute phages were acquired from Proteobacteria. Sequence-related prophages, which occasionally harboured antibiotic resistance genes, were identified in two distinct Negativicute orders (Veillonellales and Acidaminococcales), possibly suggesting horizontal cross-order phage infection between human gut commensals. Our results reveal ancient genomic signatures of phage and bacteria co-evolution despite horizontal phage mobilization.     

Evolution along the parasitism-mutualism continuum determines the genetic repertoire of prophages

Integrated into their bacterial hosts’ genomes, prophage sequences exhibit a wide diversity of length and gene content, from highly degraded cryptic sequences to intact, functional prophages that retain a full complement of lytic-function genes. We apply three approaches—bioinformatics, analytical modelling and computational simulation—to understand the diverse gene content of prophages. In the bioinformatics work, we examine the distributions of over 50,000 annotated prophage genes identified in 1384 prophage sequences, comparing the gene repertoires of intact and incomplete prophages. These data indicate that genes involved in the replication, packaging, and release of phage particles have been preferentially lost in incomplete prophages, while tail fiber, transposase and integrase genes are significantly enriched. Consistent with these results, our mathematical and computational approaches predict that genes involved in phage lytic function are preferentially lost, resulting in shorter prophages that often retain genes that benefit the host. Informed by these models, we offer novel hypotheses for the enrichment of integrase and transposase genes in cryptic prophages. Overall, we demonstrate that functional and cryptic prophages represent a diversity of genetic sequences that evolve along a parasitism-mutualism continuum.     

Increased rectal microbial richness is associated with the presence of colorectal adenomas in humans

Differences in the composition of the gut microbial community have been associated with diseases such as obesity, Crohn''s disease, ulcerative colitis and colorectal cancer (CRC). We used 454 titanium pyrosequencing of the V1–V2 region of the 16S rRNA gene to characterize adherent bacterial communities in mucosal biopsy samples from 33 subjects with adenomas and 38 subjects without adenomas (controls). Biopsy samples from subjects with adenomas had greater numbers of bacteria from 87 taxa than controls; only 5 taxa were more abundant in control samples. The magnitude of the differences in the distal gut microbiota between patients with adenomas and controls was more pronounced than that of any other clinical parameters including obesity, diet or family history of CRC. This suggests that sequence analysis of the microbiota could be used to identify patients at risk for developing adenomas.     

Evidence supporting the microbiota–gut–brain axis in a songbird

Recent research in mammals supports a link between cognitive ability and the gut microbiome, but little is known about this relationship in other taxa. In a captive population of 38 zebra finches (Taeniopygia guttata), we quantified performance on cognitive tasks measuring learning and memory. We sampled the gut microbiome via cloacal swab and quantified bacterial alpha and beta diversity. Performance on cognitive tasks related to beta diversity but not alpha diversity. We then identified differentially abundant genera influential in the beta diversity differences among cognitive performance categories. Though correlational, this study provides some of the first evidence of an avian microbiota–gut–brain axis, building foundations for future microbiome research in wild populations and during host development.     

Host-specificity among abundant and rare taxa in the sponge microbiome

Microbial communities have a key role in the physiology of the sponge host, and it is therefore essential to understand the stability and specificity of sponge–symbiont associations. Host-specific bacterial associations spanning large geographic distance are widely acknowledged in sponges. However, the full spectrum of specificity remains unclear. In particular, it is not known whether closely related sponges host similar or very different microbiota over wide bathymetric and geographic gradients, and whether specific associations extend to the rare members of the sponge microbiome. Using the ultra-deep Illumina sequencing technology, we conducted a comparison of sponge bacterial communities in seven closely related Hexadella species with a well-resolved host phylogeny, as well as of a distantly related sponge Mycale. These samples spanned unprecedentedly large bathymetric (15–960 m) gradients and varying European locations. In addition, this study included a bacterial community analysis of the local background seawater for both Mycale and the widespread deep-sea taxa Hexadella cf. dedritifera. We observed a striking diversity of microbes associated with the sponges, spanning 47 bacterial phyla. The data did not reveal any Hexadella microbiota co-speciation pattern, but confirmed sponge-specific and species-specific host–bacteria associations, even within extremely low abundant taxa. Oligotyping analysis also revealed differential enrichment preferences of closely related Nitrospira members in closely related sponges species. Overall, these results demonstrate highly diverse, remarkably specific and stable sponge–bacteria associations that extend to members of the rare biosphere at a very fine phylogenetic scale, over significant geographic and bathymetric gradients.     

Modular prophage interactions driven by capsule serotype select for capsule loss under phage predation

Klebsiella species are able to colonize a wide range of environments and include worrisome nosocomial pathogens. Here, we sought to determine the abundance and infectivity of prophages of Klebsiella to understand how the interactions between induced prophages and bacteria affect population dynamics and evolution. We identified many prophages in the species, placing these taxa among the top 5% of the most polylysogenic bacteria. We selected 35 representative strains of the Klebsiella pneumoniae species complex to establish a network of induced phage–bacteria interactions. This revealed that many prophages are able to enter the lytic cycle, and subsequently kill or lysogenize closely related Klebsiella strains. Although 60% of the tested strains could produce phages that infect at least one other strain, the interaction network of all pairwise cross-infections is very sparse and mostly organized in modules corresponding to the strains’ capsule serotypes. Accordingly, capsule mutants remain uninfected showing that the capsule is a key factor for successful infections. Surprisingly, experiments in which bacteria are predated by their own prophages result in accelerated loss of the capsule. Our results show that phage infectiousness defines interaction modules between small subsets of phages and bacteria in function of capsule serotype. This limits the role of prophages as competitive weapons because they can infect very few strains of the species complex. This should also restrict phage-driven gene flow across the species. Finally, the accelerated loss of the capsule in bacteria being predated by their own phages, suggests that phages drive serotype switch in nature.Subject terms: Bacteriophages, Microbial ecology, Environmental microbiology     

Identification of Specialists and Abundance-Occupancy Relationships among Intestinal Bacteria of Aves,Mammalia, and Actinopterygii

The coalescence of next-generation DNA sequencing methods, ecological perspectives, and bioinformatics analysis tools is rapidly advancing our understanding of the evolution and function of vertebrate-associated bacterial communities. Delineation of host-microbe associations has applied benefits ranging from clinical treatments to protecting our natural waters. Microbial communities follow some broad-scale patterns observed for macroorganisms, but it remains unclear how the specialization of intestinal vertebrate-associated communities to a particular host environment influences broad-scale patterns in microbial abundance and distribution. We analyzed the V6 region of 16S rRNA genes amplified from 106 fecal samples spanning Aves, Mammalia, and Actinopterygii (ray-finned fish). We investigated the interspecific abundance-occupancy relationship, where widespread taxa tend to be more abundant than narrowly distributed taxa, among operational taxonomic units (OTUs) within and among host species. In a separate analysis, we identified specialist OTUs that were highly abundant in a single host and rare in all other hosts by using a multinomial model without excluding undersampled OTUs a priori. We show that intestinal microbes in humans and other vertebrates display abundance-occupancy relationships, but because intestinal host-associated communities have undergone intense specialization, this trend is violated by a disproportionately large number of specialist taxa. Although it is difficult to distinguish the effects of dispersal limitations, host selection, historical contingency, and stochastic processes on community assembly, results suggest that intestinal bacteria can be shared among diverse hosts in ways that resemble the distribution of “free-living” bacteria in the extraintestinal environment.     

Comparison of the Distal Gut Microbiota from People and Animals in Africa

The gut microbiota plays a key role in the maintenance of healthy gut function as well as many other aspects of health. High-throughput sequence analyses have revealed the composition of the gut microbiota, showing that there is a core signature to the human gut microbiota, as well as variation in its composition between people. The gut microbiota of animals is also being investigated. We are interested in the relationship between bacterial taxa of the human gut microbiota and those in the gut microbiota of domestic and semi-wild animals. While it is clear that some human gut bacterial pathogens come from animals (showing that human – animal transmission occurs), the extent to which the usually non-pathogenic commensal taxa are shared between humans and animals has not been explored. To investigate this we compared the distal gut microbiota of humans, cattle and semi-captive chimpanzees in communities that are geographically sympatric in Uganda. The gut microbiotas of these three host species could be distinguished by the different proportions of bacterial taxa present. We defined multiple operational taxonomic units (OTUs) by sequence similarity and found evidence that some OTUs were common between human, cattle and chimpanzees, with the largest number of shared OTUs occurring between chimpanzees and humans, as might be expected with their close physiological similarity. These results show the potential for the sharing of usually commensal bacterial taxa between humans and other animals. This suggests that further investigation of this phenomenon is needed to fully understand how it drives the composition of human and animal gut microbiotas.     

Bioinformatics Analysis of Antimicrobial Resistance Genes and Prophages Colocalized in Human Gut Metagenomes

The constant increase of bacterial antibiotic-resistant strains is directly linked to a common use of antibiotics in medicine and animal breeding. It is suggested that the gut microbiota serves as a reservoir for antibiotic resistance genes that can be transferred from symbiotic bacteria to pathogenic ones, particularly due to phage transduction. In this study, using the PHASTER prophage predicting tool and CARD antibiotics resistance database we have searched for antibiotic resistance genes that are located within prophages in human gut microbiota. After analysing metagenomic assemblies of eight samples of antibiotic treated patients, lsaE, mdfA, and cpxR/cpxA genes were identified inside prophages. These genes confer resistance to antimicrobial peptides, pleuromutilin, lincomycins, streptogramins and also multidrug resistance. Three (0.46%) of 659 putative prophages predicted in the metagenomic assemblies contained antibiotics resistance genes in their sequences.     

Higher contribution of globally rare bacterial taxa reflects environmental transitions across the surface ocean

Microbial taxa range from being ubiquitous and abundant across space to extremely rare and endemic, depending on their ecophysiology and on different processes acting locally or regionally. However, little is known about how cosmopolitan or rare taxa combine to constitute communities and whether environmental variations promote changes in their relative abundances. Here we identified the Spatial Abundance Distribution (SpAD) of individual prokaryotic taxa (16S rDNA‐defined Operational Taxonomic Units, OTUs) across 108 globally‐distributed surface ocean stations. We grouped taxa based on their SpAD shape (“normal‐like”‐ abundant and ubiquitous; “logistic”‐ globally rare, present in few sites; and “bimodal”‐ abundant only in certain oceanic regions), and investigated how the abundance of these three categories relates to environmental gradients. Most surface assemblages were numerically dominated by a few cosmopolitan “normal‐like” OTUs, yet there was a gradual shift towards assemblages dominated by “logistic” taxa in specific areas with productivity and temperature differing the most from the average conditions in the sampled stations. When we performed the SpAD categorization including additional habitats (deeper layers and particles of varying sizes), the SpAD of many OTUs changed towards fewer “normal‐like” shapes, and OTUs categorized as globally rare in the surface ocean became abundant. This suggests that understanding the mechanisms behind microbial rarity and dominance requires expanding the context of study beyond local communities and single habitats. We show that marine bacterial communities comprise taxa displaying a continuum of SpADs, and that variations in their abundances can be linked to habitat transitions or barriers that delimit the distribution of community members.     

Diversity and dynamics of bacteriocins from human microbiome

Human commensal microbiota are an important determinant of health and disease of the host. Different human body sites harbour different bacterial microbiota, bacterial communities that maintain a stable balance. However, many of the factors influencing the stabilities of bacterial communities associated with humans remain unknown. In this study, we identified putative bacteriocins produced by human commensal microbiota. Bacteriocins are peptides or proteins with antimicrobial activity that contribute to the stability and dynamics of microbial communities. We employed bioinformatic analyses to identify putative bacteriocin sequences in metagenomic sequences obtained from different human body sites. Prevailing bacterial taxa of the putative bacteriocins producers matched the most abundant organisms in each human body site. Remarkably, we found that samples from different body sites contain different density of putative bacteriocin genes, with the highest in samples from the vagina, the airway, and the oral cavity and the lowest in those from gut. Inherent differences of different body sites thus influence the density and types of bacteriocins produced by commensal bacteria. Our results suggest that bacteriocins play important roles to allow different bacteria to occupy several human body sites, and to establish a long‐term commensal relationship with human hosts.     

Prophage-Mediated Dynamics of ‘Candidatus Liberibacter asiaticus’ Populations,the Destructive Bacterial Pathogens of Citrus Huanglongbing

Prophages are highly dynamic components in the bacterial genome and play an important role in intraspecies variations. There are at least two prophages in the chromosomes of Candidatus Liberibacter asiaticus’ (Las) Floridian isolates. Las is both unculturable and the most prevalent species of Liberibacter pathogens that cause huanglongbing (HLB), a worldwide destructive disease of citrus. In this study, seven new prophage variants resulting from two hyper-variable regions were identified by screening clone libraries of infected citrus, periwinkle and psyllids. Among them, Types A and B share highly conserved sequences and localize within the two prophages, FP1 and FP2, respectively. Although Types B and C were abundant in all three libraries, Type A was much more abundant in the libraries from the Las-infected psyllids than from the Las-infected plants, and Type D was only identified in libraries from the infected host plants but not from the infected psyllids. Sequence analysis of these variants revealed that the variations may result from recombination and rearrangement events. Conventional PCR results using type-specific molecular markers indicated that A, B, C and D are the four most abundant types in Las-infected citrus and periwinkle. However, only three types, A, B and C are abundant in Las-infected psyllids. Typing results for Las-infected citrus field samples indicated that mixed populations of Las bacteria present in Floridian isolates, but only the Type D population was correlated with the blotchy mottle symptom. Extended cloning and sequencing of the Type D region revealed a third prophage/phage in the Las genome, which may derive from the recombination of FP1 and FP2. Dramatic variations in these prophage regions were also found among the global Las isolates. These results are the first to demonstrate the prophage/phage-mediated dynamics of Las populations in plant and insect hosts, and their correlation with insect transmission and disease development.     

Geographical Range and Local Abundance of Tree Species in China

Most studies on the geographical distribution of species have utilized a few well-known taxa in Europe and North America, with little research in China and its wide range of climate and forest types. We assembled large datasets to quantify the geographic ranges of tree species in China and to test several biogeographic hypotheses: 1) whether locally abundant species tend to be geographically widespread; 2) whether species are more abundant towards their range-centers; and 3) how abundances are correlated between sites. Local abundances of 651 species were derived from four tree plots of 20–25 ha where all individuals ≥1 cm in stem diameter were mapped and identified taxonomically. Range sizes of these species across China were then estimated from over 460,000 geo-referenced records; a Bayesian approach was used, allowing careful measures of error of each range estimate. The log-transformed range sizes had a bell-shaped distribution with a median of 703,000 km², and >90% of 651 species had ranges >10⁵ km². There was no relationship between local abundance and range size, and no evidence for species being more abundant towards their range-centers. Finally, species’ abundances were positively correlated between sites. The widespread nature of most tree species in China suggests few are vulnerable to global extinction, and there is no indication of the double-peril that would result if rare species also had narrow ranges.     

Identification of Aminoglycoside and β-Lactam Resistance Genes from within an Infant Gut Functional Metagenomic Library

The infant gut microbiota develops rapidly during the first 2 years of life, acquiring microorganisms from diverse sources. During this time, significant opportunities exist for the infant to acquire antibiotic resistant bacteria, which can become established and constitute the infant gut resistome. With increased antibiotic resistance limiting our ability to treat bacterial infections, investigations into resistance reservoirs are highly pertinent. This study aimed to explore the nascent resistome in antibiotically-naïve infant gut microbiomes, using a combination of metagenomic approaches. Faecal samples from 22 six-month-old infants without previous antibiotic exposure were used to construct a pooled metagenomic library, which was functionally screened for ampicillin and gentamicin resistance. Our library of ∼220Mb contained 0.45 ampicillin resistant hits/Mb and 0.059 gentamicin resistant hits/Mb. PCR-based analysis of fosmid clones and uncloned metagenomic DNA, revealed a diverse and abundant aminoglycoside and β-lactam resistance reservoir within the infant gut, with resistance determinants exhibiting homology to those found in common gut inhabitants, including Escherichia coli, Enterococcus sp., and Clostridium difficile, as well as to genes from cryptic environmental bacteria. Notably, the genes identified differed from those revealed when a sequence-driven PCR-based screen of metagenomic DNA was employed. Carriage of these antibiotic resistance determinants conferred substantial, but varied (2–512x), increases in antibiotic resistance to their bacterial host. These data provide insights into the infant gut resistome, revealing the presence of a varied aminoglycoside and β-lactam resistance reservoir even in the absence of selective pressure, confirming the infant resistome establishes early in life, perhaps even at birth.     

A single-cell polony method reveals low levels of infected Prochlorococcus in oligotrophic waters despite high cyanophage abundances

Long-term stability of picocyanobacteria in the open oceans is maintained by a balance between synchronous division and death on daily timescales. Viruses are considered a major source of microbial mortality, however, current methods to measure infection have significant methodological limitations. Here we describe a method that pairs flow-cytometric sorting with a PCR-based polony technique to simultaneously screen thousands of taxonomically resolved individual cells for intracellular virus DNA, enabling sensitive, high-throughput, and direct quantification of infection by different virus lineages. Under controlled conditions with picocyanobacteria-cyanophage models, the method detected infection throughout the lytic cycle and discriminated between varying infection levels. In North Pacific subtropical surface waters, the method revealed that only a small percentage of Prochlorococcus (0.35–1.6%) were infected, predominantly by T4-like cyanophages, and that infection oscillated 2-fold in phase with the diel cycle. This corresponds to 0.35–4.8% of Prochlorococcus mortality daily. Cyanophages were 2–4-fold more abundant than Prochlorococcus, indicating that most encounters did not result in infection and suggesting infection is mitigated via host resistance, reduced phage infectivity and inefficient adsorption. This method will enable quantification of infection for key microbial taxa across oceanic regimes and will help determine the extent that viruses shape microbial communities and ecosystem level processes.Subject terms: Microbial ecology, Microbial biooceanography, Molecular ecology, Microbial ecology, Microbial biooceanography     

Phylogenetic analysis of the fecal microbial community in herbivorous land and marine iguanas of the Galápagos Islands using 16S rRNA-based pyrosequencing

Herbivorous reptiles depend on complex gut microbial communities to effectively degrade dietary polysaccharides. The composition of these fermentative communities may vary based on dietary differences. To explore the role of diet in shaping gut microbial communities, we evaluated the fecal samples from two related host species—the algae-consuming marine iguana (Amblyrhynchus cristatus) and land iguanas (LI) (genus Conolophus) that consume terrestrial vegetation. Marine and LI fecal samples were collected from different islands in the Galápagos archipelago. High-throughput 16S rRNA-based pyrosequencing was used to provide a comparative analysis of fecal microbial diversity. At the phylum level, the fecal microbial community in iguanas was predominated by Firmicutes (69.5±7.9%) and Bacteroidetes (6.2±2.8%), as well as unclassified Bacteria (20.6±8.6%), suggesting that a large portion of iguana fecal microbiota is novel and could be involved in currently unknown functions. Host species differed in the abundance of specific bacterial groups. Bacteroides spp., Lachnospiraceae and Clostridiaceae were significantly more abundant in the marine iguanas (MI) (P-value>1E−9). In contrast, Ruminococcaceae were present at >5-fold higher abundance in the LI than MI (P-value>6E−14). Archaea were only detected in the LI. The number of operational taxonomic units (OTUs) in the LI (356–896 OTUs) was >2-fold higher than in the MI (112–567 OTUs), and this increase in OTU diversity could be related to the complexity of the resident bacterial population and their gene repertoire required to breakdown the recalcitrant polysaccharides prevalent in terrestrial plants. Our findings suggest that dietary differences contribute to gut microbial community differentiation in herbivorous lizards. Most importantly, this study provides a better understanding of the microbial diversity in the iguana gut; therefore facilitating future efforts to discover novel bacterial-associated enzymes that can effectively breakdown a wide variety of complex polysaccharides.     

Patterns of gammaridean amphipod abundance and species composition associated with dominant subtidal macroalgae from the western Antarctic Peninsula

The communities of gammaridean amphipods associated with eight dominant macroalgal species were examined near Palmer Station, Western Antarctic Peninsula. A total of 78,415 individuals belonging to 32 amphipod taxa were identified with mean densities ranging up to 20 individuals/g algal wet wt. The most abundant amphipod taxon, Metaleptamphopus pectinatus, was found to associate predominately with the brown alga Desmarestia menziesii, while the second most common taxon, Jassa spp. occurred primarily on the red alga Gigartina skottsbergii. Non-metric multidimensional scaling analysis demonstrated that the population densities of each amphipod species and amphipod species composition were similar on the same algal species but dissimilar on different species of algae. Comparisons of amphipod communities associated with a given algal species but from different sampling sites indicated that although the structure of species-specific macroalgal-associated amphipod communities can vary across spatial scales of 3 km, 50% of the macroalgal species examined showed no significant inter-site differences in associated amphipod community structure. Spearman rank correlation analyses showed that higher abundances of amphipods occurred on the macroalgae with the highest number of branches. As many Antarctic amphipods are known consumers of macroalgae, their remarkable abundances are likely to play a significant role in mediating energy and nutrient transfer in nearshore Antarctic Peninsular macroalgal communities.     

Analysis of six prophages in Lactococcus lactis IL1403: different genetic structure of temperate and virulent phage populations

We report the genetic organisation of six prophages present in the genome of Lactococcus lactis IL1403. The three larger prophages (36–42 kb), belong to the already described P335 group of temperate phages, whereas the three smaller ones (13–15 kb) are most probably satellites relying on helper phage(s) for multiplication. These data give a new insight into the genetic structure of lactococcal phage populations. P335 temperate phages have variable genomes, sharing homology over only 10–33% of their length. In contrast, virulent phages have highly similar genomes sharing homology over >90% of their length. Further analysis of genetic structure in all known groups of phages active on other bacterial hosts such as Escherichia coli, Bacillus subtilis, Mycobacterium and Streptococcus thermophilus confirmed the existence of two types of genetic structure related to the phage way of life. This might reflect different intensities of horizontal DNA exchange: low among purely virulent phages and high among temperate phages and their lytic homologues. We suggest that the constraints on genetic exchange among purely virulent phages reflect their optimal genetic organisation, adapted to a more specialised and extreme form of parasitism than temperate/lytic phages.