Community cyberinfrastructure for Advanced Microbial Ecology Research and Analysis: the CAMERA resource

Open ocean and coastal strains of the N2-fixing cyanobacterium UCYN-A have distinct transcriptomes

Decades of research on marine N2 fixation focused on Trichodesmium, which are generally free-living cyanobacteria, but in recent years the endosymbiotic cyanobacterium Candidatus Atelocyanobacterium thalassa (UCYN-A) has received increasing attention. However, few studies have shed light on the influence of the host versus the habitat on UCYN-A N2 fixation and overall metabolism. Here we compared transcriptomes from natural populations of UCYN-A from oligotrophic open-ocean versus nutrient-rich coastal waters, using a microarray that targets the full genomes of UCYN-A1 and UCYN-A2 and known genes for UCYN-A3. We found that UCYN-A2, usually regarded as adapted to coastal environments, was transcriptionally very active in the open ocean and appeared to be less impacted by habitat change than UCYN-A1. Moreover, for genes with 24 h periodic expression we observed strong but inverse correlations among UCYN-A1, A2, and A3 to oxygen and chlorophyll, which suggests distinct host-symbiont relationships. Across habitats and sublineages, genes for N2 fixation and energy production had high transcript levels, and, intriguingly, were among the minority of genes that kept the same schedule of diel expression. This might indicate different regulatory mechanisms for genes that are critical to the symbiosis for the exchange of nitrogen for carbon from the host. Our results underscore the importance of N2 fixation in UCYN-A symbioses across habitats, with consequences for community interactions and global biogeochemical cycles.

Role of plant growth-promoting rhizobacteria in boosting the phytoremediation of stressed soils: Opportunities, challenges, and prospects

Soil is considered as a vital natural resource equivalent to air and water which supports growth of the plants and provides habitats to microorganisms. Changes in soil properties, productivity, and, inevitably contamination/stress are the result of urbanisation, industrialization, and long-term use of synthetic fertiliser. Therefore, in the recent scenario, reclamation of contaminated/stressed soils has become a potential challenge. Several customized, such as, physical, chemical, and biological technologies have been deployed so far to restore contaminated land. Among them, microbial-assisted phytoremediation is considered as an economical and greener approach. In recent decades, soil microbes have successfully been used to improve plants' ability to tolerate biotic and abiotic stress and strengthen their phytoremediation capacity. Therefore, in this context, the current review work critically explored the microbial assisted phytoremediation mechanisms to restore different types of stressed soil. The role of plant growth-promoting rhizobacteria (PGPR) and their potential mechanisms that foster plants' growth and also enhance phytoremediation capacity are focussed. Finally, this review has emphasized on the application of advanced tools and techniques to effectively characterize potent soil microbial communities and their significance in boosting the phytoremediation process of stressed soils along with prospects for future research.

Microbial Involvement in the Bioremediation of Total Petroleum Hydrocarbon Polluted Soils: Challenges and Perspectives

Nowadays, soil contamination by total petroleum hydrocarbons is still one of the most widespread forms of contamination. Intervention technologies are consolidated; however, full-scale interventions turn out to be not sustainable. Sustainability is essential not only in terms of costs, but also in terms of restoration of the soil resilience. Bioremediation has the possibility to fill the gap of sustainability with proper knowledge. Bioremediation should be optimized by the exploitation of the recent “omic” approaches to the study of hydrocarburoclastic microbiomes. To reach the goal, an extensive and deep knowledge in the study of bacterial and fungal degradative pathways, their interactions within microbiomes and of microbiomes with the soil matrix has to be gained. “Omic” approaches permits to study both the culturable and the unculturable soil microbial communities active in degradation processes, offering the instruments to identify the key organisms responsible for soil contaminant depletion and restoration of soil resilience. Tools for the investigation of both microbial communities, their degradation pathways and their interaction, will be discussed, describing the dedicated genomic and metagenomic approaches, as well as the interpretative tools of the deriving data, that are exploitable for both optimizing bio-based approaches for the treatment of total petroleum hydrocarbon contaminated soils and for the correct scaling up of the technologies at the industrial scale.

A Global Approach to Estimating the Abundance and Duplication of Polyketide Synthase Domains in Dinoflagellates

Many dinoflagellate species make toxins in a myriad of different molecular configurations but the underlying chemistry in all cases is presumably via modular synthases, primarily polyketide synthases. In many organisms modular synthases occur as discrete synthetic genes or domains within a gene that act in coordination thus forming a module that produces a particular fragment of a natural product. The modules usually occur in tandem as gene clusters with a syntenic arrangement that is often predictive of the resultant structure. Dinoflagellate genomes however are notoriously complex with individual genes present in many tandem repeats and very few synthetic modules occurring as gene clusters, unlike what has been seen in bacteria and fungi. However, modular synthesis in all organisms requires a free thiol group that acts as a carrier for sequential synthesis called a thiolation domain. We scanned 47 dinoflagellate transcriptomes for 23 modular synthase domain models and compared their abundance among 10 orders of dinoflagellates as well as their co-occurrence with thiolation domains. The total count of domain types was quite large with over thirty-thousand identified, 29 000 of which were in the core dinoflagellates. Although there were no specific trends in domain abundance associated with types of toxins, there were readily observable lineage specific differences. The Gymnodiniales, makers of long polyketide toxins such as brevetoxin and karlotoxin had a high relative abundance of thiolation domains as well as multiple thiolation domains within a single transcript. Orders such as the Gonyaulacales, makers of small polyketides such as spirolides, had fewer thiolation domains but a relative increase in the number of acyl transferases. Unique to the core dinoflagellates, however, were thiolation domains occurring alongside tetratricopeptide repeats that facilitate protein-protein interactions, especially hexa and hepta-repeats, that may explain the scaffolding required for synthetic complexes capable of making large toxins. Clustering analysis for each type of domain was also used to discern possible origins of duplication for the multitude of single domain transcripts. Single domain transcripts frequently clustered with synonymous domains from multi-domain transcripts such as the BurA and ZmaK like genes as well as the multi-ketosynthase genes, sometimes with a large degree of apparent gene duplication, while fatty acid synthesis genes formed distinct clusters. Surprisingly the acyl-transferases and ketoreductases involved in fatty acid synthesis (FabD and FabG, respectively) were found in very large clusters indicating an unprecedented degree of gene duplication for these genes. These results demonstrate a complex evolutionary history of core dinoflagellate modular synthases with domain specific duplications throughout the lineage as well as clues to how large protein complexes can be assembled to synthesize the largest natural products known.

Development of an Ocean Protein Portal for Interactive Discovery and Education

Proteins are critical in catalyzing chemical reactions, forming key cellular structures, and in regulating cellular processes. Investigation of marine microbial proteins by metaproteomics methods enables the discovery of numerous aspects of microbial biogeochemical processes. However, these datasets present big data challenges as they often involve many samples collected across broad geospatial and temporal scales, resulting in thousands of protein identifications, abundances, and corresponding annotation information. The Ocean Protein Portal (OPP) was created to enable data sharing and discovery among multiple scientific domains and serve both research and education functions. The portal focuses on three use case questions: "Where is my protein of interest?", "Who makes it?", and "How much is there?" and provides profile and section visualizations, real-time taxonomic analysis, and links to metadata, sequence analysis, and other external resources to enable connections to be made between biogeochemical and proteomics datasets.

Uncovering the genomic potential of the Amazon River microbiome to degrade rainforest organic matter

BACKGROUND

The Amazon River is one of the largest in the world and receives huge amounts of terrestrial organic matter (TeOM) from the surrounding rainforest. Despite this TeOM is typically recalcitrant (i.e. resistant to degradation), only a small fraction of it reaches the ocean, pointing to a substantial TeOM degradation by the river microbiome. Yet, microbial genes involved in TeOM degradation in the Amazon River were barely known. Here, we examined the Amazon River microbiome by analysing 106 metagenomes from 30 sampling points distributed along the river.

RESULTS

We constructed the Amazon River basin Microbial non-redundant Gene Catalogue (AMnrGC) that includes ~ 3.7 million non-redundant genes, affiliating mostly to bacteria. We found that the Amazon River microbiome contains a substantial gene-novelty compared to other relevant known environments (rivers and rainforest soil). Genes encoding for proteins potentially involved in lignin degradation pathways were correlated to tripartite tricarboxylates transporters and hemicellulose degradation machinery, pointing to a possible priming effect. Based on this, we propose a model on how the degradation of recalcitrant TeOM could be modulated by labile compounds in the Amazon River waters. Our results also suggest changes of the microbial community and its genomic potential along the river course.

