The probability of preservation of a newly arisen gene duplicate.

Newly emerging data from genome sequencing projects suggest that gene duplication, often accompanied by genetic map changes, is a common and ongoing feature of all genomes. This raises the possibility that differential expansion/contraction of various genomic sequences may be just as important a mechanism of phenotypic evolution as changes at the nucleotide level. However, the population-genetic mechanisms responsible for the success vs. failure of newly arisen gene duplicates are poorly understood. We examine the influence of various aspects of gene structure, mutation rates, degree of linkage, and population size (N) on the joint fate of a newly arisen duplicate gene and its ancestral locus. Unless there is active selection against duplicate genes, the probability of permanent establishment of such genes is usually no less than 1/(4N) (half of the neutral expectation), and it can be orders of magnitude greater if neofunctionalizing mutations are common. The probability of a map change (reassignment of a key function of an ancestral locus to a new chromosomal location) induced by a newly arisen duplicate is also generally >1/(4N) for unlinked duplicates, suggesting that recurrent gene duplication and alternative silencing may be a common mechanism for generating microchromosomal rearrangements responsible for postreproductive isolating barriers among species. Relative to subfunctionalization, neofunctionalization is expected to become a progressively more important mechanism of duplicate-gene preservation in populations with increasing size. However, even in large populations, the probability of neofunctionalization scales only with the square of the selective advantage. Tight linkage also influences the probability of duplicate-gene preservation, increasing the probability of subfunctionalization but decreasing the probability of neofunctionalization.     

Preservation of duplicate genes by complementary, degenerative mutations

The origin of organismal complexity is generally thought to be tightly coupled to the evolution of new gene functions arising subsequent to gene duplication. Under the classical model for the evolution of duplicate genes, one member of the duplicated pair usually degenerates within a few million years by accumulating deleterious mutations, while the other duplicate retains the original function. This model further predicts that on rare occasions, one duplicate may acquire a new adaptive function, resulting in the preservation of both members of the pair, one with the new function and the other retaining the old. However, empirical data suggest that a much greater proportion of gene duplicates is preserved than predicted by the classical model. Here we present a new conceptual framework for understanding the evolution of duplicate genes that may help explain this conundrum. Focusing on the regulatory complexity of eukaryotic genes, we show how complementary degenerative mutations in different regulatory elements of duplicated genes can facilitate the preservation of both duplicates, thereby increasing long-term opportunities for the evolution of new gene functions. The duplication-degeneration-complementation (DDC) model predicts that (1) degenerative mutations in regulatory elements can increase rather than reduce the probability of duplicate gene preservation and (2) the usual mechanism of duplicate gene preservation is the partitioning of ancestral functions rather than the evolution of new functions. We present several examples (including analysis of a new engrailed gene in zebrafish) that appear to be consistent with the DDC model, and we suggest several analytical and experimental approaches for determining whether the complementary loss of gene subfunctions or the acquisition of novel functions are likely to be the primary mechanisms for the preservation of gene duplicates. For a newly duplicated paralog, survival depends on the outcome of the race between entropic decay and chance acquisition of an advantageous regulatory mutation.Sidow 1996(p. 717) On one hand, it may fix an advantageous allele giving it a slightly different, and selectable, function from its original copy. This initial fixation provides substantial protection against future fixation of null mutations, allowing additional mutations to accumulate that refine functional differentiation. Alternatively, a duplicate locus can instead first fix a null allele, becoming a pseudogene.Walsh 1995 (p. 426) Duplicated genes persist only if mutations create new and essential protein functions, an event that is predicted to occur rarely.Nadeau and Sankoff 1997 (p. 1259) Thus overall, with complex metazoans, the major mechanism for retention of ancient gene duplicates would appear to have been the acquisition of novel expression sites for developmental genes, with its accompanying opportunity for new gene roles underlying the progressive extension of development itself.Cooke et al. 1997 (p. 362)     

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.     

Preservation of duplicate genes by originalization

Neofunctionalization, subfunctionalization and increasing gene dosage were proposed to be the possible ways to explain duplicate-gene preservation in previous studies. However, in some natural populations, such as yeast Saccharomyces cerevisiae, a considerable proportion of the duplicate genes originated from ancient whole genomic duplication (WGD) is preserved till now, which cannot be sufficiently explained by these mechanisms. In this article, we present another possible way to explain this conundrum—originalization, by which duplicate genes are both preserved intact at a high frequency in the population under only purifying selection. With approximate equal rates of mutation at the two duplicated loci, analytical, numerical and simulation results consistently show that the mean time to nonfunctionalization for unlinked haploinsufficient gene duplication might become markedly prolonged, which results from originalization. These theoretical results imply that originalization might be an alternative effective and temporary way of preserving duplicate genes.     

