Gene and genome duplication

Genomic sequencing projects have revealed the productivity of processes duplicating genes or entire chromosome segments. Substantial proportions of the yeast, Arabidopsis and human gene complements are made up of duplicates. This has prompted much interest in the processes of duplication, functional divergence and loss of genes, has renewed the debate on whether an early vertebrate genome was tetraploid, and has inspired mathematical models and algorithms in computational biology.     

Two rounds of whole genome duplication in the ancestral vertebrate

The hypothesis that the relatively large and complex vertebrate genome was created by two ancient, whole genome duplications has been hotly debated, but remains unresolved. We reconstructed the evolutionary relationships of all gene families from the complete gene sets of a tunicate, fish, mouse, and human, and then determined when each gene duplicated relative to the evolutionary tree of the organisms. We confirmed the results of earlier studies that there remains little signal of these events in numbers of duplicated genes, gene tree topology, or the number of genes per multigene family. However, when we plotted the genomic map positions of only the subset of paralogous genes that were duplicated prior to the fish–tetrapod split, their global physical organization provides unmistakable evidence of two distinct genome duplication events early in vertebrate evolution indicated by clear patterns of four-way paralogous regions covering a large part of the human genome. Our results highlight the potential for these large-scale genomic events to have driven the evolutionary success of the vertebrate lineage.     

Gene duplication,transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.     

Early duplication and functional diversification of the opsin gene family in insects

Recent analysis of the complete mosquito Anopheles gambiae genome has revealed a far higher number of opsin genes than for either the Drosophila melanogaster genome or any other known insect. In particular, the analysis revealed an extraordinary opsin gene content expansion, whereby half are long wavelength-sensitive (LW) opsin gene duplicates. We analyzed this genomic data in relationship to other known insect opsins to estimate the relative timing of the LW opsin gene duplications and to identify "missing" paralogs in extant species. The inferred branching patterns of the LW opsin gene family phylogeny indicate at least one early gene duplication within insects before the emergence of the orders Orthoptera, Mantodea, Hymenoptera, Lepidoptera, and Diptera. These data predict the existence of one more LW opsin gene than is currently known from most insects. We tested this prediction by using a degenerate PCR strategy to screen the hymenopteran genome for novel LW opsin genes. We isolated two LW opsin gene sequences from each of five bee species, Bombus impatiens, B. terrestris, Diadasia afflicta, D. rinconis, and Osmia rufa, including 1.1 to 1.2 kb from a known (LW Rh1) and 1 kb from a new opsin gene (LW Rh2). Phylogenetic analysis suggests that the novel hymenopteran gene is orthologous to A. gambiae GPRop7, a gene that is apparently missing from D. melanogaster. Relative rate tests show that LW Rh2 is evolving at a slower rate than LW Rh1 and, therefore, may be a useful marker for higher-level hymenopteran systematics. Site-specific rate tests indicate the presence of several amino acid sites between LW Rh1 and LW Rh2 that have undergone shifts in selective constraints after duplication. These sites and others are discussed in relationship to putative structural and functional differences between the two genes.     

Gene duplication and the uniqueness of vertebrate genomes circa 1970-1999

In this article I review research undertaken over the past 30 years into the role that gene duplication played in shaping vertebrate genomes. I discuss early karyotype studies that pointed to a relative stability of mammalian and avian genomes, the discovery and possible evolutionary significance of enormous genomes in urodele amphibians and lungfish, genome compaction in certain specialised bony fish, evidence for two rounds of total genome doubling in early vertebrate evolution and the fate of duplicated genes in polyploid fish.     

Gene duplication and spectral diversification of cone visual pigments of zebrafish

Zebrafish is becoming a powerful animal model for the study of vision but the genomic organization and variation of its visual opsins have not been fully characterized. We show here that zebrafish has two red (LWS-1 and LWS-2), four green (RH2-1, RH2-2, RH2-3, and RH2-4), and single blue (SWS2) and ultraviolet (SWS1) opsin genes in the genome, among which LWS-2, RH2-2, and RH2-3 are novel. SWS2, LWS-1, and LWS-2 are located in tandem and RH2-1, RH2-2, RH2-3, and RH2-4 form another tandem gene cluster. The peak absorption spectra (lambdamax) of the reconstituted photopigments from the opsin cDNAs differed markedly among them: 558 nm (LWS-1), 548 nm (LWS-2), 467 nm (RH2-1), 476 nm (RH2-2), 488 nm (RH2-3), 505 nm (RH2-4), 355 nm (SWS1), 416 nm (SWS2), and 501 nm (RH1, rod opsin). The quantitative RT-PCR revealed a considerable difference among the opsin genes in the expression level in the retina. The expression of the two red opsin genes and of three green opsin genes, RH2-1, RH2-3, and RH2-4, is significantly lower than that of RH2-2, SWS1, and SWS2. These findings must contribute to our comprehensive understanding of visual capabilities of zebrafish and the evolution of the fish visual system and should become a basis of further studies on expression and developmental regulation of the opsin genes.     

