Dataset Information

ENVirT: inference of ecological characteristics of viruses from metagenomic data.

ABSTRACT: BACKGROUND:Estimating the parameters that describe the ecology of viruses,particularly those that are novel, can be made possible using metagenomic approaches. However, the best-performing existing methods require databases to first estimate an average genome length of a viral community before being able to estimate other parameters, such as viral richness. Although this approach has been widely used, it can adversely skew results since the majority of viruses are yet to be catalogued in databases. RESULTS:In this paper, we present ENVirT, a method for estimating the richness of novel viral mixtures, and for the first time we also show that it is possible to simultaneously estimate the average genome length without a priori information. This is shown to be a significant improvement over database-dependent methods, since we can now robustly analyze samples that may include novel viral types under-represented in current databases. We demonstrate that the viral richness estimates produced by ENVirT are several orders of magnitude higher in accuracy than the estimates produced by existing methods named PHACCS and CatchAll when benchmarked against simulated data. We repeated the analysis of 20 metavirome samples using ENVirT, which produced results in close agreement with complementary in virto analyses. CONCLUSIONS:These insights were previously not captured by existing computational methods. As such, ENVirT is shown to be an essential tool for enhancing our understanding of novel viral populations.

SUBMITTER: Jayasundara D

PROVIDER: S-EPMC7394321 | biostudies-literature | 2019 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

ENVirT: inference of ecological characteristics of viruses from metagenomic data.

Jayasundara Duleepa D Herath Damayanthi D Senanayake Damith D Saeed Isaam I Yang Cheng-Yu CY Sun Yuan Y Chang Bill C BC Tang Sen-Lin SL Halgamuge Saman K SK

BMC bioinformatics 20190204 Suppl 13

<h4>Background</h4>Estimating the parameters that describe the ecology of viruses,particularly those that are novel, can be made possible using metagenomic approaches. However, the best-performing existing methods require databases to first estimate an average genome length of a viral community before being able to estimate other parameters, such as viral richness. Although this approach has been widely used, it can adversely skew results since the majority of viruses are yet to be catalogued in ...[more]

PMID: 30717665

Similar Datasets

Project description:Abstract Accelerometers in animal‐attached tags are powerful tools in behavioural ecology, they can be used to determine behaviour and provide proxies for movement‐based energy expenditure. Researchers are collecting and archiving data across systems, seasons and device types. However, using data repositories to draw ecological inference requires a good understanding of the error introduced according to sensor type and position on the study animal and protocols for error assessment and minimisation. Using laboratory trials, we examine the absolute accuracy of tri‐axial accelerometers and determine how inaccuracies impact measurements of dynamic body acceleration (DBA), a proxy for energy expenditure, in human participants. We then examine how tag type and placement affect the acceleration signal in birds, using pigeons Columba livia flying in a wind tunnel, with tags mounted simultaneously in two positions, and back‐ and tail‐mounted tags deployed on wild kittiwakes Rissa tridactyla. Finally, we present a case study where two generations of tag were deployed using different attachment procedures on red‐tailed tropicbirds Phaethon rubricauda foraging in different seasons. Bench tests showed that individual acceleration axes required a two‐level correction to eliminate measurement error. This resulted in DBA differences of up to 5% between calibrated and uncalibrated tags for humans walking at a range of speeds. Device position was associated with greater variation in DBA, with upper and lower back‐mounted tags varying by 9% in pigeons, and tail‐ and back‐mounted tags varying by 13% in kittiwakes. The tropicbird study highlighted the difficulties of attributing changes in signal amplitude to a single factor when confounding influences tend to covary, as DBA varied by 25% between seasons. Accelerometer accuracy, tag placement and attachment critically affect the signal amplitude and thereby the ability of the system to detect biologically meaningful phenomena. We propose a simple method to calibrate accelerometers that can be executed under field conditions. This should be used prior to deployments and archived with resulting data. We also suggest a way that researchers can assess accuracy in previously collected data, and caution that variable tag placement and attachment can increase sensor noise and even generate trends that have no biological meaning.