A Diffusion Approach to Approximating Preservation Probabilities for Gene Duplicates

Consider a haploid population and, within its genome, a gene whose presence is vital for the survival of any individual. Each copy of this gene is subject to mutations which destroy its function. Suppose one member of the population somehow acquires a duplicate copy of the gene, where the duplicate is fully linked to the original gene’s locus. Preservation is said to occur if eventually the entire population consists of individuals descended from this one which initially carried the duplicate. The system is modelled by a finite state-space Markov process which in turn is approximated by a diffusion process, whence an explicit expression for the probability of preservation is derived. The event of preservation can be compared to the fixation of a selectively neutral gene variant initially present in a single individual, the probability of which is the reciprocal of the population size. For very weak mutation, this and the probability of preservation are equal, while as mutation becomes stronger, the preservation probability tends to double this reciprocal. This is in excellent agreement with simulation studies.     

Duplication of accelerated evolution and growth hormone gene in passerine birds

We report the discovery of a duplication of the growth hormone (GH) gene in a major group of birds, the passerines (Aves: Passeriformes). Phylogenetic analysis of 1.3-kb partial DNA sequences of GH genes for 24 species of passerines and numerous outgroups indicates that the duplication occurred in the ancestral lineage of extant passerines. Both duplicates and their open-reading frames are preserved throughout the passerine clade, and both duplicates are expressed in the zebra finch brain, suggesting that both are likely to be functional. The estimated rates of amino acid evolution are more than 10-fold higher in passerine GH genes than in those of their closest nonpasserine relatives. In addition, although the 84 codons sequenced are generally highly conserved for both passerines and nonpasserines, comparisons of the nonsynonymous/synonymous substitution ratios and the rate of predicted amino acid changes indicate that the 2 gene duplicates are evolving under different selective pressures and may be functionally divergent. The evidence of differential selection, coupled with the preservation of both gene copies in all major lineages since the origin of passerines, suggests that the duplication may be of adaptive significance, with possible implications for the explosive diversification of the passerine clade.     

The evolution of mutator genes in bacterial populations: the roles of environmental change and timing

Recent studies have found high frequencies of bacteria with increased genomic rates of mutation in both clinical and laboratory populations. These observations may seem surprising in light of earlier experimental and theoretical studies. Mutator genes (genes that elevate the genomic mutation rate) are likely to induce deleterious mutations and thus suffer an indirect selective disadvantage; at the same time, bacteria carrying them can increase in frequency only by generating beneficial mutations at other loci. When clones carrying mutator genes are rare, however, these beneficial mutations are far more likely to arise in members of the much larger nonmutator population. How then can mutators become prevalent? To address this question, we develop a model of the population dynamics of bacteria confronted with ever-changing environments. Using analytical and simulation procedures, we explore the process by which initially rare mutator alleles can rise in frequency. We demonstrate that subsequent to a shift in environmental conditions, there will be relatively long periods of time during which the mutator subpopulation can produce a beneficial mutation before the ancestral subpopulations are eliminated. If the beneficial mutation arises early enough, the overall frequency of mutators will climb to a point higher than when the process began. The probability of producing a subsequent beneficial mutation will then also increase. In this manner, mutators can increase in frequency over successive selective sweeps. We discuss the implications and predictions of these theoretical results in relation to antibiotic resistance and the evolution of mutation rates.     

Species concepts and the recognition of ancestors

Whether or not ancestral species can be recognised depends on the species concept adopted. A “metaspecies”; is a species that completely lacks autapomorphies, and which might (or might not) be ancestral to other species. Such taxa have been identified among both living and fossil organisms. Under the most commonly‐used species concepts (biological, evolutionary, phenetic, phylogenetic, ecological, recognition and cohesion), “metaspecies”; can be assumed to be ancestral. Even if the known members of a metaspecies are not ancestral to anything, parsimony dictates that the (as yet unknown) ancestral lineage is identical to the metaspecies and, under these species concepts, assignable to the same species. Only the cladistic and monophyletic species concepts would deny “metaspecies”; ancestral status, but these species concepts are problematical and have never been used by practising systematists.     

Assessment of linkage disequilibrium by the decay of haplotype sharing, with application to fine-scale genetic mapping.

