Genivaldo G. Z. Silva

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.

Critical assessment of metagenome interpretation − a benchmark of computational metagenomics software

Methods for assembly, taxonomic profiling and binning are key to interpreting metagenome data, but a lack of consensus about benchmarking complicates performance assessment. The Critical Assessment of Metagenome Interpretation (CAMI) challenge has engaged the global developer community to benchmark their programs on highly complex and realistic data sets, generated from ∼700 newly sequenced microorganisms and ∼600 novel viruses and plasmids and representing common experimental setups. Assembly and genome binning programs performed well for species represented by individual genomes but were substantially affected by the presence of related strains. Taxonomic profiling and binning programs were proficient at high taxonomic ranks, with a notable performance decrease below family level. Parameter settings markedly affected performance, underscoring their importance for program reproducibility. The CAMI results highlight current challenges but also provide a roadmap for software selection to answer specific research questions.

Metagenomic Analysis of Healthy and White Plague-Affected Mussismilia braziliensis Corals

Coral health is under threat throughout the world due to regional and global stressors. White plague disease (WP) is one of the most important threats affecting the major reef builder of the Abrolhos Bank in Brazil, the endemic coral Mussismilia braziliensis. We performed a metagenomic analysis of healthy and WP-affected M. braziliensis in order to determine the types of microbes associated with this coral species. We also optimized a protocol for DNA extraction from coral tissues. Our taxonomic analysis revealed Proteobacteria, Bacteroidetes, Firmicutes, Cyanobacteria, and Actinomycetes as the main groups in all healthy and WP-affected corals. Vibrionales, members of the Cytophaga–Flavobacterium–Bacteroides complex, Rickettsiales, and Neisseriales were more abundant in the WP-affected corals. Diseased corals also had more eukaryotic metagenomic sequences identified as Alveolata and Apicomplexa. Our results suggest that WP disease in M. braziliensis is caused by a polymicrobial consortium.

Taxonomic and Functional Microbial Signatures of the Endemic Marine Sponge Arenosclera brasiliensis

The endemic marine sponge Arenosclera brasiliensis (Porifera, Demospongiae, Haplosclerida) is a known source of secondary metabolites such as arenosclerins A-C. In the present study, we established the composition of the A. brasiliensis microbiome and the metabolic pathways associated with this community. We used 454 shotgun pyrosequencing to generate approximately 640,000 high-quality sponge-derived sequences (∼150 Mb). Clustering analysis including sponge, seawater and twenty-three other metagenomes derived from marine animal microbiomes shows that A. brasiliensis contains a specific microbiome. Fourteen bacterial phyla (including Proteobacteria, Cyanobacteria, Actinobacteria, Bacteroidetes, Firmicutes and Cloroflexi) were consistently found in the A. brasiliensis metagenomes. The A. brasiliensis microbiome is enriched for Betaproteobacteria (e.g., Burkholderia) and Gammaproteobacteria (e.g., Pseudomonas and Alteromonas) compared with the surrounding planktonic microbial communities. Functional analysis based on Rapid Annotation using Subsystem Technology (RAST) indicated that the A. brasiliensis microbiome is enriched for sequences associated with membrane transport and one-carbon metabolism. In addition, there was an overrepresentation of sequences associated with aerobic and anaerobic metabolism as well as the synthesis and degradation of secondary metabolites. This study represents the first analysis of sponge-associated microbial communities via shotgun pyrosequencing, a strategy commonly applied in similar analyses in other marine invertebrate hosts, such as corals and algae. We demonstrate that A. brasiliensis has a unique microbiome that is distinct from that of the surrounding planktonic microbes and from other marine organisms, indicating a species-specific microbiome.

Oxygen minimum zones harbour novel viral communities with low diversity

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.

FOCUS: an alignment-free model to identify organisms in metagenomes using non-negative least squares.

One of the major goals in metagenomics is to identify the organisms present in a microbial community from unannotated shotgun sequencing reads. Taxonomic profiling has valuable applications in biological and medical research, including disease diagnostics. Most currently available approaches do not scale well with increasing data volumes, which is important because both the number and lengths of the reads provided by sequencing platforms keep increasing. Here we introduce FOCUS, an agile composition based approach using non-negative least squares (NNLS) to report the organisms present in metagenomic samples and profile their abundances. FOCUS was tested with simulated and real metagenomes, and the results show that our approach accurately predicts the organisms present in microbial communities. FOCUS was implemented in Python. The source code and web-sever are freely available at http://edwards.sdsu.edu/FOCUS.

Critical Assessment of Metagenome Interpretation — a benchmark of metagenomics software

Methods for assembly, taxonomic profiling and binning are key to interpreting metagenome data, but a lack of consensus about benchmarking complicates performance assessment. The Critical Assessment of Metagenome Interpretation (CAMI) challenge has engaged the global developer community to benchmark their programs on highly complex and realistic data sets, generated from ∼700 newly sequenced microorganisms and ∼600 novel viruses and plasmids and representing common experimental setups. Assembly and genome binning programs performed well for species represented by individual genomes but were substantially affected by the presence of related strains. Taxonomic profiling and binning programs were proficient at high taxonomic ranks, with a notable performance decrease below family level. Parameter settings markedly affected performance, underscoring their importance for program reproducibility. The CAMI results highlight current challenges but also provide a roadmap for software selection to answer specific research questions.

Genomic Taxonomy of the Genus Prochlorococcus

The genus Prochlorococcus is globally abundant and dominates the total phytoplankton biomass and production in the oligotrophic ocean. The single species, Prochlorococcus marinus, comprises six named ecotypes. Our aim was to analyze the taxonomic structure of the genus Prochlorococcus. We analyzed the complete genomes of 13 cultured P. marinus type and reference strains by means of several genomic taxonomy tools (i.e., multilocus sequence analysis, amino acid identity, Karlin genomic signature, and genome to genome distance). In addition, we estimated the diversity of Prochlorococcus species in over 100 marine metagenomes from all the major oceanic provinces. According to our careful taxonomic analysis, the 13 strains corresponded, in fact, to ten different Prochlorococcus species. This analysis establishes a new taxonomic framework for the genus Prochlorococcus. Further, the analysis of the metagenomic data suggests that, in total, there may only be 35 Prochlorococcus species in the worlds oceans. We propose that the dearth of species observed in this study is driven by high selective pressures that limit diversification in the global ocean.

SUPER-FOCUS: a tool for agile functional analysis of shotgun metagenomic data.

Summary: Analyzing the functional profile of a microbial community from unannotated shotgun sequencing reads is one of the important goals in metagenomics. Functional profiling has valuable applications in biological research because it identifies the abundances of the functional genes of the organisms present in the original sample, answering the question what they can do. Currently, available tools do not scale well with increasing data volumes, which is important because both the number and lengths of the reads produced by sequencing platforms keep increasing. Here, we introduce SUPER-FOCUS, SUbsystems Profile by databasE Reduction using FOCUS, an agile homology-based approach using a reduced reference database to report the subsystems present in metagenomic datasets and profile their abundances. SUPER-FOCUS was tested with over 70 real metagenomes, the results showing that it accurately predicts the subsystems present in the profiled microbial communities, and is up to 1000 times faster than other tools. Availability and implementation: SUPER-FOCUS was implemented in Python, and its source code and the tool website are freely available at https://edwards.sdsu.edu/SUPERFOCUS. Contact: redwards@mail.sdsu.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Multi-Analytical Approach Reveals Potential Microbial Indicators in Soil for Sugarcane Model Systems

This study focused on the effects of organic and inorganic amendments and straw retention on the microbial biomass (MB) and taxonomic groups of bacteria in sugarcane-cultivated soils in a greenhouse mesocosm experiment monitored for gas emissions and chemical factors. The experiment consisted of combinations of synthetic nitrogen (N), vinasse (V; a liquid waste from ethanol production), and sugarcane-straw blankets. Increases in CO2-C and N2O-N emissions were identified shortly after the addition of both N and V to the soils, thus increasing MB nitrogen (MB-N) and decreasing MB carbon (MB-C) in the N+V-amended soils and altering soil chemical factors that were correlated with the MB. Across 57 soil metagenomic datasets, Actinobacteria (31.5%), Planctomycetes (12.3%), Deltaproteobacteria (12.3%), Alphaproteobacteria (12.0%) and Betaproteobacteria (11.1%) were the most dominant bacterial groups during the experiment. Differences in relative abundance of metagenomic sequences were mainly revealed for Acidobacteria, Actinobacteria, Gammaproteobacteria and Verrucomicrobia with regard to N+V fertilization and straw retention. Differential abundances in bacterial groups were confirmed using 16S rRNA gene-targeted phylum-specific primers for real-time PCR analysis in all soil samples, whose results were in accordance with sequence data, except for Gammaproteobacteria. Actinobacteria were more responsive to straw retention with Rubrobacterales, Bifidobacteriales and Actinomycetales related to the chemical factors of N+V-amended soils. Acidobacteria subgroup 7 and Opitutae, a verrucomicrobial class, were related to the chemical factors of soils without straw retention as a surface blanket. Taken together, the results showed that MB-C and MB-N responded to changes in soil chemical factors and CO2-C and N2O-N emissions, especially for N+V-amended soils. The results also indicated that several taxonomic groups of bacteria, such as Acidobacteria, Actinobacteria and Verrucomicrobia, and their subgroups acted as early-warning indicators of N+V amendments and straw retention in sugarcane-cultivated soils, which can alter the soil chemical factors.