Efficient representation of uncertainty in multiple sequence alignments using directed acyclic graphs

Abstract

BackgroundA standard procedure in many areas of bioinformatics is to use a single multiple sequence alignment (MSA) as the basis for various types of analysis. However, downstream results may be highly sensitive to the alignment used, and neglecting the uncertainty in the alignment can lead to significant bias in the resulting inference. In recent years, a number of approaches have been developed for probabilistic sampling of alignments, rather than simply generating a single optimum. However, this type of probabilistic information is currently not widely used in the context of downstream inference, since most existing algorithms are set up to make use of a single alignment.ResultsIn this work we present a framework for representing a set of sampled alignments as a directed acyclic graph (DAG) whose nodes are alignment columns; each path through this DAG then represents a valid alignment. Since the probabilities of individual columns can be estimated from empirical frequencies, this approach enables sample-based estimation of posterior alignment probabilities. Moreover, due to conditional independencies between columns, the graph structure encodes a much larger set of alignments than the original set of sampled MSAs, such that the effective sample size is greatly increased.ConclusionsThe alignment DAG provides a natural way to represent a distribution in the space of MSAs, and allows for existing algorithms to be efficiently scaled up to operate on large sets of alignments. As an example, we show how this can be used to compute marginal probabilities for tree topologies, averaging over a very large number of MSAs. This framework can also be used to generate a statistically meaningful summary alignment; example applications show that this summary alignment is consistently more accurate than the majority of the alignment samples, leading to improvements in downstream tree inference.Implementations of the methods described in this article are available at http://statalign.github.io/WeaveAlign.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-015-0516-1) contains supplementary material, which is available to authorized users.

Scaling statistical multiple sequence alignment to large datasets

“…The impact of multiple sequence alignment on downstream analyses is known to be substantial, with errors in multiple sequence alignment producing increased error rates in phylogeny estimation, false detection of positive selection, difficulties in detecting active sites in proteins, etc. [5]. Thus, highly accurate multiple sequence alignment, especially of large datasets spanning large evolutionary distances, is one of the major outstanding bioinformatics problems [6].…”

Historian: accurate reconstruction of ancestral sequences and evolutionary rates

“…These are expected to be more accurate, but generally take much longer. MCMC samplers may be usefully supplemented by decisiontheoretic approaches to summarize a sampling run (Herman et al, 2015). Other potential ways to improve accuracy include context-dependent gap penalties, as used by Muscle (Edgar, 2004), and explicit modeling of tandem duplications (Szalkowski and Anisimova, 2013).…”

Building the Resilience of WSS Utilities to Climate Change and Other Threats

“…Applying this process will help the WSS utility manager (or decision-makers) reconcile short-and long-term priorities, and make decisions to increase the utility's resilience, despite deep uncertainties about future conditions (be these related to climate, demand, or budget). For instance, it has been used to integrate long-term infrastructure investments, near-term operational strategies, and medium-term financial instruments (Herman et al 2015).…”

“…and analysts must select those that are robust to deep uncertainties or adaptable to a changing future. A review of methodologies can be found inHerman et al (2015) andDittrich, Wreford, and Moran (2016). For instance, decision scaling(Brown et al 2012) uses a weather generator to stress-test water projects against a wide range of future climate and hydrological conditions and adapt them so that they perform well in a large number of scenarios.…”