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.

Using viromes to predict novel immune proteins in non-model organisms.

Immunity is mostly studied in a few model organisms, leaving the majority of immune systems on the planet unexplored. To characterize the immune systems of non-model organisms alternative approaches are required. Viruses manipulate host cell biology through the expression of proteins that modulate the immune response. We hypothesized that metagenomic sequencing of viral communities would be useful to identify both known and unknown host immune proteins. To test this hypothesis, a mock human virome was generated and compared to the human proteome using tBLASTn, resulting in 36 proteins known to be involved in immunity. This same pipeline was then applied to reef-building coral, a non-model organism that currently lacks traditional molecular tools like transgenic animals, gene-editing capabilities, and in vitro cell cultures. Viromes isolated from corals and compared with the predicted coral proteome resulted in 2503 coral proteins, including many proteins involved with pathogen sensing and apoptosis. There were also 159 coral proteins predicted to be involved with coral immunity but currently lacking any functional annotation. The pipeline described here provides a novel method to rapidly predict host immune components that can be applied to virtually any system with the potential to discover novel immune proteins.

Acidobacteria Subgroups and Their Metabolic Potential for Carbon Degradation in Sugarcane Soil Amended With Vinasse and Nitrogen Fertilizers.

Acidobacteria is a predominant bacterial phylum in tropical agricultural soils, including sugarcane cultivated soils. The increased need for fertilizers due to the expansion of sugarcane production is a threat to the ability of the soil to maintain its potential for self-regulation in the long term, in witch carbon degradation has essential role. In this study, a culture-independent approach based on high-throughput DNA sequencing and microarray technology was used to perform taxonomic and functional profiling of the Acidobacteria community in a tropical soil under sugarcane (Saccharum spp.) that was supplemented with nitrogen (N) combined with vinasse. These analyses were conducted to identify the subgroup-level responses to chemical changes and the carbon (C) degradation potential of the different Acidobacteria subgroups. Eighteen Acidobacteria subgroups from a total of 26 phylogenetically distinct subgroups were detected based on high-throughput DNA sequencing, and 16 gene families associated with C degradation were quantified using Acidobacteria-derived DNA microarray probes. The subgroups Gp13 and Gp18 presented the most positive correlations with the gene families associated with C degradation, especially those involved in hemicellulose degradation. However, both subgroups presented low abundance in the treatment containing vinasse. In turn, the Gp4 subgroup was the most abundant in the treatment that received vinasse, but did not present positive correlations with the gene families for C degradation analyzed in this study. The metabolic potential for C degradation of the different Acidobacteria subgroups in sugarcane soil amended with N and vinasse can be driven in part through the increase in soil nutrient availability, especially calcium (Ca), magnesium (Mg), potassium (K), aluminum (Al), boron (B) and zinc (Zn). This soil management practice reduces the abundance of Acidobacteria subgroups, including those potentially involved with C degradation in this agricultural soil.

Cyanobacterial biodiversity of semiarid public drinking water supply reservoirs assessed via next-generation DNA sequencing technology.

Next-generation DNA sequencing technology was applied to generate molecular data from semiarid reservoirs during well-defined seasons. Target sequences of 16S-23S rRNA ITS and cpcBA-IGS were used to reveal the taxonomic groups of cyanobacteria present in the samples, and genes coding for cyanotoxins such as microcystins (mcyE), saxitoxins (sxtA), and cylindrospermopsins (cyrJ) were investigated. The presence of saxitoxins in the environmental samples was evaluated using ELISA kit. Taxonomic analyses of high-throughput DNA sequencing data showed the dominance of the genus Microcystis in Mundaú reservoir. Furthermore, it was the most abundant genus in the dry season in Ingazeira reservoir. In the rainy season, 16S-23S rRNA ITS analysis revealed that Cylindrospermopsis raciborskii comprised 46.8% of the cyanobacterial community in Ingazeira reservoir, while the cpcBAIGS region revealed that C. raciborskii (31.8%) was the most abundant taxon followed by Sphaerospermopsis aphanizomenoides (17.3%) and Planktothrix zahidii (16.6%). Despite the presence of other potential toxin-producing genera, the detected sxtA gene belonged to C. raciborskii, while the mcyE gene belonged to Microcystis in both reservoirs. The detected mcyE gene had good correlation with MC content, while the amplification of the sxtA gene was related to the presence of STX. The cyrJ gene was not detected in these samples. Using DNA analyses, our results showed that the cyanobacterial composition of Mundaú reservoir was similar in successive dry seasons, and it varied between seasons in Ingazeira reservoir. In addition, our data suggest that some biases of analysis influenced the cyanobacterial communities seen in the NGS output of Ingazeira reservoir.