CONCLUSIONS

Our work contributes to expand significantly our comprehension of the world's largest river microbiome and its potential metabolism related to TeOM degradation. Furthermore, the produced gene catalogue (AMnrGC) represents an important resource for future research in tropical rivers. Video abstract.

Resources and tools for the high-throughput, multi-omic study of intestinal microbiota

The human gut microbiome impacts several aspects of human health and disease, including digestion, drug metabolism and the propensity to develop various inflammatory, autoimmune and metabolic diseases. Many of the molecular processes that play a role in the activity and dynamics of the microbiota go beyond species and genic composition and thus, their understanding requires advanced bioinformatics support. This article aims to provide an up-to-date view of the resources and software tools that are being developed and used in human gut microbiome research, in particular data integration and systems-level analysis efforts. These efforts demonstrate the power of standardized and reproducible computational workflows for integrating and analysing varied omics data and gaining deeper insights into microbe community structure and function as well as host-microbe interactions.

The diversity and abundance of fungi and bacteria on the healthy and dandruff affected human scalp

Dandruff is a skin condition that affects the scalp of up to half the world's population, it is characterised by an itchy, flaky scalp and is associated with colonisation of the skin by Malassezia spp. Management of this condition is typically via antifungal therapies, however the precise role of microbes in the aggravation of the condition are incompletely characterised. Here, a combination of 454 sequencing and qPCR techniques were used to compare the scalp microbiota of dandruff and non-dandruff affected Chinese subjects. Based on 454 sequencing of the scalp microbiome, the two most abundant bacterial genera found on the scalp surface were Cutibacterium (formerly Propionibacterium) and Staphylococcus, while Malassezia was the main fungal inhabitant. Quantitative PCR (qPCR) analysis of four scalp taxa (M. restricta, M. globosa, C. acnes and Staphylococcus spp.) believed to represent the bulk of the overall population was additionally carried out. Metataxonomic and qPCR analyses were performed on healthy and lesional buffer scrub samples to facilitate assessment of whether the scalp condition is associated with differential microbial communities on the sampled skin. Dandruff was associated with greater frequencies of M. restricta and Staphylococcus spp. compared with the healthy population (p<0.05). Analysis also revealed the presence of an unclassified fungal taxon that could represent a novel Malassezia species.

MNEMONIC: MetageNomic Experiment Mining to create an OTU Network of Inhabitant Correlations

Background

The number of publicly available metagenomic experiments in various environments has been rapidly growing, empowering the potential to identify similar shifts in species abundance between different experiments. This could be a potentially powerful way to interpret new experiments, by identifying common themes and causes behind changes in species abundance.

Results

We propose a novel framework for comparing microbial shifts between conditions. Using data from one of the largest human metagenome projects to date, the American Gut Project (AGP), we obtain differential abundance vectors for microbes using experimental condition information provided with the AGP metadata, such as patient age, dietary habits, or health status. We show it can be used to identify similar and opposing shifts in microbial species, and infer putative interactions between microbes. Our results show that groups of shifts with similar effects on microbiome can be identified and that similar dietary interventions display similar microbial abundance shifts.

Conclusions

Without comparison to prior data, it is difficult for experimentalists to know if their observed changes in species abundance have been observed by others, both in their conditions and in others they would never consider comparable. Yet, this can be a very important contextual factor in interpreting the significance of a shift. We’ve proposed and tested an algorithmic solution to this problem, which also allows for comparing the metagenomic signature shifts between conditions in the existing body of data.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2623-x) contains supplementary material, which is available to authorized users.

The Transcriptional Cycle Is Suited to Daytime N2 Fixation in the Unicellular Cyanobacterium "Candidatus Atelocyanobacterium thalassa" (UCYN-A)

The symbiotic N₂-fixing cyanobacterium UCYN-A, which is closely related to Braarudosphaera bigelowii, and its eukaryotic algal host have been shown to be globally distributed and important in open-ocean N₂ fixation. These unique cyanobacteria have reduced metabolic capabilities, even lacking genes for oxygenic photosynthesis and carbon fixation. Cyanobacteria generally use energy from photosynthesis for nitrogen fixation but require mechanisms for avoiding inactivation of the oxygen-sensitive nitrogenase enzyme by ambient oxygen (O₂) or the O₂ evolved through photosynthesis. This study showed that symbiosis between the N₂-fixing cyanobacterium UCYN-A and its eukaryotic algal host has led to adaptation of its daily gene expression pattern in order to enable daytime aerobic N₂ fixation, which is likely more energetically efficient than fixing N₂ at night, as found in other unicellular marine cyanobacteria.

Symbiosis between a marine alga and a N₂-fixing cyanobacterium (Cyanobacterium UCYN-A) is geographically widespread in the oceans and is important in the marine N cycle. UCYN-A is uncultivated and is an unusual unicellular cyanobacterium because it lacks many metabolic functions, including oxygenic photosynthesis and carbon fixation, which are typical in cyanobacteria. It is now presumed to be an obligate symbiont of haptophytes closely related to Braarudosphaera bigelowii. N₂-fixing cyanobacteria use different strategies to avoid inhibition of N₂ fixation by the oxygen evolved in photosynthesis. Most unicellular cyanobacteria temporally separate the two incompatible activities by fixing N₂ only at night, but, surprisingly, UCYN-A appears to fix N₂ during the day. The goal of this study was to determine how the unicellular UCYN-A strain coordinates N₂ fixation and general metabolism compared to other marine cyanobacteria. We found that UCYN-A has distinct daily cycles of many genes despite the fact that it lacks two of the three circadian clock genes found in most cyanobacteria. We also found that the transcription patterns in UCYN-A are more similar to those in marine cyanobacteria that are capable of aerobic N₂ fixation in the light, such as Trichodesmium and heterocyst-forming cyanobacteria, than to those in Crocosphaera or Cyanothece species, which are more closely related to unicellular marine cyanobacteria evolutionarily. Our findings suggest that the symbiotic interaction has resulted in a shift of transcriptional regulation to coordinate UCYN-A metabolism with that of the phototrophic eukaryotic host, thus allowing efficient coupling of N₂ fixation (by the cyanobacterium) to the energy obtained from photosynthesis (by the eukaryotic unicellular alga) in the light.

GLOSSary: the GLobal Ocean 16S subunit web accessible resource

BACKGROUND

Environmental metagenomics is a challenging approach that is exponentially spreading in the scientific community to investigate taxonomic diversity and possible functions of the biological components. The massive amount of sequence data produced, often endowed with rich environmental metadata, needs suitable computational tools to fully explore the embedded information. Bioinformatics plays a key role in providing methodologies to manage, process and mine molecular data, integrated with environmental metagenomics collections. One such relevant example is represented by the Tara Ocean Project.

RESULTS

We considered the Tara 16S miTAGs released by the consortium, representing raw sequences from a shotgun metagenomics approach with similarities to 16S rRNA genes. We generated assembled 16S rDNA sequences, which were classified according to their lengths, the possible presence of chimeric reads, the putative taxonomic affiliation. The dataset was included in GLOSSary (the GLobal Ocean 16S Subunit web accessible resource), a bioinformatics platform to organize environmental metagenomics data. The aims of this work were: i) to present alternative computational approaches to manage challenging metagenomics data; ii) to set up user friendly web-based platforms to allow the integration of environmental metagenomics sequences and of the associated metadata; iii) to implement an appropriate bioinformatics platform supporting the analysis of 16S rDNA sequences exploiting reference datasets, such as the SILVA database. We organized the data in a next-generation NoSQL "schema-less" database, allowing flexible organization of large amounts of data and supporting native geospatial queries. A web interface was developed to permit an interactive exploration and a visual geographical localization of the data, either raw miTAG reads or 16S contigs, from our processing pipeline. Information on unassembled sequences is also available. The taxonomic affiliations of contigs and miTAGs, and the spatial distribution of the sampling sites and their associated sequence libraries, as they are contained in the Tara metadata, can be explored by a query interface, which allows both textual and visual investigations. In addition, all the sequence data were made available for a dedicated BLAST-based web application alongside the SILVA collection.

CONCLUSIONS

GLOSSary provides an expandable bioinformatics environment, able to support the scientific community in current and forthcoming environmental metagenomics analyses.

Comparative Genomics Highlights Symbiotic Capacities and High Metabolic Flexibility of the Marine Genus Pseudovibrio

Pseudovibrio is a marine bacterial genus members of which are predominantly isolated from sessile marine animals, and particularly sponges. It has been hypothesized that Pseudovibrio spp. form mutualistic relationships with their hosts. Here, we studied Pseudovibrio phylogeny and genetic adaptations that may play a role in host colonization by comparative genomics of 31 Pseudovibrio strains, including 25 sponge isolates. All genomes were highly similar in terms of encoded core metabolic pathways, albeit with substantial differences in overall gene content. Based on gene composition, Pseudovibrio spp. clustered by geographic region, indicating geographic speciation. Furthermore, the fact that isolates from the Mediterranean Sea clustered by sponge species suggested host-specific adaptation or colonization. Genome analyses suggest that Pseudovibrio hongkongensis UST20140214-015BT is only distantly related to other Pseudovibrio spp., thereby challenging its status as typical Pseudovibrio member. All Pseudovibrio genomes were found to encode numerous proteins with SEL1 and tetratricopeptide repeats, which have been suggested to play a role in host colonization. For evasion of the host immune system, Pseudovibrio spp. may depend on type III, IV, and VI secretion systems that can inject effector molecules into eukaryotic cells. Furthermore, Pseudovibrio genomes carry on average seven secondary metabolite biosynthesis clusters, reinforcing the role of Pseudovibrio spp. as potential producers of novel bioactive compounds. Tropodithietic acid, bacteriocin, and terpene biosynthesis clusters were highly conserved within the genus, suggesting an essential role in survival, for example through growth inhibition of bacterial competitors. Taken together, these results support the hypothesis that Pseudovibrio spp. have mutualistic relations with sponges.

Maser: one-stop platform for NGS big data from analysis to visualization

A major challenge in analyzing the data from high-throughput next-generation sequencing (NGS) is how to handle the huge amounts of data and variety of NGS tools and visualize the resultant outputs. To address these issues, we developed a cloud-based data analysis platform, Maser (Management and Analysis System for Enormous Reads), and an original genome browser, Genome Explorer (GE). Maser enables users to manage up to 2 terabytes of data to conduct analyses with easy graphical user interface operations and offers analysis pipelines in which several individual tools are combined as a single pipeline for very common and standard analyses. GE automatically visualizes genome assembly and mapping results output from Maser pipelines, without requiring additional data upload. With this function, the Maser pipelines can graphically display the results output from all the embedded tools and mapping results in a web browser. Therefore Maser realized a more user-friendly analysis platform especially for beginners by improving graphical display and providing the selected standard pipelines that work with built-in genome browser. In addition, all the analyses executed on Maser are recorded in the analysis history, helping users to trace and repeat the analyses. The entire process of analysis and its histories can be shared with collaborators or opened to the public. In conclusion, our system is useful for managing, analyzing, and visualizing NGS data and achieves traceability, reproducibility, and transparency of NGS analysis.

Database URL: http://cell-innovation.nig.ac.jp/maser/

A multi-source domain annotation pipeline for quantitative metagenomic and metatranscriptomic functional profiling

BACKGROUND

Biochemical and regulatory pathways have until recently been thought and modelled within one cell type, one organism and one species. This vision is being dramatically changed by the advent of whole microbiome sequencing studies, revealing the role of symbiotic microbial populations in fundamental biochemical functions. The new landscape we face requires the reconstruction of biochemical and regulatory pathways at the community level in a given environment. In order to understand how environmental factors affect the genetic material and the dynamics of the expression from one environment to another, we want to evaluate the quantity of gene protein sequences or transcripts associated to a given pathway by precisely estimating the abundance of protein domains, their weak presence or absence in environmental samples.

RESULTS

MetaCLADE is a novel profile-based domain annotation pipeline based on a multi-source domain annotation strategy. It applies directly to reads and improves identification of the catalog of functions in microbiomes. MetaCLADE is applied to simulated data and to more than ten metagenomic and metatranscriptomic datasets from different environments where it outperforms InterProScan in the number of annotated domains. It is compared to the state-of-the-art non-profile-based and profile-based methods, UProC and HMM-GRASPx, showing complementary predictions to UProC. A combination of MetaCLADE and UProC improves even further the functional annotation of environmental samples.

CONCLUSIONS

Learning about the functional activity of environmental microbial communities is a crucial step to understand microbial interactions and large-scale environmental impact. MetaCLADE has been explicitly designed for metagenomic and metatranscriptomic data and allows for the discovery of patterns in divergent sequences, thanks to its multi-source strategy. MetaCLADE highly improves current domain annotation methods and reaches a fine degree of accuracy in annotation of very different environments such as soil and marine ecosystems, ancient metagenomes and human tissues.

Stoichio-Metagenomics of Ocean Waters: A Molecular Evolution Approach to Trace the Dynamics of Nitrogen Conservation in Natural Communities

Nitrogen is crucially limiting in ocean surface waters, and its availability varies substantially with coastal regions typically richer in nutrients than open oceans. In a biological stoichiometry framework, a parsimonious strategy of nitrogen allocation predicts nitrogen content of proteins to be lower in communities adapted to open ocean than to coastal regions. To test this hypothesis we have directly interrogated marine microbial communities, using a series of metagenomics datasets with a broad geographical distribution from the Global Ocean Sampling Expedition. Analyzing over 20 million proteins, we document a ubiquitous signal of nitrogen conservation in open ocean communities, both in membrane and non-membrane proteins. Efficient nitrogen allocation is expected to specifically target proteins that are expressed at high rate in response to nitrogen starvation. Furthermore, in order to preserve protein functional efficiency, economic nitrogen allocation is predicted to target primarily the least functionally constrained regions of proteins. Contrasting the NtcA-induced pathway, typically up-regulated in response to nitrogen starvation, with the arginine anabolic pathway, which is instead up-regulated in response to nitrogen abundance, we show how both these predictions are fulfilled. Using evolutionary rates as an informative proxy of functional constraints, we show that variation in nitrogen allocation between open ocean and coastal communities is primarily localized in the least functionally constrained regions of the genes triggered by NtcA. As expected, such a pattern is not detectable in the genes involved in the arginine anabolic pathway. These results directly link environmental nitrogen availability to different adaptive strategies of genome evolution, and emphasize the relevance of the material costs of evolutionary change in natural ecosystems.

Viruses of the Nahant Collection, characterization of 251 marine Vibrionaceae viruses

Viruses are highly discriminating in their interactions with host cells and are thought to play a major role in maintaining diversity of environmental microbes. However, large-scale ecological and genomic studies of co-occurring virus-host pairs, required to characterize the mechanistic and genomic foundations of virus-host interactions, are lacking. Here, we present the largest dataset of cultivated and sequenced co-occurring virus-host pairs that captures ecologically representative fine-scale diversity. Using the ubiquitous and ecologically diverse marine Vibrionaceae as a host platform, we isolate and sequence 251 dsDNA viruses and their hosts from three time points within a 93-day time-series study. The virus collection includes representatives of the three Caudovirales tailed virus morphotypes, a novel family of nontailed viruses, and the smallest (10,046 bp) and largest (348,911 bp) Vibrio virus genomes described. We provide general characterization and annotation of the viruses and describe read-mapping protocols to standardize genome presentation. The rich ecological and genomic contextualization of hosts and viruses make the Nahant Collection a unique platform for high-resolution studies of environmental virus-host infection networks.

Taxon Richness of "Megaviridae" Exceeds those of Bacteria and Archaea in the Ocean

Since the discovery of the giant mimivirus, evolutionarily related viruses have been isolated or identified from various environments. Phylogenetic analyses of this group of viruses, tentatively referred to as the family “Megaviridae”, suggest that it has an ancient origin that may predate the emergence of major eukaryotic lineages. Environmental genomics has since revealed that Megaviridae represents one of the most abundant and diverse groups of viruses in the ocean. In the present study, we compared the taxon richness and phylogenetic diversity of Megaviridae, Bacteria, and Archaea using DNA-dependent RNA polymerase as a common marker gene. By leveraging existing microbial metagenomic data, we found higher richness and phylogenetic diversity in this single viral family than in the two prokaryotic domains. We also obtained results showing that the evolutionary rate alone cannot account for the observed high diversity of Megaviridae lineages. These results suggest that the Megaviridae family has a deep co-evolutionary history with diverse marine protists since the early “Big-Bang” radiation of the eukaryotic tree of life.

Computational operon prediction in whole-genomes and metagenomes