Project description:BackgroundCharacterizing animal space use is critical for understanding ecological relationships. Animal telemetry technology has revolutionized the fields of ecology and conservation biology by providing high quality spatial data on animal movement. Radio-telemetry with very high frequency (VHF) radio signals continues to be a useful technology because of its low cost, miniaturization, and low battery requirements. Despite a number of statistical developments synthetically integrating animal location estimation and uncertainty with spatial process models using satellite telemetry data, we are unaware of similar developments for azimuthal telemetry data. As such, there are few statistical options to handle these unique data and no synthetic framework for modeling animal location uncertainty and accounting for it in ecological models.We developed a hierarchical modeling framework to provide robust animal location estimates from one or more intersecting or non-intersecting azimuths. We used our azimuthal telemetry model (ATM) to account for azimuthal uncertainty with covariates and propagate location uncertainty into spatial ecological models. We evaluate the ATM with commonly used estimators (Lenth (1981) maximum likelihood and M-Estimators) using simulation. We also provide illustrative empirical examples, demonstrating the impact of ignoring location uncertainty within home range and resource selection analyses. We further use simulation to better understand the relationship among location uncertainty, spatial covariate autocorrelation, and resource selection inference.ResultsWe found the ATM to have good performance in estimating locations and the only model that has appropriate measures of coverage. Ignoring animal location uncertainty when estimating resource selection or home ranges can have pernicious effects on ecological inference. Home range estimates can be overly confident and conservative when ignoring location uncertainty and resource selection coefficients can lead to incorrect inference and over confidence in the magnitude of selection. Furthermore, our simulation study clarified that incorporating location uncertainty helps reduce bias in resource selection coefficients across all levels of covariate spatial autocorrelation.ConclusionThe ATM can accommodate one or more azimuths when estimating animal locations, regardless of how they intersect; this ensures that all data collected are used for ecological inference. Our findings and model development have important implications for interpreting historical analyses using this type of data and the future design of radio-telemetry studies.

Project description:Giant viruses have remarkable genomic repertoires-blurring the line with cellular life-and act as top-down controls of eukaryotic plankton. However, to date only six cultured giant virus genomes are available from the pelagic ocean. We used at-sea flow cytometry with staining and sorting designed to target wild predatory eukaryotes, followed by DNA sequencing and assembly, to recover novel giant viruses from the Pacific Ocean. We retrieved four 'PacV' partial genomes that range from 421 to 1605 Kb, with 13 contigs on average, including the largest marine viral genomic assembly reported to date. Phylogenetic analyses indicate that three of the new viruses span a clade with deep-branching members of giant Mimiviridae, incorporating the Cafeteria roenbergensis virus, the uncultivated terrestrial Faunusvirus, one PacV from a choanoflagellate and two PacV with unclear hosts. The fourth virus, oPacV-421, is phylogenetically related to viruses that infect haptophyte algae. About half the predicted proteins in each PacV have no matches in NCBI nr (e-value < 10-5), totalling 1735 previously unknown proteins; the closest affiliations of the other proteins were evenly distributed across eukaryotes, prokaryotes and viruses of eukaryotes. The PacVs encode many translational proteins and two encode eukaryotic-like proteins from the Rh family of the ammonium transporter superfamily, likely influencing the uptake of nitrogen during infection. cPacV-1605 encodes a microbial viral rhodopsin (VirR) and the biosynthesis pathway for the required chromophore, the second finding of a choanoflagellate-associated virus that encodes these genes. In co-collected metatranscriptomes, 85% of cPacV-1605 genes were expressed, with capsids, heat shock proteins and proteases among the most highly expressed. Based on orthologue presence-absence patterns across the PacVs and other eukaryotic viruses, we posit the observed viral groupings are connected to host lifestyles as heterotrophs or phototrophs. This article is part of a discussion meeting issue 'Single cell ecology'.

Project description:BackgroundMetagenomics enables the analysis of bacterial population composition and the study of emergent population features, such as shared metabolic pathways. Recently, we have shown that metagenomics datasets can be leveraged to characterize population-wide transcriptional regulatory networks, or meta-regulons, providing insights into how bacterial populations respond collectively to specific triggers. Here we formalize a Bayesian inference framework to analyze the composition of transcriptional regulatory networks in metagenomes by determining the probability of regulation of orthologous gene sequences. We assess the performance of this approach on synthetic datasets and we validate it by analyzing the copper-homeostasis network of Firmicutes species in the human gut microbiome.ResultsAssessment on synthetic datasets shows that our method provides a robust and interpretable metric for assessing putative regulation by a transcription factor on sets of promoter sequences mapping to an orthologous gene cluster. The inference framework integrates the regulatory contribution of secondary sites and can discern false positives arising from multiple instances of a clonal sequence. Posterior probabilities for orthologous gene clusters decline sharply when less than 20 % of mapped promoters have binding sites, but we introduce a sensitivity adjustment procedure to speed up computation that enhances regulation assessment in heterogeneous ortholog clusters. Analysis of the copper-homeostasis regulon governed by CsoR in the human gut microbiome Firmicutes reveals that CsoR controls itself and copper-translocating P-type ATPases, but not CopZ-type copper chaperones. Our analysis also indicates that CsoR frequently targets promoters with dual CsoR-binding sites, suggesting that it exploits higher-order binding conformations to fine-tune its activity.ConclusionsWe introduce and validate a method for the analysis of transcriptional regulatory networks from metagenomic data that enables inference of meta-regulons in a systematic and interpretable way. Validation of this method on the CsoR meta-regulon of gut microbiome Firmicutes illustrates the usefulness of the approach, revealing novel properties of the copper-homeostasis network in poorly characterized bacterial species and putting forward evidence of new mechanisms of DNA binding for this transcriptional regulator. Our approach will enable the comparative analysis of regulatory networks across metagenomes, yielding novel insights into the evolution of transcriptional regulatory networks.

