Stability and evolution of overlapping genes

Abstract.— When the same sequence of nucleotides codes for regions of more than one functional polypeptide, this sequence contains overlapping genes. Overlap is most common in rapidly evolving genomes with high mutation rates such as viruses, bacteria, and mitochondria. Overlap is thought to be important as: (1) a means of compressing a maximum amount of information into short sequences of structural genes; and (2) as a mechanism for regulating gene expression through translational coupling of functionally related polypeptides. The stability of overlapping codes is examined in relation to the information cost of overlap and the mutation rate of the genome. The degree of overlap in a given population will tend to become monomorphic. Evolution toward partial overlap of genes is shown to depend on a convex cost function of overlap. Overlap does not evolve when expression of overlapping genes is mutually exclusive and produced by rare mutations to the wild-type genome. Assuming overlap increases coupling between functionally related genes, the conditions favoring overlap are explored in relation to the kinetics of gene activation and decay. Coupling is most effective for genes in which the gene overlapping at its 5'end (leading gene) decays rapidly, while the gene overlapping at the 3'end (induced gene) decays slowly. If gene expression can feedback on itself (autocatalysis), then high rates of activation favor overlap.     

Evolutionary coupling between the deleteriousness of gene mutations and the amount of non-coding sequences

The phenotypic effects of random mutations depend on both the architecture of the genome and the gene-trait relationships. Both levels thus play a key role in the mutational variability of the phenotype, and hence in the long-term evolutionary success of the lineage. Here, by simulating the evolution of organisms with flexible genomes, we show that the need for an appropriate phenotypic variability induces a relationship between the deleteriousness of gene mutations and the quantity of non-coding sequences maintained in the genome. The more deleterious the gene mutations, the shorter the intergenic sequences. Indeed, in a shorter genome, fewer genes are affected by rearrangements (duplications, deletions, inversions, translocations) at each replication, which compensates for the higher impact of each gene mutation. This spontaneous adjustment of genome structure allows the organisms to retain the same average fitness loss per replication, despite the higher impact of single gene mutations. These results show how evolution can generate unexpected couplings between distinct organization levels.     

Purging deleterious mutations under self fertilization: paradoxical recovery in fitness with increasing mutation rate in Caenorhabditis elegans

Background

The accumulation of deleterious mutations can drastically reduce population mean fitness. Self-fertilization is thought to be an effective means of purging deleterious mutations. However, widespread linkage disequilibrium generated and maintained by self-fertilization is predicted to reduce the efficacy of purging when mutations are present at multiple loci.

Methodology/Principal Findings

We tested the ability of self-fertilizing populations to purge deleterious mutations at multiple loci by exposing obligately self-fertilizing populations of Caenorhabditis elegans to a range of elevated mutation rates and found that mutations accumulated, as evidenced by a reduction in mean fitness, in each population. Therefore, purging in obligate selfing populations is overwhelmed by an increase in mutation rate. Surprisingly, we also found that obligate and predominantly self-fertilizing populations exposed to very high mutation rates exhibited consistently greater fitness than those subject to lesser increases in mutation rate, which contradicts the assumption that increases in mutation rate are negatively correlated with fitness. The high levels of genetic linkage inherent in self-fertilization could drive this fitness increase.

Conclusions

Compensatory mutations can be more frequent under high mutation rates and may alleviate a portion of the fitness lost due to the accumulation of deleterious mutations through epistatic interactions with deleterious mutations. The prolonged maintenance of tightly linked compensatory and deleterious mutations facilitated by self-fertilization may be responsible for the fitness increase as linkage disequilibrium between the compensatory and deleterious mutations preserves their epistatic interaction.     

Evolution by Gene Duplication and Compensatory Advantageous Mutations

