An empirical test of the concomitantly variable codon hypothesis

Further Simulations and Analyses Demonstrate Open Problems of Phylostratigraphy

Phylostratigraphy, originally designed for gene age estimation by BLAST-based protein homology searches of sequenced genomes, has been widely used for studying patterns and inferring mechanisms of gene origination and evolution. We previously showed by computer simulation that phylostratigraphy underestimates gene age for a nonnegligible fraction of genes and that the underestimation is severer for genes with certain properties such as fast evolution and short protein sequences. Consequently, many previously reported age distributions of gene properties may have been methodological artifacts rather than biological realities. Domazet-Lošo and colleagues recently argued that our simulations were flawed and that phylostratigraphic bias does not impact inferences about gene emergence and evolution. Here we discuss conceptual difficulties of phylostratigraphy, identify numerous problems in Domazet-Lošo et al.’s argument, reconfirm phylostratigraphic error using simulations suggested by Domazet-Lošo and colleagues, and demonstrate that a phylostratigraphic trend claimed to be robust to error disappears when genes likely to be error-resistant are analyzed. We conclude that extreme caution is needed in interpreting phylostratigraphic results because of the inherent biases of the method and that reanalysis using genes exhibiting no error in realistic simulations may help reduce spurious findings.

Cytolethal Distending Toxin: A Unique Variation on the AB Toxin Paradigm

Some of the most potent toxins produced by plants and bacteria are members of a large family known as the AB toxins. AB toxins are generally characterized by a heterogenous complex consisting of two protein chains arranged in various monomeric or polymeric configurations. The newest class within this superfamily is the cytolethal distending toxin (Cdt). The Cdt is represented by a subfamily of toxins produced by a group of taxonomically distinct Gram negative bacteria. Members of this subfamily have a related AB-type chain or subunit configuration and properties distinctive to the AB paradigm. In this review, the unique structural and cytotoxic properties of the Cdt subfamily, target cell specificities, intoxication pathway, modes of action, and relationship to the AB toxin superfamily are compared and contrasted.

Strong epistatic interactions within a single protein

A central question in protein molecular evolution is whether an amino acid that occurs at a given site makes an independent contribution to fitness or whether its contribution depends on other amino acid sites in the protein sequence, a phenomenon known as intragenic epistasis. In the absence of intragenic epistasis, natural selection acts on a protein mutation independent of its genetic background, and the experimentally determined fitness for a mutation should be the same across all genetic backgrounds. We tested this hypothesis by using site-directed mutagenesis to introduce a well-characterized deleterious single amino acid substitution in 56 different hepatitis C virus NS3 protease variants. The catalytic efficiency of the mutated proteases was determined and compared with the corresponding wild-type protein. Fitness effects ranged from lethality to significantly beneficial. Although primarily deleterious and lethal effects were observed (41 and 5 out of 56 tested variants, respectively), deleterious effects ranged from near neutral (-26.7% reduction of fitness) to near lethal (-98.5%). Our findings demonstrate that the introduced amino acid substitution has different fitness effects in different protein variants and provide independent support for the relevant role of intragenic epistasis in protein evolution.

Interchanging functionality among homologous elongation factors using signatures of heterotachy

Numerous models of molecular evolution have been formulated to describe the forces that shape sequence divergence among homologous proteins. These models have greatly enhanced our understanding of evolutionary processes. Rarely are such models empirically tested in the laboratory, and even more rare, are such models exploited to generate novel molecules useful for synthetic biology. Here, we experimentally demonstrate that the heterotachy model of evolution captures signatures of functional divergence among homologous elongation factors (EFs) between bacterial EF-Tu and eukaryotic eEF1A. These EFs are GTPases that participate in protein translation by presenting aminoacylated-tRNAs to the ribosome. Upon release from the ribosome, the EFs are recharged by nucleotide exchange factors EF-Ts in bacteria or eEF1B in eukaryotes. The two nucleotide exchange factors perform analogous functions despite not being homologous proteins. The heterotachy model was used to identify a set of sites in eEF1A/EF-Tu associated with eEF1B binding in eukaryotes and another reciprocal set associated with EF-Ts binding in bacteria. Introduction of bacterial EF-Tu residues at these sites into eEF1A protein efficiently disrupted binding of cognate eEF1B as well as endowed eEF1A with the novel ability to bind bacterial EF-Ts. We further demonstrate that eEF1A variants, unlike yeast wild-type, can function in a reconstituted in vitro bacterial translation system.

The contribution of coevolving residues to the stability of KDO8P synthase

Background

The evolutionary tree of 3-deoxy-D-manno-octulosonate 8-phosphate (KDO8P) synthase (KDO8PS), a bacterial enzyme that catalyzes a key step in the biosynthesis of bacterial endotoxin, is evenly divided between metal and non-metal forms, both having similar structures, but diverging in various degrees in amino acid sequence. Mutagenesis, crystallographic and computational studies have established that only a few residues determine whether or not KDO8PS requires a metal for function. The remaining divergence in the amino acid sequence of KDO8PSs is apparently unrelated to the underlying catalytic mechanism.

Methodology/Principal Findings

The multiple alignment of all known KDO8PS sequences reveals that several residue pairs coevolved, an indication of their possible linkage to a structural constraint. In this study we investigated by computational means the contribution of coevolving residues to the stability of KDO8PS. We found that about 1/4 of all strongly coevolving pairs probably originated from cycles of mutation (decreasing stability) and suppression (restoring it), while the remaining pairs are best explained by a succession of neutral or nearly neutral covarions.

