Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution

It has been proposed that two events of duplication of the entire genome occurred early in vertebrate history (2R hypothesis). Several phylogenetic studies with a few gene families (mostly Hox genes and proteins from the MHC) have tried to confirm these polyploidization events. However, data from a single locus cannot explain the evolutionary history of a complete genome. To study this 2R hypothesis, we have taken advantage of the phylogenetic position of the lamprey to study the history of gene duplications in vertebrates. We selected most gene families that contain several paralogous genes in vertebrates and for which lamprey genes and an out-group are known in databases. In addition, we isolated members of the nuclear receptor superfamily in lamprey. Hagfish genes were also analyzed and found to confirm the lamprey gene analysis. Consistent with the 2R hypothesis, the phylogenetic analysis of 33 selected gene families, dispersed through the whole genome, revealed that one period of gene duplication arose before the lamprey-gnathostome split and this was followed by a second period of gene duplication after the lamprey-gnathostome split. Nevertheless, our analysis suggests that numerous gene losses and other gene-genome duplications occurred during the evolution of the vertebrate genomes. Thus, the complexity of all the paralogy groups present in vertebrates should be explained by the contribution of genome duplications (2R hypothesis), extra gene duplications, and gene losses.     

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.     

Diversity in membrane binding sites of ankyrins. Brain ankyrin, erythrocyte ankyrin, and processed erythrocyte ankyrin associate with distinct sites in kidney microsomes

This report presents evidence for diversity in membrane binding sites between three forms of ankyrin: brain ankyrin, erythrocyte ankyrin, and a variant of erythrocyte ankyrin (protein 2.2) present in circulating human erythrocytes that is missing a regulatory domain. These ankyrins were compared with respect to binding to kidney microsomes and exhibited the following behavior. 1) Brain and erythrocyte ankyrin each bind to distinct sites. 2) Protein 2.2 is an activated ankyrin that binds to all of the sites accessible to both brain and erythrocyte ankyrin and, in addition, associates with its own specialized sites. 3) The specificity of these membrane sites for various ankyrins is not absolute but reflects 2.5-10-fold differences in relative affinities. Further evidence that binding sites of different ankyrins share some common features is that the cytoplasmic domain of the erythrocyte anion transporter associates with all three ankyrins and displaces binding of the ankyrin variants to kidney membranes. The differences between erythrocyte and brain ankyrins in association with kidney membranes are likely to have physiological relevance to kidney because immunologically related isoforms of ankyrin are expressed in this tissue: erythroid ankyrin which is restricted to the basolateral domains of two cell types and a brain-related ankyrin expressed in all cells and present on apical as well as basolateral membrane surfaces. An unanticipated observation was the discovery of a membrane-associated ankyrin protease in kidney that is specific for erythrocyte ankyrin and may selectively activate the erythroid isoform of ankyrin. The variety of binding sites within this group of ankyrin proteins supports the idea that ankyrins are capable of linking a number of different membrane proteins to the spectrin-actin skeleton.     

Positive Selection and Functional Divergence After Melanopsin Gene Duplication

A newly discovered melanopsin gene (Opn4) encodes a member of the opsins, melanopsin. Two melanopsin genes, mammalian-like Opn4m and Xenopus-like Opn4x, have been described in nonmammalian vertebrates, but the underlying evolutionary mechanisms behind the duplication of melanopsin genes remain unclear. We conducted a comprehensive evolutionary analysis within a phylogenetic framework. In our phylogenetic tree, the duplication of Opn4m and Opn4x probably occurred prior to the emergence of vertebrates, and subsequently Opn4x disappeared in the lineages leading to mammalian species. Evolutionary analyses show strong purifying selection during melanopsin evolution. We also provide evidence that Opn4x underwent positive selection after the early gene duplication events. It has been indicated that functional divergence and altered functional constraints occurred between Opn4m and Opn4x duplicates with the identification of positively selected amino acids. Our findings highlight the evolutionary malleability in vertebrate melanopsin genes and provide a genetic basis for comparative studies of functional properties of these two melanopsins.     

