Dataset Information

Unbiased Estimate of Synonymous and Nonsynonymous Substitution Rates with Nonstationary Base Composition.

ABSTRACT: The measurement of synonymous and nonsynonymous substitution rates (dS and dN) is useful for assessing selection operating on protein sequences or for investigating mutational processes affecting genomes. In particular, the ratio dNdS is expected to be a good proxy for ω, the ratio of fixation probabilities of nonsynonymous mutations relative to that of neutral mutations. Standard methods for estimating dN, dS, or ω rely on the assumption that the base composition of sequences is at the equilibrium of the evolutionary process. In many clades, this assumption of stationarity is in fact incorrect, and we show here through simulations and analyses of empirical data that nonstationarity biases the estimate of dN, dS, and ω. We show that the bias in the estimate of ω can be fixed by explicitly taking into consideration nonstationarity in the modeling of codon evolution, in a maximum likelihood framework. Moreover, we propose an exact method for estimating dN and dS on branches, based on stochastic mapping, that can take into account nonstationarity. This method can be directly applied to any kind of codon evolution model, as long as neutrality is clearly parameterized.

SUBMITTER: Gueguen L

PROVIDER: S-EPMC5850866 | biostudies-literature | 2018 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Unbiased Estimate of Synonymous and Nonsynonymous Substitution Rates with Nonstationary Base Composition.

Guéguen Laurent L Duret Laurent L

Molecular biology and evolution 20180301 3

The measurement of synonymous and nonsynonymous substitution rates (dS and dN) is useful for assessing selection operating on protein sequences or for investigating mutational processes affecting genomes. In particular, the ratio dNdS is expected to be a good proxy for ω, the ratio of fixation probabilities of nonsynonymous mutations relative to that of neutral mutations. Standard methods for estimating dN, dS, or ω rely on the assumption that the base composition of sequences is at the equilibr ...[more]

PMID: 29220511

Similar Datasets

Project description:Gene duplication leads to the formation of gene families, wherein purifying or neutral selection maintains the original gene function, while diversifying selection confers new functions onto duplicated genes. The B56 gene family is highly conserved; it is encoded by one gene in protists and fungi, and five genes in vertebrates. B56 regulates protein phosphatase 2A (PP2A), an abundant heterotrimeric serine/threonine phosphatase that functions as a tumor suppressor and consists of a scaffolding "A" and catalytic "C" subunit heterodimer bound to a regulatory "B" subunit. Individual regulatory B56 subunits confer disparate functions onto PP2A in various cell-cell signaling pathways. B56 proteins share a conserved central core domain, but have divergent N- and C-termini which play a role in isoform specificity. We carried out a nonsynonymous/synonymous substitution analysis to better understand the divergence of vertebrate B56 genes. When five B56 paralogs from ten vertebrate species were analyzed, the gene family displayed purifying selection; stronger purifying selection was revealed when individual B56 isoforms were analyzed separately. The B56 core experienced stronger purifying selection than the N- and C-termini, which correlates with the presence of several contacts between the core and the AC heterodimer. Indeed, the majority of the contact points that we analyzed between B56 and the AC heterodimer experienced strong purifying selection. B56 subfamilies showed distinct patterns of selection in their N- and C-termini. The C-terminus of the B56-1 subfamily and the N-terminus of the B56-2 subfamily exhibited strong purifying selection, suggesting that these termini carry out subfamily-specific functions, while the opposite termini exhibited diversifying selection and likely carry out isoform-specific functions. We also found reduced synonymous substitutions at the N- and C-termini when grouping B56 genes by species but not by isoform, suggesting species-specific codon bias may have a role in regulating B56 gene expression.

Project description:BackgroundIt has long been known that rates of synonymous substitutions are unusually low in mitochondrial genes of flowering and other land plants. Although two dramatic exceptions to this pattern have recently been reported, it is unclear how often major increases in substitution rates occur during plant mitochondrial evolution and what the overall magnitude of substitution rate variation is across plants.ResultsA broad survey was undertaken to evaluate synonymous substitution rates in mitochondrial genes of angiosperms and gymnosperms. Although most taxa conform to the generality that plant mitochondrial sequences evolve slowly, additional cases of highly accelerated rates were found. We explore in detail one of these new cases, within the genus Silene. A roughly 100-fold increase in synonymous substitution rate is estimated to have taken place within the last 5 million years and involves only one of ten species of Silene sampled in this study. Examples of unusually slow sequence evolution were also identified. Comparison of the fastest and slowest lineages shows that synonymous substitution rates vary by four orders of magnitude across seed plants. In other words, some plant mitochondrial lineages accumulate more synonymous change in 10,000 years than do others in 100 million years. Several perplexing cases of gene-to-gene variation in sequence divergence within a plant were uncovered. Some of these probably reflect interesting biological phenomena, such as horizontal gene transfer, mitochondrial-to-nucleus transfer, and intragenomic variation in mitochondrial substitution rates, whereas others are likely the result of various kinds of errors.ConclusionThe extremes of synonymous substitution rates measured here constitute by far the largest known range of rate variation for any group of organisms. These results highlight the utility of examining absolute substitution rates in a phylogenetic context rather than by traditional pairwise methods. Why substitution rates are generally so low in plant mitochondrial genomes yet occasionally increase dramatically remains mysterious.

Dataset Information

Unbiased Estimate of Synonymous and Nonsynonymous Substitution Rates with Nonstationary Base Composition.

Publications

Unbiased Estimate of Synonymous and Nonsynonymous Substitution Rates with Nonstationary Base Composition.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets