Search

Oxygen minimum zones (OMZs) are oceanographic features that affect ocean productivity and biodiversity, and contribute to ocean nitrogen loss and greenhouse gas emissions. Here we describe the viral communities associated with the Eastern Tropical South Pacific (ETSP) OMZ off Iquique, Chile for the first time through abundance estimates and viral metagenomic analysis. The viral-to-microbial ratio (VMR) in the ETSP OMZ fluctuated in the oxycline and declined in the anoxic core to below one on several occasions. The number of viral genotypes (unique genomes as defined by sequence assembly) ranged from 2040 at the surface to 98 in the oxycline, which is the lowest viral diversity recorded to date in the ocean. Within the ETSP OMZ viromes, only 4.95% of genotypes were shared between surface and anoxic core viromes using reciprocal BLASTn sequence comparison. ETSP virome comparison with surface marine viromes (Sargasso Sea, Gulf of Mexico, Kingman Reef, Chesapeake Bay) revealed a dissimilarity of ETSP OMZ viruses to those from other oceanic regions. From the 1.4 million non-redundant DNA sequences sampled within the altered oxygen conditions of the ETSP OMZ, more than 97.8% were novel. Of the average 3.2% of sequences that showed similarity to the SEED non-redundant database, phage sequences dominated the surface viromes, eukaryotic virus sequences dominated the oxycline viromes, and phage sequences dominated the anoxic core viromes. The viral community of the ETSP OMZ was characterized by fluctuations in abundance, taxa and diversity across the oxygen gradient. The ecological significance of these changes was difficult to predict; however, it appears that the reduction in oxygen coincides with an increased shedding of eukaryotic viruses in the oxycline, and a shift to unique viral genotypes in the anoxic core.

Yeasts, as a group of single-celled fungi, are widely distributed in nature and play important roles in biotechnological applications. However, how the types of wastewater and the treatment processes influence yeast populations is not clear. In this study, both cultivation and culture-independent methods were used to investigate the distribution and diversity of yeasts in three typical full-scale plants processing biopharmaceutical, papermaking and municipal wastewater. Cultivable yeasts were very abundant ranging from 10 2 to 10 5 cfu g −1 sludge, and highly diverse with 48 taxons belonging to 21 different genera, thus exceeding the yeast richness reported for marine and sugar-rich soil habitats. Genera Rhodotorula , Candida , Trichosporon , Pichia and some unidentified Ascomycetes were the most frequent populations cultivated. However, the compositions of yeast community structures in the plants were dissimilar and were shaped primarily by the type of wastewater treated but also by processing conditions. Culturing yeasts is a powerful way to complement culture-independent approaches to study their diversity, since significant more yeast species, especially quite a few possible novel species were recorded and isolated by cultivation method.

Summary Viruses have been detected in the different stages of wastewater treatment plants (WWTPs) at concentrations of 108–1010 ml−1 of virus-like particles (VLPs), 10–1000 times higher than in natural aquatic environments, suggesting that WWTPs can be considered as an important reservoir and source of viruses. This study revealed novel diversity and function with the DNA viral communities in the influent, activated sludge, anaerobic digester, and effluent of a domestic WWTP using metagenomics. WWTP was a very specific environment, with less than 5% of the > 936 000 metagenomic sequences obtained (∼70–119 Mbp per sample) similar to sequences present in other environmental viromes. Many viruses found in the WWTP were novel, resulting in only 82% were shared. These data suggested that VLPs of most viral types could be present between 1 and 30 days in the process before they were discharged. Viruses in WWTP and the discharged ones can have potential impacts on the functioning of the wastewater treatment system and on the dynamics of microbial community in the surrounding aquatic environments respectively.

Viral diversity and life cycles are poorly understood in the human gut and other body habitats. Phages and their encoded functions may provide informative signatures of a human microbiota and of microbial community responses to various disturbances, and may indicate whether community health or dysfunction is manifest after apparent recovery from a disease or therapeutic intervention. Here we report sequencing of the viromes (metagenomes) of virus-like particles isolated from faecal samples collected from healthy adult female monozygotic twins and their mothers at three time points over a one-year period. We compared these data sets with data sets of sequenced bacterial 16S ribosomal RNA genes and total-faecal-community DNA. Co-twins and their mothers share a significantly greater degree of similarity in their faecal bacterial communities than do unrelated individuals. In contrast, viromes are unique to individuals regardless of their degree of genetic relatedness. Despite remarkable interpersonal variations in viromes and their encoded functions, intrapersonal diversity is very low, with >95% of virotypes retained over the period surveyed, and with viromes dominated by a few temperate phages that exhibit remarkable genetic stability. These results indicate that a predatory viral-microbial dynamic, manifest in a number of other characterized environmental ecosystems, is notably absent in the very distal intestine.

