Reference-independent comparative metagenomics using cross-assembly: crAss

Solving genomic puzzles: computational methods for metagenomic binning

Metagenomics involves the study of genetic material obtained directly from communities of microorganisms living in natural environments. The field of metagenomics has provided valuable insights into the structure, diversity and ecology of microbial communities. Once an environmental sample is sequenced and processed, metagenomic binning clusters the sequences into bins representing different taxonomic groups such as species, genera, or higher levels. Several computational tools have been developed to automate the process of metagenomic binning. These tools have enabled the recovery of novel draft genomes of microorganisms allowing us to study their behaviors and functions within microbial communities. This review classifies and analyzes different approaches of metagenomic binning and different refinement, visualization, and evaluation techniques used by these methods. Furthermore, the review highlights the current challenges and areas of improvement present within the field of research.

Abolishment of morphology-based taxa and change to binomial species names: 2022 taxonomy update of the ICTV bacterial viruses subcommittee

This article summarises the activities of the Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses for the period of March 2021-March 2022. We provide an overview of the new taxa proposed in 2021, approved by the Executive Committee, and ratified by vote in 2022. Significant changes to the taxonomy of bacterial viruses were introduced: the paraphyletic morphological families Podoviridae, Siphoviridae, and Myoviridae as well as the order Caudovirales were abolished, and a binomial system of nomenclature for species was established. In addition, one order, 22 families, 30 subfamilies, 321 genera, and 862 species were newly created, promoted, or moved.

Unraveling the viral dark matter through viral metagenomics

Viruses are part of the microbiome and have essential roles in immunology, evolution, biogeochemical cycles, health, and disease progression. Viruses influence a wide variety of systems and processes, and the continued discovery of novel viruses is anticipated to reveal new mechanisms influencing the biology of diverse environments. While the identity and roles of viruses continue to be discovered and understood through viral metagenomics, most of the sequences in virome datasets cannot be attributed to known viruses or may be only distantly related to species already described in public sequence databases, at best. Such viruses are known as the viral dark matter. Ongoing discoveries from the viral dark matter have provided insights into novel viruses from a variety of environments, as well as their potential in immunological processes, virus evolution, health, disease, therapeutics, and surveillance. Increased understanding of the viral dark matter will continue with a combination of cultivation, microscopy, sequencing, and bioinformatic efforts, which are discussed in the present review.

PathoLive—Real-Time Pathogen Identification from Metagenomic Illumina Datasets

Over the past years, NGS has become a crucial workhorse for open-view pathogen diagnostics. Yet, long turnaround times result from using massively parallel high-throughput technologies as the analysis can only be performed after sequencing has finished. The interpretation of results can further be challenged by contaminations, clinically irrelevant sequences, and the sheer amount and complexity of the data. We implemented PathoLive, a real-time diagnostics pipeline for the detection of pathogens from clinical samples hours before sequencing has finished. Based on real-time alignment with HiLive2, mappings are scored with respect to common contaminations, low-entropy areas, and sequences of widespread, non-pathogenic organisms. The results are visualized using an interactive taxonomic tree that provides an easily interpretable overview of the relevance of hits. For a human plasma sample that was spiked in vitro with six pathogenic viruses, all agents were clearly detected after only 40 of 200 sequencing cycles. For a real-world sample from Sudan, the results correctly indicated the presence of Crimean-Congo hemorrhagic fever virus. In a second real-world dataset from the 2019 SARS-CoV-2 outbreak in Wuhan, we found the presence of a SARS coronavirus as the most relevant hit without the novel virus reference genome being included in the database. For all samples, clinically irrelevant hits were correctly de-emphasized. Our approach is valuable to obtain fast and accurate NGS-based pathogen identifications and correctly prioritize and visualize them based on their clinical significance: PathoLive is open source and available on GitLab and BioConda.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Metagenome SNP calling via read-colored de Bruijn graphs

MOTIVATION

Metagenomics refers to the study of complex samples containing of genetic contents of multiple individual organisms and, thus, has been used to elucidate the microbiome and resistome of a complex sample. The microbiome refers to all microbial organisms in a sample, and the resistome refers to all of the antimicrobial resistance (AMR) genes in pathogenic and non-pathogenic bacteria. Single-nucleotide polymorphisms (SNPs) can be effectively used to 'fingerprint' specific organisms and genes within the microbiome and resistome and trace their movement across various samples. However, to effectively use these SNPs for this traceability, a scalable and accurate metagenomics SNP caller is needed. Moreover, such an SNP caller should not be reliant on reference genomes since 95% of microbial species is unculturable, making the determination of a reference genome extremely challenging. In this article, we address this need.

RESULTS

We present LueVari, a reference-free SNP caller based on the read-colored de Bruijn graph, an extension of the traditional de Bruijn graph that allows repeated regions longer than the k-mer length and shorter than the read length to be identified unambiguously. LueVari is able to identify SNPs in both AMR genes and chromosomal DNA from shotgun metagenomics data with reliable sensitivity (between 91% and 99%) and precision (between 71% and 99%) as the performance of competing methods varies widely. Furthermore, we show that LueVari constructs sequences containing the variation, which span up to 97.8% of genes in datasets, which can be helpful in detecting distinct AMR genes in large metagenomic datasets.

AVAILABILITY AND IMPLEMENTATION

Code and datasets are publicly available at https://github.com/baharpan/cosmo/tree/LueVari.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Impact of phosphate dosing on the microbial ecology of drinking water distribution systems: Fieldwork studies in chlorinated networks

Phosphate is routinely dosed to ensure regulatory compliance for lead in drinking water distribution systems. Little is known about the impact of the phosphate dose on the microbial ecology in these systems and in particular the endemic biofilms. Disturbance of the biofilms and embedded material in distribution can cause regulatory failures for turbidity and metals. To investigate the impact of phosphate on developing biofilms, pipe wall material from four independent pipe sections was mobilised and collected using two twin-flushing operations a year apart in a chlorinated UK network pre- and post-phosphate dosing. Intensive monitoring was undertaken, including turbidity and water physico-chemistry, traditional microbial culture-based indicators, and microbial community structure via sequencing the 16S rRNA gene for bacteria and the ITS2 gene for fungi. Whole metagenome sequencing was used to study shifts in functional characteristics following the addition of phosphate. As an operational consequence, turbidity responses from the phosphate-enriched water were increased, particularly from cast iron pipes. Differences in the taxonomic composition of both bacteria and fungi were also observed, emphasising a community shift towards microorganisms able to use or metabolise phosphate. Phosphate increased the relative abundance of bacteria such as Pseudomonas, Paenibacillus, Massilia, Acinetobacter and the fungi Cadophora, Rhizophagus and Eupenicillium. Whole metagenome sequencing showed with phosphate a favouring of sequences related to Gram-negative bacterium type cell wall function, virions and thylakoids, but a reduction in the number of sequences associated to vitamin binding, methanogenesis and toxin biosynthesis. With current faecal indicator tests only providing risk detection in bulk water samples, this work improves understanding of how network changes effect microbial ecology and highlights the potential for new approaches to inform future monitoring or control strategies to protect drinking water quality.