Molecular evolution of ankyrin: gain of function in vertebrates by acquisition of an obscurin/titin-binding-related domain

Ankyrins form a family of modular adaptor proteins that link between integral membrane proteins and the cytoskeleton. They evolved within the Metazoa as an adaptation for organizing membrane microstructure and directing membrane traffic. Molecular cloning has identified one Caenorhabditis elegans (unc-44), two Drosophila (Dank1, Dank2), and three mammalian (Ank1, Ank2, Ank3) genes. We have previously identified a 76-amino acid (aa) alternatively spliced sequence that is present in muscle polypeptides encoded by the rat Ank3 gene. A closely related sequence in a muscle Ank1 product binds the cytoskeletal muscle proteins obscurin and titin. This obscurin/titin-binding-related domain (OTBD) contains repeated modules of 18 aa: three are encoded by Ank1 and Ank2, two by Ank3; this pattern is conserved throughout vertebrate ankyrin genes. The C. elegans ankyrin, UNC-44, contains one 18-aa module as does the ankyrin gene in the urochordate Ciona intestinalis, but the insect ankyrins contain none. Our data indicate that an ancestral ankyrin acquired an 18-aa module which was preserved in the Ecdysozoa/deuterostome divide, but it was subsequently lost from arthropods. Successive duplications of the module led to a gain of function in vertebrates as it acquired obscurin/titin-binding activity. We suggest that the OTBD represents an adaptation of the cytoskeleton that confers muscle cells with resilience to the forces associated with vertebrate life.     

Novel interactions of ankyrins-G at the costameres: the muscle-specific Obscurin/Titin-Binding-related Domain (OTBD) binds plectin and filamin C

Ankyrins, the adapters of the spectrin skeleton, are involved in local accumulation and stabilization of integral proteins to the appropriate membrane domains. In striated muscle, tissue-dependent alternative splicing generates unique Ank3 gene products (ankyrins-G); they share the Obscurin/Titin-Binding-related Domain (OTBD), a muscle-specific insert of the C-terminal domain which is highly conserved among ankyrin genes, and binds obscurin and titin to Ank1 gene products. We previously proposed that OTBD sequences constitute a novel domain of protein–protein interactions which confers ankyrins with specific cellular functions in muscle. Here we searched for muscle proteins binding to ankyrin-G OTBD by yeast two hybrid assay, and we found plectin and filamin C, two organizing elements of the cytoskeleton with essential roles in myogenesis, muscle cell cytoarchitecture, and muscle disease. The three proteins coimmunoprecipitate from skeletal muscle extracts and colocalize at costameres in adult muscle fibers. During in vitro myogenesis, muscle ankyrins-G are first expressed in postmitotic myocytes undergoing fusion to myotubes. In western blots of subcellular fractions from C2C12 cells, the majority of muscle ankyrins-G appear associated with membrane compartments. Occasional but not extensive co-localization at nascent costameres suggested that ankyrin-G interactions with plectin and filamin C are not involved in costamere assembly; they would rather reinforce stability and/or modulate molecular interactions in sarcolemma microdomains by establishing novel links between muscle-specific ankyrins-G and the two costameric dystrophin-associated glycoprotein and integrin-based protein complexes. These results report the first protein–protein interactions involving the ankyrin-G OTBD domain and support the hypothesis that OTBD sequences confer ankyrins with a gain of function in vertebrates, bringing further consolidation and resilience of the linkage between sarcomeres and sarcolemma.     

基因倍增和脊椎动物起源