Microbial diversity in unique environmental settings enables abrupt responses catalysed by altering the gene regulation and formation of gene clusters called operons. Operons increases bacterial adaptability, which in turn increases their survival. This review article presents the emergence of computational operon prediction methods for whole microbial genomes and metagenomes, and discusses their strengths and limitations. Most of the whole-genome operon prediction methods struggle to generalize on unrelated genomes. The applicability of universal whole-genome operon prediction methods to metagenomic data is an interesting yet less investigated question. We have evaluated the potential of various operon prediction features for genomic and metagenomic data. Most of operon prediction methods with high accuracy have been compiled into databases. Despite of the high predictive performance, the data among many databases are not completely consistent for similar species. We performed a correlation analysis between the computationally predicted operon databases and experimentally validated data for Escherichia coli, Bacillus subtilis and Mycobacterium tuberculosis. Operon prediction for most of the less characterized microbes cannot be verified due to absence of experimentally validated operons. The generation of validated information for other microbes would test the authenticity of operon databases for other less annotated microbes as well. Advances in sequencing technologies and development of better analysis methods will help researchers to overcome the technological hurdles (such as long sequencing reads and improved contig size) and further improve operon predictions and better utilize operonic information.

Flexible metagenome analysis using the MGX framework

BACKGROUND

The characterization of microbial communities based on sequencing and analysis of their genetic information has become a popular approach also referred to as metagenomics; in particular, the recent advances in sequencing technologies have enabled researchers to study even the most complex communities. Metagenome analysis, the assignment of sequences to taxonomic and functional entities, however, remains a tedious task: large amounts of data need to be processed. There are a number of approaches addressing particular aspects, but scientific questions are often too specific to be answered by a general-purpose method.

RESULTS

We present MGX, a flexible and extensible client/server-framework for the management and analysis of metagenomic datasets; MGX features a comprehensive set of adaptable workflows required for taxonomic and functional metagenome analysis, combined with an intuitive and easy-to-use graphical user interface offering customizable result visualizations. At the same time, MGX allows to include own data sources and devise custom analysis pipelines, thus enabling researchers to perform basic as well as highly specific analyses within a single application.

CONCLUSIONS

With MGX, we provide a novel metagenome analysis platform giving researchers access to the most recent analysis tools. MGX covers taxonomic and functional metagenome analysis, statistical evaluation, and a wide range of visualizations easing data interpretation. Its default taxonomic classification pipeline provides equivalent or superior results in comparison to existing tools.

The Resistome of Low-Impacted Marine Environments Is Composed by Distant Metallo-β-Lactamases Homologs

The worldwide dispersion and sudden emergence of new antibiotic resistance genes (ARGs) determined the need in uncovering which environment participate most as their source and reservoir. ARGs closely related to those currently found in human pathogens occur in the resistome of anthropogenic impacted environments. However, the role of pristine environment as the origin and source of ARGs remains underexplored and controversy, particularly, the marine environments represented by the oceans. Here, due to the ocean nature, we hypothesized that the resistome of this pristine/low-impacted marine environment is represented by distant ARG homologs. To test this hypothesis we performed an in silico analysis on the Global Ocean Sampling (GOS) metagenomic project dataset focusing on the metallo-β-lactamases (MβLs) as the ARG model. MβLs have been a challenge to public health, since they hydrolyze the carbapenems, one of the last therapeutic choice in clinics. Using Hidden Markov Model (HMM) profiles, we were successful in identifying a high diversity of distant MβL homologs, related to the B1, B2, and B3 subclasses. The majority of them were distributed across the Atlantic, Indian, and Pacific Oceans being related to the chromosomally encoded MβL GOB present in Elizabethkingia genus. It was observed only a reduced number of metagenomic sequence homologs related to the acquired MβL enzymes (VIM, SPM-1, and AIM-1) that currently have impact in clinics. Therefore, low antibiotic impacted marine environment, as the ocean, are unlikely the source of ARGs that have been causing enormous threat to the public health.

Multispecies reconstructions uncover widespread conservation, and lineage-specific elaborations in eukaryotic mRNA metabolism

The degree of conservation and evolution of cytoplasmic mRNA metabolism pathways across the eukaryotes remains incompletely resolved. In this study, we describe a comprehensive genome and transcriptome-wide analysis of proteins involved in mRNA maturation, translation, and mRNA decay across representative organisms from the six eukaryotic super-groups. We demonstrate that eukaryotes share common pathways for mRNA metabolism that were almost certainly present in the last eukaryotic common ancestor, and show for the first time a correlation between intron density and a selective absence of some Exon Junction Complex (EJC) components in eukaryotes. In addition, we identify pathways that have diversified in individual lineages, with a specific focus on the unique gene gains and losses in members of the Excavata and SAR groups that contribute to their unique gene expression pathways compared to other organisms.

Bacterial SBP56 identified as a Cu-dependent methanethiol oxidase widely distributed in the biosphere

Oxidation of methanethiol (MT) is a significant step in the sulfur cycle. MT is an intermediate of metabolism of globally significant organosulfur compounds including dimethylsulfoniopropionate (DMSP) and dimethylsulfide (DMS), which have key roles in marine carbon and sulfur cycling. In aerobic bacteria, MT is degraded by a MT oxidase (MTO). The enzymatic and genetic basis of MT oxidation have remained poorly characterized. Here, we identify for the first time the MTO enzyme and its encoding gene (mtoX) in the DMS-degrading bacterium Hyphomicrobium sp. VS. We show that MTO is a homotetrameric metalloenzyme that requires Cu for enzyme activity. MTO is predicted to be a soluble periplasmic enzyme and a member of a distinct clade of the Selenium-binding protein (SBP56) family for which no function has been reported. Genes orthologous to mtoX exist in many bacteria able to degrade DMS, other one-carbon compounds or DMSP, notably in the marine model organism Ruegeria pomeroyi DSS-3, a member of the Rhodobacteraceae family that is abundant in marine environments. Marker exchange mutagenesis of mtoX disrupted the ability of R. pomeroyi to metabolize MT confirming its function in this DMSP-degrading bacterium. In R. pomeroyi, transcription of mtoX was enhanced by DMSP, methylmercaptopropionate and MT. Rates of MT degradation increased after pre-incubation of the wild-type strain with MT. The detection of mtoX orthologs in diverse bacteria, environmental samples and its abundance in a range of metagenomic data sets point to this enzyme being widely distributed in the environment and having a key role in global sulfur cycling.

Novel archaeal thermostable cellulases from an oil reservoir metagenome

Microbial assemblages were sampled from an offshore deep sub-surface petroleum reservoir 2.5 km below the ocean floor off the coast of Norway, providing conditions of high temperature and pressure, to identify new thermostable enzymes. In this study, we used DNA sequences obtained directly from the sample metagenome and from a derived fosmid library to survey the functional diversity of this extreme habitat. The metagenomic fosmid library containing 11,520 clones was screened using function- and sequence-based methods to identify recombinant clones expressing carbohydrate-degrading enzymes. Open reading frames (ORFs) encoding carbohydrate-degrading enzymes were predicted by BLAST against the CAZy database, and many fosmid clones expressing carbohydrate-degrading activities were discovered by functional screening using Escherichia coli as a heterologous host. Each complete ORF predicted to encode a cellulase identified from sequence- or function-based screening was subcloned in an expression vector. Five subclones was found to have significant activity using a fluorescent cellulose model substrate, and three of these were observed to be highly thermostable. Based on phylogenetic analyses, the thermostable cellulases were derived from thermophilic Archaea and are distinct from known cellulases. Cellulase F1, obtained from function-based screening, contains two distinct cellulase modules, perhaps resulting from fusion of two archaeal cellulases and with a novel protein structure that may result in enhanced activity and thermostability. This enzyme was found to exhibit exocellulase function and to have a remarkably high activity compared to commercially available enzymes. Results from this study highlight the complementarity of hybrid approaches for enzyme discovery, combining sequence- and function-based screening.

Characterisation of the canine faecal virome in healthy dogs and dogs with acute diarrhoea using shotgun metagenomics