Evolutionary Study of the Crassphage Virus at Gene Level

crAss-like viruses are a putative family of bacteriophages recently discovered. The eponym of the clade, crAssphage, is an enteric bacteriophage estimated to be present in at least half of the human population and it constitutes up to 90% of the sequences in some human fecal viral metagenomic datasets. We focused on the evolutionary dynamics of the genes encoded on the crAssphage genome. By investigating the conservation of the genes, a consistent variation in the evolutionary rates across the different functional groups was found. Gene duplications in crAss-like genomes were detected. By exploring the differences among the functional categories of the genes, we confirmed that the genes encoding capsid proteins were the most ubiquitous, despite their overall low sequence conservation. It was possible to identify a core of proteins whose evolutionary trees strongly correlate with each other, suggesting their genetic interaction. This group includes the capsid proteins, which are thus established as extremely suitable for rebuilding the phylogenetic tree of this viral clade. A negative correlation between the ubiquity and the conservation of viral protein sequences was shown. Together, this study provides an in-depth picture of the evolution of different genes in crAss-like viruses.

Shining Light on Human Gut Bacteriophages

The human gut is a complex environment that contains a multitude of microorganisms that are collectively termed the microbiome. Multiple factors have a role to play in driving the composition of human gut bacterial communities either toward homeostasis or the instability that is associated with many disease states. One of the most important forces are likely to be bacteriophages, bacteria-infecting viruses that constitute by far the largest portion of the human gut virome. Despite this, bacteriophages (phages) are the one of the least studied residents of the gut. This is largely due to the challenges associated with studying these difficult to culture entities. Modern high throughput sequencing technologies have played an important role in improving our understanding of the human gut phageome but much of the generated sequencing data remains uncharacterised. Overcoming this requires database-independent bioinformatic pipelines and even those phages that are successfully characterized only provide limited insight into their associated biological properties, and thus most viral sequences have been characterized as “viral dark matter.” Fundamental to understanding the role of phages in shaping the human gut microbiome, and in turn perhaps influencing human health, is how they interact with their bacterial hosts. An essential aspect is the isolation of novel phage-bacteria host pairs by direct isolation through various screening methods, which can transform in silico phages into a biological reality. However, this is also beset with multiple challenges including culturing difficulties and the use of traditional methods, such as plaquing, which may bias which phage-host pairs that can be successfully isolated. Phage-bacteria interactions may be influenced by many aspects of complex human gut biology which can be difficult to reproduce under laboratory conditions. Here we discuss some of the main findings associated with the human gut phageome to date including composition, our understanding of phage-host interactions, particularly the observed persistence of virulent phages and their hosts, as well as factors that may influence these highly intricate relationships. We also discuss current methodologies and bottlenecks hindering progression in this field and identify potential steps that may be useful in overcoming these hurdles.

Assessing human-specific CrAssphage recovery after acidification-filtration concentrating method in environmental water

Pinpointing water pollution sources using host-specific gastrointestinal microbes, known as microbial source tracking (MST), have significant benefits for countries with water quality management issues related to pollution. A recently discovered bacteriophage, crAssphage, shows promise as a human-specific MST marker. However, loss of genetic materials during the recovery and the detection processes could alter the ability to measure virus quantities in a water sample. This study determined the crAssphage recovery efficiencies in water sources, including seawater, freshwater, and influent and effluent from a wastewater treatment plant, by spiking natural crAssphage concentrates prior to DNA extraction and quantitative PCR analysis. The results showed that river and seawater with no or low crAssphage background experienced no recovery loss. Evaluating recovery efficiencies in samples with high crAssphage backgrounds posed a challenge due to the inability to prepare high crAssphage titers. This study highlights the importance of intra-laboratory assessment of recovery efficiency in environmental samples for retrieving absolute crAssphage quantification with correction of bias among water samples and increase in data accuracy. PRACTITIONER POINTS: In laboratory assessment of recovery efficiency is crucial for bias correction and data accuracy for absolute crAssphage quantification in water samples. No loss in crAssphage recovery was observed in river and seawater that contained no or low crAssphage backgrounds. Inability to prepare high crAssphage spike concentrations remains the major limitation for evaluating recovery in samples with high crAssphage backgrounds. The results underline the importance of evaluating method recovery in real environmental samples that reflect actual matrix effect. Absolute crAssphage quantification, as human-specific pollution marker, could be used for prioritizing water quality restoration and area-based management plan.

Gut Bacteriophage: Current Understanding and Challenges

The gut microbiome is widely accepted to have a significant impact on human health yet, despite years of research on this complex ecosystem, the contributions of different forces driving microbial population structure remain to be fully elucidated. The viral component of the human gut microbiome is dominated by bacteriophage, which are known to play crucial roles in shaping microbial composition, driving bacterial diversity, and facilitating horizontal gene transfer. Bacteriophage are also one of the most poorly understood components of the human gut microbiome, with the vast majority of viral sequences sharing little to no homology to reference databases. If we are to understand the dynamics of bacteriophage populations, their interaction with the human microbiome and ultimately their influence on human health, we will depend heavily on sequence based approaches and in silico tools. This is complicated by the fact that, as with any research field in its infancy, methods of analyses vary and this can impede our ability to compare the outputs of different studies. Here, we discuss the major findings to date regarding the human virome and reflect on our current understanding of how gut bacteriophage shape the microbiome. We consider whether or not the virome field is built on unstable foundations and if so, how can we provide a solid basis for future experimentation. The virome is a challenging yet crucial piece of the human microbiome puzzle. In order to develop our understanding, we will discuss the need to underpin future studies with robust research methods and suggest some solutions to existing challenges.

