Dataset Information

Genetic distance for a general non-stationary markov substitution process.

ABSTRACT: The genetic distance between biological sequences is a fundamental quantity in molecular evolution. It pertains to questions of rates of evolution, existence of a molecular clock, and phylogenetic inference. Under the class of continuous-time substitution models, the distance is commonly defined as the expected number of substitutions at any site in the sequence. We eschew the almost ubiquitous assumptions of evolution under stationarity and time-reversible conditions and extend the concept of the expected number of substitutions to nonstationary Markov models where the only remaining constraint is of time homogeneity between nodes in the tree. Our measure of genetic distance reduces to the standard formulation if the data in question are consistent with the stationarity assumption. We apply this general model to samples from across the tree of life to compare distances so obtained with those from the general time-reversible model, with and without rate heterogeneity across sites, and the paralinear distance, an empirical pairwise method explicitly designed to address nonstationarity. We discover that estimates from both variants of the general time-reversible model and the paralinear distance systematically overestimate genetic distance and departure from the molecular clock. The magnitude of the distance bias is proportional to departure from stationarity, which we demonstrate to be associated with longer edge lengths. The marked improvement in consistency between the general nonstationary Markov model and sequence alignments leads us to conclude that analyses of evolutionary rates and phylogenies will be substantively improved by application of this model.

SUBMITTER: Kaehler BD

PROVIDER: S-EPMC4380038 | biostudies-literature | 2015 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Genetic distance for a general non-stationary markov substitution process.

Kaehler Benjamin D BD Yap Von Bing VB Zhang Rongli R Huttley Gavin A GA

Systematic biology 20141209 2

The genetic distance between biological sequences is a fundamental quantity in molecular evolution. It pertains to questions of rates of evolution, existence of a molecular clock, and phylogenetic inference. Under the class of continuous-time substitution models, the distance is commonly defined as the expected number of substitutions at any site in the sequence. We eschew the almost ubiquitous assumptions of evolution under stationarity and time-reversible conditions and extend the concept of t ...[more]

PMID: 25503772

Similar Datasets

Project description:INTRODUCTION:Allergic rhino-conjunctivitis (ARC) is an IgE-mediated disease that occurs after exposure to indoor or outdoor allergens, or to non-specific triggers. Effective treatment options for seasonal ARC are available, but the economic aspects and burden of these therapies are not of secondary importance, also considered that the prevalence of ARC has been estimated at 23% in Europe. For these reasons, we propose a novel flexible cost-effectiveness analysis (CEA) model, intended to provide healthcare professionals and policymakers with useful information aimed at cost-effective interventions for grass-pollen induced allergic rhino-conjunctivitis (ARC). METHODS:Treatments compared are: 1. no AIT, first-line symptomatic drug-therapy with no allergoid immunotherapy (AIT). 2. SCIT, subcutaneous immunotherapy. 3. SLIT, sublingual immunotherapy. The proposed model is a non-stationary Markovian model, that is flexible enough to reflect those treatment-related problems often encountered in real-life and clinical practice, but that cannot be adequately represented in randomized clinical trials (RCTs). At the same time, we described in detail all the structural elements of the model as well as its input parameters, in order to minimize any issue of transparency and facilitate the reproducibility and circulation of the results among researchers. RESULTS:Using the no AIT strategy as a comparator, and the Incremental Cost Effectiveness Ratio (ICER) as a statistic to summarize the cost-effectiveness of a health care intervention, we could conclude that: SCIT systematically outperforms SLIT, except when a full societal perspective is considered. For example, for T = 9 and a pollen season of 60 days, we have ICER = €16,729 for SCIT vs. ICER = €15,116 for SLIT (in the full societal perspective).For longer pollen seasons or longer follow-up duration the ICER decreases, because each patient experiences a greater clinical benefit over a larger time span, and Quality-adjusted Life Year (QALYs) gained per cycle increase accordingly.Assuming that no clinical benefit is achieved after premature discontinuation, and that at least three years of immunotherapy are required to improve clinical manifestations and perceiving a better quality of life, ICERs become far greater than €30,000.If the immunotherapy is effective only at the peak of the pollen season, the relative ICERs rise sharply. For example, in the scenario where no clinical benefit is present after premature discontinuation of immunotherapy, we have ICER = €74,770 for SCIT vs. ICER = €152,110 for SLIT.The distance between SCIT and SLIT strongly depends on under which model the interventions are meta-analyzed. CONCLUSIONS:Even though there is a considerable evidence that SCIT outperforms SLIT, we could not state that both SCIT and SLIT (or only one of these two) can be considered cost-effective for ARC, as a reliable threshold value for cost-effectiveness set by national regulatory agencies for pharmaceutical products is missing. Moreover, the impact of model input parameters uncertainty on the reliability of our conclusions needs to be investigated further.

Project description:Intrinsically disordered proteins (IDPs) or proteins with disordered regions (IDRs) do not have a well-defined tertiary structure, but perform a multitude of functions, often relying on their native disorder to achieve the binding flexibility through changing to alternative conformations. Intrinsic disorder is frequently found in all three kingdoms of life, and may occur in short stretches or span whole proteins. To date most studies contrasting the differences between ordered and disordered proteins focused on simple summary statistics. Here, we propose an evolutionary approach to study IDPs, and contrast patterns specific to ordered protein regions and the corresponding IDRs.Two empirical Markov models of amino acid substitutions were estimated, based on a large set of multiple sequence alignments with experimentally verified annotations of disordered regions from the DisProt database of IDPs. We applied new methods to detect differences in Markovian evolution and evolutionary rates between IDRs and the corresponding ordered protein regions. Further, we investigated the distribution of IDPs among functional categories, biochemical pathways and their preponderance to contain tandem repeats.We find significant differences in the evolution between ordered and disordered regions of proteins. Most importantly we find that disorder promoting amino acids are more conserved in IDRs, indicating that in some cases not only amino acid composition but the specific sequence is important for function. This conjecture is also reinforced by the observation that for of our data set IDRs evolve more slowly than the ordered parts of the proteins, while we still support the common view that IDRs in general evolve more quickly. The improvement in model fit indicates a possible improvement for various types of analyses e.g. de novo disorder prediction using a phylogenetic Hidden Markov Model based on our matrices showed a performance similar to other disorder predictors.

Dataset Information

Genetic distance for a general non-stationary markov substitution process.

Publications

Genetic distance for a general non-stationary markov substitution process.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets