Beltran Rodriguez-Mueller

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.

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 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.

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.

Metagenomic islands of hyperhalophiles: the case of Salinibacter ruber

BackgroundSaturated 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.ResultsThree 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.ConclusionThe 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.

Bacteriophage-encoded shiga toxin gene in atypical bacterial host.

BackgroundContamination 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.ResultsTwenty 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.ConclusionsThe 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.

Microbial community function and biomarker discovery in the human microbiome

Microbial communities carry out the majority of the biochemical activity on the planet, and they play integral roles in processes such as metabolism and immune homeostasis in the human microbiome. Whole genome shotgun sequencing of such communities’ metagenomes is becoming an increasingly feasible complement to obtaining organismal information from taxonomic markers. However, the resultant dataset typically comprises short reads from hundreds of different organisms, making it challenging to assemble and functionally annotate these sequences in the standard manner for single-organism genomes. We describe an alternative to this approach to infer the functional and metabolic potential of a microbial community metagenome by determining whether gene families and pathways are present or absent, as well as their relative abundances, directly from short sequence reads. We validated this methodology using synthetic metagenomes, recovering the presence and abundance of large pathways and of small functional modules with high accuracy. We subsequently applied this approach to the microbial communities of 649 metagenomes drawn from 7 primary body sites on 102 individuals as part of the Human Microbiome Project (HMP), demonstrating the scalability of our methodology and the critical importance of microbial metabolism in the human microbiota. This provided a framework in which to define functional diversity in comparison to organismal ecology, including an example of microbial metabolism linked to specific organisms and to host phenotype (vaginal pH) in the posterior fornix. We provide profiles of 168 functional modules and 196 metabolic pathways that were determined to be specific to one or more niches within the human microbiome, including details of glycosaminoglycan degradation in the gut. Understanding how and why these biomolecular activities differ among environmental conditions or disease phenotypes is, more broadly, one of the central questions addressed by high-throughput biology. We have thus developed the linear discriminant analysis (LDA) effect size algorithm (LEfSe) to discover and explain microbial and functional biomarkers in the human microbiota and other microbiomes. We demonstrate this method to be effective for mining human microbiomes for metagenomic biomarkers associated with mucosal tissues and with different levels of oxygen availability. Similarly, when applied to 16S rRNA gene data from a murine ulcerative colitis gut community, LEfSe confirms the key role played by Bifidobacterium in this disease and suggests the involvement of additional clades, including the Clostridia and Metascardovia. A quantitative validation of LEfSe highlights a lower false positive rate, consistent ranking of biomarker relevance, and concise representations of taxonomic and functional shifts in microbial communities associated with environmental conditions or disease phenotypes. Implementations of both methodologies are available at the Huttenhower laboratory’s website [1,2]. Together, they provide a way to accurately and efficiently characterize microbial metabolic pathways and functional modules directly from high-throughput sequencing reads and, subsequently, to identify organisms, genes or pathways that consistently explain the differences between two or more microbial communities. This has allowed the determination of community roles in the HMP cohort, as well as their niche and population specificity, which we anticipate will be applicable to future metagenomic studies.