Role of gene duplication in evolution

It is now known that many multigene and supergene families exist in eukaryote genomes: multigene families with uniform copy members like genes for ribosomal RNA, those with variable members like immunoglobulin genes, and supergene families such as those for various growth factor and hormone receptors. Many such examples indicate that gene duplication and subsequent differentiation are extremely important for organismal evolution. In particular, gene duplication could well have been the primary mechanism for the evolution of complexity in higher organisms. Population genetic models for the origin of gene families with diverse functions are presented, in which natural selection favors those genomes with more useful mutants in duplicated genes. Since any gene has a certain probability of degenerating by mutation, success versus failure in acquiring a new gene by duplication may be expressed as the ratio of probabilities of spreading of useful versus detrimental mutations in redundant gene copies. Also examined are the effects of gene duplication on evolution by compensatory advantageous mutations. Results of the analyses show that both natural selection and random drift are important for the origin of gene families. In addition, interaction between molecular mechanisms such as unequal crossing-over and gene conversion, and selection or drift is found to have a large effect on evolution by gene duplication.     

Evolutionary mechanisms underlying the functional divergence of duplicate genes involved in vertebrates' circadian rhythm pathway

Circadian rhythms, that are governed physiologically and behaviorally by endogenous clock, have been described in many species. Living organisms use this endogenous circadian clock to anticipate environmental transitions, perform activities at biologically advantageous times during the day, and undergo characteristic seasonal responses. Gene duplication is one of the most important mechanisms in the evolution of gene diversity. After duplication, one or both of duplicates can accumulate amino acid changes, thereby promoting functional divergence through the action of natural selection. The circadian system, like many other multigene families, has undergone this genetic revolution, and so circadian genes that are found in single copies in insects are duplicated in vertebrates. We analyzed six groups of genes involved in vertebrates' circadian rhythm pathway to find signatures of molecular evolutionary processes such as gene duplication, natural selection, recombination, and functional divergence. The obtained results, then, were used to determine what evolutionary forces have influenced the fates of duplicated genes of each group. We showed in this research that recombination has not been widespread during the evolution of circadian genes and that purifying selection has been the prominent natural pressure operating on circadian genes. We also showed that the evolution of circadian genes has been depended on gene duplication and functional divergence. Finally, we put forward models best describing the evolutionary fates of circadian duplicates.     

Multigene families and the evolution of complexity

Summary Higher organisms are complex, and their developmental processes are controlled by the sequential expression of genes that often form multigene families. Facts are surveyed on how functional diversity of genes is related to duplication of genes or segments of genes, by emphasizing that diversity is often enhanced by alternate splicing and proteolytic cleavage involving duplicated genes or gene segments. Analyses of a population genetics model for the origin of gene families suggest that positive Darwinian selection is needed for acquiring gene families with desirable functions. Based on these considerations, examples that show acceleration of amino acid substitution relative to synonymous change during evolutionary processes are surveyed. Some of such examples strongly suggest that positive selection has worked. In other cases it is difficult to judge whether or not acceleration is caused by positive Darwinian selection. As a general pattern, acceleration of amino acid substitution is often found to be related to gene duplication. It is thought that complexity and diversity of gene function have been advantageous in the long evolutionary course of higher organisms.     

FURTHER SIMULATION STUDIES ON EVOLUTION BY GENE DUPLICATION

In order to understand the origin of multigene families, Monte Carlo simulations were performed to see how a genetic system evolves under unequal crossing-over, mutation, random genetic drift and natural selection, starting from a single gene copy. Both haploid and diploid models were examined. Beneficial, neutral, and detrimental mutations were incorporated, and “positive” selection favors those chromosomes (haploid) or individuals (diploid) with more beneficial mutations than others. The same model for haploids was previously investigated with special reference to the evolution of gene organization, and the ratio of the numbers of beneficial genes to pseudogenes was found to be a rough indicator of the relative strengths of positive and negative (against deleterious alleles) natural selection (Ohta, 1987b). In the present paper, the evolution of gene organization and of sequence divergence among genes in the multigene family is examined. It is shown that positive selection accelerates the accumulation of arrays containing different beneficial mutations, but that total divergence including both neutral and beneficial mutations is not very sensitive to positive selection, under this model. The proportion of beneficial mutations in the total mutations accumulated is a better indicator of positive selection than is the total divergence. It is pointed out that various observed examples in which amino-acid substitutions are accelerated, as compared with synonymous substitutions in duplicated genes (Li, 1985), may reflect the effect of selection similar to the present scheme. The diploid model is shown to be more efficient for accumulating beneficial mutations in duplicated genes than the haploid one, and the relevance of this finding to the advantage of sexual reproduction is discussed.     