Relaxation of selective constraint is thought to play an important role for evolution by gene duplication, in connection with compensatory advantageous mutant substitutions. Models were investigated by incorporating gene duplication by unequal crossing over, selection, mutation and random genetic drift into Monte Carlo simulations. Compensatory advantageous mutations were introduced, and simulations were carried out with and without relaxation, when genes are redundant on chromosomes. Relaxation was introduced by assuming that deleterious mutants have no effect on fitness, so long as one or more genes free of such mutations remain in the array. Compensatory mutations are characterized by the intermediate deleterious step of their substitutions, and therefore relaxation by gene redundancy is important. Through extensive Monte Carlo simulations, it was found that compensatory mutant substitutions require relaxation in addition to gene duplication, when mutant effects are large. However when mutant effects are small, such that the product of selection coefficient and population size is around unity, evolution by compensatory mutation is enhanced by gene duplication even without relaxation.     

The fate of duplicated genes: loss or new function?

Gene duplication events are important sources of novel gene functions. However, more often than not, a duplicate gene may lose its function and become a pseudogene. What is the relative frequency of these two scenarios: functional divergence versus gene loss? Given that most non-neutral mutations are deleterious, gene loss should be far more frequent than divergence. However, a recent empirical study suggests that about 50% of all gene duplications will lead to functional divergence. The study infers the frequency of functional divergence from the size distribution of gene families produced by two successive genome duplications early in vertebrate evolution. Reasons for this unexpectedly high frequency of functional divergence are discussed.     

Compensatory mutations are repeatable and clustered within proteins

Compensatory mutations improve fitness in genotypes that contain deleterious mutations but have no beneficial effects otherwise. As such, compensatory mutations represent a very specific form of epistasis. We show that intragenic compensatory mutations occur non-randomly over gene sequence. Compensatory mutations are more likely to appear at some sites than others. Moreover, the sites of compensatory mutations are more likely than expected by chance to be near the site of the original deleterious mutation. Furthermore, compensatory mutations tend to occur more commonly in certain regions of the protein even when controlling for clustering around the site of the deleterious mutation. These results suggest that compensatory evolution at the protein level is partially predictable and may be convergent.     

The genetic basis of inbreeding depression

Data on the effects of inbreeding on fitness components are reviewed in the light of population genetic models of the possible genetic causes of inbreeding depression. Deleterious mutations probably play a major role in causing inbreeding depression. Putting together the different kinds of quantitative genetic data, it is difficult to account for the very large effects of inbreeding on fitness in Drosophila and outcrossing plants without a significant contribution from variability maintained by selection. Overdominant effects of alleles on fitness components seem not to be important in most cases. Recessive or partially recessive deleterious effects of alleles, some maintained by mutation pressure and some by balancing selection, thus seem to be the most important source of inbreeding depression. Possible experimental approaches to resolving outstanding questions are discussed.     

Engineering gene overlaps to sustain genetic constructs in vivo

Evolution is often an obstacle to the engineering of stable biological systems due to the selection of mutations inactivating costly gene circuits. Gene overlaps induce important constraints on sequences and their evolution. We show that these constraints can be harnessed to increase the stability of costly genes by purging loss-of-function mutations. We combine computational and synthetic biology approaches to rationally design an overlapping reading frame expressing an essential gene within an existing gene to protect. Our algorithm succeeded in creating overlapping reading frames in 80% of E. coli genes. Experimentally, scoring mutations in both genes of such overlapping construct, we found that a significant fraction of mutations impacting the gene to protect have a deleterious effect on the essential gene. Such an overlap thus protects a costly gene from removal by natural selection by associating the benefit of this removal with a larger or even lethal cost. In our synthetic constructs, the overlap converts many of the possible mutants into evolutionary dead-ends, reducing the evolutionary potential of the system and thus increasing its stability over time.     

Historical Contingency Causes Divergence in Adaptive Expression of the lac Operon

