Global patterns in the biogeography of bacterial taxa.

Bacteria control major nutrient cycles and directly influence plant, animal and human health. However, we know relatively little about the forces shaping their large-scale ecological ranges. Here, we reveal patterns in the distribution of individual bacterial taxa at multiple levels of phylogenetic resolution within and between Earth's major habitat types. Our analyses suggest that while macro-scale habitats structure bacterial distribution to some degree, abundant bacteria (i.e. detectable using 16S rRNA gene sequencing methods) are confined to single assemblages. Additionally, we show that the most cosmopolitan taxa are also the most abundant in individual assemblages. These results add to the growing body of data that support that the diversity of the overall bacterial metagenome is tremendous. The mechanisms governing microbial distribution remain poorly understood, but our analyses provide a framework with which to test the importance of macro-ecological environmental gradients, relative abundance, neutral processes and the ecological strategies of individual taxa in structuring microbial communities.

Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data.

Next-generation sequencing techniques, and PhyloChip, have made simultaneous phylogenetic analyses of hundreds of microbial communities possible. Insight into community structure has been limited by the inability to integrate and visualize such vast datasets. Fast UniFrac overcomes these issues, allowing integration of larger numbers of sequences and samples into a single analysis. Its new array-based implementation offers orders of magnitude improvements over the original version. New 3D visualization of principal coordinates analysis results, with the option to view multiple coordinate axes simultaneously, provides a powerful way to quickly identify patterns that relate vast numbers of microbial communities. We show the potential of Fast UniFrac using examples from three data types: Sanger-sequencing studies of diverse free-living and animal-associated bacterial assemblages and from the gut of obese humans as they diet, pyrosequencing data integrated from studies of the human hand and gut, and PhyloChip data from a study of citrus pathogens. We show that a Fast UniFrac analysis using a reference tree recaptures patterns that could not be detected without considering phylogenetic relationships and that Fast UniFrac, coupled with BLAST-based sequence assignment, can be used to quickly analyze pyrosequencing runs containing hundreds of thousands of sequences, showing patterns relating human and gut samples. Finally, we show that the application of Fast UniFrac to PhyloChip data could identify well-defined subcategories associated with infection. Together, these case studies point the way toward a broad range of applications and show some of the new features of Fast UniFrac.

Bacterial community variation in human body habitats across space and time.

Elucidating the biogeography of bacterial communities on the human body is critical for establishing healthy baselines from which to detect differences associated with diseases. To obtain an integrated view of the spatial and temporal distribution of the human microbiota, we surveyed bacteria from up to 27 sites in seven to nine healthy adults on four occasions. We found that community composition was determined primarily by body habitat. Within habitats, interpersonal variability was high, whereas individuals exhibited minimal temporal variability. Several skin locations harbored more diverse communities than the gut and mouth, and skin locations differed in their community assembly patterns. These results indicate that our microbiota, although personalized, varies systematically across body habitats and time; such trends may ultimately reveal how microbiome changes cause or prevent disease.

High-fat diet determines the composition of the murine gut microbiome independently of obesity.

BACKGROUND & AIMS: The composition of the gut microbiome is affected by host phenotype, genotype, immune function, and diet. Here, we used the phenotype of RELMbeta knockout (KO) mice to assess the influence of these factors. METHODS: Both wild-type and RELMbeta KO mice were lean on a standard chow diet, but, upon switching to a high-fat diet, wild-type mice became obese, whereas RELMbeta KO mice remained comparatively lean. To investigate the influence of diet, genotype, and obesity on microbiome composition, we used deep sequencing to characterize 25,790 16S rDNA sequences from uncultured bacterial communities from both genotypes on both diets. RESULTS: We found large alterations associated with switching to the high-fat diet, including a decrease in Bacteroidetes and an increase in both Firmicutes and Proteobacteria. This was seen for both genotypes (ie, in the presence and absence of obesity), indicating that the high-fat diet itself, and not the obese state, mainly accounted for the observed changes in the gut microbiota. The RELMbeta genotype also modestly influenced microbiome composition independently of diet. Metagenomic analysis of 537,604 sequence reads documented extensive changes in gene content because of a high-fat diet, including an increase in transporters and 2-component sensor responders as well as a general decrease in metabolic genes. Unexpectedly, we found a substantial amount of murine DNA in our samples that increased in proportion on a high-fat diet. CONCLUSIONS: These results demonstrate the importance of diet as a determinant of gut microbiome composition and suggest the need to control for dietary variation when evaluating the composition of the human gut microbiome.

Characterization of airborne microbial communities at a high-elevation site and their potential to act as atmospheric ice nuclei.