A Phage Protein Aids Bacterial Symbionts in Eukaryote Immune Evasion

Phages are increasingly recognized as important members of host-associated microbiomes, with a vast genomic diversity. The new frontier is to understand how phages may affect higher order processes, such as in the context of host-microbe interactions. Here, we use marine sponges as a model to investigate the interplay between phages, bacterial symbionts, and eukaryotic hosts. Using viral metagenomics, we find that sponges, although massively filtering seawater, harbor species-specific and even individually unique viral signatures that are taxonomically distinct from other environments. We further discover a symbiont phage-encoded ankyrin-domain-containing protein, which is widely spread in phages of many host-associated contexts including human. We confirm in macrophage infection assays that the ankyrin protein (ANKp) modulates the eukaryotic host immune response against bacteria. We predict that the role of ANKp in nature is to facilitate coexistence in the tripartite interplay between phages, symbionts, and sponges and possibly many other host-microbe associations.

On the Role of Bioinformatics and Data Science in Industrial Microbiome Applications

Advances in sequencing and computational biology have drastically increased our capability to explore the taxonomic and functional compositions of microbial communities that play crucial roles in industrial processes. Correspondingly, commercial interest has risen for applications where microbial communities make important contributions. These include food production, probiotics, cosmetics, and enzyme discovery. Other commercial applications include software that takes the user's gut microbiome data as one of its inputs and outputs evidence-based, automated, and personalized diet recommendations for balanced blood sugar levels. These applications pose several bioinformatic and data science challenges that range from requiring strain-level resolution in community profiles to the integration of large datasets for predictive machine learning purposes. In this perspective, we provide our insights on such challenges by touching upon several industrial areas, and briefly discuss advances and future directions of bioinformatics and data science in microbiome research.

Recentrifuge: Robust comparative analysis and contamination removal for metagenomics

Metagenomic sequencing is becoming widespread in biomedical and environmental research, and the pace is increasing even more thanks to nanopore sequencing. With a rising number of samples and data per sample, the challenge of efficiently comparing results within a specimen and between specimens arises. Reagents, laboratory, and host related contaminants complicate such analysis. Contamination is particularly critical in low microbial biomass body sites and environments, where it can comprise most of a sample if not all. Recentrifuge implements a robust method for the removal of negative-control and crossover taxa from the rest of samples. With Recentrifuge, researchers can analyze results from taxonomic classifiers using interactive charts with emphasis on the confidence level of the classifications. In addition to contamination-subtracted samples, Recentrifuge provides shared and exclusive taxa per sample, thus enabling robust contamination removal and comparative analysis in environmental and clinical metagenomics. Regarding the first area, Recentrifuge's novel approach has already demonstrated its benefits showing that microbiomes of Arctic and Antarctic solar panels display similar taxonomic profiles. In the clinical field, to confirm Recentrifuge's ability to analyze complex metagenomes, we challenged it with data coming from a metagenomic investigation of RNA in plasma that suffered from critical contamination to the point of preventing any positive conclusion. Recentrifuge provided results that yielded new biological insight into the problem, supporting the growing evidence of a blood microbiota even in healthy individuals, mostly translocated from the gut, the oral cavity, and the genitourinary tract. We also developed a synthetic dataset carefully designed to rate the robust contamination removal algorithm, which demonstrated a significant improvement in specificity while retaining a high sensitivity even in the presence of cross-contaminants. Recentrifuge's official website is www.recentrifuge.org. The data and source code are anonymously and freely available on GitHub and PyPI. The computing code is licensed under the AGPLv3. The Recentrifuge Wiki is the most extensive and continually-updated source of documentation for Recentrifuge, covering installation, use cases, testing, and other useful topics.

Taking Advantage of the Genomics Revolution for Monitoring and Conservation of Chondrichthyan Populations

Chondrichthyes (sharks, rays, skates and chimaeras) are among the oldest extant predators and are vital to top-down regulation of oceanic ecosystems. They are an ecologically diverse group occupying a wide range of habitats and are thus, exploited by coastal, pelagic and deep-water fishing industries. Chondrichthyes are among the most data deficient vertebrate species groups making design and implementation of regulatory and conservation measures challenging. High-throughput sequencing technologies have significantly propelled ecological investigations and understanding of marine and terrestrial species’ populations, but there remains a paucity of NGS based research on chondrichthyan populations. We present a brief review of current methods to access genomic and metagenomic data from Chondrichthyes and discuss applications of these datasets to increase our understanding of chondrichthyan taxonomy, evolution, ecology and population structures. Last, we consider opportunities and challenges offered by genomic studies for conservation and management of chondrichthyan populations.

The gill-associated microbiome is the main source of wood plant polysaccharide hydrolases and secondary metabolite gene clusters in the mangrove shipworm Neoteredo reynei

Teredinidae are a family of highly adapted wood-feeding and wood-boring bivalves, commonly known as shipworms, whose evolution is linked to the acquisition of cellulolytic gammaproteobacterial symbionts harbored in bacteriocytes within the gills. In the present work we applied metagenomics to characterize microbiomes of the gills and digestive tract of Neoteredo reynei, a mangrove-adapted shipworm species found over a large range of the Brazilian coast. Comparative metagenomics grouped the gill symbiont community of different N. reynei specimens, indicating closely related bacterial types are shared. Similarly, the intestine and digestive gland communities were related, yet were more diverse than and showed no overlap with the gill community. Annotation of assembled metagenomic contigs revealed that the gill symbiotic community of N. reynei encodes a plethora of plant cell wall polysaccharides degrading glycoside hydrolase encoding genes, and Biosynthetic Gene Clusters (BGCs). In contrast, the digestive tract microbiomes seem to play little role in wood digestion and secondary metabolites biosynthesis. Metagenome binning recovered the nearly complete genome sequences of two symbiotic Teredinibacter strains from the gills, a representative of Teredinibacter turnerae “clade I” strain, and a yet to be cultivated Teredinibacter sp. type. These Teredinibacter genomes, as well as un-binned gill-derived gammaproteobacteria contigs, also include an endo-β-1,4-xylanase/acetylxylan esterase multi-catalytic carbohydrate-active enzyme, and a trans-acyltransferase polyketide synthase (trans-AT PKS) gene cluster with the gene cassette for generating β-branching on complex polyketides. Finally, we use multivariate analyses to show that the secondary metabolome from the genomes of Teredinibacter representatives, including genomes binned from N. reynei gills’ metagenomes presented herein, stands out within the Cellvibrionaceae family by size, and enrichments for polyketide, nonribosomal peptide and hybrid BGCs. Results presented here add to the growing characterization of shipworm symbiotic microbiomes and indicate that the N. reynei gill gammaproteobacterial community is a prolific source of biotechnologically relevant enzymes for wood-digestion and bioactive compounds production.

