Bacteriophages have key roles in microbial communities, to a large extent shaping the taxonomic and functional composition of the microbiome, but data on the connections between phage diversity and the composition of communities are scarce. Using taxon-specific marker genes, we identified and monitored 20 viral taxa in 252 human gut metagenomic samples, mostly at the level of genera. On average, five phage taxa were identified in each sample, with up to three of these being highly abundant. The abundances of most phage taxa vary by up to four orders of magnitude between the samples, and several taxa that are highly abundant in some samples are absent in others. Significant correlations exist between the abundances of some phage taxa and human host metadata: for example, 'Group 936 lactococcal phages' are more prevalent and abundant in Danish samples than in samples from Spain or the United States of America. Quantification of phages that exist as integrated prophages revealed that the abundance profiles of prophages are highly individual-specific and remain unique to an individual over a 1-year time period, and prediction of prophage lysis across the samples identified hundreds of prophages that are apparently active in the gut and vary across the samples, in terms of presence and lytic state. Finally, a prophage-host network of the human gut was established and includes numerous novel host-phage associations.

Viruses are the most abundant biological entities on earth and encompass a vast amount of genetic diversity. The recent rapid increase in the number of sequenced viral genomes has created unprecedented opportunities for gaining new insight into the structure and evolution of the virosphere. Here we present an update of the Phage Orthologous Groups (POGs), a collection of 4,542 clusters of orthologous genes from bacteriophages that now also includes viruses infecting archaea and encompasses more than 1,000 distinct virus genomes. Analysis of this expanded dataset shows that the number of POGs keeps growing without saturation and that a substantial majority of the POGs remain specific to viruses, lacking homologues in prokaryotic cells, outside of known proviruses. Thus, the great majority of virus genes apparently remains to be discovered. A complementary observation is that numerous viral genomes remain poorly if at all covered by POGs. The genome coverage by POGs is expected to increase as more genomes are sequenced. Taxon-specific, single-copy signature genes, that are not observed in prokaryotic genomes outside of detected proviruses, were identified for two-thirds of the 57 taxa (those with genomes available from at least 3 distinct viruses, with half of these present in all members of the respective taxon. These signatures can be used to specifically identify the presence and quantify the abundance of viruses from particular taxa in metagenomic samples and thus gain new insights into the ecology and evolution of viruses in relation to their hosts.

Whereas large-scale efforts have rapidly advanced the understanding and practical impact of human genomic variation, the practical impact of variation is largely unexplored in the human microbiome. We therefore developed a framework for metagenomic variation analysis and applied it to 252 faecal metagenomes of 207 individuals from Europe and North America. Using 7.4 billion reads aligned to 101 reference species, we detected 10.3 million single nucleotide polymorphisms (SNPs), 107,991 short insertions/deletions, and 1,051 structural variants. The average ratio of non-synonymous to synonymous polymorphism rates of 0.11 was more variable between gut microbial species than across human hosts. Subjects sampled at varying time intervals exhibited individuality and temporal stability of SNP variation patterns, despite considerable composition changes of their gut microbiota. This indicates that individual-specific strains are not easily replaced and that an individual might have a unique metagenomic genotype, which may be exploitable for personalized diet or drug intake.

The breakthrough of next generation sequencing-technologies has enabled large-scale studies of natural microbial communities and the 16S rRNA genes have been widely used as a phylogenetic marker to study community structure. However, major limitations of this approach are that neither strain-level resolution nor genomic context of microorganisms can be provided. This information, however, is crucial to answer fundamental questions about the temporal stability and distinctiveness of natural microbial communities.

Despite the high abundance of Archaea in the global ocean, their metabolism and biogeochemical roles remain largely unresolved. We investigated the population dynamics and metabolic activity of Thaumarchaeota in polar environments, where these microorganisms are particularly abundant and exhibit seasonal growth. Thaumarchaeota were more abundant in deep Arctic and Antarctic waters and grew throughout the winter at surface and deeper Arctic halocline waters. However, in situ single-cell activity measurements revealed a low activity of this group in the uptake of both leucine and bicarbonate (<5% Thaumarchaeota cells active), which is inconsistent with known heterotrophic and autotrophic thaumarchaeal lifestyles. These results suggested the existence of alternative sources of carbon and energy. Our analysis of an environmental metagenome from the Arctic winter revealed that Thaumarchaeota had pathways for ammonia oxidation and, unexpectedly, an abundance of genes involved in urea transport and degradation. Quantitative PCR analysis confirmed that most polar Thaumarchaeota had the potential to oxidize ammonia, and a large fraction of them had urease genes, enabling the use of urea to fuel nitrification. Thaumarchaeota from Arctic deep waters had a higher abundance of urease genes than those near the surface suggesting genetic differences between closely related archaeal populations. In situ measurements of urea uptake and concentration in Arctic waters showed that small-sized prokaryotes incorporated the carbon from urea, and the availability of urea was often higher than that of ammonium. Therefore, the degradation of urea may be a relevant pathway for Thaumarchaeota and other microorganisms exposed to the low-energy conditions of dark polar waters.

Due to the complexity of the protocols and a limited knowledge of the nature of microbial communities, simulating metagenomic sequences plays an important role in testing the performance of existing tools and data analysis methods with metagenomic data. We developed metagenomic read simulators with platform-specific (Sanger, pyrosequencing, Illumina) base-error models, and simulated metagenomes of differing community complexities. We first evaluated the effect of rigorous quality control on Illumina data. Although quality filtering removed a large proportion of the data, it greatly improved the accuracy and contig lengths of resulting assemblies. We then compared the quality-trimmed Illumina assemblies to those from Sanger and pyrosequencing. For the simple community (10 genomes) all sequencing technologies assembled a similar amount and accurately represented the expected functional composition. For the more complex community (100 genomes) Illumina produced the best assemblies and more correctly resembled the expected functional composition. For the most complex community (400 genomes) there was very little assembly of reads from any sequencing technology. However, due to the longer read length the Sanger reads still represented the overall functional composition reasonably well. We further examined the effect of scaffolding of contigs using paired-end Illumina reads. It dramatically increased contig lengths of the simple community and yielded minor improvements to the more complex communities. Although the increase in contig length was accompanied by increased chimericity, it resulted in more complete genes and a better characterization of the functional repertoire. The metagenomic simulators developed for this research are freely available.

Chlorinated solvents are among the most prevalent groundwater contaminants in the industrialized world. Biodegradation with Dehalococcoides-containing mixed cultures is an effective remediation technology. To elucidate transcribed genes in a Dehalococcoides-containing mixed culture, a shotgun metagenome microarray was created and used to investigate gene transcription during vinyl chloride (VC) dechlorination and during starvation (no chlorinated compounds) by a microbial enrichment culture called KB-1. In both treatment conditions, methanol was amended as an electron donor. Subsequently, spots were sequenced that contained the genes most differentially transcribed between the VC-degrading and methanol-only conditions, as well as spots with the highest intensities. Sequencing revealed that during VC degradation Dehalococcoides genes involved in transcription, translation, metabolic energy generation, and amino acid and lipid metabolism and transport were overrepresented in the transcripts compared to the average Dehalococcoides genome. KB-1 rdhA14 (vcrA) was the only reductive dehalogenase homologous (RDH) gene with higher transcript levels during VC degradation, while multiple RDH genes had higher transcript levels in the absence of VC. Numerous hypothetical genes from Dehalococcoides also had higher transcript levels in methanol-only treatments, indicating that many uncharacterized proteins are involved in cell maintenance in the absence of chlorinated substrates. In addition, microarray results prompted biological experiments confirming that electron acceptor limiting conditions activated a Dehalococcoides prophage. Transcripts from Spirochaetes, Chloroflexi, Geobacter, and methanogens demonstrate the importance of non-Dehalococcoides organisms to the culture, and sequencing of identified shotgun clones of interest provided information for follow-on targeted studies.

Despite the current wealth of sequencing data, one-third of all biochemically characterized metabolic enzymes lack a corresponding gene or protein sequence, and as such can be considered orphan enzymes. They represent a major gap between our molecular and biochemical knowledge, and consequently are not amenable to modern systemic analyses. As 555 of these orphan enzymes have metabolic pathway neighbours, we developed a global framework that utilizes the pathway and (meta)genomic neighbour information to assign candidate sequences to orphan enzymes. For 131 orphan enzymes (37% of those for which (meta)genomic neighbours are available), we associate sequences to them using scoring parameters with an estimated accuracy of 70%, implying functional annotation of 16 345 gene sequences in numerous (meta)genomes. As a case in point, two of these candidate sequences were experimentally validated to encode the predicted activity. In addition, we augmented the currently available genome-scale metabolic models with these new sequence-function associations and were able to expand the models by on average 8%, with a considerable change in the flux connectivity patterns and improved essentiality prediction.

Vinyl chloride (VC) is a carcinogenic contaminant commonly found in groundwater. Much research has focused on anaerobic reductive dechlorination of VC, and recently on aerobic VC degradation. In this study, the stable carbon isotope enrichment factor associated with aerobic VC assimilation was determined for Mycobacterium sp. strains JS60, JS61, and JS617 and Nocardioides sp. strain JS614. The enrichment factors ranged from -8.2+/-0.1 to -7.0+/-0.3 % and did not change as a function of biomass concentration. The measured enrichment factors for aerobic VC degradation were smaller than those reported for anaerobic VC degradation. Enrichment factors can also be expressed in terms of kinetic isotope effects (KIEs), 12k/13k, which result from the difference in reaction rates of bonds containing light and heavy isotopes. The KIEs for aerobic VC degradation (1.01+/-0.001) were smaller than those for anaerobic VC degradation (1.03+/-0.007). From the perspective of bond breakage during a chemical reaction, the larger KIE associated with anaerobic VC degradation as compared to aerobic VC degradation agrees with KIE theory. This theory predicts that larger fractionations can be expected in reactions where heavier atoms are involved (i.e., C-Cl bond for anaerobic versus C=C for aerobic) and in reactions involving large changes in vibrational frequencies of the molecule between its ground state and transition state (i.e., C-Cl cleavage versus C=C epoxidation). The significant fractionation observed during aerobic VC degradation suggests that stable carbon isotope measurements may be used as a tool to distinguish between biodegraded and nonbiodegraded VC.

Degenerate primers were used to amplify 14 distinct reductive-dehalogenase-homologous (RDH) genes from the Dehalococcoides-containing mixed culture KB1. Most of the corresponding predicted proteins were highly similar (97 to >99% amino acid identity) to previously reported Dehalococcoides reductive dehalogenases. To examine the differential transcription of these RDH genes, KB1 was split into five subcultures amended with either trichloroethene, cis-1,2-dichloroethene, vinyl chloride, 1,2-dichlorethane, or no chlorinated electron acceptor. Total RNA was extracted following the onset of reductive dechlorination, and RDH transcripts were reverse transcribed and amplified using degenerate primers. The results indicate that the transcription of RDH genes requires the presence of a chlorinated electron acceptor, and for all treatments, multiple RDH genes were simultaneously transcribed, with transcripts of two of the genes being present under all four electron-accepting conditions. Two of the transcribed sequences were highly similar to reported vinyl chloride reductase genes, namely, vcrA from Dehalococcoides sp. strain VS and bvcA from Dehalococcoides sp. strain BAV1. These findings suggest that multiple RDH genes are induced by a single chlorinated substrate and that multiple reductive dehalogenases contribute to chloroethene degradation in KB1.

An extensive microcosm survey of perchlorate-contaminated sites was undertaken to assess the ability of indigenous microorganisms to degrade perchlorate. Samples from 12 contaminated sites and from one pristine location were analysed. Perchlorate was degraded to below detection limit in all electron donor-amended microcosms. Perchlorate-reducing microorganisms (PRMs) were numerous at most of these sites. Sixteen distinct PRMs were isolated that were phylogenetically related to either Dechloromonas in the Beta Proteobacteria (9/16 isolates) or to Azospirillum in the Alpha Proteobacteria (7/16 isolates). The majority of previously isolated PRMs are in the Beta Proteobacteria related to Dechloromonas or Dechlorosoma. This study indicates that PRMs of the genus Azospirillum may be more prevalent at contaminated sites than the current record of isolates suggests. Cell yields, electron donor to perchlorate ratios and maximum specific growth rates were similar among the isolates and similar to the few previously published values. However, the Monod half-saturation constants for perchlorate for the two Azospirillum isolates characterized were lower than those measured for other genera, suggesting that they may be more effective at low concentrations of perchlorate. These results extend the current understanding of PRMs from diverse environments and provide added confidence that microbial perchlorate reduction is ubiquitous, even at highly contaminated sites, and can be harnessed effectively for bioremediation.