Bacteria and fungi are ubiquitous in the atmosphere. The diversity and abundance of airborne microbes may be strongly influenced by atmospheric conditions or even influence atmospheric conditions themselves by acting as ice nucleators. However, few comprehensive studies have described the diversity and dynamics of airborne bacteria and fungi based on culture-independent techniques. We document atmospheric microbial abundance, community composition, and ice nucleation at a high-elevation site in northwestern Colorado. We used a standard small-subunit rRNA gene Sanger sequencing approach for total microbial community analysis and a bacteria-specific 16S rRNA bar-coded pyrosequencing approach (4,864 sequences total). During the 2-week collection period, total microbial abundances were relatively constant, ranging from 9.6 x 10(5) to 6.6 x 10(6) cells m(-3) of air, and the diversity and composition of the airborne microbial communities were also relatively static. Bacteria and fungi were nearly equivalent, and members of the proteobacterial groups Burkholderiales and Moraxellaceae (particularly the genus Psychrobacter) were dominant. These taxa were not always the most abundant in freshly fallen snow samples collected at this site. Although there was minimal variability in microbial abundances and composition within the atmosphere, the number of biological ice nuclei increased significantly during periods of high relative humidity. However, these changes in ice nuclei numbers were not associated with changes in the relative abundances of the most commonly studied ice-nucleating bacteria.

Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale.

Soils harbor enormously diverse bacterial populations, and soil bacterial communities can vary greatly in composition across space. However, our understanding of the specific changes in soil bacterial community structure that occur across larger spatial scales is limited because most previous work has focused on either surveying a relatively small number of soils in detail or analyzing a larger number of soils with techniques that provide little detail about the phylogenetic structure of the bacterial communities. Here we used a bar-coded pyrosequencing technique to characterize bacterial communities in 88 soils from across North and South America, obtaining an average of 1,501 sequences per soil. We found that overall bacterial community composition, as measured by pairwise UniFrac distances, was significantly correlated with differences in soil pH (r = 0.79), largely driven by changes in the relative abundances of Acidobacteria, Actinobacteria, and Bacteroidetes across the range of soil pHs. In addition, soil pH explains a significant portion of the variability associated with observed changes in the phylogenetic structure within each dominant lineage. The overall phylogenetic diversity of the bacterial communities was also correlated with soil pH (R(2) = 0.50), with peak diversity in soils with near-neutral pHs. Together, these results suggest that the structure of soil bacterial communities is predictable, to some degree, across larger spatial scales, and the effect of soil pH on bacterial community composition is evident at even relatively coarse levels of taxonomic resolution.

Regulation of myocardial ketone body metabolism by the gut microbiota during nutrient deprivation.

Studies in mice indicate that the gut microbiota promotes energy harvest and storage from components of the diet when these components are plentiful. Here we examine how the microbiota shapes host metabolic and physiologic adaptations to periods of nutrient deprivation. Germ-free (GF) mice and mice who had received a gut microbiota transplant from conventionally raised donors were compared in the fed and fasted states by using functional genomic, biochemical, and physiologic assays. A 24-h fast produces a marked change in gut microbial ecology. Short-chain fatty acids generated from microbial fermentation of available glycans are maintained at higher levels compared with GF controls. During fasting, a microbiota-dependent, Ppar alpha-regulated increase in hepatic ketogenesis occurs, and myocardial metabolism is directed to ketone body utilization. Analyses of heart rate, hydraulic work, and output, mitochondrial morphology, number, and respiration, plus ketone body, fatty acid, and glucose oxidation in isolated perfused working hearts from GF and colonized animals (combined with in vivo assessments of myocardial physiology) revealed that the fasted GF heart is able to sustain its performance by increasing glucose utilization, but heart weight, measured echocardiographically or as wet mass and normalized to tibial length or lean body weight, is significantly reduced in both fasted and fed mice. This myocardial-mass phenotype is completely reversed in GF mice by consumption of a ketogenic diet. Together, these results illustrate benefits provided by the gut microbiota during periods of nutrient deprivation, and emphasize the importance of further exploring the relationship between gut microbes and cardiovascular health.

CodonExplorer: an interactive online database for the analysis of codon usage and sequence composition.

The analysis of DNA composition and codon usage reveals many factors that influence the evolution of genes and genomes. In this chapter, we show how to use CodonExplorer, a web tool and interactive database that contains millions of genes, to better understand the principles governing evolution at the single gene and whole-genome level. We present principles and practical procedures for using analyses of GC content and codon usage frequency to identify highly expressed or horizontally transferred genes and to study the relative contribution of different types of mutation to gene and genome composition. CodonExplorer's combination of a user-friendly web interface and a comprehensive genomic database makes these diverse analyses fast and straightforward to perform. CodonExplorer is thus a powerful tool that facilitates and automates a wide range of compositional analyses.

Microbial community profiling for human microbiome projects: Tools, techniques, and challenges.