The human oropharynx is a reservoir for many potential pathogens, including streptococcal species that cause endocarditis. Although oropharyngeal microbes have been well described, viral communities are essentially uncharacterized. We conducted a metagenomic study to determine the composition of oropharyngeal DNA viral communities (both phage and eukaryotic viruses) in healthy individuals and to evaluate oropharyngeal swabs as a rapid method for viral detection. Viral DNA was extracted from 19 pooled oropharyngeal swabs and sequenced. Viral communities consisted almost exclusively of phage, and complete genomes of several phage were recovered, including Escherichia coli phage T3, Propionibacterium acnes phage PA6, and Streptococcus mitis phage SM1. Phage relative abundances changed dramatically depending on whether samples were chloroform treated or filtered to remove microbial contamination. pblA and pblB genes of phage SM1 were detected in the metagenomes. pblA and pblB mediate the attachment of S. mitis to platelets and play a significant role in S. mitis virulence in the endocardium, but have never previously been detected in the oral cavity. These genes were also identified in salivary metagenomes from three individuals at three time points and in individual saliva samples by PCR. Additionally, we demonstrate that phage SM1 can be induced by commonly ingested substances. Our results indicate that the oral cavity is a reservoir for pblA and pblB genes and for phage SM1 itself. Further studies will determine the association between pblA and pblB genes in the oral cavity and the risk of endocarditis.

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.

[Background] Viruses have unique properties, small genome and regions of high similarity, whose effects on metagenomic assemblies have not been characterized so far. This study uses diverse in silico simulated viromes to evaluate how extensively genomes can be assembled using different sequencing platforms and assemblers. Further, it investigates the suitability of different methods to estimate viral diversity in metagenomes., [Results] We created in silico metagenomes mimicking various platforms at different sequencing depths. The CLC assembler revealed subpar compared to IDBA_UD and CAMERA, which are metagenomic-specific. Up to a saturation point, Illumina platforms proved more capable of reconstructing large portions of viral genomes compared to 454. Read length was an important factor for limiting chimericity, while scaffolding marginally improved contig length and accuracy. The genome length of the various viruses in the metagenomes did not significantly affect genome reconstruction, but the co-existence of highly similar genomes was detrimental. When evaluating diversity estimation tools, we found that PHACCS results were more accurate than those from CatchAll and clustering, which were both orders of magnitude above expected., [Conclusions] Assemblers designed specifically for the analysis of metagenomes should be used to facilitate the creation of high-quality long contigs. Despite the high coverage possible, scientists should not expect to always obtain complete genomes, because their reconstruction may be hindered by co-existing species bearing highly similar genomic regions. Further development of metagenomics-oriented assemblers may help bypass these limitations in future studies. Meanwhile, the lack of fully reconstructed communities keeps methods to estimate viral diversity relevant. While none of the three methods tested had absolute precision, only PHACCS was deemed suitable for comparative studies. © 2015 Aguirre de Cárcer et al.. licensee BioMed Central Ltd., This work was funded in part by the Spanish Ministry of Science and Innovation grant CTM2011-15091-E/ANT. Daniel Aguirre de Cárcer was supported by the Marie Curie International Incoming Fellow grant PIIF-GA- 2012-328287. Florent Angly was supported by the Australian Research Council’s Discovery Early Career Research Award DE120101213, We acknowledge support by the CSIC Open Access Publication Initiative through its Unit of Information Resources for Research (URICI)

Increasing atmospheric partial pressure of CO2 (pCO2) is a major threat to coral reefs, but some argue that the threat is mitigated by factors such as the variability in the response of coral calcification to acidification, differences in bleaching susceptibility, and the potential for rapid adaptation to anthropogenic warming. However the evidence for these mitigating factors tends to involve experimental studies on corals, as opposed to coral reefs, and rarely includes the influence of multiple variables (e.g., temperature and acidification) within regimes that include diurnal and seasonal variability. Here, we demonstrate that the inclusion of all these factors results in the decalcification of patch-reefs under business-as-usual scenarios and reduced, although positive, calcification under reduced-emission scenarios. Primary productivity was found to remain constant across all scenarios, despite significant bleaching and coral mortality under both future scenarios. Daylight calcification decreased and nocturnal decalcification increased sharply from the preindustrial and control conditions to the future scenarios of low (reduced emissions) and high (business-as-usual) increases in pCO2. These changes coincided with deeply negative carbonate budgets, a shift toward smaller carbonate sediments, and an increase in the abundance of sediment microbes under the business-as-usual emission scenario. Experimental coral reefs demonstrated highest net calcification rates and lowest rates of coral mortality under preindustrial conditions, suggesting that reef processes may not have been able to keep pace with the relatively minor environmental changes that have occurred during the last century. Taken together, our results have serious implications for the future of coral reefs under business-as-usual environmental changes projected for the coming decades and century.

Diuron is a herbicide commonly used in agricultural areas where excess application causes it to leach into rivers, reach sensitive marine environments like the Great Barrier Reef (GBR) lagoon and pose risks to marine life. To investigate the impact of diuron on whole prokaryotic communities that underpin the marine food web and are integral to coral reef health, GBR lagoon water was incubated with diuron at environmentally-relevant concentration (8 µg/L), and sequenced at specific time points over the following year. 16S rRNA gene amplicon profiling revealed no significant short- or long-term effect of diuron on microbiome structure. The relative abundance of prokaryotic phototrophs was not significantly altered by diuron, which suggests that they were largely tolerant at this concentration. Assembly of a metagenome derived from waters sampled at a similar location in the GBR lagoon did not reveal the presence of mutations in the cyanobacterial photosystem that could explain diuron tolerance. However, resident phages displayed several variants of this gene and could potentially play a role in tolerance acquisition. Slow biodegradation of diuron was reported in the incubation flasks, but no correlation with the relative abundance of heterotrophs was evident. Analysis of metagenomic reads supports the hypothesis that previously uncharacterized hydrolases carried by low-abundance species may mediate herbicide degradation in the GBR lagoon. Overall, this study offers evidence that pelagic phototrophs of the GBR lagoon may be more tolerant of diuron than other tropical organisms, and that heterotrophs in the microbial seed bank may have the potential to degrade diuron and alleviate local anthropogenic stresses to inshore GBR ecosystems.

The role of microorganisms in maintaining coral reef health is increasingly recognized. Riverine floodwater containing herbicides and excess nutrients from fertilizers compromises water quality in the inshore Great Barrier Reef (GBR), with unknown consequences for planktonic marine microbial communities and thus coral reefs. In this baseline study, inshore GBR microbial communities were monitored along a 124 km long transect between 2011 and 2013 using 16S rRNA gene amplicon sequencing. Members of the bacterial orders Rickettsiales (e.g., Pelagibacteraceae) and Synechococcales (e.g.,Prochlorococcus), and of the archaeal class Marine Group II were prevalent in all samples, exhibiting a clear seasonal dynamics. Microbial communities near the Tully river mouth included a mixture of taxa from offshore marine sites and from the river system. The environmental parameters collected could be summarized into four groups, represented by salinity, rainfall, temperature and water quality, that drove the composition of microbial communities. During the wet season, lower salinity and a lower water quality index resulting from higher river discharge corresponded to increases in riverine taxa at sites near the river mouth. Particularly large, transient changes in microbial community structure were seen during the extreme wet season 2010–11, and may be partially attributed to the effects of wind and waves, which resuspend sediments and homogenize the water column in shallow near-shore regions. This work shows that anthropogenic floodwaters and other environmental parameters work in conjunction to drive the spatial distribution of microorganisms in the GBR lagoon, as well as their seasonal and daily dynamics.