Metagenomics of Coral Reefs Under Phase Shift and High Hydrodynamics

Local and global stressors have affected coral reef ecosystems worldwide. Switches from coral to algal dominance states and microbialization are the major processes underlying the global decline of coral reefs. However, most of the knowledge concerning microbialization has not considered physical disturbances (e.g., typhoons, waves, and currents). Southern Japan reef systems have developed under extreme physical disturbances. Here, we present analyses of a three-year investigation on the coral reefs of Ishigaki Island that comprised benthic and fish surveys, water quality analyses, metagenomics and microbial abundance data. At the four studied sites, inorganic nutrient concentrations were high and exceeded eutrophication thresholds. The dissolved organic carbon (DOC) concentration (up to 233.3 μM) and microbial abundance (up to 2.5 × 105 cell/mL) values were relatively high. The highest vibrio counts coincided with the highest turf cover (∼55-85%) and the lowest coral cover (∼4.4-10.2%) and fish biomass (0.06 individuals/m2). Microbiome compositions were similar among all sites and were dominated by heterotrophs. Our data suggest that a synergic effect among several regional stressors are driving coral decline. In a high hydrodynamics reef environment, high algal/turf cover, stimulated by eutrophication and low fish abundance due to overfishing, promote microbialization. Together with crown-of-thorns starfish (COTS) outbreaks and possible of climate changes impacts, theses coral reefs are likely to collapse.

Novel phages of healthy skin metaviromes from South Africa

Recent skin metagenomic studies have investigated the harbored viral diversity and its possible influence on healthy skin microbial populations, and tried to establish global patterns of skin-phage evolution. However, the detail associated with the phages that potentially play a role in skin health has not been investigated. While skin metagenome and -metavirome studies have indicated that the skin virome is highly site specific and shows marked interpersonal variation, they have not assessed the presence/absence of individual phages. Here, we took a semi-culture independent approach (metaviromic) to better understand the composition of phage communities on skin from South African study participants. Our data set adds over 130 new phage species of the skin to existing databases. We demonstrated that identical phages were present on different individuals and in different body sites, and we conducted a detailed analysis of the structural organization of these phages. We further found that a bacteriophage related to the Staphylococcus capitis phage Stb20 may be a common skin commensal virus potentially regulating its host and its activities on the skin.

The Human Virome Protein Cluster Database (HVPC): A Human Viral Metagenomic Database for Diversity and Function Annotation

Human virome, including those of bacteria (bacteriophages) have received an increasing attention recently, owing to the rapid developments in human microbiome research and the awareness of the far-reaching influence of microbiomes on health and disease. Nevertheless, human viromes are still underrepresented in literature making viruses a virtually untapped resource of diversity, functional and physiological information. Here we present the human virome protein cluster database as an effort to improve functional annotation and characterization of human viromes. The database was built out of hundreds of virome datasets from six different body sites. We also show the utility of this database through its use for the characterization of three bronchoalveolar lavage (BAL) viromes from one healthy control in addition to one moderate and one severe chronic obstructive pulmonary disease (COPD) patients. The use of the database allowed for a better functional annotation, which were otherwise poorly characterized when limited to annotation using sequences from full-length viral genomes. In addition, our BAL samples gave a first insight into viral communities of COPD patients and confirm a state of dysbiosis for viruses that increases with disease progression. Moreover, they shed light on the potential role of phages in the horizontal gene transfer of bacterial virulence factors, a phenomenon that highlights a possible contribution of phages to etiopathology.

Optimizing and evaluating the reconstruction of Metagenome-assembled microbial genomes

Background

Microbiome/host interactions describe characteristics that affect the host's health. Shotgun metagenomics includes sequencing a random subset of the microbiome to analyze its taxonomic and metabolic potential. Reconstruction of DNA fragments into genomes from metagenomes (called metagenome-assembled genomes) assigns unknown fragments to taxa/function and facilitates discovery of novel organisms. Genome reconstruction incorporates sequence assembly and sorting of assembled sequences into bins, characteristic of a genome. However, the microbial community composition, including taxonomic and phylogenetic diversity may influence genome reconstruction. We determine the optimal reconstruction method for four microbiome projects that had variable sequencing platforms (IonTorrent and Illumina), diversity (high or low), and environment (coral reefs and kelp forests), using a set of parameters to select for optimal assembly and binning tools.

Methods

We tested the effects of the assembly and binning processes on population genome reconstruction using 105 marine metagenomes from 4 projects. Reconstructed genomes were obtained from each project using 3 assemblers (IDBA, MetaVelvet, and SPAdes) and 2 binning tools (GroopM and MetaBat). We assessed the efficiency of assemblers using statistics that including contig continuity and contig chimerism and the effectiveness of binning tools using genome completeness and taxonomic identification.

Results

We concluded that SPAdes, assembled more contigs (143,718 ± 124 contigs) of longer length (N50 = 1632 ± 108 bp), and incorporated the most sequences (sequences-assembled = 19.65%). The microbial richness and evenness were maintained across the assembly, suggesting low contig chimeras. SPAdes assembly was responsive to the biological and technological variations within the project, compared with other assemblers. Among binning tools, we conclude that MetaBat produced bins with less variation in GC content (average standard deviation: 1.49), low species richness (4.91 ± 0.66), and higher genome completeness (40.92 ± 1.75) across all projects. MetaBat extracted 115 bins from the 4 projects of which 66 bins were identified as reconstructed metagenome-assembled genomes with sequences belonging to a specific genus. We identified 13 novel genomes, some of which were 100% complete, but show low similarity to genomes within databases.

Conclusions

In conclusion, we present a set of biologically relevant parameters for evaluation to select for optimal assembly and binning tools. For the tools we tested, SPAdes assembler and MetaBat binning tools reconstructed quality metagenome-assembled genomes for the four projects. We also conclude that metagenomes from microbial communities that have high coverage of phylogenetically distinct, and low taxonomic diversity results in highest quality metagenome-assembled genomes.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4294-1) contains supplementary material, which is available to authorized users.