Populations of Escherichia coli selected in constant and fluctuating environments containing lactose often adapt by substituting mutations in the lacI repressor that cause constitutive expression of the lac operon. These mutations occur at a high rate and provide a significant benefit. Despite this, eight of 24 populations evolved for 8,000 generations in environments containing lactose contained no detectable repressor mutations. We report here on the basis of this observation. We find that, given relevant mutation rates, repressor mutations are expected to have fixed in all evolved populations if they had maintained the same fitness effect they confer when introduced to the ancestor. In fact, reconstruction experiments demonstrate that repressor mutations have become neutral or deleterious in those populations in which they were not detectable. Populations not fixing repressor mutations nevertheless reached the same fitness as those that did fix them, indicating that they followed an alternative evolutionary path that made redundant the potential benefit of the repressor mutation, but involved unique mutations of equivalent benefit. We identify a mutation occurring in the promoter region of the uspB gene as a candidate for influencing the selective choice between these paths. Our results detail an example of historical contingency leading to divergent evolutionary outcomes.     

Comparable Rates of Gene Loss and Functional Divergence after Genome Duplications Early in Vertebrate Evolution

Duplicated genes are an important source of new protein functions and novel developmental and physiological pathways. Whereas most models for fate of duplicated genes show that they tend to be rapidly lost, models for pathway evolution suggest that many duplicated genes rapidly acquire novel functions. Little empirical evidence is available, however, for the relative rates of gene loss vs. divergence to help resolve these contradictory expectations. Gene families resulting from genome duplications provide an opportunity to address this apparent contradiction. With genome duplication, the number of duplicated genes in a gene family is at most 2(n), where n is the number of duplications. The size of each gene family, e.g., 1, 2, 3, . . . , 2(n), reflects the patterns of gene loss vs. functional divergence after duplication. We focused on gene families in humans and mice that arose from genome duplications in early vertebrate evolution and we analyzed the frequency distribution of gene family size, i.e., the number of families with two, three or four members. All the models that we evaluated showed that duplicated genes are almost as likely to acquire a new and essential function as to be lost through acquisition of mutations that compromise protein function. An explanation for the unexpectedly high rate of functional divergence is that duplication allows genes to accumulate more neutral than disadvantageous mutations, thereby providing more opportunities to acquire diversified functions and pathways.     

Simulation Study of a Multigene Family, with Special Reference to the Evolution of Compensatory Advantageous Mutations

We investigate the evolution of a multigene family incorporating the forces of drift, mutation, gene conversion, unequal crossing over and selection. The use of simulation studies is required due to the complexity of the model. Selection is modeled in two modes: positive selection as a function of the number of different beneficial alleles and negative selection against deleterious alleles. We assume that gene conversion is unbiased, and that all mutations are initially deleterious. Compensation between mutants creates beneficial and neutral alleles, and allowances are made for compensatory mutations either within or between the members of a multigene family. We find that gene conversion can enhance the rate of acquisition of compensatory advantageous mutations when genes are redundant.     

Epistasis and the adaptability of an RNA virus

We have explored the patterns of fitness recovery in the vesicular stomatitis RNA virus. We show that, in our experimental setting, reversions to the wild-type genotype were rare and fitness recovery was at least partially driven by compensatory mutations. We compared compensatory adaptation for genotypes carrying (1) mutations with varying deleterious fitness effects, (2) one or two deleterious mutations, and (3) pairs of mutations showing differences in the strength and sign of epistasis. In all cases, we found that the rate of fitness recovery and the proportion of reversions were positively affected by population size. Additionally, we observed that mutations with large fitness effect were always compensated faster than mutations with small fitness effect. Similarly, compensatory evolution was faster for genotypes carrying a single deleterious mutation than for those carrying pairs of mutations. Finally, for genotypes carrying two deleterious mutations, we found evidence of a negative correlation between the epistastic effect and the rate of compensatory evolution.     

PURGING THE GENOME WITH SEXUAL SELECTION: REDUCING MUTATION LOAD THROUGH SELECTION ON MALES