Viruses have been detected in the different stages of wastewater treatment plants (WWTPs) at concentrations of 10(8) -10(10) ml(-1) of virus-like particles (VLPs), 10-1000 times higher than in natural aquatic environments, suggesting that WWTPs can be considered as an important reservoir and source of viruses. This study revealed novel diversity and function with the DNA viral communities in the influent, activated sludge, anaerobic digester, and effluent of a domestic WWTP using metagenomics. WWTP was a very specific environment, with less than 5% of the 936 000 metagenomic sequences obtained (∼70-119 Mbp per sample) similar to sequences present in other environmental viromes. Many viruses found in the WWTP were novel, resulting in only 5-20% of the reads being phylogenetically or functionally assigned. DNA metabolism was observed as the most abundant function with DNA methylase detected at levels 4.2-fold higher than other published viromes, while carbohydrate and amino acids metabolisms were 3.7- and 4.2-fold less abundant respectively. These specific aspects of the WWTP community functions are likely due to high biomass concentration, turnover rate and microbial activity in WWTPs, and likely include mechanisms that help viruses increase their infectivity. Among ∼500 genotypes estimated in individual WWTP viromes, 82% were shared. These data suggested that VLPs of most viral types could be present between 1 and 30 days in the process before they were discharged. Viruses in WWTP and the discharged ones can have potential impacts on the functioning of the wastewater treatment system and on the dynamics of microbial community in the surrounding aquatic environments respectively.

Viruses are the most abundant biological entities on the planet and play an important role in balancing microbes within an ecosystem and facilitating horizontal gene transfer. Although bacteriophages are abundant in rumen environments, little is known about the types of viruses present or their interaction with the rumen microbiome. We undertook random pyrosequencing of virus-enriched metagenomes (viromes) isolated from bovine rumen fluid and analysed the resulting data using comparative metagenomics. A high level of diversity was observed with up to 28,000 different viral genotypes obtained from each environment. The majority (~78%) of sequences did not match any previously described virus. Prophages outnumbered lytic phages approximately 2:1 with the most abundant bacteriophage and prophage types being associated with members of the dominant rumen phyla (Firmicutes and Proteobacteria). Metabolic profiling based on SEED subsystems revealed an enrichment of sequences with putative functional roles in DNA and protein metabolism, but a surprisingly low proportion of sequences assigned to carbohydrate and amino acid metabolism. We expanded our analysis to include previously described metagenomic data and 14 reference genomes. Clustered regularly interspaced short palindromic repeats (CRISPR) were detected in most of the microbial genomes, suggesting previous interactions between viral and microbial communities.

A Modular Open-Source Assembler (AMOS) was designed to offer a modular approach to genome assembly. AMOS includes a wide range of tools for assembly, including the lightweight de novo assemblers Minimus and Minimo, and Bambus 2, a robust scaffolder able to handle metagenomic and polymorphic data. This protocol describes how to configure and use AMOS for the assembly of Next Generation sequence data. Additionally, we provide three tutorial examples that include bacterial, viral, and metagenomic datasets with specific tips for improving assembly quality.

Summary: The development of bioinformatic solutions for microbial ecology in Perl is limited by the lack of modules to represent and manipulate microbial community profiles from amplicon and meta-omics studies. Here we introduce Bio-Community, an open-source, collaborative toolkit that extends BioPerl. Bio-Community interfaces with commonly used programs using various file formats, including BIOM, and provides operations such as rarefaction and taxonomic summaries. Bio-Community will help bioinformaticians to quickly piece together custom analysis pipelines and develop novel software. Availability an implementation: Bio-Community is cross-platform Perl code available from http://search.cpan.org/dist/Bio-Community under the Perl license. A readme file describes software installation and how to contribute. Contact: f.angly@uq.edu.au Supplementary information: Supplementary data are available at Bioinformatics online