有机体基因复制导致基因复杂性增加及其和脊椎动物起源的关系已经成为进化生物学研究的热点。20世纪70年代由Ohno提出后经Holland等修正的原始脊索动物经两轮基因组复制产生脊椎动物的假设目前已被广泛接受。脊椎动物起源和进化过程中发生过两轮基因组复制的主要证据有三点：（1）据估计脊椎动物基因组内编码基因数目大约相当于果蝇、海鞘等无脊椎动物的4倍;原口动物如果蝇和后口动物如头索动物文昌鱼的基因组大都只有单拷贝的基因,而脊椎动物的基因组则通常有4个同属于一个家族的基因。（2）无脊椎动物如节肢动物、海胆和头索动物文昌鱼都只有一个Hox基因簇,而脊椎动物除鱼类外,有7个具有Hox基因簇,其余都具有4个Hox基因簇。（3）基因作图证明,不但在鱼类和哺乳动物染色体广大片段上基因顺序相似,而且有证据显示哺乳动物基因组不同染色体之间存在相似性。据认为第一次基因倍增发生在脊椎动物与头索动物分开之后,第二次基因倍增发生在有颌类脊椎动物和无颌类脊椎动物分开以后。但是,基因是逐个发生倍增,还是通过基因组内某些DNA片段抑或整个基因组的加倍而实现的,目前还颇有争议。     

Molecular evolution of the vertebrate hexokinase gene family: Identification of a conserved fifth vertebrate hexokinase gene

Hexokinases (HK) phosphorylate sugar immediately upon its entry into cells allowing these sugars to be metabolized. A total of four hexokinases have been characterized in a diversity of vertebrates—HKI, HKII, HKIII, and HKIV. HKIV is often called glucokinase (GCK) and has half the molecular weight of the other hexokinases, as it only has one hexokinase domain, while other vertebrate HKs have two. Differing hypothesis has been proposed to explain the diversification of the hexokinase gene family. We used a genomic approach to characterize hexokinase genes in a diverse array of vertebrate species and close relatives. Surprisingly we identified a fifth hexokinase-like gene, HKDC1 that exists and is expressed in diverse vertebrates. Analysis of the amino acid sequence of HKDC1 suggests that it may function as a hexokinase. To understand the evolution of the vertebrate hexokinase gene family we established a phylogeny of the hexokinase domain in all of the vertebrate hexokinase genes, as well as hexokinase genes from close relatives of the vertebrates. Our phylogeny demonstrates that duplication of the hexokinase domain, yielding a HK with two hexokinase domains, occurred prior to the diversification of the hexokinase gene family. We also establish that GCK evolved from a two hexokinase domain-containing gene, but has lost its N-terminal hexokinase domain. We also show that parallel changes in enzymatic function of HKI and HKIII have occurred.     

Evolution of the ATP-binding-cassette transmembrane transporters of vertebrates

The ATP-binding-cassette transmembrane transporters (ABC transporters) known from vertebrates belong to four major subfamilies: (1) the P- glycoproteins (Pgp); (2) the cystic fibrosis transmembrane conductance regulators (CFTR); (3) the Tap proteins encoded with the major histocompatibility complex of mammals; and (4) the peroxisomal membrane proteins. Both Pgp and CFTR have a structure suggesting a past internal gene duplication; a phylogenetic analysis indicated that these duplications occurred independently, while an independent tandem gene duplication occurred in the case of the Tap family. Both the Pgp and Tap proteins show evidence of relationship to bacterial ABC transporters lacking internal duplication, and both are significantly more closely related to the HlyB and MsbA families of transporters from purple bacteria than they are to ABC transporters from nonpurple bacteria. The simplest hypothesis to explain this observation is that eukaryotic Pgp and Tap genes are descended from a mitochondrial gene or genes that were subsequently translocated to the nuclear genome. The Pgp genes of eukaryotes are characterized by a remarkable degree of convergent evolution between the ATP-binding cassettes of their N- terminal and C-terminal halves, whereas no such convergence is seen between the two halves of CFTR genes or between the duplicated Tap genes. Exon 13 of the CFTR gene, which encodes a putative regulatory domain not found in other ABC transporters apart from CFTR, showed high levels of both synonymous and nonsynonymous difference in comparisons among different mammalian species, suggesting that this region is a mutational hot spot.      

Kinesin-related genes from diplomonad,sponge, amphioxus,and cyclostomes: divergence pattern of kinesin family and evolution of giardial membrane-bounded organella

