An integrated metagenomics pipeline for strain profiling reveals novel patterns of bacterial transmission and biogeography.

We present the Metagenomic Intra-species Diversity Analysis System (MIDAS), which is an integrated computational pipeline for quantifying bacterial species abundance and strain-level genomic variation, including gene content and single-nucleotide polymorphisms (SNPs), from shotgun metagenomes. Our method leverages a database of more than 30,000 bacterial reference genomes that we clustered into species groups. These cover the majority of abundant species in the human microbiome but only a small proportion of microbes in other environments, including soil and seawater. We applied MIDAS to stool metagenomes from 98 Swedish mothers and their infants over one year and used rare SNPs to track strains between hosts. Using this approach, we found that although species compositions of mothers and infants converged over time, strain-level similarity diverged. Specifically, early colonizing bacteria were often transmitted from an infant's mother, while late colonizing bacteria were often transmitted from other sources in the environment and were enriched for spore-formation genes. We also applied MIDAS to 198 globally distributed marine metagenomes and used gene content to show that many prevalent bacterial species have population structure that correlates with geographic location. Strain-level genetic variants present in metagenomes clearly reveal extensive structure and dynamics that are obscured when data are analyzed at a coarser taxonomic resolution.

Metabolic reconstruction for metagenomic data and its application to the human microbiome.

Microbial communities carry out the majority of the biochemical activity on the planet, and they play integral roles in processes including metabolism and immune homeostasis in the human microbiome. Shotgun sequencing of such communities' metagenomes provides information complementary to organismal abundances from taxonomic markers, but the resulting data typically comprise short reads from hundreds of different organisms and are at best challenging to assemble comparably to single-organism genomes. Here, we describe an alternative approach to infer the functional and metabolic potential of a microbial community metagenome. We determined the gene families and pathways present or absent within a community, as well as their relative abundances, directly from short sequence reads. We validated this methodology using a collection of synthetic metagenomes, recovering the presence and abundance both of large pathways and of small functional modules with high accuracy. We subsequently applied this method, HUMAnN, to the microbial communities of 649 metagenomes drawn from seven primary body sites on 102 individuals as part of the Human Microbiome Project (HMP). This provided a means to compare functional diversity and organismal ecology in the human microbiome, and we determined a core of 24 ubiquitously present modules. Core pathways were often implemented by different enzyme families within different body sites, and 168 functional modules and 196 metabolic pathways varied in metagenomic abundance specifically to one or more niches within the microbiome. These included glycosaminoglycan degradation in the gut, as well as phosphate and amino acid transport linked to host phenotype (vaginal pH) in the posterior fornix. An implementation of our methodology is available at http://huttenhower.sph.harvard.edu/humann. This provides a means to accurately and efficiently characterize microbial metabolic pathways and functional modules directly from high-throughput sequencing reads, enabling the determination of community roles in the HMP cohort and in future metagenomic studies.

Microbial to reef scale interactions between the reef-building coral Montastraea annularis and benthic algae.

Competition between reef-building corals and benthic algae is of key importance for reef dynamics. These interactions occur on many spatial scales, ranging from chemical to regional. Using microprobes, 16S rDNA pyrosequencing and underwater surveys, we examined the interactions between the reef-building coral Montastraea annularis and four types of benthic algae. The macroalgae Dictyota bartayresiana and Halimeda opuntia, as well as a mixed consortium of turf algae, caused hypoxia on the adjacent coral tissue. Turf algae were also associated with major shifts in the bacterial communities at the interaction zones, including more pathogens and virulence genes. In contrast to turf algae, interactions with crustose coralline algae (CCA) and M. annularis did not appear to be antagonistic at any scale. These zones were not hypoxic, the microbes were not pathogen-like and the abundance of coral-CCA interactions was positively correlated with per cent coral cover. We propose a model in which fleshy algae (i.e. some species of turf and fleshy macroalgae) alter benthic competition dynamics by stimulating bacterial respiration and promoting invasion of virulent bacteria on corals. This gives fleshy algae a competitive advantage over corals when human activities, such as overfishing and eutrophication, remove controls on algal abundance. Together, these results demonstrate the intricate connections and mechanisms that structure coral reefs.

Bacteriophage-encoded shiga toxin gene in atypical bacterial host.

BACKGROUND: Contamination from fecal bacteria in recreational waters is a major health concern since bacteria capable of causing human disease can be found in animal feces. The Dog Beach area of Ocean Beach in San Diego, California is a beach prone to closures due to high levels of fecal indicator bacteria (FIB). A potential source of these FIB could be the canine feces left behind by owners who do not clean up after their pets. We tested this hypothesis by screening the DNA isolated from canine feces for the bacteriophage-encoded stx gene normally found in the virulent strains of the fecal bacterium Escherichia coli. RESULTS: Twenty canine fecal samples were collected, processed for total and bacterial fraction DNA, and screened by PCR for the stx gene. The stx gene was detected in the total and bacterial fraction DNA of one fecal sample. Bacterial isolates were then cultivated from the stx-positive fecal sample. Eighty nine of these canine fecal bacterial isolates were screened by PCR for the stx gene. The stx gene was detected in five of these isolates. Sequencing and phylogenetic analyses of 16S rRNA gene PCR products from the canine fecal bacterial isolates indicated that they were Enterococcus and not E. coli. CONCLUSIONS: The bacteriophage-encoded stx gene was found in multiple species of bacteria cultivated from canine fecal samples gathered at the shoreline of the Dog Beach area of Ocean Beach in San Diego, California. The canine fecal bacteria carrying the stx gene were not the typical E. coli host and were instead identified through phylogenetic analyses as Enterococcus. This suggests a large degree of horizontal gene transfer of exotoxin genes in recreational waters.