Linkage disequilibrium (LD) is of great interest for gene mapping and the study of population history. We propose a multilocus model for LD, based on the decay of haplotype sharing (DHS). The DHS model is most appropriate when the LD in which one is interested is due to the introduction of a variant on an ancestral haplotype, with recombinations in succeeding generations resulting in preservation of only a small region of the ancestral haplotype around the variant. This is generally the scenario of interest for gene mapping by LD. The DHS parameter is a measure of LD that can be interpreted as the expected genetic distance to which the ancestral haplotype is preserved, or, equivalently, 1/(time in generations to the ancestral haplotype). The method allows for multiple origins of alleles and for mutations, and it takes into account missing observations and ambiguities in haplotype determination, via a hidden Markov model. Whereas most commonly used measures of LD apply to pairs of loci, the DHS measure is designed for application to the densely mapped haplotype data that are increasingly available. The DHS method explicitly models the dependence among multiple tightly linked loci on a chromosome. When the assumptions about population structure are sufficiently tractable, the estimate of LD is obtained by maximum likelihood. For more-complicated models of population history, we find means and covariances based on the model and solve a quasi-score estimating equation. Simulations show that this approach works extremely well both for estimation of LD and for fine mapping. We apply the DHS method to published data sets for cystic fibrosis and progressive myoclonus epilepsy.     

Molecular trajectories leading to the alternative fates of duplicate genes

Gene duplication generates extra gene copies in which mutations can accumulate without risking the function of pre-existing genes. Such mutations modify duplicates and contribute to evolutionary novelties. However, the vast majority of duplicates appear to be short-lived and experience duplicate silencing within a few million years. Little is known about the molecular mechanisms leading to these alternative fates. Here we delineate differing molecular trajectories of a relatively recent duplication event between humans and chimpanzees by investigating molecular properties of a single duplicate: DNA sequences, gene expression and promoter activities. The inverted duplication of the Glutathione S-transferase Theta 2 (GSTT2) gene had occurred at least 7 million years ago in the common ancestor of African great apes and is preserved in chimpanzees (Pan troglodytes), whereas a deletion polymorphism is prevalent in humans. The alternative fates are associated with expression divergence between these species, and reduced expression in humans is regulated by silencing mutations that have been propagated between duplicates by gene conversion. In contrast, selective constraint preserved duplicate divergence in chimpanzees. The difference in evolutionary processes left a unique DNA footprint in which dying duplicates are significantly more similar to each other (99.4%) than preserved ones. Such molecular trajectories could provide insights for the mechanisms underlying duplicate life and death in extant genomes.     

Inference from gene trees in a subdivided population

This paper studies gene trees in subdivided populations which are constructed as perfect phylogenies from the pattern of mutations in a sample of DNA sequences and presents a new recursion for the probability distribution of such gene trees. The underlying evolutionary model is the coalescent process in a subdivided population. The infinitely-many-sites model of mutation is assumed. Ancestral inference questions that are discussed are maximum likelihood estimation of migration and mutation rates; detection of population growth by likelihood techniques; determining the distribution of the time to the most recent common ancestor of a sample of sequences; determining the distribution of the age of the mutations on the gene tree; determining in which subpopulation the most recent common ancestor of all the sequences was; determining subpopulation ancestors, where they were, and times to them; and determining in which subpopulations mutations occurred. A computational technique of Griffiths and Tavaré used is a computer intensive Markov chain simulation, which simulates gene trees conditional on their topology implied by the mutation pattern in the sample of DNA sequences. The software GENETREE, which implements these ancestral inference techniques, is available.     

An ancient frame-shifting event in the highly conserved KPNA gene family has undergone extensive compensation by natural selection in vertebrates

One of the prevailing arguments in evolutionary theory is that the duplicates of genes can acquire novel functionality. This is because only one of the paralogs need maintain the ancestral function, leaving room for natural experimentation due to a respite in purifying selection. Although many duplicates can subsequently become disabled by nullifying mutations, a few may also go on to diverge along a novel evolutionary trajectory. Here, evidence is provided that demonstrates how this scenario may not always be true. Rather, in the case of the highly conserved KPNA importin family, an initial relaxation in selection induced a frameshift that was later suppressed and heavily compensated for as part of a reparative and optimizing process. Despite a resulting divergence, there remains a distinct preservation of both sequence and functionality among the paralogs. This would indicate that duplicates can be retained by selection for reasons related to their redundant functionality. It also shows that, even when positive selection is inferred in duplicate genes, this may be of a compensatory nature rather than one representing any biochemical innovation. Generally, this development would perhaps be a more common outcome for gene duplication than is currently maintained.     

Analysing gene function after duplication

After gene duplication, mutations cause the gene copies to diverge. The classical model predicts that these mutations will generally lead to the loss of function of one gene copy; rarely, new functions will be created and both duplicate genes are conserved. In contrast, under the subfunctionalization model both duplicates are preserved due to the partition of different functions between the duplicates. A recent study provides support for the subfunctionalization model, identifying several expressed gene duplicates common to humans and mice that contain regions conserved in one duplicate but variable in the other (and vice versa). We discuss both the methodology used in this study and also how gene phylogeny may lead to additional evidence for the importance of subfunctionalization in the evolution of new genes.     