High-throughput sequencing studies and new software tools are revolutionizing microbial community analyses, yet the variety of experimental and computational methods can be daunting. In this review, we discuss some of the different approaches to community profiling, highlighting strengths and weaknesses of various experimental approaches, sequencing methodologies, and analytical methods. We also address one key question emerging from various Human Microbiome Projects: Is there a substantial core of abundant organisms or lineages that we all share? It appears that in some human body habitats, such as the hand and the gut, the diversity among individuals is so great that we can rule out the possibility that any species is at high abundance in all individuals: It is possible that the focus should instead be on higher-level taxa or on functional genes instead.

CodonExplorer: an online tool for analyzing codon usage and sequence composition, scaling from genes to genomes.

DNA composition in general, and codon usage in particular, is crucial for understanding gene function and evolution. CodonExplorer, available online at http://bmf.colorado.edu/codonexplorer/, is an online tool and interactive database that contains millions of genes, allowing rapid exploration of the factors governing gene and genome compositional evolution and exploiting GC content and codon usage frequency to identify genes with composition suggesting high levels of expression or horizontal transfer.

A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses.

Acidobacteria are ubiquitous and abundant members of soil bacterial communities. However, an ecological understanding of this important phylum has remained elusive because its members have been difficult to culture and few molecular investigations have focused exclusively on this group. We generated an unprecedented number of acidobacterial DNA sequence data using pyrosequencing and clone libraries (39,707 and 1787 sequences, respectively) to characterize the relative abundance, diversity and composition of acidobacterial communities across a range of soil types. To gain insight into the ecological characteristics of acidobacterial taxa, we investigated the large-scale biogeographic patterns exhibited by acidobacterial communities, and related soil and site characteristics to acidobacterial community assemblage patterns. The 87 soils analyzed by pyrosequencing contained more than 8600 unique acidobacterial phylotypes (at the 97% sequence similarity level). One phylotype belonging to Acidobacteria subgroup 1, but not closely related to any cultured representatives, was particularly abundant, accounting for 7.4% of bacterial sequences and 17.6% of acidobacterial sequences, on average, across the soils. The abundance of Acidobacteria relative to other bacterial taxa was highly variable across the soils examined, but correlated strongly with soil pH (R=-0.80, P<0.001). Soil pH was also the best predictor of acidobacterial community composition, regardless of how the communities were characterized, and the relative abundances of the dominant Acidobacteria subgroups were readily predictable. Acidobacterial communities were more phylogenetically clustered as soil pH departed from neutrality, suggesting that pH is an effective habitat filter, restricting community membership to progressively more narrowly defined lineages as pH deviates from neutrality.

A core gut microbiome in obese and lean twins.

The human distal gut harbours a vast ensemble of microbes (the microbiota) that provide important metabolic capabilities, including the ability to extract energy from otherwise indigestible dietary polysaccharides. Studies of a few unrelated, healthy adults have revealed substantial diversity in their gut communities, as measured by sequencing 16S rRNA genes, yet how this diversity relates to function and to the rest of the genes in the collective genomes of the microbiota (the gut microbiome) remains obscure. Studies of lean and obese mice suggest that the gut microbiota affects energy balance by influencing the efficiency of calorie harvest from the diet, and how this harvested energy is used and stored. Here we characterize the faecal microbial communities of adult female monozygotic and dizygotic twin pairs concordant for leanness or obesity, and their mothers, to address how host genotype, environmental exposure and host adiposity influence the gut microbiome. Analysis of 154 individuals yielded 9,920 near full-length and 1,937,461 partial bacterial 16S rRNA sequences, plus 2.14 gigabases from their microbiomes. The results reveal that the human gut microbiome is shared among family members, but that each person's gut microbial community varies in the specific bacterial lineages present, with a comparable degree of co-variation between adult monozygotic and dizygotic twin pairs. However, there was a wide array of shared microbial genes among sampled individuals, comprising an extensive, identifiable 'core microbiome' at the gene, rather than at the organismal lineage, level. Obesity is associated with phylum-level changes in the microbiota, reduced bacterial diversity and altered representation of bacterial genes and metabolic pathways. These results demonstrate that a diversity of organismal assemblages can nonetheless yield a core microbiome at a functional level, and that deviations from this core are associated with different physiological states (obese compared with lean).

The influence of sex, handedness, and washing on the diversity of hand surface bacteria.

