Maximum principle and mutation thresholds for four-letter sequence evolution

A four-state mutation-selection model for the evolution of populations of DNA-sequences is investigated with particular interest in the phenomenon of error thresholds. The mutation model considered is the Kimura 3ST mutation scheme, fitness functions, which determine the selection process, come from the permutation-invariant class. Error thresholds can be found for various fitness functions, the phase diagrams are more interesting than for equivalent two-state models. Results for (small) finite sequence lengths are compared with those for infinite sequence length, obtained via a maximum principle that is equivalent to the principle of minimal free energy in physics.Comment: 25 pages, 16 figure

Developing and applying heterogeneous phylogenetic models with XRate

Modeling sequence evolution on phylogenetic trees is a useful technique in computational biology. Especially powerful are models which take account of the heterogeneous nature of sequence evolution according to the "grammar" of the encoded gene features. However, beyond a modest level of model complexity, manual coding of models becomes prohibitively labor-intensive. We demonstrate, via a set of case studies, the new built-in model-prototyping capabilities of XRate (macros and Scheme extensions). These features allow rapid implementation of phylogenetic models which would have previously been far more labor-intensive. XRate's new capabilities for lineage-specific models, ancestral sequence reconstruction, and improved annotation output are also discussed. XRate's flexible model-specification capabilities and computational efficiency make it well-suited to developing and prototyping phylogenetic grammar models. XRate is available as part of the DART software package: http://biowiki.org/DART .Comment: 34 pages, 3 figures, glossary of XRate model terminolog

Stops making sense: translational trade-offs and stop codon reassignment

Background Efficient gene expression involves a trade-off between (i) premature termination of protein synthesis; and (ii) readthrough, where the ribosome fails to dissociate at the terminal stop. Sense codons that are similar in sequence to stop codons are more susceptible to nonsense mutation, and are also likely to be more susceptible to transcriptional or translational errors causing premature termination. We therefore expect this trade-off to be influenced by the number of stop codons in the genetic code. Although genetic codes are highly constrained, stop codon number appears to be their most volatile feature. Results In the human genome, codons readily mutable to stops are underrepresented in coding sequences. We construct a simple mathematical model based on the relative likelihoods of premature termination and readthrough. When readthrough occurs, the resultant protein has a tail of amino acid residues incorrectly added to the C-terminus. Our results depend strongly on the number of stop codons in the genetic code. When the code has more stop codons, premature termination is relatively more likely, particularly for longer genes. When the code has fewer stop codons, the length of the tail added by readthrough will, on average, be longer, and thus more deleterious. Comparative analysis of taxa with a range of stop codon numbers suggests that genomes whose code includes more stop codons have shorter coding sequences. Conclusions We suggest that the differing trade-offs presented by alternative genetic codes may result in differences in genome structure. More speculatively, multiple stop codons may mitigate readthrough, counteracting the disadvantage of a higher rate of nonsense mutation. This could help explain the puzzling overrepresentation of stop codons in the canonical genetic code and most variants

Actinopolyspora algeriensis sp. nov., a novel halophilic actinomycete isolated from a Saharan soil

A halophilic actinomycete strain designated H19T, was isolated from a Saharan soil in the Bamendil region (Ouargla province, South Algeria) and was characterized taxonomically by using a polyphasic approach. The morphological and chemotaxonomic characteristics of the strain were consistent with those of members of the genus Actinopolyspora, and 16S rRNA gene sequence analysis confirmed that strain H19T was a novel species of the genus Actinopolyspora. DNA–DNA hybridization value between strain H19T and the nearest Actinopolyspora species, A. halophila, was clearly below the 70 % threshold. The genotypic and phenotypic data showed that the organism represents a novel species of the genus Actinopolyspora for which the name Actinopolyspora algeriensis sp. nov. is proposed, with the type strain H19T (= DSM 45476T = CCUG 62415T)

Nutrition issues in Codex: health claims, nutrient reference values and WTO agreements: a conference report

BACKGROUND: Codex documents may be used as educational and consensus materials for member governments. Also, the WTO SPS Agreement recognizes Codex as the presumptive international authority on food issues. Nutrient bioavailability is a critical factor in determining the ability of nutrients to provide beneficial effects. Bioavailability also influences the quantitative dietary requirements that are the basis of nutrient intake recommendations and NRVs. HEALTH CLAIMS: Codex, EFSA and some national regulatory authorities have established guidelines or regulations that will permit several types of health claims. The scientific basis for claims has been established by the US FDA and EFSA, but not yet by Codex. Evidence-based nutrition differs from evidence-based medicine, but the differences are only recently gaining recognition. Health claims on foods may provide useful information to consumers, but many will interpret the information to mean that they can rely upon the food or nutrient to eliminate a disease risk. NUTRIENT REFERENCE VALUES: NRVs are designed to provide a quantitative basis for comparing the nutritive values of foods, helping to illustrate how specific foods fit into the overall diet. The INL-98 and the mean of adult male and female values provide NRVs that are sufficient when used as targets for individual intakes by most adults. WORLD TRADE ORGANIZATION AGREEMENTS: WTO recognizes Codex as the primary international authority on food issues. Current regulatory schemes based on recommended dietary allowances are trade restrictive. A substantial number of decisions by the EFSA could lead to violation of WTO agreements