The species composition and metabolic potential of microbial and viral communities are predictable and stable for most ecosystems. This apparent stability contradicts theoretical models as well as the viral-microbial dynamics observed in simple ecosystems, both of which show Kill-the-Winner behavior causing cycling of the dominant taxa. Microbial and viral metagenomes were obtained from four human-controlled aquatic environments at various time points separated by one day to >1 year. These environments were maintained within narrow geochemical bounds and had characteristic species composition and metabolic potentials at all time points. However, underlying this stability were rapid changes at the fine-grained level of viral genotypes and microbial strains. These results suggest a model wherein functionally redundant microbial and viral taxa are cycling at the level of viral genotypes and virus-sensitive microbial strains. Microbial taxa, viral taxa, and metabolic function persist over time in stable ecosystems and both communities fluctuate in a Kill-the-Winner manner at the level of viral genotypes and microbial strains. The ISME Journal (2010) 4, 739-751; doi:10.1038/ismej.2010.1; published online 11 February 2010

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., Author Summary Metagenomics uses DNA or RNA sequences isolated directly from the environment to determine what viruses or microorganisms exist in natural communities and what metabolic activities they encode. Typically, metagenomic sequences are compared to annotated sequences in public databases using the BLAST search tool. Our methods, implemented in the Genome relative Abundance and Average Size (GAAS) software, improve the way BLAST searches are processed to estimate the taxonomic composition of communities and their average genome length. GAAS provides a more accurate picture of community composition by correcting for a systematic sampling bias towards larger genomes, and is useful in situations where organisms with small genomes are abundant, such as disease outbreaks caused by small RNA viruses. Microbial average genome length relates to environmental complexity and the distribution of genome lengths describes community diversity. A study of the average genome length of viruses and microorganisms in four different biomes using GAAS on 169 metagenomes showed significantly different average genome sizes between biomes, and large variability within biomes as well. This also revealed that microbial and viral average genome sizes in the same environment are independent of each other, which reflects the different ways that microorganisms and viruses respond to stress and environmental conditions.

Roseophage SIO1 is a lytic marine phage that infects Roseobacter SIO67, a member of the Roseobacter clade of near-shore alphaproteobacteria. Roseophage SIO1 was first isolated in 1989 and sequenced in 2000. We have re-sequenced and re-annotated the original isolate. Our current annotation could only assign functions to seven additional open reading frames, indicating that, despite the advances in bioinformatics tools and increased genomic resources, we are still far from being able to translate phage genomic sequences into biological functions. In 2001, we isolated four new strains of Roseophage SIO1 from California near-shore locations. The genomes of all four were sequenced and compared against the original Roseophage SIO1 isolated in 1989. A high degree of conservation was evident across all five genomes; comparisons at the nucleotide level yielded an average 97% identity. The observed differences were clustered in protein-encoding regions and were mostly synonymous. The one strain that was found to possess an expanded host range also showed notable changes in putative tail protein-coding regions. Despite the possibly rapid evolution of phage and the mostly uncharacterized diversity found in viral metagenomic data sets, these findings indicate that viral genomes such as the genome of SIO1-like Roseophages can be stably maintained over ecologically significant time and distance (i.e. over a decade and approximately 50 km).

The human respiratory tract is constantly exposed to a wide variety of viruses, microbes and inorganic particulates from environmental air, water and food. Physical characteristics of inhaled particles and airway mucosal immunity determine which viruses and microbes will persist in the airways. Here we present the first metagenomic study of DNA viral communities in the airways of diseased and non-diseased individuals. We obtained sequences from sputum DNA viral communities in 5 individuals with cystic fibrosis (CF) and 5 individuals without the disease. Overall, diversity of viruses in the airways was low, with an average richness of 175 distinct viral genotypes. The majority of viral diversity was uncharacterized. CF phage communities were highly similar to each other, whereas Non-CF individuals had more distinct phage communities, which may reflect organisms in inhaled air. CF eukaryotic viral communities were dominated by a few viruses, including human herpesviruses and retroviruses. Functional metagenomics showed that all Non-CF viromes were similar, and that CF viromes were enriched in aromatic amino acid metabolism. The CF metagenomes occupied two different metabolic states, probably reflecting different disease states. There was one outlying CF virome which was characterized by an over-representation of Guanosine-5'-triphosphate,3'-diphosphate pyrophosphatase, an enzyme involved in the bacterial stringent response. Unique environments like the CF airway can drive functional adaptations, leading to shifts in metabolic profiles. These results have important clinical implications for CF, indicating that therapeutic measures may be more effective if used to change the respiratory environment, as opposed to shifting the taxonomic composition of resident microbiota.

P>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 sequencing of DNA viruses from the feces of a healthy week-old infant revealed a viral community with extremely low diversity. The identifiable sequences were dominated by phages, which likely influence the diversity and abundance of co-occurring microbes. The most abundant fecal viral sequences did not originate from breast milk or formula, suggesting a non-dietary initial source of viruses. Certain sequences were stable in the infant's gut over the first 3 months of life, but microarray experiments demonstrated that the overall viral community composition changed dramatically between 1 and 2 weeks of age.