Bacteria thrive on and within the human body. One of the largest human-associated microbial habitats is the skin surface, which harbors large numbers of bacteria that can have important effects on health. We examined the palmar surfaces of the dominant and nondominant hands of 51 healthy young adult volunteers to characterize bacterial diversity on hands and to assess its variability within and between individuals. We used a novel pyrosequencing-based method that allowed us to survey hand surface bacterial communities at an unprecedented level of detail. The diversity of skin-associated bacterial communities was surprisingly high; a typical hand surface harbored >150 unique species-level bacterial phylotypes, and we identified a total of 4,742 unique phylotypes across all of the hands examined. Although there was a core set of bacterial taxa commonly found on the palm surface, we observed pronounced intra- and interpersonal variation in bacterial community composition: hands from the same individual shared only 17% of their phylotypes, with different individuals sharing only 13%. Women had significantly higher diversity than men, and community composition was significantly affected by handedness, time since last hand washing, and an individual's sex. The variation within and between individuals in microbial ecology illustrated by this study emphasizes the challenges inherent in defining what constitutes a "healthy" bacterial community; addressing these challenges will be critical for the International Human Microbiome Project.

The convergence of carbohydrate active gene repertoires in human gut microbes.

The extreme variation in gene content among phylogenetically related microorganisms suggests that gene acquisition, expansion, and loss are important evolutionary forces for adaptation to new environments. Accordingly, phylogenetically disparate organisms that share a habitat may converge in gene content as they adapt to confront shared challenges. This response should be especially pronounced for functional genes that are important for survival in a particular habitat. We illustrate this principle by showing that the repertoires of two different types of carbohydrate-active enzymes, glycoside hydrolases and glycosyltransferases, have converged in bacteria and archaea that live in the human gut and that this convergence is largely due to horizontal gene transfer rather than gene family expansion. We also identify gut microbes that may have more similar dietary niches in the human gut than would be expected based on phylogeny. The techniques used to obtain these results should be broadly applicable to understanding the functional genes and evolutionary processes important for adaptation in many environments and useful for interpreting the large number of reference microbial genome sequences being generated for the International Human Microbiome Project.

Worlds within worlds: evolution of the vertebrate gut microbiota.

In this Analysis we use published 16S ribosomal RNA gene sequences to compare the bacterial assemblages that are associated with humans and other mammals, metazoa and free-living microbial communities that span a range of environments. The composition of the vertebrate gut microbiota is influenced by diet, host morphology and phylogeny, and in this respect the human gut bacterial community is typical of an omnivorous primate. However, the vertebrate gut microbiota is different from free-living communities that are not associated with animal body habitats. We propose that the recently initiated international Human Microbiome Project should strive to include a broad representation of humans, as well as other mammalian and environmental samples, as comparative analyses of microbiotas and their microbiomes are a powerful way to explore the evolutionary history of the biosphere.

MotifCluster: an interactive online tool for clustering and visualizing sequences using shared motifs.

MotifCluster finds related motifs in a set of sequences, and clusters the sequences into families using the motifs they contain. MotifCluster, at http://bmf.colorado.edu/motifcluster, lets users test whether proteins are related, cluster sequences by shared conserved motifs, and visualize motifs mapped onto trees, sequences and three-dimensional structures. We demonstrate MotifCluster's accuracy using gold-standard protein superfamilies; using recommended settings, families were assigned to the correct superfamilies with 0.17% false positive and no false negative assignments.

Bioinformatics assessment of beta-myosin mutations reveals myosin's high sensitivity to mutations.

More than 200 mutations in the beta-myosin gene (MYH7) that cause clinically distinct cardiac and/or skeletal myopathies have been reported, but to date, no comprehensive statistical analysis of these mutations has been performed. As a part of this review, we developed a new interactive database and research tool called MyoMAPR (Myopathic Mutation Analysis Profiler and Repository). We report that the distribution of mutations along the beta-myosin gene is not homogeneous, and that myosin is a highly constrained molecule with an uncommon sensitivity to amino acid substitutions. Increasing knowledge of the characteristics of MH7 mutations may provide a valuable resource for scientists and clinicians studying diagnosis, risk stratification, and treatment of disease associated with these mutations.

Evolution of mammals and their gut microbes.

Mammals are metagenomic in that they are composed of not only their own gene complements but also those of all of their associated microbes. To understand the coevolution of the mammals and their indigenous microbial communities, we conducted a network-based analysis of bacterial 16S ribosomal RNA gene sequences from the fecal microbiota of humans and 59 other mammalian species living in two zoos and in the wild. The results indicate that host diet and phylogeny both influence bacterial diversity, which increases from carnivory to omnivory to herbivory; that bacterial communities codiversified with their hosts; and that the gut microbiota of humans living a modern life-style is typical of omnivorous primates.

Error-correcting barcoded primers for pyrosequencing hundreds of samples in multiplex.

We constructed error-correcting DNA barcodes that allow one run of a massively parallel pyrosequencer to process up to 1,544 samples simultaneously. Using these barcodes we processed bacterial 16S rRNA gene sequences representing microbial communities in 286 environmental samples, corrected 92% of sample assignment errors, and thus characterized nearly as many 16S rRNA genes as have been sequenced to date by Sanger sequencing.