Healthy males are likely to have higher mating success than unhealthy males because of differential expression of condition‐dependent traits such as mate searching intensity, fighting ability, display vigor, and some types of exaggerated morphological characters. We therefore expect that most new mutations that are deleterious for overall fitness may also be deleterious for male mating success. From this perspective, sexual selection is not limited to influencing those genes directly involved in exaggerated morphological traits but rather affects most, if not all, genes in the genome. If true, sexual selection can be an important force acting to reduce the frequency of deleterious mutations and, as a result, mutation load. We review the literature and find various forms of indirect evidence that sexual selection helps to eliminate deleterious mutations. However, direct evidence is scant, and there are almost no data available to address a key issue: is selection in males stronger than selection in females? In addition, the total effect of sexual selection on mutation load is complicated by possible increases in mutation rate that may be attributable to sexual selection. Finally, sexual selection affects population fitness not only through mutation load but also through sexual conflict, making it difficult to empirically measure how sexual selection affects load. Several lines of enquiry are suggested to better fill large gaps in our understanding of sexual selection and its effect on genetic load.     

Deleterious Mutations, Apparent Stabilizing Selection and the Maintenance of Quantitative Variation

Apparent stabilizing selection on a quantitative trait that is not causally connected to fitness can result from the pleiotropic effects of unconditionally deleterious mutations, because as N. Barton noted, "...individuals with extreme values of the trait will tend to carry more deleterious alleles...." We use a simple model to investigate the dependence of this apparent selection on the genomic deleterious mutation rate, U; the equilibrium distribution of K, the number of deleterious mutations per genome; and the parameters describing directional selection against deleterious mutations. Unlike previous analyses, we allow for epistatic selection against deleterious alleles. For various selection functions and realistic parameter values, the distribution of K, the distribution of breeding values for a pleiotropically affected trait, and the apparent stabilizing selection function are all nearly Gaussian. The additive genetic variance for the quantitative trait is kQa2, where k is the average number of deleterious mutations per genome, Q is the proportion of deleterious mutations that affect the trait, and a2 is the variance of pleiotropic effects for individual mutations that do affect the trait. In contrast, when the trait is measured in units of its additive standard deviation, the apparent fitness function is essentially independent of Q and a2; and beta, the intensity of selection, measured as the ratio of additive genetic variance to the "variance" of the fitness curve, is very close to s = U/k, the selection coefficient against individual deleterious mutations at equilibrium. Therefore, this model predicts appreciable apparent stabilizing selection if s exceeds about 0.03, which is consistent with various data. However, the model also predicts that beta must equal Vm/VG, the ratio of new additive variance for the trait introduced each generation by mutation to the standing additive variance. Most, although not all, estimates of this ratio imply apparent stabilizing selection weaker than generally observed. A qualitative argument suggests that even when direct selection is responsible for most of the selection observed on a character, it may be essentially irrelevant to the maintenance of variation for the character by mutation-selection balance. Simple experiments can indicate the fraction of observed stabilizing selection attributable to the pleiotropic effects of deleterious mutations.     

A Pleiotropic Nonadditive Model of Variation in Quantitative Traits