To understand the question of whether divergence of eukaryotic genes by gene duplications and domain shufflings proceeded gradually or intermittently during evolution, we have cloned and sequenced Giardia lamblia cDNAs encoding kinesins and kinesin-related proteins and have obtained 13 kinesin-related cDNAs, some of which are likely homologs of vertebrate kinesins involved in vesicle transfer to ER, Golgi, and plasma membrane. A phylogenetic tree of the kinesin family revealed that most gene duplications that gave rise to different kinesin subfamilies with distinct functions have been completed before the earliest divergence of extant eukaryotes. This suggests that the complex endomembrane system has arisen very early in eukaryotic evolution, and the diminutive ER and Golgi apparatus recognized in the giardial cells, together with the absence of mitochondria, might be characters acquired secondarily during the evolution of parasitism. To understand the divergence pattern of the kinesin family in the lineage leading to vertebrates, seven more Unc104-related cDNAs have been cloned from sponge, amphioxus, hagfish, and lamprey. The divergence pattern of the animal Unc104/KIF1 subfamily is characterized by two active periods in gene duplication interrupted by a considerably long period of silence, instead of proceeding gradually: animals underwent extensive gene duplications before the parazoan-eumetazoan split. In the early evolution of vertebrates around the cyclostome-gnathostome split, further gene duplications occurred, by which a variety of genes with similar structures over the entire regions were generated. This pattern of divergence is similar to those of animal genes involved in cell-cell communication and developmental control.     

Evolution and diversity of axon guidance Robo receptor family genes

The diversity of axon guidance (AG) receptors reflects gains in complexity of the animal nervous system during evolution. Members of the Roundabout (Robo) family of receptors interact with Slit proteins and play important roles in many developmental processes, including AG and neural crest cell migration. There are four members of the Robo gene family. However, the evolutionary history of Robo family genes remain obscure. We analyzed the distribution of Robo family members in metazoan species ranging in complexity from hydras to humans. We undertook a phylogenetic analysis in metazoans, synteny analysis, and ancestral chromosome mapping in vertebrates, and detected selection pressure and functional divergence among four mammalian Robo paralogs. Based on our analysis, we proposed that the ancestral Robo gene could have undergone a tandem duplication in the vertebrate ancestor; then one round of whole genome duplication events occurred before the divergence of ancestral lamprey and gnathostome, generating four paralogs in early vertebrates. Robo4 paralog underwent segmental loss in the following evolutionary process. Our results showed that Robo3 paralog is under more powerful purifying selection pressure compared with other three paralogs, which could correlate with its unique expression pattern and function. Furthermore, we found four sites under positive selection pressure on the Ig1‐2 domains of Robo4 that might interfere with its binding to Slits ligand. Diverge analysis at the amino acid level showed that Robo4 paralog have relatively greater functional diversifications than other Robo paralogs. This coincides with the fact that Robo4 predominantly functions in vascular endothelial cells but not the nervous system.     

Evolutionary history of the vertebrate mitogen activated protein kinases family

Background

The mitogen activated protein kinases (MAPK) family pathway is implicated in diverse cellular processes and pathways essential to most organisms. Its evolution is conserved throughout the eukaryotic kingdoms. However, the detailed evolutionary history of the vertebrate MAPK family is largely unclear.

Methodology/Principal Findings

The MAPK family members were collected from literatures or by searching the genomes of several vertebrates and invertebrates with the known MAPK sequences as queries. We found that vertebrates had significantly more MAPK family members than invertebrates, and the vertebrate MAPK family originated from 3 progenitors, suggesting that a burst of gene duplication events had occurred after the divergence of vertebrates from invertebrates. Conservation of evolutionary synteny was observed in the vertebrate MAPK subfamilies 4, 6, 7, and 11 to 14. Based on synteny and phylogenetic relationships, MAPK12 appeared to have arisen from a tandem duplication of MAPK11 and the MAPK13-MAPK14 gene unit was from a segmental duplication of the MAPK11-MAPK12 gene unit. Adaptive evolution analyses reveal that purifying selection drove the evolution of MAPK family, implying strong functional constraints of MAPK genes. Intriguingly, however, intron losses were specifically observed in the MAPK4 and MAPK7 genes, but not in their flanking genes, during the evolution from teleosts to amphibians and mammals. The specific occurrence of intron losses in the MAPK4 and MAPK7 subfamilies might be associated with adaptive evolution of the vertebrates by enhancing the gene expression level of both MAPK genes.