The virome has been increasingly investigated in numerous animal species and in different sites of the body, facilitating the identification and discovery of a variety of viruses. In spite of this, the faecal virome of healthy dogs has not been investigated. In this study we describe the faecal virome of healthy dogs and dogs with acute diarrhoea in Australia, using a shotgun metagenomic approach. Viral sequences from a range of different virus families, including both RNA and DNA families, and known pathogens implicated in enteric disease were documented. Twelve viral families were identified, of which four were bacteriophages. Eight eukaryotic viral families were detected: Astroviridae, Coronaviridae, Reoviridae, Picornaviridae, Caliciviridae, Parvoviridae, Adenoviridae and Papillomaviridae. Families Astroviridae, Picornaviridae and Caliciviridae were found only in dogs with acute diarrhoea, with Astroviridae being the most common family identified in this group. Due to its prevalence, characterisation the complete genome of a canine astrovirus was performed. These studies indicate that metagenomic analyses are useful for the investigation of viral populations in the faeces of dogs. Further studies to elucidate the epidemiological and biological relevance of these findings are warranted.

Metabolomic and Metagenomic Analysis of Two Crude Oil Production Pipelines Experiencing Differential Rates of Corrosion

Corrosion processes in two North Sea oil production pipelines were studied by analyzing pig envelope samples via metagenomic and metabolomic techniques. Both production systems have similar physico-chemical properties and injection waters are treated with nitrate, but one pipeline experiences severe corrosion and the other does not. Early and late pigging material was collected to gain insight into the potential causes for differential corrosion rates. Metabolites were extracted and analyzed via ultra-high performance liquid chromatography/high-resolution mass spectrometry with electrospray ionization (ESI) in both positive and negative ion modes. Metabolites were analyzed by comparison with standards indicative of aerobic and anaerobic hydrocarbon metabolism and by comparison to predicted masses for KEGG metabolites. Microbial community structure was analyzed via 16S rRNA gene qPCR, sequencing of 16S PCR products, and MySeq Illumina shotgun sequencing of community DNA. Metagenomic data were used to reconstruct the full length 16S rRNA genes and genomes of dominant microorganisms. Sequence data were also interrogated via KEGG annotation and for the presence of genes related to terminal electron accepting (TEA) processes as well as aerobic and anaerobic hydrocarbon degradation. Significant and distinct differences were observed when comparing the ‘high corrosion’ (HC) and the ‘low corrosion’ (LC) pipeline systems, especially with respect to the TEA utilization potential. The HC samples were dominated by sulfate-reducing bacteria (SRB) and archaea known for their ability to utilize simple carbon substrates, whereas LC samples were dominated by pseudomonads with the genetic potential for denitrification and aerobic hydrocarbon degradation. The frequency of aerobic hydrocarbon degradation genes was low in the HC system, and anaerobic hydrocarbon degradation genes were not detected in either pipeline. This is in contrast with metabolite analysis, which demonstrated the presence of several succinic acids in HC samples that are diagnostic of anaerobic hydrocarbon metabolism. Identifiable aerobic metabolites were confined to the LC samples, consistent with the metagenomic data. Overall, these data suggest that corrosion management might benefit from a more refined understanding of microbial community resilience in the face of disturbances such as nitrate treatment or pigging, which frequently prove insufficient to alter community structure toward a stable, less-corrosive assemblage.

Viral Metagenomic Analysis Displays the Co-Infection Situation in Healthy and PMWS Affected Pigs

The development of high-throughput sequencing technologies have allowed the possibility to investigate and characterise the entire microbiome of individuals, providing better insight to the complex interaction between different microorganisms. This will help to understand how the microbiome influence the susceptibility of secondary agents and development of disease. We have applied viral metagenomics to investigate the virome of lymph nodes from Swedish pigs suffering from the multifactorial disease postweaning multisystemic wasting syndrome (PMWS) as well as from healthy pigs. The aim is to increase knowledge of potential viruses, apart from porcine circovirus type 2 (PCV2), involved in PMWS development as well as to increase knowledge on the virome of healthy individuals. In healthy individuals, a diverse viral flora was seen with several different viruses present simultaneously. The majority of the identified viruses were small linear and circular DNA viruses, such as different circoviruses, anelloviruses and bocaviruses. In the pigs suffering from PMWS, PCV2 sequences were, as expected, detected to a high extent but other viruses were also identified in the background of PCV2. Apart from DNA viruses also RNA viruses were identified, among them were a porcine pestivirus showing high similarity to a recently (in 2015) discovered atypical porcine pestivirus in the US. Majority of the viruses identified in the background of PCV2 in PMWS pigs could also be identified in the healthy pigs. PCV2 sequences were also identified in the healthy pigs but to a much lower extent than in PMWS affected pigs. Although the method used here is not quantitative the very clear difference in amount of PCV2 sequences in PMWS affected pigs and healthy pigs most likely reflect the very strong replication of PCV2 known to be a hallmark of PMWS. Taken together, these findings illustrate that pigs appear to have a considerable viral flora consisting to a large extent of small single-stranded and circular DNA viruses. Future research on these types of viruses will help to better understand the role that these ubiquitous viruses may have on health and disease of pigs. We also demonstrate for the first time, in Europe, the presence of a novel porcine pestivirus.

Distribution and Classification of Serine β-Lactamases in Brazilian Hospital Sewage and Other Environmental Metagenomes Deposited in Public Databases

β-lactam is the most used antibiotic class in the clinical area and it acts on blocking the bacteria cell wall synthesis, causing cell death. However, some bacteria have evolved resistance to these antibiotics mainly due the production of enzymes known as β-lactamases. Hospital sewage is an important source of dispersion of multidrug-resistant bacteria in rivers and oceans. In this work, we used next-generation DNA sequencing to explore the diversity and dissemination of serine β-lactamases in two hospital sewage from Rio de Janeiro, Brazil (South Zone, SZ and North Zone, NZ), presenting different profiles, and to compare them with public environmental data available. Also, we propose a Hidden-Markov-Model approach to screen potential serine β-lactamases genes (in public environments samples and generated hospital sewage data), exploring its evolutionary relationships. Due to the high variability in β-lactamases, we used a position-specific scoring matrix search method (RPS-BLAST) against conserved domain database profiles (CDD, Pfam, and COG) followed by visual inspection to detect conserved motifs, to increase the reliability of the results and remove possible false positives. We were able to identify novel β-lactamases from Brazilian hospital sewage and to estimate relative abundance of its types. The highest relative abundance found in SZ was the Class A (50%), while Class D is predominant in NZ (55%). CfxA (65%) and ACC (47%) types were the most abundant genes detected in SZ, while in NZ the most frequent were OXA-10 (32%), CfxA (28%), ACC (21%), CEPA (20%), and FOX (19%). Phylogenetic analysis revealed β-lactamases from Brazilian hospital sewage grouped in the same clade and close to sequences belonging to Firmicutes and Bacteroidetes groups, but distant from potential β-lactamases screened from public environmental data, that grouped closer to β-lactamases of Proteobacteria. Our results demonstrated that HMM-based approach identified homologs of serine β-lactamases, indicating the specificity and high sensitivity of this approach in large datasets, contributing for the identification and classification of a large number of homologous genes, comprising possible new ones. Phylogenetic analysis revealed the potential reservoir of β-lactam resistance genes in the environment, contributing to understanding the evolution and dissemination of these genes.

A metagenomic analysis displays the diverse microbial community of a vermicomposting system in Uganda

Background

Vermicomposting is a mesophilic process using earthworms to efficiently and at low cost process large volumes of organic waste. It has been suggested to not only increase soil fertility but also increase biomass of beneficial bacteria while reducing harmful bacteria. The aim of this study was to set up a strategy to investigate and characterise the viral as well as the bacterial composition of a vermicomposting system.

Material and methods

The vermicomposting unit used in this study was placed at the Makerere University Agricultural Research Institute Kabanyolo on the outskirts of Kampala, Uganda, and was fed with 80% cattle manure and 20% food waste. On Day 172, the compost was terminated and compost samples were collected from three layers of the unit: the top, the middle and the bottom layer. A metagenomic approach was then applied to characterise the viral and bacterial composition of the vermicomposting system.

Results and discussion

A high abundance and diversity of bacteria were identified. Proteobacteria was the largest phyla in the compost (mainly Alpha-, Gamma- and Betaproteobacteria), constituting almost 65% of the bacterial reads in the data sets. DNA samples from several possible pathogenic bacteria, such as Salmonella spp., Escherichia coli, Enterobacter spp., Enterococcus spp. and Clostridium spp, were detected in the vermicompost, suggesting that there might still be harmful bacteria in the vermicast. Phages constituted the main viral group; apart from phages, mainly insect viruses were identified. The only animal or human virus identified was kobuvirus. In summary, metagenomic analysis was shown to be an efficient technology to characterise the microbial composition of vermicast. The data from this study contribute to a better understanding of the microbes present in this kind of composting system and can help determine measures necessary for safe manure handling.