The skin microbiome of the common thresher shark (Alopias vulpinus) has low taxonomic and gene function β-diversity

The health of sharks, like all organisms, is linked to their microbiome. At the skin interface, sharks have dermal denticles that protrude above the mucus, which may affect the types of microbes that occur here. We characterized the microbiome from the skin of the common thresher shark (Alopias vulpinus) to investigate the structure and composition of the skin microbiome. On average 618 812 (80.9% ± S.D. 0.44%) reads per metagenomic library contained open reading frames; of those, between 7.6% and 12.8% matched known protein sequences. Genera distinguishing the A. vulpinus microbiome from the water column included, Pseudoalteromonas (12.8% ± 4.7 of sequences), Erythrobacter (5. 3% ± 0.5) and Idiomarina (4.2% ± 1.2) and distinguishing gene pathways included, cobalt, zinc and cadmium resistance (2.2% ± 0.1); iron acquisition (1.2% ± 0.1) and ton/tol transport (1.3% ± 0.08). Taxonomic community overlap (100 - dissimilarity index) was greater in the skin microbiome (77.6), relative to the water column microbiome (70.6) and a reference host-associated microbiome (algae: 71.5). We conclude the A. vulpinus skin microbiome is influenced by filtering processes, including biochemical and biophysical components of the shark skin and result in a structured microbiome.

Environmental drivers of viral community composition in Antarctic soils identified by viromics

BACKGROUND

The Antarctic continent is considered the coldest and driest place on earth with simple ecosystems, devoid of higher plants. Soils in the ice-free regions of Antarctica are known to harbor a wide range of microorganisms from primary producers to grazers, yet their ecology and particularly the role of viruses is poorly understood. In this study, we examined the virus community structures of 14 soil samples from the Mackay Glacier region.

METHODS

Viral communities were extracted from soil and the dsDNA was extracted, amplified using single-primer amplification, and sequenced using the Ion Torrent Proton platform. Metadata on soil physico-chemistry was collected from all sites. Both read and contig datasets were analyzed with reference-independent and reference-dependent methods to assess viral community structures and the influence of environmental parameters on their distribution.

RESULTS

We observed a high heterogeneity in virus signatures, independent of geographical proximity. Tailed bacteriophages were dominant in all samples, but the incidences of the affiliated families Siphoviridae and Myoviridae were inversely correlated, suggesting direct competition for hosts. Viruses of the families Phycodnaviridae and Mimiviridae were present at significant levels in high-diversity soil samples and were found to co-occur, implying little competition between them. Combinations of soil factors, including pH, calcium content, and site altitude, were found to be the main drivers of viral community structure.

CONCLUSIONS

The pattern of viral community structure with higher levels of diversity at lower altitude and pH, and co-occurring viral families, suggests that these cold desert soil viruses interact with each other, the host, and the environment in an intricate manner, playing a potentially crucial role in maintaining host diversity and functioning of the microbial ecosystem in the extreme environments of Antarctic soil.

Metagenomic Approaches to Assess Bacteriophages in Various Environmental Niches

Bacteriophages are ubiquitous and numerous parasites of bacteria and play a critical evolutionary role in virtually every ecosystem, yet our understanding of the extent of the diversity and role of phages remains inadequate for many ecological niches, particularly in cases in which the host is unculturable. During the past 15 years, the emergence of the field of viral metagenomics has drastically enhanced our ability to analyse the so-called viral ‘dark matter’ of the biosphere. Here, we review the evolution of viral metagenomic methodologies, as well as providing an overview of some of the most significant applications and findings in this field of research.

Pharmacometabolomics Informs Viromics toward Precision Medicine

Nowadays, we are experiencing the big data era with the emerging challenge of single data interpretation. Although the advent of high-throughput technologies as well as chemo- and bio-informatics tools presents pan-omics data as the way forward to precision medicine, personalized health care and tailored-made therapeutics can be only envisaged when interindividual variability in response to/toxicity of xenobiotics can be interpreted and thus, predicted. We know that such variability is the net outcome of genetics (host and microbiota) and environmental factors (diet, lifestyle, polypharmacy, and microbiota) and for this, tremendous efforts have been made to clarify key-molecules from correlation to causality to clinical significance. Herein, we focus on the host–microbiome interplay and its direct and indirect impact on efficacy and toxicity of xenobiotics and we inevitably wonder about the role of viruses, as the least acknowledged ones. We present the emerging discipline of pharmacometabolomics-informed viromics, in which pre-dose metabotypes can assist modeling and prediction of interindividual response to/toxicity of xenobiotics. Such features, either alone or in combination with host genetics, can power biomarker discovery so long as the features are variable among patients, stable enough to be of predictive value, and better than pre-existing tools for predicting therapeutic efficacy/toxicity.

Freshwater Metaviromics and Bacteriophages: A Current Assessment of the State of the Art in Relation to Bioinformatic Challenges

Advances in bioinformatics and sequencing technologies have allowed for the analysis of complex microbial communities at an unprecedented rate. While much focus is often placed on the cellular members of these communities, viruses play a pivotal role, particularly bacteria-infecting viruses (bacteriophages); phages mediate global biogeochemical processes and drive microbial evolution through bacterial grazing and horizontal gene transfer. Despite their importance and ubiquity in nature, very little is known about the diversity and structure of viral communities. Though the need for culture-based methods for viral identification has been somewhat circumvented through metagenomic techniques, the analysis of metaviromic data is marred with many unique issues. In this review, we examine the current bioinformatic approaches for metavirome analyses and the inherent challenges facing the field as illustrated by the ongoing efforts in the exploration of freshwater phage populations.

MetaFast: fast reference-free graph-based comparison of shotgun metagenomic data

MOTIVATION

High-throughput metagenomic sequencing has revolutionized our view on the structure and metabolic potential of microbial communities. However, analysis of metagenomic composition is often complicated by the high complexity of the community and the lack of related reference genomic sequences. As a start point for comparative metagenomic analysis, the researchers require efficient means for assessing pairwise similarity of the metagenomes (beta-diversity). A number of approaches were used to address this task, however, most of them have inherent disadvantages that limit their scope of applicability. For instance, the reference-based methods poorly perform on metagenomes from previously unstudied niches, while composition-based methods appear to be too abstract for straightforward interpretation and do not allow to identify the differentially abundant features.