Conclusions/Significance

These results provide valuable insight into the evolutionary history of the vertebrate MAPK family.     

Small, Membrane-bound, Alternatively Spliced Forms of Ankyrin 1 Associated with the Sarcoplasmic Reticulum of Mammalian Skeletal Muscle

We have recently found that the erythroid ankyrin gene, Ank1, expresses isoforms in mouse skeletal muscle, several of which share COOH-terminal sequence with previously known Ank1 isoforms but have a novel, highly hydrophobic 72–amino acid segment at their NH₂ termini. Here, through the use of domainspecific peptide antibodies, we report the presence of the small ankyrins in rat and rabbit skeletal muscle and demonstrate their selective association with the sarcoplasmic reticulum. In frozen sections of rat skeletal muscle, antibodies to the spectrin-binding domain (anti-p65) react only with a 210-kD Ank1 and label the sarcolemma and nuclei, while antibodies to the COOH terminus of the small ankyrin (anti-p6) react with peptides of 20 to 26 kD on immunoblots and decorate the myoplasm in a reticular pattern. Mice homozygous for the normoblastosis mutation (gene symbol nb) are deficient in the 210-kD ankyrin but contain normal levels of the small ankyrins in the myoplasm. In nb/nb skeletal muscle, anti-p65 label is absent from the sarcolemma, whereas anti-p6 label shows the same distribution as in control skeletal muscle. In normal skeletal muscle of the rat, anti-p6 decorates Z lines, as defined by antidesmin distribution, and is also present at M lines where it surrounds the thick myosin filaments. Immunoblots of the proteins isolated with rabbit sarcoplasmic reticulum indicate that the small ankyrins are highly enriched in this fraction. When expressed in transfected HEK 293 cells, the small ankyrins are distributed in a reticular pattern resembling the ER if the NH₂-terminal hydrophobic domain is present, but they are uniformly distributed in the cytosol if this domain is absent. These results suggest that the small ankyrins are integral membrane proteins of the sarcoplasmic reticulum. We propose that, unlike the 210-kD form of Ank1, previously localized to the sarcolemma and believed to be a part of the supporting cytoskeleton, the small Ank1 isoforms may stabilize the sarcoplasmic reticulum by linking it to the contractile apparatus.     

Evolution of the proteasome components

 A phylogenetic analysis of proteasome subunits revealed two major families (α and β) which originated by an ancient gene duplication prior to the divergence of archaebacteria and eukaryotes. Numerous gene duplications have subsequently occurred in eukaryotes; at least nine of these duplications were shown to have occurred prior to the divergence of animals and fungi. In mammals, two genes encoding proteasome subunits (LMP2 and LMP7) are located in the major histocompatibility complex (MHC) region and play a specific role in generation of peptides for presentation by class I MHC molecules. Phylogenetic analysis of LMP7 and related sequences from mammals and lower vertebrates indicated that this locus arose by gene duplication prior to the divergence of jawed and jawless vertebrates; the time of this duplication was estimated to have been about 600 million years ago. The evolutionary history of the proteasome subunits provides support for a model of the evolution of new gene function postulating that, after gene duplication, the proteins encoded by daughter loci can adapt to specialized functions previously performed by the product of a single generalized ancestral locus. Received: 19 August 1996 / Revised: 24 December 1996     

Phylogenetic and evolutionary analysis of the septin protein family in metazoan