Bacterial Human Virulence Genes across Diverse Habitats As Assessed by In silico Analysis of Environmental Metagenomes

The occurrence and distribution of clinically relevant bacterial virulence genes across natural (non-human) environments is not well understood. We aimed to investigate the occurrence of homologs to bacterial human virulence genes in a variety of ecological niches to better understand the role of natural environments in the evolution of bacterial virulence. Twenty four bacterial virulence genes were analyzed in 46 diverse environmental metagenomic datasets, representing various soils, seawater, freshwater, marine sediments, hot springs, the deep-sea, hypersaline mats, microbialites, gutless worms and glacial ice. Homologs to 16 bacterial human virulence genes, involved in urinary tract infections, gastrointestinal diseases, skin diseases, and wound and systemic infections, showed global ubiquity. A principal component analysis did not demonstrate clear trends across the metagenomes with respect to occurrence and frequency of observed gene homologs. Full-length (>95%) homologs of several virulence genes were identified, and translated sequences of the environmental and clinical genes were up to 50-100% identical. Furthermore, phylogenetic analyses indicated deep branching positions of some of the environmental gene homologs, suggesting that they represent ancient lineages in the phylogeny of the clinical genes. Fifteen virulence gene homologs were detected in metatranscriptomes, providing evidence of environmental expression. The ubiquitous presence and transcription of the virulence gene homologs in non-human environments point to an important ecological role of the genes for the activity and survival of environmental bacteria. Furthermore, the high degree of sequence conservation between several of the environmental and clinical genes suggests common ancestral origins.

Comparative genomics of Paracoccus sp. SM22M-07 isolated from coral mucus: insights into bacteria-host interactions

One of the main goals of coral microbiology is to understand the ways in which coral-bacteria associations are established and maintained. This work describes the sequencing of the genome of Paracoccus sp. SM22M-07 isolated from the mucus of the endemic Brazilian coral species Mussismilia hispida. Comparative analysis was used to identify unique genomic features of SM22M-07 that might be involved in its adaptation to the marine ecosystem and the nutrient-rich environment provided by coral mucus, as well as in the establishment and strengthening of the interaction with the host. These features included genes related to the type IV protein secretion system, erythritol catabolism, and succinoglycan biosynthesis. We experimentally confirmed the production of succinoglycan by Paracoccus sp. SM22M-07 and we hypothesize that it may be involved in the association of the bacterium with coral surfaces.

Direct 16S rRNA-seq from bacterial communities: a PCR-independent approach to simultaneously assess microbial diversity and functional activity potential of each taxon

The analysis of environmental microbial communities has largely relied on a PCR-dependent amplification of genes entailing species identity as 16S rRNA. This approach is susceptible to biases depending on the level of primer matching in different species. Moreover, possible yet-to-discover taxa whose rRNA could differ enough from known ones would not be revealed. DNA-based methods moreover do not provide information on the actual physiological relevance of each taxon within an environment and are affected by the variable number of rRNA operons in different genomes. To overcome these drawbacks we propose an approach of direct sequencing of 16S ribosomal RNA without any primer- or PCR-dependent step. The method was tested on a microbial community developing in an anammox bioreactor sampled at different time-points. A conventional PCR-based amplicon pyrosequencing was run in parallel. The community resulting from direct rRNA sequencing was highly consistent with the known biochemical processes operative in the reactor. As direct rRNA-seq is based not only on taxon abundance but also on physiological activity, no comparison between its results and those from PCR-based approaches can be applied. The novel principle is in this respect proposed not as an alternative but rather as a complementary methodology in microbial community studies.

HoloVir: A Workflow for Investigating the Diversity and Function of Viruses in Invertebrate Holobionts

Abundant bioinformatics resources are available for the study of complex microbial metagenomes, however their utility in viral metagenomics is limited. HoloVir is a robust and flexible data analysis pipeline that provides an optimized and validated workflow for taxonomic and functional characterization of viral metagenomes derived from invertebrate holobionts. Simulated viral metagenomes comprising varying levels of viral diversity and abundance were used to determine the optimal assembly and gene prediction strategy, and multiple sequence assembly methods and gene prediction tools were tested in order to optimize our analysis workflow. HoloVir performs pairwise comparisons of single read and predicted gene datasets against the viral RefSeq database to assign taxonomy and additional comparison to phage-specific and cellular markers is undertaken to support the taxonomic assignments and identify potential cellular contamination. Broad functional classification of the predicted genes is provided by assignment of COG microbial functional category classifications using EggNOG and higher resolution functional analysis is achieved by searching for enrichment of specific Swiss-Prot keywords within the viral metagenome. Application of HoloVir to viral metagenomes from the coral Pocillopora damicornis and the sponge Rhopaloeides odorabile demonstrated that HoloVir provides a valuable tool to characterize holobiont viral communities across species, environments, or experiments.

"Every Gene Is Everywhere but the Environment Selects": Global Geolocalization of Gene Sharing in Environmental Samples through Network Analysis

The spatial distribution of microbes on our planet is famously formulated in the Baas Becking hypothesis as “everything is everywhere but the environment selects.” While this hypothesis does not strictly rule out patterns caused by geographical effects on ecology and historical founder effects, it does propose that the remarkable dispersal potential of microbes leads to distributions generally shaped by environmental factors rather than geographical distance. By constructing sequence similarity networks from uncultured environmental samples, we show that microbial gene pool distributions are not influenced nearly as much by geography as ecology, thus extending the Bass Becking hypothesis from whole organisms to microbial genes. We find that gene pools are shaped by their broad ecological niche (such as sea water, fresh water, host, and airborne). We find that freshwater habitats act as a gene exchange bridge between otherwise disconnected habitats. Finally, certain antibiotic resistance genes deviate from the general trend of habitat specificity by exhibiting a high degree of cross-habitat mobility. The strong cross-habitat mobility of antibiotic resistance genes is a cause for concern and provides a paradigmatic example of the rate by which genes colonize new habitats when new selective forces emerge.

New biosynthetic pathway for pink pigments from uncultured oceanic viruses

The pink open-chain tetrapyrrole pigment phycoerythrobilin (PEB) is employed by marine cyanobacteria, red algae and cryptophytes as a light-harvesting chromophore in phycobiliproteins. Genes encoding biosynthesis proteins for PEB have also been discovered in cyanophages, viruses that infect cyanobacteria, and mimic host pigment biosynthesis with the exception of PebS which combines the enzymatic activities of two host enzymes. In this study, we have identified novel members of the PEB biosynthetic enzyme families, heme oxygenases and ferredoxin-dependent bilin reductases. Encoding genes were found in metagenomic datasets and could be traced back to bacteriophage but not cyanophage origin. While the heme oxygenase exhibited standard activity, a new bilin reductase with highest homology to the teal pigment producing enzyme PcyA revealed PEB biosynthetic activity. Although PcyX possesses PebS-like activity both enzymes share only 9% sequence identity and likely catalyze the reaction via two independent mechanisms. Our data point towards the presence of phycobilin biosynthetic genes in phages that probably infect alphaproteobacteria and, therefore, further support a role of phycobilins outside oxygenic phototrophs.

Discovery Strategies of Bioactive Compounds Synthesized by Nonribosomal Peptide Synthetases and Type-I Polyketide Synthases Derived from Marine Microbiomes

Considering that 70% of our planet's surface is covered by oceans, it is likely that undiscovered biodiversity is still enormous. A large portion of marine biodiversity consists of microbiomes. They are very attractive targets of bioprospecting because they are able to produce a vast repertoire of secondary metabolites in order to adapt in diverse environments. In many cases secondary metabolites of pharmaceutical and biotechnological interest such as nonribosomal peptides (NRPs) and polyketides (PKs) are synthesized by multimodular enzymes named nonribosomal peptide synthetases (NRPSes) and type-I polyketide synthases (PKSes-I), respectively. Novel findings regarding the mechanisms underlying NRPS and PKS evolution demonstrate how microorganisms could leverage their metabolic potential. Moreover, these findings could facilitate synthetic biology approaches leading to novel bioactive compounds. Ongoing advances in bioinformatics and next-generation sequencing (NGS) technologies are driving the discovery of NRPs and PKs derived from marine microbiomes mainly through two strategies: genome-mining and metagenomics. Microbial genomes are now sequenced at an unprecedented rate and this vast quantity of biological information can be analyzed through genome mining in order to identify gene clusters encoding NRPSes and PKSes of interest. On the other hand, metagenomics is a fast-growing research field which directly studies microbial genomes and their products present in marine environments using culture-independent approaches. The aim of this review is to examine recent developments regarding discovery strategies of bioactive compounds synthesized by NRPS and type-I PKS derived from marine microbiomes and to highlight the vast diversity of NRPSes and PKSes present in marine environments by giving examples of recently discovered bioactive compounds.