RESULTS

We developed MetaFast, an approach that allows to represent a shotgun metagenome from an arbitrary environment as a modified de Bruijn graph consisting of simplified components. For multiple metagenomes, the resulting representation is used to obtain a pairwise similarity matrix. The dimensional structure of the metagenomic components preserved in our algorithm reflects the inherent subspecies-level diversity of microbiota. The method is computationally efficient and especially promising for an analysis of metagenomes from novel environmental niches.

AVAILABILITY AND IMPLEMENTATION

Source code and binaries are freely available for download at https://github.com/ctlab/metafast The code is written in Java and is platform independent (tested on Linux and Windows x86_64).

CONTACT

ulyantsev@rain.ifmo.ru

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Metagenomic analysis of lacustrine viral diversity along a latitudinal transect of the Antarctic Peninsula

Environmental viruses constitute the most abundant biological entities on earth, and harbor an enormous genetic diversity. While their strong influence on the ecosystem is widely acknowledged, current knowledge about their diversity and distribution remains limited. Here we present the metagenomic study of viral communities from freshwater bodies located along a transect of the Antarctic Peninsula. These ecosystems were chosen on the basis of environmental and biogeographical variation. The results obtained indicate that the virus assemblages were diverse, and that the larger fraction represented viruses with no close relatives in the databases. Comparisons to existing metaviromes showed that the communities studied were dissimilar to other freshwater viromes including those from the Arctic. Finally, we observed no indication of there being a reduction in either viral richness or diversity estimates with increasing latitude along the studied transect, further adding to the controversy regarding the possible existence of latitudinal gradients of diversity in the microbial world.

Lytic to temperate switching of viral communities

Microbial viruses can control host abundances via density-dependent lytic predator-prey dynamics. Less clear is how temperate viruses, which coexist and replicate with their host, influence microbial communities. Here we show that virus-like particles are relatively less abundant at high host densities. This suggests suppressed lysis where established models predict lytic dynamics are favoured. Meta-analysis of published viral and microbial densities showed that this trend was widespread in diverse ecosystems ranging from soil to freshwater to human lungs. Experimental manipulations showed viral densities more consistent with temperate than lytic life cycles at increasing microbial abundance. An analysis of 24 coral reef viromes showed a relative increase in the abundance of hallmark genes encoded by temperate viruses with increased microbial abundance. Based on these four lines of evidence, we propose the Piggyback-the-Winner model wherein temperate dynamics become increasingly important in ecosystems with high microbial densities; thus 'more microbes, fewer viruses'.

Analysis of dsDNA and RNA viromes in methanogenic digesters reveals novel viral genetic diversity

Although viruses are not the key players of the anaerobic digestion process, they may affect the dynamics of bacterial and archaeal populations involved in biogas production. Until now viruses have received very little attention in this specific habitat; therefore, as a first step towards their characterization, we optimized a virus filtration protocol from anaerobic sludge. Afterwards, to assess dsDNA and RNA viral diversity in sludge samples from nine different reactors fed either with waste water, agricultural residues or solid municipal waste plus agro‐food residues, we performed metagenomic analyses. As a result we showed that, while the dsDNA viromes (21 assigned families in total) were dominated by dsDNA phages of the order Caudovirales, RNA viruses (14 assigned families in total) were less diverse and were for the main part plant‐infecting viruses. Interestingly, less than 2% of annotated contigs were assigned as putative human and animal pathogens. Our study greatly extends the existing view of viral genetic diversity in methanogenic reactors and shows that these viral assemblages are distinct not only among the reactor types but also from nearly 30 other environments already studied, including the human gut, fermented food, deep sea sediments and other aquatic habitats.

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

Metagenomic analysis of planktonic microbial consortia from a non-tidal urban-impacted segment of James River

Knowledge of the diversity and ecological function of the microbial consortia of James River in Virginia, USA, is essential to developing a more complete understanding of the ecology of this model river system. Metagenomic analysis of James River's planktonic microbial community was performed for the first time using an unamplified genomic library and a 16S rDNA amplicon library prepared and sequenced by Ion PGM and MiSeq, respectively. From the 0.46-Gb WGS library (GenBank:SRR1146621; MG-RAST:4532156.3), 4 × 10(6) reads revealed >3 × 10(6) genes, 240 families of prokaryotes, and 155 families of eukaryotes. From the 0.68-Gb 16S library (GenBank:SRR2124995; MG-RAST:4631271.3; EMB:2184), 4 × 10(6) reads revealed 259 families of eubacteria. Results of the WGS and 16S analyses were highly consistent and indicated that more than half of the bacterial sequences were Proteobacteria, predominantly Comamonadaceae. The most numerous genera in this group were Acidovorax (including iron oxidizers, nitrotolulene degraders, and plant pathogens), which accounted for 10 % of assigned bacterial reads. Polaromonas were another 6 % of all bacterial reads, with many assignments to groups capable of degrading polycyclic aromatic hydrocarbons. Albidiferax (iron reducers) and Variovorax (biodegraders of a variety of natural biogenic compounds as well as anthropogenic contaminants such as polycyclic aromatic hydrocarbons and endocrine disruptors) each accounted for an additional 3 % of bacterial reads. Comparison of these data to other publically-available aquatic metagenomes revealed that this stretch of James River is highly similar to the upper Mississippi River, and that these river systems are more similar to aquaculture and sludge ecosystems than they are to lakes or to a pristine section of the upper Amazon River. Taken together, these analyses exposed previously unknown aspects of microbial biodiversity, documented the ecological responses of microbes to urban effects, and revealed the noteworthy presence of 22 human-pathogenic bacterial genera (e.g., Enterobacteriaceae, pathogenic Pseudomonadaceae, and 'Vibrionales') and 6 pathogenic eukaryotic genera (e.g., Trypanosomatidae and Vahlkampfiidae). This information about pathogen diversity may be used to promote human epidemiological studies, enhance existing water quality monitoring efforts, and increase awareness of the possible health risks associated with recreational use of James River.

The human gut virome: a multifaceted majority

