Models of compensatory molecular evolution: Effects of back mutation

Effects of gene duplication, epistasis, recombination and gene conversion on the fixation time of compensatory mutations

Gene duplication is one of the major mechanisms of molecular evolution. Gene duplication enables copies of a gene to accumulate mutations through functional redundancy. If a gene encodes a specific protein that interacts with other proteins, RNA, or DNA, the relaxation of selective constraints caused by gene duplication might contribute to the fixation of compensatory mutations that occur at the interacting sites. In this study, we investigate the effect of gene duplication, epistasis among the duplicated copies and gene conversion on the fixation time of compensatory mutations by extending the original model of compensatory evolution proposed by Kimura in 1985. Our simulation results reveal that the time to fixation of compensatory mutations can be decreased remarkably by gene duplication if one of the duplicated loci can completely mask the deleterious effects of a mutation that occurs at the other locus. Conversely, the fixation time can be increased by gene duplication if such functional compensation is weak. We also show that the combination of the degree of functional compensation and the rate of gene conversion between duplicate loci have contrasting effects on the time to fixation of compensatory mutations.

Secondary structure conservation of the stem-loop IV sub-domain of internal ribosomal entry sites in human rhinovirus clinical isolates

OBJECTIVES

The aim of this study was to investigate the genetic diversity in the stem-loop (SL) IV sub-domain of the human rhinovirus (HRV) internal ribosomal entry site (IRES), which plays key roles in the initiation of viral translation by host protein interaction.

METHODS

The primary SL-IV sequences of 194 HRVs, consisting of 97 reference strains and 97 clinical isolates, including the IRES sub-domains SL-IVa, SL-IVb, SL-IVc, and SL-IVd, were analyzed using Lasergene, MEGA 4, and WebLogo. Additionally, secondary structures of SL-IV were predicted and classified by RNAfold and CentroidHomfold-LAST.

RESULTS

The predicted secondary structures of SL-IV showed variations in the position of bulbs, size of the loop, and length of stems. SL-IVc had the most highly conserved nucleotide sequence, with structures classified into two groups by the location of the poly(C) loop. Of the SL-IV sequences analyzed, 74 (79.56%) were classified in the major group and 19 (20.44%) in the minor group. Thirteen compensatory substitution pairs of SL-IVc contributed to maintaining the stem structure.

CONCLUSIONS

This study showed that the IRES secondary structures of a large number of reference and clinical HRVs were highly conserved, with several compensatory substitutions. It is expected that these results will facilitate investigations into HRV function based on IRES secondary structures.

A model of compensatory molecular evolution involving multiple sites in RNA molecules

Consider two sites under compensatory fitness interaction, such as a Watson-Crick base pair in an RNA helix or two interacting residues in a protein. A mutation at any one of these two sites may reduce the fitness of an individual. However, fitness may be restored by the occurrence of a second mutation at the other site. Kimura modeled this process using a two-locus haploid fitness scheme with two alleles at each locus. He predicted that compensatory evolution following this model is very rare unless selection against the deleterious single mutations is weak and linkage between the interacting sites is tight. Here we investigate the question whether the rate of compensatory evolution increases if we take the context of the two directly interacting sites into account. By "context", we mean the effect of neighboring sites in an RNA helix. Interaction between the focal pair of sites under consideration and the context may lead to so-called indirect compensation. Thus, extending Kimura's classical model of compensatory evolution, we study the effects of both direct and indirect compensation on the rate of compensatory evolution. It is shown that the effects of indirect compensation are very strong. We find that recombination does not slow down the rate of compensatory evolution as predicted by the classical model. Instead, compensatory substitutions may be relatively frequent, even if linkage between the focal interacting sites is loose, selection against deleterious mutations is strong, and mutation rate is low. We compare our theoretical results with data on RNA secondary structures from vertebrate introns.

The dynamics of alternative pathways to compensatory substitution

The role of epistatic interactions among loci is a central question in evolutionary biology and is increasingly relevant in the genomic age. While the population genetics of compensatory substitution have received considerable attention, most studies have focused on the case when natural selection is very strong against deleterious intermediates. In the biologically-plausible scenario of weak to moderate selection there exist two alternate pathways for compensatory substitution. In one pathway, a deleterious mutation becomes fixed prior to occurrence of the compensatory mutation. In the other, the two loci are simultaneously polymorphic. The rates of compensatory substitution along these two pathways and their relative probabilities are functions of the population size, selection strength, mutation rate, and recombination rate. In this paper these rates and path probabilities are derived analytically and verified using population genetic simulations. The expected time durations of these two paths are similar when selection is moderate, but not when selection is weak. The effect of recombination on the dynamics of the substitution process are explored using simulation. Using the derived rates, a phylogenetic substitution model of the compensatory evolution process is presented that could be used for inference of population genetic parameters from interspecific data.