Microbial activities shape the biogeochemistry of the planet and macroorganism health. Determining the metabolic processes performed by microbes is important both for understanding and for manipulating ecosystems (for example, disruption of key processes that lead to disease, conservation of environmental services, and so on). Describing microbial function is hampered by the inability to culture most microbes and by high levels of genomic plasticity. Metagenomic approaches analyse microbial communities to determine the metabolic processes that are important for growth and survival in any given environment. Here we conduct a metagenomic comparison of almost 15 million sequences from 45 distinct microbiomes and, for the first time, 42 distinct viromes and show that there are strongly discriminatory metabolic profiles across environments. Most of the functional diversity was maintained in all of the communities, but the relative occurrence of metabolisms varied, and the differences between metagenomes predicted the biogeochemical conditions of each environment. The magnitude of the microbial metabolic capabilities encoded by the viromes was extensive, suggesting that they serve as a repository for storing and sharing genes among their microbial hosts and influence global evolutionary and metabolic processes.

Metagenomic analyses suggest that the rank–abundance curve for marine phage communities follows a power law distribution. A new type of power law dependence based on a simple model in which a modified version of Lotka–Volterra predator–prey dynamics is sampled uniformly in time is presented. Biologically, the model embodies a kill the winner hypothesis and a neutral evolution hypothesis. The model can match observed power law distributions and uses very few parameters that are readily identifiable and characterize phage ecosystems. The model makes new untested predictions: (1) it is unlikely that the most abundant phage genotype will be the same at different time points and (2) the long-term decay of isolated phage populations follows a power law.

We introduce Grinder (http://sourceforge.net/ projects/biogrinder/), an open-source bioinformatic tool to simulate amplicon and shotgun (genomic, metagenomic, transcriptomic and metatranscriptomic) datasets from reference sequences. This is the first tool to simulate amplicon datasets (e.g. 16S rRNA) widely used by microbial ecologists. Grinder can create sequence libraries with a specific community structure, α and β diversities and experimental biases (e.g. chimeras, gene copy number variation) for commonly used sequencing platforms. This versatility allows the creation of simple to complex read datasets necessary for hypothesis testing when developing bioinformatic software, benchmarking existing tools or designing sequence-based experiments. Grinder is particularly useful for simulating clinical or environmental microbial communities and complements the use of in vitro mock communities.

Nature 452, 629–632 (2008) In this Letter, functional diversity and evenness were calculated using log10, but it has been drawn to our attention that the calculations should have been made with natural logs (Table 2). Recalculation does not change the relative levels of diversity; however, this reanalysis showed that a mistake was made in the original evenness calculation.

BackgroundThe complex microbiome of the ceca of chickens plays an important role in nutrient utilization, growth and well-being of these animals. Since we have a very limited understanding of the capabilities of most species present in the cecum, we investigated the role of the microbiome by comparative analyses of both the microbial community structure and functional gene content using random sample pyrosequencing. The overall goal of this study was to characterize the chicken cecal microbiome using a pathogen-free chicken and one that had been challenged with Campylobacter jejuni.Methodology/principal findingsComparative metagenomic pyrosequencing was used to generate 55,364,266 bases of random sampled pyrosequence data from two chicken cecal samples. SSU rDNA gene tags and environmental gene tags (EGTs) were identified using SEED subsystems-based annotations. The distribution of phylotypes and EGTs detected within each cecal sample were primarily from the Firmicutes, Bacteroidetes and Proteobacteria, consistent with previous SSU rDNA libraries of the chicken cecum. Carbohydrate metabolism and virulence genes are major components of the EGT content of both of these microbiomes. A comparison of the twelve major pathways in the SEED Virulence Subsystem (metavirulome) represented in the chicken cecum, mouse cecum and human fecal microbiomes showed that the metavirulomes differed between these microbiomes and the metavirulomes clustered by host environment. The chicken cecum microbiomes had the broadest range of EGTs within the SEED Conjugative Transposon Subsystem, however the mouse cecum microbiomes showed a greater abundance of EGTs in this subsystem. Gene assemblies (32 contigs) from one microbiome sample were predominately from the Bacteroidetes, and seven of these showed sequence similarity to transposases, whereas the remaining sequences were most similar to those from catabolic gene families.Conclusion/significanceThis analysis has demonstrated that mobile DNA elements are a major functional component of cecal microbiomes, thus contributing to horizontal gene transfer and functional microbiome evolution. Moreover, the metavirulomes of these microbiomes appear to associate by host environment. These data have implications for defining core and variable microbiome content in a host species. Furthermore, this suggests that the evolution of host specific metavirulomes is a contributing factor in disease resistance to zoonotic pathogens.

Phages, viruses that infect prokaryotes, are the most abundant microbes in the world. A major limitation to studying these viruses is the difficulty of cultivating the appropriate prokaryotic hosts. One way around this limitation is to directly clone and sequence shotgun libraries of uncultured viral communities (i.e., metagenomic analyses). PHACCS http://phage.sdsu.edu/phaccs , Phage Communities from Contig Spectrum, is an online bioinformatic tool to assess the biodiversity of uncultured viral communities. PHACCS uses the contig spectrum from shotgun DNA sequence assemblies to mathematically model the structure of viral communities and make predictions about diversity. PHACCS builds models of possible community structure using a modified Lander-Waterman algorithm to predict the underlying contig spectrum. PHACCS finds the most appropriate structure model by optimizing the model parameters until the predicted contig spectrum is as close as possible to the experimental one. This model is the basis for making estimates of uncultured viral community richness, evenness, diversity index and abundance of the most abundant genotype. PHACCS analysis of four different environmental phage communities suggests that the power law is an important rank-abundance form to describe uncultured viral community structure. The estimates support the fact that the four phage communities were extremely diverse and that phage community biodiversity and structure may be correlated with that of their hosts.

Viruses, and more particularly phages (viruses that infect bacteria), represent one of the most abundant living entities in aquatic and terrestrial environments. The biogeography of phages has only recently been investigated and so far reveals a cosmopolitan distribution of phage genetic material (or genotypes). Here we address this cosmopolitan distribution through the analysis of phage communities in modern microbialites, the living representatives of one of the most ancient life forms on Earth. On the basis of a comparative metagenomic analysis of viral communities associated with marine (Highborne Cay, Bahamas) and freshwater (Pozas Azules II and Rio Mesquites, Mexico) microbialites, we show that some phage genotypes are geographically restricted. The high percentage of unknown sequences recovered from the three metagenomes (>97%), the low percentage similarities with sequences from other environmental viral (n = 42) and microbial (n = 36) metagenomes, and the absence of viral genotypes shared among microbialites indicate that viruses are genetically unique in these environments. Identifiable sequences in the Highborne Cay metagenome were dominated by single-stranded DNA microphages that were not detected in any other samples examined, including sea water, fresh water, sediment, terrestrial, extreme, metazoan-associated and marine microbial mats. Finally, a marine signature was present in the phage community of the Pozas Azules II microbialites, even though this environment has not been in contact with the ocean for tens of millions of years. Taken together, these results prove that viruses in modern microbialites display biogeographical variability and suggest that they may be derived from an ancient community.

Background: Culture-independent molecular surveys targeting conserved marker genes, most notably 16S rRNA, to assess microbial diversity remain semi-quantitative due to variations in the number of gene copies between species. Results: Based on 2,900 sequenced reference genomes, we show that 16S rRNA gene copy number (GCN) is strongly linked to microbial phylogenetic taxonomy, potentially under-representing Archaea in amplicon microbial profiles. Using this relationship, we inferred the GCN of all bacterial and archaeal lineages in the Greengenes database within a phylogenetic framework. We created CopyRighter, new software which uses these estimates to correct 16S rRNA amplicon microbial profiles and associated quantitative (q)PCR total abundance. CopyRighter parses microbial profiles and, because GCN estimates are pre-computed for all taxa in the reference taxonomy, rapidly corrects GCN bias. Software validation with in silico and in vitro mock communities indicated that GCN correction results in more accurate estimates of microbial relative abundance and improves the agreement between metagenomic and amplicon profiles. Analyses of human-associated and anaerobic digester microbiomes illustrate that correction makes tangible changes to estimates of qPCR total abundance, α and β diversity, and can significantly change biological interpretation. For example, human gut microbiomes from twins were reclassified into three rather than two enterotypes after GCN correction. Conclusions: The CopyRighter bioinformatic tools permits rapid correction of GCN in microbial surveys, resulting in improved estimates of microbial abundance, α and β diversity.