Radiation and functional diversification of alpha keratins during early vertebrate evolution

The conquest of land was arguably one of the most fundamental ecological transitions in vertebrates and entailed significant changes in skin structure and appendages to cope with the new environment. In extant tetrapods, the rigidity of the integument is largely created by type I and type II keratins, which are structural proteins essential in forming a strong cytoplasmic network. It is expected that such proteins have undergone fundamental changes in both stem and crown tetrapods. Here, we integrate genomic, phylogenetic, and expression data in a comprehensive study on the early evolution and functional diversification of tetrapod keratins. Our analyses reveal that all type I and type II tetrapod keratins evolved from only two genes that were present in the ancestor of extant vertebrates. Subsequently, the water-to-land transition in the stem lineage of tetrapods was associated with a major radiation and functional diversification of keratin genes. These duplications acquired functions that serve rigidity in integumental hard structures and were the prime for subsequent independent keratin diversification in tetrapod lineages.     

Genome duplication, extinction and vertebrate evolution

Vertebrate evolution has been punctuated by three episodes of widespread gene or genome duplication, which have been linked with the origin of vertebrates, gnathostomes and teleosts, respectively. These three events coincide with bursts of character acquisition and increases in phenotypic complexity, and many researchers have suggested a causal relationship between the two. However, this pattern is derived from data for living taxa only; we argue here that, when fossils are taken into account, bursts of character acquisition disappear and gen(om)e duplication in vertebrate phylogeny can no longer be correlated with the origin of body plans. If patterns of character acquisition or morphological gaps between higher taxa are a reflection of phenotypic complexity, then more inclusive data sets incorporating fossil taxa provide no support for hypotheses linking gen(om)e duplications and the evolution of complexity in vertebrates.     

Teleost fish with specific genome duplication as unique models of vertebrate evolution

Whole-genome duplication (WGD) is believed to be one of the major evolutionary events that shaped the genome organization of vertebrates. Here, we review recent research on vertebrate genome evolution, specifically on WGD and its consequences for gene and genome evolution in teleost fishes. Recent genome analyses confirmed that all vertebrates experienced two rounds of WGD early in their evolution, and that teleosts experienced a subsequent additional third-round (3R)-WGD. The 3R-WGD was estimated to have occurred 320–400 million years ago in a teleost ancestor, but after its divergence from a common ancestor with living non-teleost actinopterygians (Bichir, Sturgeon, Bowfin, and Gar) based on the analyses of teleost-specific duplicate genes. This 3R-WGD was confirmed by synteny analysis and ancestral karyotype inference using the genome sequences of Tetraodon and medaka. Most of the tetrapods, on the other hand, have not experienced an additional WGD; however, they have experienced repeated chromosomal rearrangements throughout the whole genome. Therefore, different types of chromosomal events have characterized the genomes of teleosts and tetrapods, respectively. The 3R-WGD is useful to investigate the consequences of WGD because it is an evolutionarily recent WGD and thus teleost genomes retain many more WGD-derived duplicates and “traces” of their evolution. In addition, the remarkable morphological, physiological, and ecological diversity of teleosts may facilitate understanding of macrophenotypic evolution on the basis of genetic/genomic information. We highlight the teleosts with 3R-WGD as unique models for future studies on ecology and evolution taking advantage of emerging genomics technologies and systems biology environments.     

The evolution of pepsinogen C genes in vertebrates: duplication, loss and functional diversification

Background

Aspartic proteases comprise a large group of enzymes involved in peptide proteolysis. This collection includes prominent enzymes globally categorized as pepsins, which are derived from pepsinogen precursors. Pepsins are involved in gastric digestion, a hallmark of vertebrate physiology. An important member among the pepsinogens is pepsinogen C (Pgc). A particular aspect of Pgc is its apparent single copy status, which contrasts with the numerous gene copies found for example in pepsinogen A (Pga). Although gene sequences with similarity to Pgc have been described in some vertebrate groups, no exhaustive evolutionary framework has been considered so far.

Methodology/Principal Findings