Project description:Rapid, sensitive, and specific virus detection is an important component of clinical diagnostics. Massively parallel sequencing enables new diagnostic opportunities that complement traditional serological and PCR based techniques. While massively parallel sequencing promises the benefits of being more comprehensive and less biased than traditional approaches, it presents new analytical challenges, especially with respect to detection of pathogen sequences in metagenomic contexts. To a first approximation, the initial detection of viruses can be achieved simply through alignment of sequence reads or assembled contigs to a reference database of pathogen genomes with tools such as BLAST. However, recognition of highly divergent viral sequences is problematic, and may be further complicated by the inherently high mutation rates of some viral types, especially RNA viruses. In these cases, increased sensitivity may be achieved by leveraging position-specific information during the alignment process. Here, we constructed HMMER3-compatible profile hidden Markov models (profile HMMs) from all the virally annotated proteins in RefSeq in an automated fashion using a custom-built bioinformatic pipeline. We then tested the ability of these viral profile HMMs ("vFams") to accurately classify sequences as viral or non-viral. Cross-validation experiments with full-length gene sequences showed that the vFams were able to recall 91% of left-out viral test sequences without erroneously classifying any non-viral sequences into viral protein clusters. Thorough reanalysis of previously published metagenomic datasets with a set of the best-performing vFams showed that they were more sensitive than BLAST for detecting sequences originating from more distant relatives of known viruses. To facilitate the use of the vFams for rapid detection of remote viral homologs in metagenomic data, we provide two sets of vFams, comprising more than 4,000 vFams each, in the HMMER3 format. We also provide the software necessary to build custom profile HMMs or update the vFams as more viruses are discovered (http://derisilab.ucsf.edu/software/vFam).

Project description:IntroductionBiological communities present in soil are essential to sustainable and productive agricultural practices; however, an accurate determination of the ecological status of agricultural soils remains to date an elusive task. An ideal indicator should be pervasive, play a relevant role in the ecosystem, show a rapid and proportional answer to external perturbations and be easily and economically measurable. Rhizobacteria play a major role in determining soil properties, becoming an attractive candidate for the detection of ecological indicators. The application of massive sequencing technologies to metagenomic analysis is providing an increasingly more precise view of the structure and composition of soil communities. In this work, we analyse soil rhizobacterial composition under various stress levels to search for potential ecological indicators.General biodiversity indicatorsOur results suggest that the Shannon index requires observation of a relatively large number of individuals to be representative of the true population diversity, and that the Simpson index may underestimate rare taxa in rhizobacterial environments.Taxonomical classification methodsDetection of indicator taxa requires comparison of taxonomical classification of sequences. We have compared RDP classifier, RTAX and similarity-based taxonomical classification and selected the latter for taxonomical assignment because it provides larger detail.Taxonomy-based ecological indicatorsThe study of significant variations in common, clearly identified, taxa, using paired datasets allows minimization of non-treatment effects and avoidance of false positives. We have identified taxa associated to specific perturbations as well as taxa generally affected in treated soils. Changes in these taxa, or combinations of them, may be used as ecological indicators of soil health. The overall number and magnitude of changes detected in taxonomic groups does also increase with stress. These changes constitute an alternative indicator to measuring specific taxa, although their determination requires large sample sizes, better obtained by massive sequencing.SummaryThe main ecological indicators available are the Shannon index, OTU counts and estimators, overall detection of the number and proportion of changes, and changes of specific indicator taxa. Massive sequencing remains the most accurate tool to measure rhizobacterial ecological indicators. When massive sequencing is not an option, various cultivable taxonomic groups, such as specific groups in the Actinobacteria tree, are attractive as potential indicators of large disruptions to the rhizobiome.

Dataset Information

ENVirT: inference of ecological characteristics of viruses from metagenomic data.

Publications

ENVirT: inference of ecological characteristics of viruses from metagenomic data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure