Accounting for alignment uncertainty in phylogenomics.
Ontology highlight
ABSTRACT:
SUBMITTER: Wu M
PROVIDER: S-EPMC3260272 | biostudies-literature | 2012
REPOSITORIES: biostudies-literature
Publications
Accounting for alignment uncertainty in phylogenomics.
Wu Martin M Chatterji Sourav S Eisen Jonathan A JA
PloS one 20120117 1
Uncertainty in multiple sequence alignments has a large impact on phylogenetic analyses. Little has been done to evaluate the quality of individual positions in protein sequence alignments, which directly impact the accuracy of phylogenetic trees. Here we describe ZORRO, a probabilistic masking program that accounts for alignment uncertainty by assigning confidence scores to each alignment position. Using the BALIBASE database and in simulation studies, we demonstrate that masking by ZORRO signi ...[more]
Similar Datasets
Project description:Inference of multiple sequence alignments (MSAs) is a critical part of phylogenetic and comparative genomics studies. However, from the same set of sequences different MSAs are often inferred, depending on the methodologies used and the assumed parameters. Much effort has recently been devoted to improving the ability to identify unreliable alignment regions. Detecting such unreliable regions was previously shown to be important for downstream analyses relying on MSAs, such as the detection of positive selection. Here we developed GUIDANCE2, a new integrative methodology that accounts for: (i) uncertainty in the process of indel formation, (ii) uncertainty in the assumed guide tree and (iii) co-optimal solutions in the pairwise alignments, used as building blocks in progressive alignment algorithms. We compared GUIDANCE2 with seven methodologies to detect unreliable MSA regions using extensive simulations and empirical benchmarks. We show that GUIDANCE2 outperforms all previously developed methodologies. Furthermore, GUIDANCE2 also provides a set of alternative MSAs which can be useful for downstream analyses. The novel algorithm is implemented as a web-server, available at: http://guidance.tau.ac.il.
| S-EPMC4489236 | biostudies-other
Project description:As in many other areas of science, systems biology makes extensive use of statistical association and significance estimates in contingency tables, a type of categorical data analysis known in this field as enrichment (also over-representation or enhancement) analysis. In spite of efforts to create probabilistic annotations, especially in the Gene Ontology context, or to deal with uncertainty in high throughput-based datasets, current enrichment methods largely ignore this probabilistic information since they are mainly based on variants of the Fisher Exact Test.We developed an open-source R-based software to deal with probabilistic categorical data analysis, ProbCD, that does not require a static contingency table. The contingency table for the enrichment problem is built using the expectation of a Bernoulli Scheme stochastic process given the categorization probabilities. An on-line interface was created to allow usage by non-programmers and is available at: http://xerad.systemsbiology.net/ProbCD/.We present an analysis framework and software tools to address the issue of uncertainty in categorical data analysis. In particular, concerning the enrichment analysis, ProbCD can accommodate: (i) the stochastic nature of the high-throughput experimental techniques and (ii) probabilistic gene annotation.
| S-EPMC2169266 | biostudies-literature
Project description:Toxicologists and pharmacologists often describe toxicity of a chemical using parameters of a nonlinear regression model. Thus estimation of parameters of a nonlinear regression model is an important problem. The estimates of the parameters and their uncertainty estimates depend upon the underlying error variance structure in the model. Typically, a priori the researcher would know if the error variances are homoscedastic (i.e., constant across dose) or if they are heteroscedastic (i.e., the variance is a function of dose). Motivated by this concern, in this article we introduce an estimation procedure based on preliminary test which selects an appropriate estimation procedure accounting for the underlying error variance structure. Since outliers and influential observations are common in toxicological data, the proposed methodology uses M-estimators. The asymptotic properties of the preliminary test estimator are investigated; in particular its asymptotic covariance matrix is derived. The performance of the proposed estimator is compared with several standard estimators using simulation studies. The proposed methodology is also illustrated using a data set obtained from the National Toxicology Program.
| S-EPMC3278194 | biostudies-literature
Project description:Sensorimotor coordination is thought to rely on cerebellar-based internal models for state estimation, but the underlying neural mechanisms and specific contribution of the cerebellar components is unknown. A central aspect of any inferential process is the representation of uncertainty or conversely precision characterizing the ensuing estimates. Here, we discuss the possible contribution of inhibition to the encoding of precision of neural representations in the granular layer of the cerebellar cortex. Within this layer, Golgi cells influence excitatory granule cells, and their action is critical in shaping information transmission downstream to Purkinje cells. In this review, we equate the ensuing excitation-inhibition balance in the granular layer with the outcome of a precision-weighted inferential process, and highlight the physiological characteristics of Golgi cell inhibition that are consistent with such computations.
| S-EPMC8059656 | biostudies-literature
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.
| S-EPMC6058391 | biostudies-literature
Project description:Advancements in next-generation sequencing techniques have led to a substantial increase in the genomic information available for analyses in evolutionary biology. As such, this data requires the exponential growth in bioinformatic methods and expertise required to understand such vast quantities of genomic data. Alignment-free phylogenomics offer an alternative approach for large-scale analyses that may have the potential to address these challenges. The evolutionary relationships between various species within the trypanosomatid family, specifically members belonging to the genera Leishmania and Trypanosoma have been extensively studies over the last 30 years. However, there is a need for a more exhaustive analysis of the Trypanosomatidae, summarising the evolutionary patterns amongst the entire family of these important protists. The mitochondrial DNA of the trypanosomatids, better known as the kinetoplast, represents a valuable taxonomic marker given its unique presence across all kinetoplastid protozoans. The aim of this study was to validate the reliability and robustness of alignment-free approaches for phylogenomic analyses and its applicability to reconstruct the evolutionary relationships between the trypanosomatid family. In the present study, alignment-free analyses demonstrated the strength of these methods, particularly when dealing with large datasets compared to the traditional phylogenetic approaches. We present a maxicircle genome phylogeny of 46 species spanning the trypanosomatid family, demonstrating the superiority of the maxicircle for the analysis and taxonomic resolution of the Trypanosomatidae.
| S-EPMC6789588 | biostudies-literature
Project description:Until now the most efficient solution to align nucleotide sequences containing open reading frames was to use indirect procedures that align amino acid translation before reporting the inferred gap positions at the codon level. There are two important pitfalls with this approach. Firstly, any premature stop codon impedes using such a strategy. Secondly, each sequence is translated with the same reading frame from beginning to end, so that the presence of a single additional nucleotide leads to both aberrant translation and alignment.We present an algorithm that has the same space and time complexity as the classical Needleman-Wunsch algorithm while accommodating sequencing errors and other biological deviations from the coding frame. The resulting pairwise coding sequence alignment method was extended to a multiple sequence alignment (MSA) algorithm implemented in a program called MACSE (Multiple Alignment of Coding SEquences accounting for frameshifts and stop codons). MACSE is the first automatic solution to align protein-coding gene datasets containing non-functional sequences (pseudogenes) without disrupting the underlying codon structure. It has also proved useful in detecting undocumented frameshifts in public database sequences and in aligning next-generation sequencing reads/contigs against a reference coding sequence.MACSE is distributed as an open-source java file executable with freely available source code and can be used via a web interface at: http://mbb.univ-montp2.fr/macse.
| S-EPMC3174933 | biostudies-literature
Project description:BACKGROUND:Despite evidence suggesting that air pollution-related health effects differ by emissions source, epidemiologic studies on fine particulate matter (PM2.5) infrequently differentiate between particles from different sources. Those that do rarely account for the uncertainty of source apportionment methods. METHODS:For each day in a 12-year period (1998-2010) in Atlanta, GA, we estimated daily PM2.5 source contributions from a Bayesian ensemble model that combined four source apportionment methods including chemical transport and receptor-based models. We fit Poisson generalized linear models to estimate associations between source-specific PM2.5 concentrations and cardiorespiratory emergency department visits (n = 1,598,117). We propagated uncertainty in the source contribution estimates through analyses using multiple imputation. RESULTS:Respiratory emergency department visits were positively associated with biomass burning and secondary organic carbon. For a 1 µg/m increase in PM2.5 from biomass burning during the past 3 days, the rate of visits for all respiratory outcomes increased by 0.4% (95% CI 0.0%, 0.7%). There was less evidence for associations between PM2.5 sources and cardiovascular outcomes, with the exception of ischemic stroke, which was positively associated with most PM2.5 sources. Accounting for the uncertainty of source apportionment estimates resulted, on average, in an 18% increase in the standard error for rate ratio estimates for all respiratory and cardiovascular emergency department visits, but inflation varied across specific sources and outcomes, ranging from 2% to 39%. CONCLUSIONS:This study provides evidence of associations between PM2.5 sources and some cardiorespiratory outcomes and quantifies the impact of accounting for variability in source apportionment approaches.
| S-EPMC6768727 | biostudies-literature
Project description:Alternative splicing is a biological process during gene expression that allows a single gene to code for multiple proteins. However, splicing patterns can be altered in some conditions or diseases. Here, we present BANDITS, a R/Bioconductor package to perform differential splicing, at both gene and transcript level, based on RNA-seq data. BANDITS uses a Bayesian hierarchical structure to explicitly model the variability between samples and treats the transcript allocation of reads as latent variables. We perform an extensive benchmark across both simulated and experimental RNA-seq datasets, where BANDITS has extremely favourable performance with respect to the competitors considered.
| S-EPMC7075019 | biostudies-literature
Project description:BackgroundA common approach to the application of epidemiological models is to determine a single (point estimate) parameterisation using the information available in the literature. However, in many cases there is considerable uncertainty about parameter values, reflecting both the incomplete nature of current knowledge and natural variation, for example between farms. Furthermore model outcomes may be highly sensitive to different parameter values. Paratuberculosis is an infection for which many of the key parameter values are poorly understood and highly variable, and for such infections there is a need to develop and apply statistical techniques which make maximal use of available data.ResultsA technique based on Latin hypercube sampling combined with a novel reweighting method was developed which enables parameter uncertainty and variability to be incorporated into a model-based framework for estimation of prevalence. The method was evaluated by applying it to a simulation of paratuberculosis in dairy herds which combines a continuous time stochastic algorithm with model features such as within herd variability in disease development and shedding, which have not been previously explored in paratuberculosis models. Generated sample parameter combinations were assigned a weight, determined by quantifying the model's resultant ability to reproduce prevalence data. Once these weights are generated the model can be used to evaluate other scenarios such as control options. To illustrate the utility of this approach these reweighted model outputs were used to compare standard test and cull control strategies both individually and in combination with simple husbandry practices that aim to reduce infection rates.ConclusionsThe technique developed has been shown to be applicable to a complex model incorporating realistic control options. For models where parameters are not well known or subject to significant variability, the reweighting scheme allowed estimated distributions of parameter values to be combined with additional sources of information, such as that available from prevalence distributions, resulting in outputs which implicitly handle variation and uncertainty. This methodology allows for more robust predictions from modelling approaches by allowing for parameter uncertainty and combining different sources of information, and is thus expected to be useful in application to a large number of disease systems.
| S-EPMC3544565 | biostudies-literature