The yjdF riboswitch candidate regulates gene expression by binding diverse azaaromatic compounds

The yjdF motif RNA is an orphan riboswitch candidate that almost exclusively associates with the yjdF protein-coding gene in many bacteria. The function of the YjdF protein is unknown, which has made speculation regarding the natural ligand for this putative riboswitch unusually challenging. By using a structure-probing assay for ligand binding, we found that a surprisingly broad diversity of nitrogen-containing aromatic heterocycles, or "azaaromatics," trigger near-identical changes in the structures adopted by representative yjdF motif RNAs. Regions of the RNA that undergo ligand-induced structural modulation reside primarily in portions of the putative aptamer region that are highly conserved in nucleotide sequence, as is typical for riboswitches. Some azaaromatic molecules are bound by the RNA with nanomolar dissociation constants, and a subset of these ligands activate riboswitch-mediated gene expression in cells. Furthermore, genetic elements most commonly adjacent to the yjdF motif RNA or to the yjdF protein-coding region are homologous to protein regulators implicated in mitigating the toxic effects of diverse phenolic acids or polycyclic compounds. Although the precise type of natural ligand sensed by yjdF motif RNAs remains unknown, our findings suggest that this riboswitch class might serve as part of a genetic response system to toxic or signaling compounds with chemical structures similar to azaaromatics.

Deciphering ocean carbon in a changing world

Dissolved organic matter (DOM) in the oceans is one of the largest pools of reduced carbon on Earth, comparable in size to the atmospheric CO2 reservoir. A vast number of compounds are present in DOM, and they play important roles in all major element cycles, contribute to the storage of atmospheric CO2 in the ocean, support marine ecosystems, and facilitate interactions between organisms. At the heart of the DOM cycle lie molecular-level relationships between the individual compounds in DOM and the members of the ocean microbiome that produce and consume them. In the past, these connections have eluded clear definition because of the sheer numerical complexity of both DOM molecules and microorganisms. Emerging tools in analytical chemistry, microbiology, and informatics are breaking down the barriers to a fuller appreciation of these connections. Here we highlight questions being addressed using recent methodological and technological developments in those fields and consider how these advances are transforming our understanding of some of the most important reactions of the marine carbon cycle.

Analysis of metagenomic data reveals common features of halophilic viral communities across continents

Microbial communities from hypersaline ponds, dominated by halophilic archaea, are considered specific of such extreme conditions. The associated viral communities have accordingly been shown to display specific features, such as similar morphologies among different sites. However, little is known about the genetic diversity of these halophilic viral communities across the Earth. Here, we studied viral communities in hypersaline ponds sampled on the coast of Senegal (8-36% of salinity) using metagenomics approach, and compared them with hypersaline viromes from Australia and Spain. The specificity of hyperhalophilic viruses could first be demonstrated at a community scale, salinity being a strong discriminating factor between communities. For the major viral group detected in all samples (Caudovirales), only a limited number of halophilic Caudovirales clades were highlighted. These clades gather viruses from different continents and display consistent genetic composition, indicating that they represent related lineages with a worldwide distribution. Non-tailed hyperhalophilic viruses display a greater rate of gene transfer and recombination, with uncharacterized genes conserved across different kind of viruses and plasmids. Thus, hypersaline viral communities around the world appear to form a genetically consistent community that are likely to harbour new genes coding for enzymes specifically adapted to these environments.

Web Resources for Metagenomics Studies

The development of next-generation sequencing (NGS) platforms spawned an enormous volume of data. This explosion in data has unearthed new scalability challenges for existing bioinformatics tools. The analysis of metagenomic sequences using bioinformatics pipelines is complicated by the substantial complexity of these data. In this article, we review several commonly-used online tools for metagenomics data analysis with respect to their quality and detail of analysis using simulated metagenomics data. There are at least a dozen such software tools presently available in the public domain. Among them, MGRAST, IMG/M, and METAVIR are the most well-known tools according to the number of citations by peer-reviewed scientific media up to mid-2015. Here, we describe 12 online tools with respect to their web link, annotation pipelines, clustering methods, online user support, and availability of data storage. We have also done the rating for each tool to screen more potential and preferential tools and evaluated five best tools using synthetic metagenome. The article comprehensively deals with the contemporary problems and the prospects of metagenomics from a bioinformatics viewpoint.

Metagenomics: Retrospect and Prospects in High Throughput Age

In recent years, metagenomics has emerged as a powerful tool for mining of hidden microbial treasure in a culture independent manner. In the last two decades, metagenomics has been applied extensively to exploit concealed potential of microbial communities from almost all sorts of habitats. A brief historic progress made over the period is discussed in terms of origin of metagenomics to its current state and also the discovery of novel biological functions of commercial importance from metagenomes of diverse habitats. The present review also highlights the paradigm shift of metagenomics from basic study of community composition to insight into the microbial community dynamics for harnessing the full potential of uncultured microbes with more emphasis on the implication of breakthrough developments, namely, Next Generation Sequencing, advanced bioinformatics tools, and systems biology.

A novel group of diverse Polinton-like viruses discovered by metagenome analysis

Background

The rapidly growing metagenomic databases provide increasing opportunities for computational discovery of new groups of organisms. Identification of new viruses is particularly straightforward given the comparatively small size of viral genomes, although fast evolution of viruses complicates the analysis of novel sequences. Here we report the metagenomic discovery of a distinct group of diverse viruses that are distantly related to the eukaryotic virus-like transposons of the Polinton superfamily.

Results

The sequence of the putative major capsid protein (MCP) of the unusual linear virophage associated with Phaeocystis globosa virus (PgVV) was used as a bait to identify potential related viruses in metagenomic databases. Assembly of the contigs encoding the PgVV MCP homologs followed by comprehensive sequence analysis of the proteins encoded in these contigs resulted in the identification of a large group of Polinton-like viruses (PLV) that resemble Polintons (polintoviruses) and virophages in genome size, and share with them a conserved minimal morphogenetic module that consists of major and minor capsid proteins and the packaging ATPase. With a single exception, the PLV lack the retrovirus-type integrase that is encoded in the genomes of all Polintons and the Mavirus group of virophages. However, some PLV encode a newly identified tyrosine recombinase-integrase that is common in bacteria and bacteriophages and is also found in the Organic Lake virophage group. Although several PLV genomes and individual genes are integrated into algal genomes, it appears likely that most of the PLV are viruses. Given the absence of protease and retrovirus-type integrase, the PLV could resemble the ancestral polintoviruses that evolved from bacterial tectiviruses. Apart from the conserved minimal morphogenetic module, the PLV widely differ in their genome complements but share a gene network with Polintons and virophages, suggestive of multiple gene exchanges within a shared gene pool.

Conclusions

The discovery of PLV substantially expands the emerging class of eukaryotic viruses and transposons that also includes Polintons and virophages. This class of selfish elements is extremely widespread and might have been a hotbed of eukaryotic virus, transposon and plasmid evolution. New families of these elements are expected to be discovered.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0207-4) contains supplementary material, which is available to authorized users.

COGNIZER: A Framework for Functional Annotation of Metagenomic Datasets

BACKGROUND

Recent advances in sequencing technologies have resulted in an unprecedented increase in the number of metagenomes that are being sequenced world-wide. Given their volume, functional annotation of metagenomic sequence datasets requires specialized computational tools/techniques. In spite of having high accuracy, existing stand-alone functional annotation tools necessitate end-users to perform compute-intensive homology searches of metagenomic datasets against "multiple" databases prior to functional analysis. Although, web-based functional annotation servers address to some extent the problem of availability of compute resources, uploading and analyzing huge volumes of sequence data on a shared public web-service has its own set of limitations. In this study, we present COGNIZER, a comprehensive stand-alone annotation framework which enables end-users to functionally annotate sequences constituting metagenomic datasets. The COGNIZER framework provides multiple workflow options. A subset of these options employs a novel directed-search strategy which helps in reducing the overall compute requirements for end-users. The COGNIZER framework includes a cross-mapping database that enables end-users to simultaneously derive/infer KEGG, Pfam, GO, and SEED subsystem information from the COG annotations.

RESULTS