Diel population and functional synchrony of microbial communities on coral reefs.

On coral reefs, microorganisms are essential for recycling nutrients to primary producers through the remineralization of benthic-derived organic matter. Diel investigations of reef processes are required to holistically understand the functional roles of microbial players in these ecosystems. Here we report a metagenomic analysis characterizing microbial communities in the water column overlying 16 remote forereef sites over a diel cycle. Our results show that microbial community composition is more dissimilar between day and night samples collected from the same site than between day or night samples collected across geographically distant reefs. Diel community differentiation is largely driven by the flux of Psychrobacter sp., which is two-orders of magnitude more abundant during the day. Nighttime communities are enriched with species of Roseobacter, Halomonas, and Alteromonas encoding a greater variety of pathways for carbohydrate catabolism, further illustrating temporal patterns of energetic provisioning between different marine microbes. Dynamic diel fluctuations of microbial populations could also support the efficient trophic transfer of energy posited in coral reef food webs.

An Agile Functional Analysis of Metagenomic Data Using SUPER-FOCUS.

One of the main goals in metagenomics is to identify the functional profile of a microbial community from unannotated shotgun sequencing reads. Functional annotation is important in biological research because it enables researchers to identify the abundance of functional genes of the organisms present in the sample, answering the question, "What can the organisms in the sample do?" Most currently available approaches do not scale with increasing data volumes, which is important because both the number and lengths of the reads provided by sequencing platforms keep increasing. Here, we present 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 real metagenomes, and the results show that it accurately predicts the subsystems present in the profiled microbial communities, is computationally efficient, and up to 1000 times faster than other tools. SUPER-FOCUS is freely available at http://edwards.sdsu.edu/SUPERFOCUS .

Draft Genome Sequence of Cylindrospermopsis raciborskii (Cyanobacteria) Strain ITEP-A1 Isolated from a Brazilian Semiarid Freshwater Body: Evidence of Saxitoxin and Cylindrospermopsin Synthetase Genes.

Cylindrospermopsis raciborskii ITEP-A1 is a saxitoxin-producing cyanobacterium. We report the draft genome sequence of ITEP-A1, which comprised 195 contigs that were assembled with SPAdes and annotated with Rapid Annotation using Subsystem Technology. The identified genome sequence had 3,605,836 bp, 40.1% G+C, and predicted 3,553 coding sequences (including the synthetase genes).

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.

Sequencing at sea: challenges and experiences in Ion Torrent PGM sequencing during the 2013 Southern Line Islands Research Expedition.

Genomics and metagenomics have revolutionized our understanding of marine microbial ecology and the importance of microbes in global geochemical cycles. However, the process of DNA sequencing has always been an abstract extension of the research expedition, completed once the samples were returned to the laboratory. During the 2013 Southern Line Islands Research Expedition, we started the first effort to bring next generation sequencing to some of the most remote locations on our planet. We successfully sequenced twenty six marine microbial genomes, and two marine microbial metagenomes using the Ion Torrent PGM platform on the Merchant Yacht Hanse Explorer. Onboard sequence assembly, annotation, and analysis enabled us to investigate the role of the microbes in the coral reef ecology of these islands and atolls. This analysis identified phosphonate as an important phosphorous source for microbes growing in the Line Islands and reinforced the importance of L-serine in marine microbial ecosystems. Sequencing in the field allowed us to propose hypotheses and conduct experiments and further sampling based on the sequences generated. By eliminating the delay between sampling and sequencing, we enhanced the productivity of the research expedition. By overcoming the hurdles associated with sequencing on a boat in the middle of the Pacific Ocean we proved the flexibility of the sequencing, annotation, and analysis pipelines.