Here, we outline our current understanding of the human gut virome, in particular the phage component of this ecosystem, highlighting progress, and challenges in viral discovery in this arena. We reveal how developments in high-throughput sequencing technologies and associated data analysis methodologies are helping to illuminate this abundant 'biological dark matter.' Current evidence suggests that the human gut virome is a highly individual but temporally stable collective, dominated by phages exhibiting a temperate lifestyle. This viral community also appears to encode a surprisingly rich functional repertoire that confers a range of attributes to their bacterial hosts, ranging from bacterial virulence and pathogenesis to maintaining host-microbiome stability and community resilience. Despite the significant advances in our understanding of the gut virome in recent years, it is clear that we remain in a period of discovery and revelation, as new methods and technologies begin to provide deeper understanding of the inherent ecological characteristics of this viral ecosystem. As our understanding increases, the nature of the multi-partite interactions occurring between host and microbiome will become clearer, helping us to more rationally define the concepts and principles that will underpin approaches to using human gut virome components for medical or biotechnological applications.

From cultured to uncultured genome sequences: metagenomics and modeling microbial ecosystems

Microorganisms and the viruses that infect them are the most numerous biological entities on Earth and enclose its greatest biodiversity and genetic reservoir. With strength in their numbers, these microscopic organisms are major players in the cycles of energy and matter that sustain all life. Scientists have only scratched the surface of this vast microbial world through culture-dependent methods. Recent developments in generating metagenomes, large random samples of nucleic acid sequences isolated directly from the environment, are providing comprehensive portraits of the composition, structure, and functioning of microbial communities. Moreover, advances in metagenomic analysis have created the possibility of obtaining complete or nearly complete genome sequences from uncultured microorganisms, providing important means to study their biology, ecology, and evolution. Here we review some of the recent developments in the field of metagenomics, focusing on the discovery of genetic novelty and on methods for obtaining uncultured genome sequences, including through the recycling of previously published datasets. Moreover we discuss how metagenomics has become a core scientific tool to characterize eco-evolutionary patterns of microbial ecosystems, thus allowing us to simultaneously discover new microbes and study their natural communities. We conclude by discussing general guidelines and challenges for modeling the interactions between uncultured microorganisms and viruses based on the information contained in their genome sequences. These models will significantly advance our understanding of the functioning of microbial ecosystems and the roles of microbes in the environment.

Combining metagenomics, metatranscriptomics and viromics to explore novel microbial interactions: towards a systems-level understanding of human microbiome

The advances in experimental methods and the development of high performance bioinformatic tools have substantially improved our understanding of microbial communities associated with human niches. Many studies have documented that changes in microbial abundance and composition of the human microbiome is associated with human health and diseased state. The majority of research on human microbiome is typically focused in the analysis of one level of biological information, i.e., metagenomics or metatranscriptomics. In this review, we describe some of the different experimental and bioinformatic strategies applied to analyze the 16S rRNA gene profiling and shotgun sequencing data of the human microbiome. We also discuss how some of the recent insights in the combination of metagenomics, metatranscriptomics and viromics can provide more detailed description on the interactions between microorganisms and viruses in oral and gut microbiomes. Recent studies on viromics have begun to gain importance due to the potential involvement of viruses in microbial dysbiosis. In addition, metatranscriptomic combined with metagenomic analysis have shown that a substantial fraction of microbial transcripts can be differentially regulated relative to their microbial genomic abundances. Thus, understanding the molecular interactions in the microbiome using the combination of metagenomics, metatranscriptomics and viromics is one of the main challenges towards a system level understanding of human microbiome.

Biodiversity and distribution of polar freshwater DNA viruses

Viruses constitute the most abundant biological entities and a large reservoir of genetic diversity on Earth. Despite the recent surge in their study, our knowledge on their actual biodiversity and distribution remains sparse. We report the first metagenomic analysis of Arctic freshwater viral DNA communities and a comparative analysis with other freshwater environments. Arctic viromes are dominated by unknown and single-stranded DNA viruses with no close relatives in the database. These unique viral DNA communities mostly relate to each other and present some minor genetic overlap with other environments studied, including an Arctic Ocean virome. Despite common environmental conditions in polar ecosystems, the Arctic and Antarctic DNA viromes differ at the fine-grain genetic level while sharing a similar taxonomic composition. The study uncovers some viral lineages with a bipolar distribution, suggesting a global dispersal capacity for viruses, and seemingly indicates that viruses do not follow the latitudinal diversity gradient known for macroorganisms. Our study sheds light into the global biogeography and connectivity of viral communities.

The effects of variable sample biomass on comparative metagenomics

Longitudinal studies that integrate samples with variable biomass are essential to understand microbial community dynamics across space or time. Shotgun metagenomics is widely used to investigate these communities at the functional level, but little is known about the effects of combining low and high biomass samples on downstream analysis. We investigated the interacting effects of DNA input and library amplification by polymerase chain reaction on comparative metagenomic analysis using dilutions of a single complex template from an Arabidopsis thaliana-associated microbial community. We modified the Illumina Nextera kit to generate high-quality large-insert (680 bp) paired-end libraries using a range of 50 pg to 50 ng of input DNA. Using assembly-based metagenomic analysis, we demonstrate that DNA input level has a significant impact on community structure due to overrepresentation of low-GC genomic regions following library amplification. In our system, these differences were largely superseded by variations between biological replicates, but our results advocate verifying the influence of library amplification on a case-by-case basis. Overall, this study provides recommendations for quality filtering and de-replication prior to analysis, as well as a practical framework to address the issue of low biomass or biomass heterogeneity in longitudinal metagenomic surveys.

Purifying the impure: sequencing metagenomes and metatranscriptomes from complex animal-associated samples

The accessibility of high-throughput sequencing has revolutionized many fields of biology. In order to better understand host-associated viral and microbial communities, a comprehensive workflow for DNA and RNA extraction was developed. The workflow concurrently generates viral and microbial metagenomes, as well as metatranscriptomes, from a single sample for next-generation sequencing. The coupling of these approaches provides an overview of both the taxonomical characteristics and the community encoded functions. The presented methods use Cystic Fibrosis (CF) sputum, a problematic sample type, because it is exceptionally viscous and contains high amount of mucins, free neutrophil DNA, and other unknown contaminants. The protocols described here target these problems and successfully recover viral and microbial DNA with minimal human DNA contamination. To complement the metagenomics studies, a metatranscriptomics protocol was optimized to recover both microbial and host mRNA that contains relatively few ribosomal RNA (rRNA) sequences. An overview of the data characteristics is presented to serve as a reference for assessing the success of the methods. Additional CF sputum samples were also collected to (i) evaluate the consistency of the microbiome profiles across seven consecutive days within a single patient, and (ii) compare the consistency of metagenomic approach to a 16S ribosomal RNA gene-based sequencing. The results showed that daily fluctuation of microbial profiles without antibiotic perturbation was minimal and the taxonomy profiles of the common CF-associated bacteria were highly similar between the 16S rDNA libraries and metagenomes generated from the hypotonic lysis (HL)-derived DNA. However, the differences between 16S rDNA taxonomical profiles generated from total DNA and HL-derived DNA suggest that hypotonic lysis and the washing steps benefit in not only removing the human-derived DNA, but also microbial-derived extracellular DNA that may misrepresent the actual microbial profiles.

Comparative metagenomic analysis of coral microbial communities using a reference-independent approach

By comparing the SEED and Pfam functional profiles of metagenomes of two Brazilian coral species with 29 datasets that are publicly available, we were able to identify some functions, such as protein secretion systems, that are overrepresented in the metagenomes of corals and may play a role in the establishment and maintenance of bacteria-coral associations. However, only a small percentage of the reads of these metagenomes could be annotated by these reference databases, which may lead to a strong bias in the comparative studies. For this reason, we have searched for identical sequences (99% of nucleotide identity) among these metagenomes in order to perform a reference-independent comparative analysis, and we were able to identify groups of microbial communities that may be under similar selective pressures. The identification of sequences shared among the metagenomes was found to be even better for the identification of groups of communities with similar niche requirements than the traditional analysis of functional profiles. This approach is not only helpful for the investigation of similarities between microbial communities with high proportion of unknown reads, but also enables an indirect overview of gene exchange between communities.

Species-specific viromes in the ancestral holobiont Hydra

Recent evidence showing host specificity of colonizing bacteria supports the view that multicellular organisms are holobionts comprised of the macroscopic host in synergistic interdependence with a heterogeneous and host-specific microbial community. Whereas host-bacteria interactions have been extensively investigated, comparatively little is known about host-virus interactions and viral contribution to the holobiont. We sought to determine the viral communities associating with different Hydra species, whether these viral communities were altered with environmental stress, and whether these viruses affect the Hydra-associated holobiont. Here we show that each species of Hydra harbors a diverse host-associated virome. Primary viral families associated with Hydra are Myoviridae, Siphoviridae, Inoviridae, and Herpesviridae. Most Hydra-associated viruses are bacteriophages, a reflection of their involvement in the holobiont. Changes in environmental conditions alter the associated virome, increase viral diversity, and affect the metabolism of the holobiont. The specificity and dynamics of the virome point to potential viral involvement in regulating microbial associations in the Hydra holobiont. While viruses are generally regarded as pathogenic agents, our study suggests an evolutionary conserved ability of viruses to function as holobiont regulators and, therefore, constitutes an emerging paradigm shift in host-microbe interactions.

A highly abundant bacteriophage discovered in the unknown sequences of human faecal metagenomes

Metagenomics, or sequencing of the genetic material from a complete microbial community, is a promising tool to discover novel microbes and viruses. Viral metagenomes typically contain many unknown sequences. Here we describe the discovery of a previously unidentified bacteriophage present in the majority of published human faecal metagenomes, which we refer to as crAssphage. Its ~97 kbp genome is six times more abundant in publicly available metagenomes than all other known phages together; it comprises up to 90% and 22% of all reads in virus-like particle (VLP)-derived metagenomes and total community metagenomes, respectively; and it totals 1.68% of all human faecal metagenomic sequencing reads in the public databases. The majority of crAssphage-encoded proteins match no known sequences in the database, which is why it was not detected before. Using a new co-occurrence profiling approach, we predict a Bacteroides host for this phage, consistent with Bacteroides-related protein homologues and a unique carbohydrate-binding domain encoded in the phage genome.

Metagenomic studies of microbial communities often report DNA sequences from unidentified viruses. Here, Dutilh et al. analyse metagenomic data to reveal the complete genome of an abundant, ubiquitous virus from human faeces, and predict that the virus infects bacteria of the Bacteroides group.

Exploring variation-aware contig graphs for (comparative) metagenomics using MaryGold

MOTIVATION

Although many tools are available to study variation and its impact in single genomes, there is a lack of algorithms for finding such variation in metagenomes. This hampers the interpretation of metagenomics sequencing datasets, which are increasingly acquired in research on the (human) microbiome, in environmental studies and in the study of processes in the production of foods and beverages. Existing algorithms often depend on the use of reference genomes, which pose a problem when a metagenome of a priori unknown strain composition is studied. In this article, we develop a method to perform reference-free detection and visual exploration of genomic variation, both within a single metagenome and between metagenomes.

RESULTS

We present the MaryGold algorithm and its implementation, which efficiently detects bubble structures in contig graphs using graph decomposition. These bubbles represent variable genomic regions in closely related strains in metagenomic samples. The variation found is presented in a condensed Circos-based visualization, which allows for easy exploration and interpretation of the found variation. We validated the algorithm on two simulated datasets containing three respectively seven Escherichia coli genomes and showed that finding allelic variation in these genomes improves assemblies. Additionally, we applied MaryGold to publicly available real metagenomic datasets, enabling us to find within-sample genomic variation in the metagenomes of a kimchi fermentation process, the microbiome of a premature infant and in microbial communities living on acid mine drainage. Moreover, we used MaryGold for between-sample variation detection and exploration by comparing sequencing data sampled at different time points for both of these datasets.

AVAILABILITY

MaryGold has been written in C++ and Python and can be downloaded from http://bioinformatics.tudelft.nl/software

Going viral: next-generation sequencing applied to phage populations in the human gut

Over the past decade, researchers have begun to characterize viral diversity using metagenomic methods. These studies have shown that viruses, the majority of which infect bacteria, are probably the most genetically diverse components of the biosphere. Here, we briefly review the incipient rise of a phage biology renaissance, which has been catalysed by advances in next-generation sequencing. We explore how work characterizing phage diversity and lifestyles in the human gut is changing our view of ourselves as supra-organisms. Finally, we discuss how a renewed appreciation of phage dynamics may yield new applications for phage therapies designed to manipulate the structure and functions of our gut microbiomes.