Validation experiments performed with real-world metagenomes and metatranscriptomes, generated using diverse sequencing technologies, indicate that the novel directed-search strategy employed in COGNIZER helps in reducing the compute requirements without significant loss in annotation accuracy. A comparison of COGNIZER's results with pre-computed benchmark values indicate the reliability of the cross-mapping database employed in COGNIZER.

CONCLUSION

The COGNIZER framework is capable of comprehensively annotating any metagenomic or metatranscriptomic dataset from varied sequencing platforms in functional terms. Multiple search options in COGNIZER provide end-users the flexibility of choosing a homology search protocol based on available compute resources. The cross-mapping database in COGNIZER is of high utility since it enables end-users to directly infer/derive KEGG, Pfam, GO, and SEED subsystem annotations from COG categorizations. Furthermore, availability of COGNIZER as a stand-alone scalable implementation is expected to make it a valuable annotation tool in the field of metagenomic research.

AVAILABILITY AND IMPLEMENTATION

A Linux implementation of COGNIZER is freely available for download from the following links: http://metagenomics.atc.tcs.com/cognizer, https://metagenomics.atc.tcs.com/function/cognizer.

Evaluation of shotgun metagenomics sequence classification methods using in silico and in vitro simulated communities

Background

The field of metagenomics (study of genetic material recovered directly from an environment) has grown rapidly, with many bioinformatics analysis methods being developed. To ensure appropriate use of such methods, robust comparative evaluation of their accuracy and features is needed. For taxonomic classification of sequence reads, such evaluation should include use of clade exclusion, which better evaluates a method’s accuracy when identical sequences are not present in any reference database, as is common in metagenomic analysis. To date, relatively small evaluations have been performed, with evaluation approaches like clade exclusion limited to assessment of new methods by the authors of the given method. What is needed is a rigorous, independent comparison between multiple major methods, using the same in silico and in vitro test datasets, with and without approaches like clade exclusion, to better characterize accuracy under different conditions.

Results

An overview of the features of 38 bioinformatics methods is provided, evaluating accuracy with a focus on 11 programs that have reference databases that can be modified and therefore most robustly evaluated with clade exclusion. Taxonomic classification of sequence reads was evaluated using both in silico and in vitro mock bacterial communities. Clade exclusion was used at taxonomic levels from species to class—identifying how well methods perform in progressively more difficult scenarios. A wide range of variability was found in the sensitivity, precision, overall accuracy, and computational demand for the programs evaluated. In experiments where distilled water was spiked with only 11 bacterial species, frequently dozens to hundreds of species were falsely predicted by the most popular programs. The different features of each method (forces predictions or not, etc.) are summarized, and additional analysis considerations discussed.

Conclusions

The accuracy of shotgun metagenomics classification methods varies widely. No one program clearly outperformed others in all evaluation scenarios; rather, the results illustrate the strengths of different methods for different purposes. Researchers must appreciate method differences, choosing the program best suited for their particular analysis to avoid very misleading results. Use of standardized datasets for method comparisons is encouraged, as is use of mock microbial community controls suitable for a particular metagenomic analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0788-5) contains supplementary material, which is available to authorized users.

The role of biofilms as environmental reservoirs of antibiotic resistance

Antibiotic resistance has become a significant and growing threat to public and environmental health. To face this problem both at local and global scales, a better understanding of the sources and mechanisms that contribute to the emergence and spread of antibiotic resistance is required. Recent studies demonstrate that aquatic ecosystems are reservoirs of resistant bacteria and antibiotic resistance genes as well as potential conduits for their transmission to human pathogens. Despite the wealth of information about antibiotic pollution and its effect on the aquatic microbial resistome, the contribution of environmental biofilms to the acquisition and spread of antibiotic resistance has not been fully explored in aquatic systems. Biofilms are structured multicellular communities embedded in a self-produced extracellular matrix that acts as a barrier to antibiotic diffusion. High population densities and proximity of cells in biofilms also increases the chances for genetic exchange among bacterial species converting biofilms in hot spots of antibiotic resistance. This review focuses on the potential effect of antibiotic pollution on biofilm microbial communities, with special emphasis on ecological and evolutionary processes underlying acquired resistance to these compounds.

Big Data Provenance: Challenges, State of the Art and Opportunities

Ability to track provenance is a key feature of scientific workflows to support data lineage and reproducibility. The challenges that are introduced by the volume, variety and velocity of Big Data, also pose related challenges for provenance and quality of Big Data, defined as veracity. The increasing size and variety of distributed Big Data provenance information bring new technical challenges and opportunities throughout the provenance lifecycle including recording, querying, sharing and utilization. This paper discusses the challenges and opportunities of Big Data provenance related to the veracity of the datasets themselves and the provenance of the analytical processes that analyze these datasets. It also explains our current efforts towards tracking and utilizing Big Data provenance using workflows as a programming model to analyze Big Data.

Provirophages in the Bigelowiella genome bear testimony to past encounters with giant viruses

Virophages are recently discovered double-stranded DNA virus satellites that prey on giant viruses (nucleocytoplasmic large DNA viruses; NCLDVs), which are themselves parasites of unicellular eukaryotes. This coupled parasitism can result in the indirect control of eukaryotic cell mortality by virophages. However, the details of such tripartite relationships remain largely unexplored. We have discovered ∼300 predicted genes of putative virophage origin in the nuclear genome of the unicellular alga Bigelowiella natans. Physical clustering of these genes indicates that virophage genomes are integrated into the B. natans genome. Virophage inserts show high levels of similarity and synteny between each other, indicating that they are closely related. Virophage genes are transcribed not only in the sequenced B. natans strain but also in other Bigelowiella isolates, suggesting that transcriptionally active virophage inserts are widespread in Bigelowiella populations. Evidence that B. natans is also a host to NCLDV members is provided by the identification of NCLDV inserts in its genome. These putative large DNA viruses may be infected by B. natans virophages. We also identify four repeated elements sharing structural and genetic similarities with transpovirons--a class of mobile elements first discovered in giant viruses--that were probably independently inserted in the B. natans genome. We argue that endogenized provirophages may be beneficial to both the virophage and B. natans by (i) increasing the chances for the virophage to coinfect the host cell with an NCLDV prey and (ii) defending the host cell against fatal NCLDV infections.

Microbial Consortium Associated with the Antarctic Marine Ciliate Euplotes focardii: An Investigation from Genomic Sequences

We report the characterization of the bacterial consortium associated to Euplotes focardii, a strictly psychrophilic marine ciliate that was maintained in laboratory cultures at 4 °C after its first isolation from Terra Nova Bay, in Antarctica. By Illumina genome analyser, we obtained 11,179 contigs of potential prokaryotic origin and classified them according to the NCBI's prokaryotic attributes table. The majority of these sequences correspond to either Bacteroidetes (16 %) or Proteobacteria (78 %). The latter were dominated by gamma- (39 %, including sequences related to the pathogenic genus Francisella), and alpha-proteobacterial (30 %) sequences. Analysis of the Pfam domain family and Gene Ontology term variation revealed that the most frequent terms that appear unique to this consortium correspond to proteins involved in "transmembrane transporter activity" and "oxidoreductase activity". Furthermore, we identified genes that encode for enzymes involved in the catabolism of complex substance for energy reserves. We also characterized members of the transposase and integrase superfamilies, whose role in bacterial evolution is well documented, as well as putative antifreeze proteins. Antibiotic treatments of E. focardii cultures delayed the cell division of the ciliate. To conclude, our results indicate that this consortium is largely represented by bacteria derived from the original Antarctic sample and may contribute to the survival of E. focardii in laboratory condition. Furthermore, our results suggest that these bacteria may have a more general role in E. focardii survival in its natural cold and oxidative environment.

New classes of self-cleaving ribozymes revealed by comparative genomics analysis

Enzymes made of RNA catalyze reactions that are essential for protein synthesis and RNA processing. However, such natural ribozymes are exceedingly rare, as evident by the fact that the discovery rate for new classes has dropped to one per decade from about one per year during the 1980s. Indeed, only 11 distinct ribozyme classes have been experimentally validated to date. Recently, we recognized that self-cleaving ribozymes frequently associate with certain types of genes from bacteria. Herein this synteny was exploited to identify divergent architectures for two previously known ribozyme classes and to discover additional noncoding RNA motifs that are self-cleaving RNA candidates. Three new self-cleaving classes, named twister sister, pistol and hatchet, have been identified from this collection, suggesting that even more ribozymes remain hidden in modern cells.