Metagenomic islands of hyperhalophiles: the case of Salinibacter ruber.

BACKGROUND: Saturated brines are extreme environments of low diversity. Salinibacter ruber is the only bacterium that inhabits this environment in significant numbers. In order to establish the extent of genetic diversity in natural populations of this microbe, the genomic sequence of reference strain DSM 13855 was compared to metagenomic fragments recovered from climax saltern crystallizers and obtained with 454 sequencing technology. This kind of analysis reveals the presence of metagenomic islands, i.e. highly variable regions among the different lineages in the population. RESULTS: Three regions of the sequenced isolate were scarcely represented in the metagenome thus appearing to vary among co-occurring S. ruber cells. These metagenomic islands showed evidence of extensive genomic corruption with atypically low GC content, low coding density, high numbers of pseudogenes and short hypothetical proteins. A detailed analysis of island gene content showed that the genes in metagenomic island 1 code for cell surface polysaccharides. The strain-specific genes of metagenomic island 2 were found to be involved in biosynthesis of cell wall polysaccharide components. Finally, metagenomic island 3 was rich in DNA related enzymes. CONCLUSION: The genomic organisation of S. ruber variable genomic regions showed a number of convergences with genomic islands of marine microbes studied, being largely involved in variable cell surface traits. This variation at the level of cell envelopes in an environment devoid of grazing pressure probably reflects a global strategy of bacteria to escape phage predation.

The GAAS metagenomic tool and its estimations of viral and microbial average genome size in four major biomes.

Metagenomic studies characterize both the composition and diversity of uncultured viral and microbial communities. BLAST-based comparisons have typically been used for such analyses; however, sampling biases, high percentages of unknown sequences, and the use of arbitrary thresholds to find significant similarities can decrease the accuracy and validity of estimates. Here, we present Genome relative Abundance and Average Size (GAAS), a complete software package that provides improved estimates of community composition and average genome length for metagenomes in both textual and graphical formats. GAAS implements a novel methodology to control for sampling bias via length normalization, to adjust for multiple BLAST similarities by similarity weighting, and to select significant similarities using relative alignment lengths. In benchmark tests, the GAAS method was robust to both high percentages of unknown sequences and to variations in metagenomic sequence read lengths. Re-analysis of the Sargasso Sea virome using GAAS indicated that standard methodologies for metagenomic analysis may dramatically underestimate the abundance and importance of organisms with small genomes in environmental systems. Using GAAS, we conducted a meta-analysis of microbial and viral average genome lengths in over 150 metagenomes from four biomes to determine whether genome lengths vary consistently between and within biomes, and between microbial and viral communities from the same environment. Significant differences between biomes and within aquatic sub-biomes (oceans, hypersaline systems, freshwater, and microbialites) suggested that average genome length is a fundamental property of environments driven by factors at the sub-biome level. The behavior of paired viral and microbial metagenomes from the same environment indicated that microbial and viral average genome sizes are independent of each other, but indicative of community responses to stressors and environmental conditions.

Metagenomic analysis of stressed coral holobionts.

The coral holobiont is the community of metazoans, protists and microbes associated with scleractinian corals. Disruptions in these associations have been correlated with coral disease, but little is known about the series of events involved in the shift from mutualism to pathogenesis. To evaluate structural and functional changes in coral microbial communities, Porites compressa was exposed to four stressors: increased temperature, elevated nutrients, dissolved organic carbon loading and reduced pH. Microbial metagenomic samples were collected and pyrosequenced. Functional gene analysis demonstrated that stressors increased the abundance of microbial genes involved in virulence, stress resistance, sulfur and nitrogen metabolism, motility and chemotaxis, fatty acid and lipid utilization, and secondary metabolism. Relative changes in taxonomy also demonstrated that coral-associated microbiota (Archaea, Bacteria, protists) shifted from a healthy-associated coral community (e.g. Cyanobacteria, Proteobacteria and the zooxanthellae Symbiodinium) to a community (e.g. Bacteriodetes, Fusobacteria and Fungi) of microbes often found on diseased corals. Additionally, low-abundance Vibrio spp. were found to significantly alter microbiome metabolism, suggesting that the contribution of a just a few members of a community can profoundly shift the health status of the coral holobiont.

Metagenomic analysis indicates that stressors induce production of herpes-like viruses in the coral Porites compressa.

During the last several decades corals have been in decline and at least one-third of all coral species are now threatened with extinction. Coral disease has been a major contributor to this threat, but little is known about the responsible pathogens. To date most research has focused on bacterial and fungal diseases; however, viruses may also be important for coral health. Using a combination of empirical viral metagenomics and real-time PCR, we show that Porites compressa corals contain a suite of eukaryotic viruses, many related to the Herpesviridae. This coral-associated viral consortium was found to shift in response to abiotic stressors. In particular, when exposed to reduced pH, elevated nutrients, and thermal stress, the abundance of herpes-like viral sequences rapidly increased in 2 separate experiments. Herpes-like viral sequences were rarely detected in apparently healthy corals, but were abundant in a majority of stressed samples. In addition, surveys of the Nematostella and Hydra genomic projects demonstrate that even distantly related Cnidarians contain numerous herpes-like viral genes, likely as a result of latent or endogenous viral infection. These data support the hypotheses that corals experience viral infections, which are exacerbated by stress, and that herpes-like viruses are common in Cnidarians.