Adaptive evolution of animal toxin multigene families

Animal toxins comprise a diverse array of proteins that have a variety of biochemical and pharmacological functions. A large number of animal toxins are encoded by multigene families. From studies of several toxin multigene families at the gene level the picture is emerging that most have been functionally diversified by gene duplication and adaptive evolution. The number of pharmacological activities in most toxin multigene families results from their adaptive evolution. The molecular evolution of animal toxins has been analysed in some multigene families, at both the intraspecies and interspecies levels. In most toxin multigene families, the rate of non-synonymous to synonymous substitutions (dN/dS) is higher than one. Thus natural selection has acted to diversify coding sequences and consequently the toxin functions. The selection pressure for the rapid adaptive evolution of animal toxins is the need for quick immobilization of the prey in classical predator and prey interactions. Currently available evidence for adaptive evolution in animal toxin multigene families will be considered in this review.     

Domain loss facilitates accelerated evolution and neofunctionalization of duplicate snake venom metalloproteinase toxin genes

Gene duplication is a key mechanism for the adaptive evolution and neofunctionalization of gene families. Large multigene families often exhibit complex evolutionary histories as a result of frequent gene duplication acting in concordance with positive selection pressures. Alterations in the domain structure of genes, causing changes in the molecular scaffold of proteins, can also result in a complex evolutionary history and has been observed in functionally diverse multigene toxin families. Here, we investigate the role alterations in domain structure have on the tempo of evolution and neofunctionalization of multigene families using the snake venom metalloproteinases (SVMPs) as a model system. Our results reveal that the evolutionary history of viperid (Serpentes: Viperidae) SVMPs is repeatedly punctuated by domain loss, with the single loss of the cysteine-rich domain, facilitating the formation of P-II class SVMPs, occurring prior to the convergent loss of the disintegrin domain to form multiple P-I SVMP structures. Notably, the majority of phylogenetic branches where domain loss was inferred to have occurred exhibited highly significant evidence of positive selection in surface-exposed amino acid residues, resulting in the neofunctionalization of P-II and P-I SVMP classes. These results provide a valuable insight into the mechanisms by which complex gene families evolve and detail how the loss of domain structures can catalyze the accelerated evolution of novel gene paralogues. The ensuing generation of differing molecular scaffolds encoded by the same multigene family facilitates gene neofunctionalization while presenting an evolutionary advantage through the retention of multiple genes capable of encoding functionally distinct proteins.     

Selectionism and neutralism in molecular evolution

Charles Darwin proposed that evolution occurs primarily by natural selection, but this view has been controversial from the beginning. Two of the major opposing views have been mutationism and neutralism. Early molecular studies suggested that most amino acid substitutions in proteins are neutral or nearly neutral and the functional change of proteins occurs by a few key amino acid substitutions. This suggestion generated an intense controversy over selectionism and neutralism. This controversy is partially caused by Kimura's definition of neutrality, which was too strict (|2Ns|< or =1). If we define neutral mutations as the mutations that do not change the function of gene products appreciably, many controversies disappear because slightly deleterious and slightly advantageous mutations are engulfed by neutral mutations. The ratio of the rate of nonsynonymous nucleotide substitution to that of synonymous substitution is a useful quantity to study positive Darwinian selection operating at highly variable genetic loci, but it does not necessarily detect adaptively important codons. Previously, multigene families were thought to evolve following the model of concerted evolution, but new evidence indicates that most of them evolve by a birth-and-death process of duplicate genes. It is now clear that most phenotypic characters or genetic systems such as the adaptive immune system in vertebrates are controlled by the interaction of a number of multigene families, which are often evolutionarily related and are subject to birth-and-death evolution. Therefore, it is important to study the mechanisms of gene family interaction for understanding phenotypic evolution. Because gene duplication occurs more or less at random, phenotypic evolution contains some fortuitous elements, though the environmental factors also play an important role. The randomness of phenotypic evolution is qualitatively different from allele frequency changes by random genetic drift. However, there is some similarity between phenotypic and molecular evolution with respect to functional or environmental constraints and evolutionary rate. It appears that mutation (including gene duplication and other DNA changes) is the driving force of evolution at both the genic and the phenotypic levels.     

Reactive site polymorphism in the murine protease inhibitor gene family is delineated using a modification of the PCR reaction (PCR + 1).

Murine protease inhibitor (alpha 1-PI) proteins are encoded by a multigene family which has undergone recent duplication. It has been suggested that the evolution of diversity within this gene family may be driven by unusual selection for novel function at the reactive site of the duplicated members (1,2,3). In an attempt to use polymerase chain reaction (PCR) to generate and sequence clones spanning the polymorphic reactive site region, a PCR artifact was identified and determined to result from heteroduplex formation during the co-amplification of the related sequences in this multigene system. This artifact results in sequences which are combinatorial mosaics of the template sequences. We present a simple and general method (PCR + 1) for overcoming this artifact and demonstrate its application in delineating five distinct alpha 1-PI reactive site sequences in C57BL/6 mice, thus providing sequence information to generate gene-specific probes. The significance of the reactive site diversity in this protease inhibitor gene family is discussed as well as the general applications and limitations of the PCR + 1 technique.     

Sequence analysis of the ly49 cluster in C57BL/6 mice: a rapidly evolving multigene family in the immune system

The cytotoxic activity of murine natural killer cells is controlled in part through the action of genes belonging to the Ly49 family. Members of this multigene family are found in a region on mouse chromosome 6 termed the natural killer gene complex. Using data available through public databases, we performed sequence analysis of a 620-kb region in C57Bl/6 (B6) mice that contains the Ly49 genes. The contiguous genomic sequence has allowed us to describe the complete B6 Ly49 gene repertoire, which includes two recently described genes as well as three partial genes. We have shown that the genes in the cluster have evolved through a series of large duplication events involving units of one or more genes and we have attempted to characterize the nature of the duplication end points. Finally, we have used information regarding gene sequence relationships and insertion of repetitive elements to construct a model for the evolution of the gene cluster. Our study illustrates that the Ly49 cluster represents an example of a rapidly evolving gene family, and continued analysis of this region in other strains will undoubtedly provide further insight into mechanisms for generating genomic diversity.     

Simulating Evolution by Gene Duplication

By considering the recent finding that unequal crossing over and other molecular interactions are contributing to the evolution of multigene families, a model of the origin of repetitive genes was studied by Monte Carlo simulations. Starting from a single gene copy, how genetic systems evolve was examined under unequal crossing over, random drift and natural selection. Both beneficial and deteriorating mutations were incorporated, and the latter were assumed to occur ten times more frequently than the former. Positive natural selection favors those chromosomes with more beneficial mutations in redundant copies than others in the population, but accumulation of deteriorating mutations (pseudogenes) have no effect on fitness so long as there remains a functional gene. The results imply the following: Positive natural selection is needed in order to acquire gene families with new functions. Without it, too many pseudogenes accumulate before attaining a functional gene family. There is a large fluctuation in the outcome even if parameters are the same. When unequal crossing over occurs more frequently, the system evolves more rapidly. It was also shown, under realistic values of parameters, that the genetic load for acquiring a new gene is not as large as J.B.S. Haldane suggested, but not so small as in a model in which a system for selection started from already redundant genes.     

Population genetics theory of concerted evolution and its application to the immunoglobulin V gene tree

Summary The previous simple model for treating concerted evolution of multigene families has been revised to be compatible with various new observations on the immunoglobulin variable region family and other families. In the previous model, gene conversion and unequal crossing-over were considered, and it was assumed that genes are randomly arranged on the chromosome; neither subdivision nor correlation of gene identity and chromosomal distance were considered. Although this model satisfactorily explains the observed amino acid diversity within and between species, it fails to predict the very ancient branching of the mouse immunoglobulin heavy chain V-gene family. By incorporating subdivided structure and genetic correlation with chromosomal distance into the simple model, the data of divergence may be satisfactorily explained, as well as the rate of nucleotide substitution and the amino acid diversity. The rate at which a V-gene is duplicated or deleted by conversion or by unequal crossing-over is estimated by the new model to be on the order of 10^–6 per year. The model may be applicable to other multigene families, such as those coding for silkmoth chorion or mammalian kallikrein.Contribution no. 1560 from the National Institute of Genetics, Mishima, 411 Japan     