Conclusions/Significance

Both sequence conservation and coevolution are involved in the preservation of the core structure of KDO8PS, but the contribution of coevolving residues is, in proportion, smaller. This is because small stability gains or losses associated with selection of certain residues in some regions of the stability landscape of KDO8PS are easily offset by a large number of possible changes in other regions. While this effect increases the tolerance of KDO8PS to deleterious mutations, it also decreases the probability that specific pairs of residues could have a strong contribution to the thermodynamic stability of the protein.

Mutual information is critically dependent on prior assumptions: would the correct estimate of mutual information please identify itself?

MOTIVATION

Mutual information (MI) is a quantity that measures the dependence between two arbitrary random variables and has been repeatedly used to solve a wide variety of bioinformatic problems. Recently, when attempting to quantify the effects of sampling variance on computed values of MI in proteins, we encountered striking differences among various novel estimates of MI. These differences revealed that estimating the 'true' value of MI is not a straightforward procedure, and minor variations of assumptions yielded remarkably different estimates.

RESULTS

We describe four formally equivalent estimates of MI, three of which explicitly account for sampling variance, that yield non-equal values of MI given exact frequencies. These MI estimates are essentially non-predictive of each other, converging only in the limit of implausibly large datasets. Lastly, we show that all four estimates are biologically reasonable estimates of MI, despite their disparity, since each is actually the Kullback-Leibler divergence between random variables conditioned on equally plausible hypotheses.

CONCLUSIONS

For sparse contingency tables of the type universally observed in protein coevolution studies, our results show that estimates of MI, and hence inferences about physical phenomena such as coevolution, are critically dependent on at least three prior assumptions. These assumptions are: (i) how observation counts relate to expected frequencies; (ii) the relationship between joint and marginal frequencies; and (iii) how non-observed categories are interpreted. In any biologically relevant data, these assumptions will affect the MI estimate as much or more-so than observed data, and are independent of uncertainty in frequency parameters.

Functionally compensating coevolving positions are neither homoplasic nor conserved in clades

We demonstrated that a pair of positions in phosphoglycerate kinase that score highly by three nonparametric covariation measures are important for function even though the positions can be occupied by aliphatic, aromatic, or charged residues. Examination of these pairs suggested that the majority of the covariation scores could be explained by within-clade conservation. However, an analysis of diversity showed that the conservation within clades of covarying pairs was indistinguishable from pairs of positions that do not covary, thus ruling out both clade conservation and extensive homoplasy as means to identify covarying positions. Mutagenesis showed that the residues in the covarying pair were epistatic, with the type of epistasis being dependent on the initial pair. The results show that nonconserved covarying positions that affect protein function can be identified with high precision.

PROCOV: maximum likelihood estimation of protein phylogeny under covarion models and site-specific covarion pattern analysis

Background

The covarion hypothesis of molecular evolution holds that selective pressures on a given amino acid or nucleotide site are dependent on the identity of other sites in the molecule that change throughout time, resulting in changes of evolutionary rates of sites along the branches of a phylogenetic tree. At the sequence level, covarion-like evolution at a site manifests as conservation of nucleotide or amino acid states among some homologs where the states are not conserved in other homologs (or groups of homologs). Covarion-like evolution has been shown to relate to changes in functions at sites in different clades, and, if ignored, can adversely affect the accuracy of phylogenetic inference.

Results

PROCOV (protein covarion analysis) is a software tool that implements a number of previously proposed covarion models of protein evolution for phylogenetic inference in a maximum likelihood framework. Several algorithmic and implementation improvements in this tool over previous versions make computationally expensive tree searches with covarion models more efficient and analyses of large phylogenomic data sets tractable. PROCOV can be used to identify covarion sites by comparing the site likelihoods under the covarion process to the corresponding site likelihoods under a rates-across-sites (RAS) process. Those sites with the greatest log-likelihood difference between a 'covarion' and an RAS process were found to be of functional or structural significance in a dataset of bacterial and eukaryotic elongation factors.

Conclusion

Covarion models implemented in PROCOV may be especially useful for phylogenetic estimation when ancient divergences between sequences have occurred and rates of evolution at sites are likely to have changed over the tree. It can also be used to study lineage-specific functional shifts in protein families that result in changes in the patterns of site variability among subtrees.

Multicopy suppression underpins metabolic evolvability

Our understanding of the origins of new metabolic functions is based upon anecdotal genetic and biochemical evidence. Some auxotrophies can be suppressed by overexpressing substrate-ambiguous enzymes (i.e., those that catalyze the same chemical transformation on different substrates). Other enzymes exhibit weak but detectable catalytic promiscuity in vitro (i.e., they catalyze different transformations on similar substrates). Cells adapt to novel environments through the evolution of these secondary activities, but neither their chemical natures nor their frequencies of occurrence have been characterized en bloc. Here, we systematically identified multifunctional genes within the Escherichia coli genome. We screened 104 single-gene knockout strains and discovered that many (20%) of these auxotrophs were rescued by the overexpression of at least one noncognate E. coli gene. The deleted gene and its suppressor were generally unrelated, suggesting that promiscuity is a product of contingency. This genome-wide survey demonstrates that multifunctional genes are common and illustrates the mechanistic diversity by which their products enhance metabolic robustness and evolvability.