Predominant gain of promoter TATA box after gene duplication associated with stress responses

TATA box, the core promoter element, exists in a broad range of eukaryotes, and the expression of TATA-containing genes usually responds to various environmental stresses. Hence, the evolution of TATA-box in duplicate genes may provide some clues for the interrelationship among environmental stress, expression differentiation, and duplicate gene preservation. In the present study, we observed that the TATA box is significantly overrepresented in duplicate genes compared with singletons in human, worm, Arabidopsis, and yeast genomes. We then conducted an extensive functional genomic analysis to investigate the evolution of TATA box along over 700 yeast gene family phylogenies. After reconstructing the ancestral TATA-box states (presence or absence), we found that significantly higher numbers of TATA box gain events than loss events had occurred after yeast gene duplications-the overall gain-loss ratio is about 3-4 to 1. Interestingly, these TATA-gain duplicate genes on average have experienced greater expression divergence from the ancestral expression states than their most closely related TATA-less duplicate partners, but only under environmental stress conditions (asymmetric evolution); indeed, under normal physiological conditions, they have similar expression divergence (symmetric evolution). Moreover, we showed that TATA-gain duplicates are enriched in stress-associated functional categories but that is not the case for TATA-ancestral duplicates (those inherited from their ancestors prior to duplication). Together, we conclude that after the gene duplication, gain of the TATA box in duplicate promoters may have played an important role in yeast duplicate preservation by accelerating expression divergence that may facilitate the adaptive evolution of the organism in response to environmental changes.     

Fixation biases affecting human SNPs

Under neutrality all classes of mutation have an equal probability of becoming fixed in a population. In this article, we describe our analysis of the frequency distributions of >5000 human SNPs and provide evident of biases in the process of fixation of certain classes of point mutation that are most likely to be attributable to biased gene conversion. The results indicate an increased fixation probability of mutations that result in the incorporation of a GC base pair. Furthermore, in transcribed regions this process exhibits strand asymmetry, and is biased towards preserving a G base on the coding strand. Biased gene conversion has the potential to explain both existence of isochores and the compositional asymmetry in mammalian transcribed regions.     

GPCR genes are preferentially retained after whole genome duplication

One of the most interesting questions in biology is whether certain pathways have been favored during evolution, and if so, what properties could cause such a preference. Due to the lack of experimental evidence, whether select gene families have been preferentially retained over time after duplication in metazoan organisms remains unclear. Here, by syntenic mapping of nonchemosensory G protein-coupled receptor genes (nGPCRs which represent half the receptome for transmembrane signaling) in the vertebrate genomes, we found that, as opposed to the 8-15% retention rate for whole genome duplication (WGD)-derived gene duplicates in the entire genome of pufferfish, greater than 27.8% of WGD-derived nGPCRs which interact with a nonpeptide ligand were retained after WGD in pufferfish Tetraodon nigroviridis. In addition, we show that concurrent duplication of cognate ligand genes by WGD could impose selection of nGPCRs that interact with a polypeptide ligand. Against less than 2.25% probability for parallel retention of a pair of WGD-derived ligands and a pair of cognate receptor duplicates, we found a more than 8.9% retention of WGD-derived ligand-nGPCR pairs--threefold greater than one would surmise. These results demonstrate that gene retention is not uniform after WGD in vertebrates, and suggest a Darwinian selection of GPCR-mediated intercellular communication in metazoan organisms.     

Toward localization of the Werner syndrome gene by linkage disequilibrium and ancestral haplotyping: lessons learned from analysis of 35 chromosome 8p11.1-21.1 markers.

Werner syndrome (WS) is an autosomal recessive disorder characterized by premature onset of a number of age-related diseases. The gene for WS, WRN, has been mapped to the 8p 11.1-21.1 region with further localization through linkage disequilibrium mapping. Here we present the results of linkage disequilibrium and ancestral haplotype analyses of 35 markers to further refine the location of WRN. We identified an interval in this region in which 14 of 18 markers tested show significant evidence of linkage disequilibrium in at least one of the two populations tested. Analysis of extended and partial haplotypes covering 21 of the markers studied supports the existence of both obligate and probable ancestral recombinant events which localize WRN almost certainly to the interval between D8S2196 and D8S2186, and most likely to the narrower interval between D8S2168 and D8S2186. These haplotype analyses also suggest that there are multiple WRN mutations in each of the two populations under study. We also present a comparison of approaches to performing disequilibrium tests with multiallelic markers, and show that some commonly used approximations for such tests perform poorly in comparison to exact probability tests. Finally, we discuss some of the difficulties introduced by the high mutation rate at microsatellite markers which influence our ability to use ancestral haplotype analysis to localize disease genes.     

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.