By combining phylogenetics and genomic analysis, we find an unexpected Pgc diversity in the vertebrate sub-phylum. We were able to reconstruct gene duplication timings relative to the divergence of major vertebrate clades. Before tetrapod divergence, a single Pgc gene tandemly expanded to produce two gene lineages (Pgbc and Pgc2). These have been differentially retained in various classes. Accordingly, we find Pgc2 in sauropsids, amphibians and marsupials, but not in eutherian mammals. Pgbc was retained in amphibians, but duplicated in the ancestor of amniotes giving rise to Pgb and Pgc1. The latter was retained in mammals and probably in reptiles and marsupials but not in birds. Pgb was kept in all of the amniote clade with independent episodes of loss in some mammalian species. Lineage specific expansions of Pgc2 and Pgbc have also occurred in marsupials and amphibians respectively. We find that teleost and tetrapod Pgc genes reside in distinct genomic regions hinting at a possible translocation.

Conclusions

We conclude that the repertoire of Pgc genes is larger than previously reported, and that tandem duplications have modelled the history of Pgc genes. We hypothesize that gene expansion lead to functional divergence in tetrapods, coincident with the invasion of terrestrial habitats.     

Sequence divergence and diversity suggests ongoing functional diversification of vertebrate NAD metabolism

NAD is not only an important cofactor in redox reactions but has also received attention in recent years because of its physiological importance in metabolic regulation, DNA repair and signaling. In contrast to the redox reactions, these regulatory processes involve degradation of NAD and therefore necessitate a constant replenishment of its cellular pool. NAD biosynthetic enzymes are common to almost all species in all clades, but the number of NAD degrading enzymes varies substantially across taxa. In particular, vertebrates, including humans, have a manifold of NAD degrading enzymes which require a high turnover of NAD. As there is currently a lack of a systematic study of how natural selection has shaped enzymes involved in NAD metabolism we conducted a comprehensive evolutionary analysis based on intraspecific variation and interspecific divergence. We compare NAD biosynthetic and degrading enzymes in four eukaryotic model species and subsequently focus on human NAD metabolic enzymes and their orthologs in other vertebrates. We find that the majority of enzymes involved in NAD metabolism are subject to varying levels of purifying selection. While NAD biosynthetic enzymes appear to experience a rather high level of evolutionary constraint, there is evidence for positive selection among enzymes mediating NAD-dependent signaling. This is particularly evident for members of the PARP family, a diverse protein family involved in DNA damage repair and programmed cell death. Based on haplotype information and substitution rate analysis we pinpoint sites that are potential targets of positive selection. We also link our findings to a three dimensional structure, which suggests that positive selection occurs in domains responsible for DNA binding and polymerization rather than the NAD catalytic domain. Taken together, our results indicate that vertebrate NAD metabolism is still undergoing functional diversification.     

Functional diversification of sea lamprey globins in evolution and development

Agnathans have a globin repertoire that markedly differs from that of jawed (gnathostome) vertebrates. The sea lamprey (Petromyzon marinus) harbors at least 18 hemoglobin, two myoglobin, two globin X, and one cytoglobin genes. However, agnathan hemoglobins and myoglobins are not orthologous to their cognates in jawed vertebrates. Thus, blood-based O₂ transport and muscle-based O₂ storage proteins emerged twice in vertebrates from a tissue-globin ancestor. Notably, the sea lamprey displays three switches in hemoglobin expression in its life cycle, analogous to hemoglobin switching in vertebrates. To study the functional changes associated with the evolution and ontogenesis of distinct globin types, we determined O₂ binding equilibria, type of quaternary assembly, and nitrite reductase enzymatic activities of one adult (aHb5a) and one embryonic/larval hemoglobin (aHb6), myoglobin (aMb1) and cytoglobin (Cygb) of the sea lamprey. We found clear functional differentiation among globin types expressed at different developmental stages and in different tissues. Cygb and aMb1 have high O₂ affinity and nitrite reductase activity, while the two hemoglobins display low O₂ affinity and nitrite reductase activity. Cygb and aHb6 but not aHb5a show cooperative O₂ binding, correlating with increased stability of dimers, as shown by gel filtration and molecular modeling. The high O₂-affinity and the lack of cooperativity confirm the identity of the sea lamprey aMb1 as O₂ storage protein of the muscle. The dimeric structure and O₂-binding properties of sea lamprey and mammalian Cygb were very similar, suggesting a conservation of function since their divergence around 500 million years ago.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification

Whole genome duplication (WGD) is often considered to be mechanistically associated with species diversification. Such ideas have been anecdotally attached to a WGD at the stem of the salmonid fish family, but remain untested. Here, we characterized an extensive set of gene paralogues retained from the salmonid WGD, in species covering the major lineages (subfamilies Salmoninae, Thymallinae and Coregoninae). By combining the data in calibrated relaxed molecular clock analyses, we provide the first well-constrained and direct estimate for the timing of the salmonid WGD. Our results suggest that the event occurred no later in time than 88 Ma and that 40–50 Myr passed subsequently until the subfamilies diverged. We also recovered a Thymallinae–Coregoninae sister relationship with maximal support. Comparative phylogenetic tests demonstrated that salmonid diversification patterns are closely allied in time with the continuous climatic cooling that followed the Eocene–Oligocene transition, with the highest diversification rates coinciding with recent ice ages. Further tests revealed considerably higher speciation rates in lineages that evolved anadromy—the physiological capacity to migrate between fresh and seawater—than in sister groups that retained the ancestral state of freshwater residency. Anadromy, which probably evolved in response to climatic cooling, is an established catalyst of genetic isolation, particularly during environmental perturbations (for example, glaciation cycles). We thus conclude that climate-linked ecophysiological factors, rather than WGD, were the primary drivers of salmonid diversification.     

Consequences of genome duplication

Polyploidy has been widely appreciated as an important force in the evolution of plant genomes, but now it is recognized as a common phenomenon throughout eukaryotic evolution. Insight into this process has been gained by analyzing the plant, animal, fungal, and recently protozoan genomes that show evidence of whole genome duplication (a transient doubling of the entire gene repertoire of an organism). Moreover, comparative analyses are revealing the evolutionary processes that occur as multiple related genomes diverge from a shared polyploid ancestor, and in individual genomes that underwent several successive rounds of duplication. Recent research including laboratory studies on synthetic polyploids indicates that genome content and gene expression can change quickly after whole genome duplication and that cross-genome regulatory interactions are important. We have a growing understanding of the relationship between whole genome duplication and speciation. Further, recent studies are providing insights into why some gene pairs survive in duplicate, whereas others do not.     

Evolutionary diversification of the vertebrate transferrin multi-gene family

In a phylogenetic analysis of vertebrate transferrins (TFs), six major clades (subfamilies) were identified: (a) S, the mammalian serotransferrins; (b) ICA, the mammalian inhibitor of carbonic anhydrase (ICA) homologs; (c) L, the mammalian lactoferrins; (d) O, the ovotransferrins of birds and reptiles; (e) M, the melanotransferrins of bony fishes, amphibians, reptiles, birds, and mammals; and (f) M-like, a newly identified TF subfamily found in bony fishes, amphibians, reptiles, and birds. A phylogenetic tree based on the joint alignment of N-lobes and C-lobes supported the hypothesis that three separate events of internal duplication occurred in vertebrate TFs: (a) in the common ancestor of the M subfamily, (b) in the common ancestor of the M-like subfamily, and (c) in the common ancestor of other vertebrate TFs. The S, ICA, and L subfamilies were found only in placental mammals, and the phylogenetic analysis supported the hypothesis that these three subfamilies arose by gene duplication after the divergence of placental mammals from marsupials. The M-like subfamily was unusual in several respects, including the presence of a uniquely high proportion of clade-specific conserved residues, including distinctive but conserved residues in the sites homologous to those functioning in carbonate binding of human serotransferrin. The M-like family also showed an unusually high proportion of cationic residues in the positively charged region corresponding to human lactoferrampin, suggesting a distinctive role of this region in the M-like subfamily, perhaps in antimicrobial defense.     

Gene expression intensity shapes evolutionary rates of the proteins encoded by the vertebrate genome

Natural selection leaves its footprints on protein-coding sequences by modulating their silent and replacement evolutionary rates. In highly expressed genes in invertebrates, these footprints are seen in the higher codon usage bias and lower synonymous divergence. In mammals, the highly expressed genes have a shorter gene length in the genome and the breadth of expression is known to constrain the rate of protein evolution. Here we have examined how the rates of evolution of proteins encoded by the vertebrate genomes are modulated by the amount (intensity) of gene expression. To understand how natural selection operates on proteins that appear to have arisen in earlier and later phases of animal evolution, we have contrasted patterns of mouse proteins that have homologs in invertebrate and protist genomes (Precambrian genes) with those that do not have such detectable homologs (vertebrate-specific genes). We find that the intensity of gene expression relates inversely to the rate of protein sequence evolution on a genomic scale. The most highly expressed genes actually show the lowest total number of substitutions per polypeptide, consistent with cumulative effects of purifying selection on individual amino acid replacements. Precambrian genes exhibit a more pronounced difference in protein evolutionary rates (up to three times) between the genes with high and low expression levels as compared to the vertebrate-specific genes, which appears to be due to the narrower breadth of expression of the vertebrate-specific genes. These results provide insights into the differential relationship and effect of the increasing complexity of animal body form on evolutionary rates of proteins.