Septins, a conserved family of cytoskeletal GTP-binding proteins, were presented in diverse eukaryotes. Here, a comprehensive phylogenetic and evolutionary analysis for septin proteins in metazoan was carried out. First, we demonstrated that all septin proteins in metazoan could be clustered into four subgroups, and the representative homologue of every subgroup was presented in the non-vertebrate chordate Ciona intestinalis, indicating that the emergence of the four septin subgroups should have occurred prior to divergence of vertebrates and invertebrates, and the expansion of the septin gene number in vertebrates was mainly by the duplication of pre-existing genes rather than by the appearance of new septin subgroup. Second, the direct orthologues of most human septins existed in zebrafish, which suggested that human septin gene repertoire was mainly formed by as far as before the split between fishes and land vertebrates. Third, we found that the evolutionary rate within septin family in mammalian lineage varies significantly, human SEPT1, SEPT 10, SEPT 12, and SEPT 14 displayed a relative elevated evolutionary rate compared with other septin members. Our data will provide new insights for the further function study of this protein family.     

Isolation and characterization of cDNAs encoding human brain ankyrins reveal a family of alternatively spliced genes

Ankyrins are a family of membrane-associated proteins that can be divided into two immunologically distinct groups: (a) erythrocyte-related isoforms (ankyrinR) that have polarized distributions in particular cell types; and (b) brain-related isoforms (ankyrinB) that display a broader distribution. In this paper, we report the isolation and sequences of cDNAs related to two ankyrinB isoforms, human brain ankyrin 1 and 2, and show that these isoforms are produced from alternatively spliced mRNAs of a single gene. Human brain ankyrin 1 and 2 share a common NH2-terminus that is similar to human erythrocyte ankyrins, with the most striking conservation occurring between areas composed of a repeated 33-amino acid motif and between areas corresponding to the central portion of the spectrin-binding domain. In contrast, COOH-terminal sequences of brain ankyrin 1 and 2 are distinct from one another and from human erythrocyte ankyrins, and thus are candidates to mediate protein interactions that distinguish these isoforms. The brain ankyrin 2 cDNA sequence includes a stop codon and encodes a polypeptide with a predicted molecular mass of 202 kD, which is similar to the Mr of the major form of ankyrin in adult bovine brain membranes. Moreover, an antibody raised against the conserved NH2-terminal domain of brain ankyrin cross-reacts with a single Mr = 220 kD polypeptide in adult human brain. These results strongly suggest that the amino acid sequence of brain ankyrin 2 determined in this report represents the complete coding sequence of the major form of ankyrin in adult human brain. In contrast, the brain ankyrin 1 cDNAs encode only part of a larger isoform. An immunoreactive polypeptide of Mr = 440 kD, which is evident in brain tissue of young rats, is a candidate to be encoded by brain ankyrin 1 mRNA. The COOH-terminal portion of brain ankyrin 1 includes 15 contiguous copies of a novel 12-amino acid repeat. Analysis of DNA from a panel of human/rodent cell hybrids linked this human brain ankyrin gene to chromosome 4. This result, coupled with previous reports assigning the human erythrocyte ankyrin gene to chromosome 8, demonstrates that human brain and erythrocyte ankyrins are encoded by distinct members of a multigene family.     

Navigation inside a protease: substrate selection and product exit in the tricorn protease from Thermoplasma acidophilum

Chaperonins are multi-subunit double-ring complexes that mediate the folding of nascent or denatured proteins. Gene duplication has been a potent force in the evolution of chaperonins in Archaea. Here we show that gene conversion has also been an important factor. We utilized a novel maximum likehood-based phylogenetic method for scanning DNA sequence alignments for regions of anomalous phylogenetic signal, such as those affected by gene conversion. Our results suggest that in crenarchaeotes, where an ancient gene duplication producing alpha and beta subunits took place in the common ancestor of the Pyrodictium, Aeropyrum, Pyrobaculum and Sulfolobus lineages, multiple independent gene conversions have occurred between the alpha and beta genes independently in each of these groups. Significantly, the conversions have repeatedly homogenized the region of the gene encoding the substrate-binding domain. This suggests that while the alpha and beta subunits in crenarchaeotes share only 50-60% overall amino acid sequence identity, they do not possess distinct roles in the binding of substrate. Cryptic gene conversion between distantly related paralogs may be more common than is currently appreciated, and could be a significant factor in slowing the functional differentiation of proteins encoded by duplicate genes long after their duplication.