A model of mutation-selection-drift balance incorporating pleiotropic and dominance effects of new mutations on quantitative traits and fitness is investigated and used to predict the amount and nature of genetic variation maintained in segregating populations. The model is based on recent information on the joint distribution of mutant effects on bristle traits and fitness in Drosophila melanogaster from experiments on the accumulation of spontaneous and P element-induced mutations. These experiments suggest a leptokurtic distribution of effects with an intermediate correlation between effects on the trait and fitness. Mutants of large effect tend to be partially recessive while those with smaller effect are on average additive, but apparently with very variable gene action. The model is parameterized with two different sets of information derived from P element insertion and spontaneous mutation data, though the latter are not fully known. They differ in the number of mutations per generation which is assumed to affect the trait. Predictions of the variance maintained for bristle number assuming parameters derived from effects of P element insertions, in which the proportion of mutations with an effect on the trait is small, fit reasonably well with experimental observations. The equilibrium genetic variance is nearly independent of the degree of dominance of new mutations. Heritabilities of between 0.4 and 0.6 are predicted with population sizes from 10(4) to 10(6), and most of the variance for the metric trait in segregating populations is due to a small proportion of mutations (about 1% of the total number) with neutral or nearly neutral effects on fitness and intermediate effects on the trait (0.1-0.5σ(P)). Much of the genetic variance is contributed by recessive or partially recessive mutants, but only a small proportion (about 10%) of the genetic variance is dominance variance. The amount of apparent selection on the trait itself generated by the model is very small. If a model is assumed in which all mutation events have an effect on the quantitative trait, the majority of the genetic variance is contributed by deleterious mutations with tiny effects on the trait. If such a model is assumed for viability, the heritability is about 0.1, independent of the population size.     

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, “marking” one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.     

Evolution of microbial genomes: sequence acquisition and loss

We present models describing the acquisition and deletion of novel sequences in populations of microorganisms. We infer that most novel sequences are neutral. Thus, sequence duplications and gene transfer between organisms sharing the same environment are rarely expected to generate adaptive functions. Two classes of models are considered: (1) a homogeneous population with constant size, and (2) an island model in which the population is subdivided into patches that are in contact through slow migration. Distributions of gene frequencies are derived in a Moran model with overlapping generations. We find that novel, neutral or near-neutral coding sequences in microorganisms will not be fixed globally because they offer large target sizes for mutations and because the populations are so large. At most, such genes may have a transient presence in only a small fraction of the population. Consequently, a microbial population is expected to have a very large diversity of transient neutral gene content. Only sequences that are under strong selection, globally or in individual patches, can be expected to persist. We suggest that genome size is maintained in microorganisms by a quasi-steady state mechanism in which random fluctuations in the effective acquisition and deletion rates result in genome sizes that vary from patch to patch. We assign the genomic identity of a global population to those genes that are required for the participation of patches in the genetic sweeps that maintain the genomic coherence of the population. In contrast, we stress the influence of sequence loss on the isolation and the divergence (speciation) of novel patches from a global population.     

Role of selection in fixation of gene duplications

New genes commonly appear through complete or partial duplications of pre-existing genes. Duplications of long DNA segments are constantly produced by rare mutations, may become fixed in a population by selection or random drift, and are subject to divergent evolution of the paralogous sequences after fixation, although gene conversion can impede this process. New data shed some light on each of these processes. Mutations which involve duplications can occur through at least two different mechanisms, backward strand slippage during DNA replication and unequal crossing-over. The background rate of duplication of a complete gene in humans is 10(-9)-10(-10) per generation, although many genes located within hot-spots of large-scale mutation are duplicated much more often. Many gene duplications affect fitness strongly, and are responsible, through gene dosage effects, for a number of genetic diseases. However, high levels of intrapopulation polymorphism caused by presence or absence of long, gene-containing DNA segments imply that some duplications are not under strong selection. The polymorphism to fixation ratios appear to be approximately the same for gene duplications and for presumably selectively neutral nucleotide substitutions, which, according to the McDonald-Kreitman test, is consistent with selective neutrality of duplications. However, this pattern can also be due to negative selection against most of segregating duplications and positive selection for at least some duplications which become fixed. Patterns in post-fixation evolution of duplicated genes do not easily reveal the causes of fixations. Many gene duplications which became fixed recently in a variety of organisms were positively selected because the increased expression of the corresponding genes was beneficial. The effects of gene dosage provide a unified framework for studying all phases of the life history of a gene duplication. Application of well-known methods of evolutionary genetics to accumulating data on new, polymorphic, and fixed duplication will enhance our understanding of the role of natural selection in the evolution by gene duplication.