Theoretical population genetics of repeated genes forming a multigene family

The evolution of repeated genes forming a multigene family in a finite population is studied with special reference to the probability of gene identity, i.e., the identity probability of two gene units chosen from the gene family. This quantity is called clonality and is defined as the sum of squares of the frequencies of gene lineages in the family. The multigene family is undergoing continuous unequal somatic crossing over, ordinary interchromosomal crossing over, mutation and random frequency drift. Two measures of clonality are used: clonality within one chromosome and that between two different chromosomes. The equilibrium properties of the means, the variances and the covariance of the two measures of clonality are investigated by using the diffusion equation method under the assumption of constant number of gene units in the multigene family. Some models of natural selection based on clonality are considered. The possible significance of the variance and covariance of clonality among the chromosomes on the adaptive differentiation of gene families such as those producing antibodies is discussed.     

An integrated approach for the comparative analysis of a multigene family: The nicotianamine synthase genes of barley

Recent genomic projects reveal that about half of the gene repertoire in plant genomes is made up by multigene families. In this paper, a set of structural and phylogenetic analyses have been applied to compare the differently sized nicotianamine synthase (NAS) gene families in barley and rice. Nicotianamine acts as a chelator of iron and other heavy metals and plays a key role in uptake, phloem transport and cytoplasmic distribution of iron, challenging efforts for the breeding of iron-efficient crop plants. Nine barley NAS genes have been mapped, and co-linearity of flanking genes in barley and rice was determined. The combined analyses reveal that the NAS multigene family members in barley originated through at least one duplication event that occurred before the divergence of rice and barley. Additional duplications appear to have occurred within each of the species. Although we detected no evidence for positive selection of recently duplicated genes within species, codon-based tests revealed evidence for positive selection having contributed to the divergence of some amino acids. The integrated comparative and phylogenetic analysis improved our current view of NAS gene family evolution, might facilitate the functional characterization of individual members and is applicable to other multigene families. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Concerted evolution of the immunoglobulin VH gene family

With the aim of understanding the concerted evolution of the immunoglobulin VH multigene family, a phylogenetic tree for the DNA sequences of 16 mouse and five human germ line genes was constructed. This tree indicates that all genes in this family have undergone substantial evolutionary divergence. The most closely related genes so far identified in the mouse genome seem to have diverged about 6 million years (MY) ago, whereas the most distantly related genes diverged about 300 MY ago. This suggests that gene duplication caused by unequal crossing-over or gene conversion occurs very slowly in this gene family. The rate of occurrence of gene duplication in the VH gene family has been estimated to be 5 x 10(-7) per gene per year, which seems to be at least about 100 times lower than that for the rRNA gene family. This low rate of concerted evolution in the VH gene family helps retain intergenic genetic variability that in turn contributes to antibody diversity. Because of accumulation of destructive mutations, however, about one-third of the mouse and human VH genes seem to have become nonfunctional. Many of these pseudogenes have apparently originated recently, but some of them seem to have existed in the genome for more than 10 MY. The rate of nucleotide substitution for the complementarity-determining regions (CDRs) is as high as that of pseudogenes. This suggests that there is virtually no purifying selection operating in the CDRs and that germ line mutations are effectively used for generating antibody diversity.      

Genome-level evolution of resistance genes in Arabidopsis thaliana

Pathogen resistance genes represent some of the most abundant and diverse gene families found within plant genomes. However, evolutionary mechanisms generating resistance gene diversity at the genome level are not well understood. We used the complete Arabidopsis thaliana genome sequence to show that most duplication of individual NBS-LRR sequences occurs at close physical proximity to the parent sequence and generates clusters of closely related NBS-LRR sequences. Deploying the statistical strength of phylogeographic approaches and using chromosomal location as a proxy for spatial location, we show that apparent duplication of NBS-LRR genes to ectopic chromosomal locations is largely the consequence of segmental chromosome duplication and rearrangement, rather than the independent duplication of individual sequences. Although accounting for a smaller fraction of NBS-LRR gene duplications, segmental chromosome duplication and rearrangement events have a large impact on the evolution of this multigene family. Intergenic exchange is dramatically lower between NBS-LRR sequences located in different chromosome regions as compared to exchange between sequences within the same chromosome region. Consequently, once translocated to new chromosome locations, NBS-LRR gene copies have a greater likelihood of escaping intergenic exchange and adopting new functions than do gene copies located within the same chromosomal region. We propose an evolutionary model that relates processes of genome evolution to mechanisms of evolution for the large, diverse, NBS-LRR gene family.     

