While large-scale efforts have rapidly advanced the understanding and practical impact of human genomic variation, the latter is largely unexplored in the human microbiome. We therefore developed a framework for metagenomic variation analysis and applied it to 252 fecal 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 indels, 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 implies 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.

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.

ABSTRACT. Cheese represents vast share in consumption of dairy products as one third of collected milk is used by cheese-makers. This is one of the rare dairy products, which per capita consumption increases even in developed countries. The paper investigates changes in per capita consumption and links them with projections for Polish cheese-making industry and consumers.

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.

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