A Novel 14C-Postlabeling Assay Using Accelerator Mass Spectrometry For the Detection of O6-Methyldeoxyguanosine Adducts

Accelerator mass spectrometry (AMS) is currently one of the most sensitive methods available for the trace detection of DNA adducts and is particularly valuable for measuring adducts in humans or animal models. However, the standard approach requires administration of a radiolabeled compound. As an alternative, we have developed a preliminary {sup 14}C-postlabeling assay for detection of the highly mutagenic O{sup 6}-MedG, by AMS. Procedures were developed for derivatizing O{sup 6}-MedG using unlabeled acetic anhydride. Using conventional LC-MS analysis, the limit of detection for the major product, triacetylated O{sup 6}-MedG, was 10 fmoles. On reaction with {sup 14}C-acetic anhydride, using a specially designed enclosed system, the predominant product was {sup 14}C-di-acetyl O{sup 6}-MedG. This change in reaction profile was due to a modification of the reaction procedure, introduced as a necessary safety precaution. The limit of detection for {sup 14}C-diacetyl O{sup 6}-MedG by AMS was determined as 79 attomoles, {approx}18,000 fold lower than that achievable by LSC. Although the assay has so far only been carried out with labeled standards, the degree of sensitivity obtained illustrates the potential of this assay for measuring O{sup 6}-MedG levels in humans

Sigma-2: Multiple sequence alignment of non-coding DNA via an evolutionary model

<p>Abstract</p> <p>Background</p> <p>While most multiple sequence alignment programs expect that all or most of their input is known to be homologous, and penalise insertions and deletions, this is not a reasonable assumption for non-coding DNA, which is much less strongly conserved than protein-coding genes. Arguing that the goal of sequence alignment should be the detection of <it>homology </it>and not <it>similarity</it>, we incorporate an evolutionary model into a previously published multiple sequence alignment program for non-coding DNA, Sigma, as a sensitive likelihood-based way to assess the significance of alignments. Version 1 of Sigma was successful in eliminating spurious alignments but exhibited relatively poor sensitivity on synthetic data. Sigma 1 used a <it>p</it>-value (the probability under the "null hypothesis" of non-homology) to assess the significance of alignments, and, optionally, a background model that captured short-range genomic correlations. Sigma version 2, described here, retains these features, but calculates the <it>p</it>-value using a sophisticated evolutionary model that we describe here, and also allows for a transition matrix for different substitution rates from and to different nucleotides. Our evolutionary model takes separate account of mutation and fixation, and can be extended to allow for locally differing functional constraints on sequence.</p> <p>Results</p> <p>We demonstrate that, on real and synthetic data, Sigma-2 significantly outperforms other programs in specificity to genuine homology (that is, it minimises alignment of spuriously similar regions that do not have a common ancestry) while it is now as sensitive as the best current programs.</p> <p>Conclusions</p> <p>Comparing these results with an extrapolation of the best results from other available programs, we suggest that conservation rates in intergenic DNA are often significantly over-estimated. It is increasingly important to align non-coding DNA correctly, in regulatory genomics and in the context of whole-genome alignment, and Sigma-2 is an important step in that direction.</p

Probabilistic Phylogenetic Inference with Insertions and Deletions

A fundamental task in sequence analysis is to calculate the probability of a multiple alignment given a phylogenetic tree relating the sequences and an evolutionary model describing how sequences change over time. However, the most widely used phylogenetic models only account for residue substitution events. We describe a probabilistic model of a multiple sequence alignment that accounts for insertion and deletion events in addition to substitutions, given a phylogenetic tree, using a rate matrix augmented by the gap character. Starting from a continuous Markov process, we construct a non-reversible generative (birth–death) evolutionary model for insertions and deletions. The model assumes that insertion and deletion events occur one residue at a time. We apply this model to phylogenetic tree inference by extending the program dnaml in phylip. Using standard benchmarking methods on simulated data and a new “concordance test” benchmark on real ribosomal RNA alignments, we show that the extended program dnamlε improves accuracy relative to the usual approach of ignoring gaps, while retaining the computational efficiency of the Felsenstein peeling algorithm

Evaluation of intra- and interspecific divergence of satellite DNA sequences by nucleotide frequency calculation and pairwise sequence comparison

Satellite DNA sequences are known to be highly variable and to have been subjected to concerted evolution that homogenizes member sequences within species. We have analyzed the mode of evolution of satellite DNA sequences in four fishes from the genus Diplodus by calculating the nucleotide frequency of the sequence array and the phylogenetic distances between member sequences. Calculation of nucleotide frequency and pairwise sequence comparison enabled us to characterize the divergence among member sequences in this satellite DNA family. The results suggest that the evolutionary rate of satellite DNA in D. bellottii is about two-fold greater than the average of the other three fishes, and that the sequence homogenization event occurred in D. puntazzo more recently than in the others. The procedures described here are effective to characterize mode of evolution of satellite DNA