Extensive and continuous duplication facilitates rapid evolution and diversification of gene families

The origin of novel gene functions through gene duplication, mutation, and natural selection represents one of the mechanisms by which organisms diversify and one of the possible paths leading to adaptation. Nonetheless, the extent, role, and consequences of duplications in the origins of ecological adaptations, especially in the context of species interactions, remain unclear. To explore the evolution of a gene family that is likely linked to species associations, we investigated the evolutionary history of the A-superfamily of conotoxin genes of predatory marine cone snails (Conus species). Members of this gene family are expressed in the venoms of Conus species and are presumably involved in predator-prey associations because of their utility in prey capture. We recovered sequences of this gene family from genomic DNA of four closely related species of Conus and reconstructed the evolutionary history of these genes. Our study is the first to directly recover conotoxin genes from Conus genomes to investigate the evolution of conotoxin gene families. Our results revealed a phenomenon of rapid and continuous gene turnover that is coupled with heightened rates of evolution. This continuous duplication pattern has not been observed previously, and the rate of gene turnover is at least two times higher than estimates from other multigene families. Conotoxin genes are among the most rapidly evolving protein-coding genes in metazoans, a phenomenon that may be facilitated by extensive gene duplications and have driven changes in conotoxin functions through neofunctionalization. Together these mechanisms led to dramatically divergent arrangements of A-superfamily conotoxin genes among closely related species of Conus. Our findings suggest that extensive and continuous gene duplication facilitates rapid evolution and drastic divergence in venom compositions among species, processes that may be associated with evolutionary responses to predator-prey interactions.     

Birth-and-death evolution of the internalin multigene family in Listeria

The birth-and-death model of multigene family evolution describes patterns of gene origination, diversification and loss within multigene families. Since it was first developed in the 1990s, the model has been found to characterize a large number of eukaryotic multigene families. In this paper, we report for the first time a bacterial multigene family that undergoes birth-and-death evolution. By analyzing the evolutionary relationships among internalins, a relatively large and diverse family of genes associated with key virulence functions in Listeria, we demonstrate the importance of birth-and-death evolution in the diversification of this important bacterial pathogen. We also detected two instances of lateral gene transfer within the internalins, but the estimated frequency would have been much higher had it not been analyzed within the context of birth-and-death evolutionary dynamics and a phenomenon that we term "paralog-sorting", which involves the unequal transmittal of gene duplicates during or subsequent to the speciation process. As such, in addition to providing the first demonstration of birth-and-death evolution within a bacterial multigene family, our results indicate that the extent of lateral transfer in bacterial multigene families should be re-examined in the light of birth-and-death evolution.     

Family business: the multidrug-resistance related protein (MRP) ABC transporter genes in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Despite the completion of the sequencing of the entire genome of Arabidopsis thaliana (L.) Heynh., the exact determination of each single gene and its function remains an open question. This is especially true for multigene families. An approach that combines analysis of genomic structure, expression data and functional genomics to ascertain the role of the members of the multidrug-resistance-related protein ( MRP) gene family, a subfamily of the ATP-binding cassette (ABC) transporters from Arabidopsis is presented. We used cDNA sequencing and alignment-based re-annotation of genomic sequences to define the exact genic structure of all known AtMRP genes. Analysis of promoter regions suggested different induction conditions even for closely related genes. Expression analysis for the entire gene family confirmed these assumptions. Phylogenetic analysis and determination of segmental duplication in the regions of AtMRP genes revealed that the evolution of the extraordinarily high number of ABC transporter genes in plants cannot solely be explained by polyploidisation during the evolution of the Arabidopsis genome. Interestingly MRP genes from Oryza sativa L. (rice; OsMRP) show very similar genomic structures to those from Arabidopsis. Screening of large populations of T-DNA-mutagenised lines of A. thaliana resulted in the isolation of AtMRP insertion mutants. This work opens the way for the defined analysis of a multigene family of important membrane transporters whose broad variety of functions expands their traditional role as cellular detoxifiers.     

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high‐amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model‐averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.