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.     

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.     

Functional diversification of B MADS-box homeotic regulators of flower development: Adaptive evolution in protein-protein interaction domains after major gene duplication events

B-class MADS-box genes have been shown to be the key regulators of petal and stamen specification in several eudicot model species such as Arabidopsis thaliana, Antirrhinum majus, and Petunia hybrida. Orthologs of these genes have been found across angiosperms and gymnosperms, and it is thought that the basic regulatory function of B proteins is conserved in seed plant lineages. The evolution of B genes is characterized by numerous duplications that might represent key elements fostering the functional diversification of duplicates with a deep impact on their role in the evolution of the floral developmental program. To evaluate this, we performed a rigorous statistical analysis with B gene sequences. Using maximum likelihood and Bayesian methods, we estimated molecular substitution rates and determined the selective regimes operating at each residue of B proteins. We implemented tests that rely on phylogenetic hypotheses and codon substitution models to detect significant differences in substitution rates (DSRs) and sites under positive adaptive selection (PS) in specific lineages before and after duplication events. With these methods, we identified several protein residues fixed by PS shortly after the origin of PISTILLATA-like and APETALA3-like lineages in angiosperms and shortly after the origin of the euAP3-like lineage in core eudicots, the 2 main B gene duplications. The residues inferred to have been fixed by positive selection lie mostly within the K domain of the protein, which is key to promote heterodimerization. Additionally, we used a likelihood method that accommodates DSRs among lineages to estimate duplication dates for AP3-PI and euAP3-TM6, calibrating with data from the fossil record. The dates obtained are consistent with angiosperm origins and diversification of core eudicots. Our results strongly suggest that novel multimer formation with other MADS proteins could have been crucial for the functional divergence of B MADS-box genes. We thus propose a mechanism of functional diversification and persistence of gene duplicates by the appearance of novel multimerization capabilities after duplications. Multimer formation in different combinations of regulatory proteins can be a mechanistic basis for the origin of novel regulatory functions and a gene regulatory mechanism for the appearance of morphological innovations.     

Seeing double: gene duplication and diversification in plant secondary metabolism

Gene duplications drive the recruitment of genes for secondary metabolism. Gene copies are gradually modified to create genes with specificities and expression patterns adapted to the needs of the new pathway in which they are involved. Duplicated genes are often in tandem repeats, forming clusters within the plant genome. However, in some cases, clusters of nonhomologous genes have also been identified as forming a functional unit. The selective forces that have caused the establishment of new pathways are far from understood and might have changed repeatedly during evolution owing to the continuously changing environment. Recent data show that the way several classes of secondary compounds are scattered among species is attributable to independent recruitment and the inactivation of biosynthetic enzymes.     

Adaptive evolution after gene duplication

One of the two ribonuclease genes in a leaf-eating monkey has adapted to a role in the digestion of bacterial RNA. Following duplication of the ancestral ribonuclease gene, adaptation occurred through a series of changes in the amino acid sequence of the protein it encodes. This example is a good illustration of how specialization of protein function after gene duplication can be as source of novel protein functions.     

Functional diversification of lepidopteran opsins following gene duplication.

A comparative approach was taken for identifying amino acid substitutions that may be under positive Darwinian selection and are correlated with spectral shifts among orthologous and paralogous lepidopteran long wavelength-sensitive (LW) opsins. Four novel LW opsin fragments were isolated, cloned, and sequenced from eye-specific cDNAs from two butterflies, Vanessa cardui (Nymphalidae) and Precis coenia (Nymphalidae), and two moths, Spodoptera exigua (Noctuidae) and Galleria mellonella (Pyralidae). These opsins were sampled because they encode visual pigments having a naturally occurring range of lambda(max) values (510-530 nm), which in combination with previously characterized lepidopteran opsins, provide a complete range of known spectral sensitivities (510-575 nm) among lepidopteran LW opsins. Two recent opsin gene duplication events were found within the papilionid but not within the nymphalid butterfly families through neighbor-joining, maximum parsimony, and maximum likelihood phylogenetic analyses of 13 lepidopteran opsin sequences. An elevated rate of evolution was detected in the red-shifted Papilio Rh3 branch following gene duplication, because of an increase in the amino acid substitution rate in the transmembrane domain of the protein, a region that forms the chromophore-binding pocket of the visual pigment. A maximum likelihood approach was used to estimate omega, the ratio of nonsynonymous to synonymous substitutions per site. Branch-specific tests of selection (free-ratio) identified one branch with omega = 2.1044, but the small number of substitutions involved was not significantly different from the expected number of changes under the neutral expectation of omega = 1. Ancestral sequences were reconstructed with a high degree of certainty from these data. Reconstructed ancestral sequences revealed several instances of convergence to the same amino acid between butterfly and vertebrate cone pigments, and between independent branches of the butterfly opsin tree that are correlated with spectral shifts.     

Gene duplication,genome duplication,and the functional diversification of vertebrate globins

The functional diversification of the vertebrate globin gene superfamily provides an especially vivid illustration of the role of gene duplication and whole-genome duplication in promoting evolutionary innovation. For example, key globin proteins that evolved specialized functions in various aspects of oxidative metabolism and oxygen signaling pathways (hemoglobin [Hb], myoglobin [Mb], and cytoglobin [Cygb]) trace their origins to two whole-genome duplication events in the stem lineage of vertebrates. The retention of the proto-Hb and Mb genes in the ancestor of jawed vertebrates permitted a physiological division of labor between the oxygen-carrier function of Hb and the oxygen-storage function of Mb. In the Hb gene lineage, a subsequent tandem gene duplication gave rise to the proto α- and β-globin genes, which permitted the formation of multimeric Hbs composed of unlike subunits (α₂β₂). The evolution of this heteromeric quaternary structure was central to the emergence of Hb as a specialized oxygen-transport protein because it provided a mechanism for cooperative oxygen-binding and allosteric regulatory control. Subsequent rounds of duplication and divergence have produced diverse repertoires of α- and β-like globin genes that are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different stages of prenatal development and postnatal life. In the ancestor of jawless fishes, the proto Mb and Hb genes appear to have been secondarily lost, and the Cygb homolog evolved a specialized respiratory function in blood-oxygen transport. Phylogenetic and comparative genomic analyses of the vertebrate globin gene superfamily have revealed numerous instances in which paralogous globins have convergently evolved similar expression patterns and/or similar functional specializations in different organismal lineages.     

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.     

Opsin gene duplication and diversification in the guppy, a model for sexual selection

Identification of genes that control variation in adaptive characters is a prerequisite for understanding the processes that drive sexual and natural selection. Male coloration and female colour perception play important roles in mate choice in the guppy (Poecilia reticulata), a model organism for studies of natural and sexual selection. We examined a potential source for the known variation in colour perception, by analysing genomic and complementary DNA sequences of genes that code for visual pigment proteins. We find high sequence variability, both within and between populations, and expanded copy number for long-wave sensitive (LWS) opsin genes. Alleles with non-synonymous changes that suggest dissimilar spectral tuning properties occur in the same population and even in the same individual, and the high frequency of non-synonymous substitutions argues for diversifying selection acting on these proteins. Therefore, variability in tuning amino acids is partitioned within individuals and populations of the guppy, in contrast to variability for LWS at higher taxonomic levels in cichlids, a second model system for differentiation owing to sexual selection. Since opsin variability parallels the extreme male colour polymorphism within guppy populations, we suggest that mate choice has been a major factor driving the coevolution of opsins and male ornaments in this species.     

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.     

Functional diversification of kir7.1 in cichlids accelerated by gene duplication

Mutation in the inward rectifier potassium channel gene, kir7.1, was previously identified as being responsible for the broader stripe zebrafish skin pattern mutant, jaguar/obelix. An amino acid substitution in this channel causes a broader stripe pattern than that of wild type zebrafish. In this study we analyzed cichlid homologs of the zebrafish kir7.1 gene. We identified two kinds of homologous genes in cichlids and named them cikir7.1 and cikir7.2. Southern hybridization using cichlid genome revealed that cichlids from the African Great Lakes, South America and Madagascar have two copies of the gene. Cichlids from Sri Lanka, however, showed only one band in this experiment. Database analysis revealed that only one copy of the kir7.1 gene exists in the genomes of the teleosts zebrafish, tetraodon, takifugu, medaka and stickleback. The deduced amino acid sequence of cikir7.1 is highly conserved among African cichlids, whereas that of cikir7.2 has several amino acid substitutions even in conserved transmembrane domains. Gene expression analysis revealed that cikir7.1 is expressed specifically in brain and eye, and cikir7.2 in testis and ovary; zebrafish kir7.1, however, is expressed in brain, eye, skin, caudal fin, testis and ovary. These results suggest that gene duplication of the cichlid kir7.1 occurred in a common ancestor of the family Cichlidae, that the function of parental kir7.1 was then divided into two genes, cikir7.1 and cikir7.2, and that the evolutionary rate of cikir7.2 might have been accelerated, thereby effecting functional diversification in the cichlid lineage. Thus, the evolution of kir7.1 genes in cichlids provides a typical example of gene duplication--one gene is conserved while the other becomes specialized for a novel function.     

Patterns of gene duplication and functional evolution during the diversification of the AGAMOUS subfamily of MADS box genes in angiosperms

Members of the AGAMOUS (AG) subfamily of MIKC-type MADS-box genes appear to control the development of reproductive organs in both gymnosperms and angiosperms. To understand the evolution of this subfamily in the flowering plants, we have identified 26 new AG-like genes from 15 diverse angiosperm species. Phylogenetic analyses of these genes within a large data set of AG-like sequences show that ancient gene duplications were critical in shaping the evolution of the subfamily. Before the radiation of extant angiosperms, one event produced the ovule-specific D lineage and the well-characterized C lineage, whose members typically promote stamen and carpel identity as well as floral meristem determinacy. Subsequent duplications in the C lineage resulted in independent instances of paralog subfunctionalization and maintained functional redundancy. Most notably, the functional homologs AG from Arabidopsis and PLENA (PLE) from Antirrhinum are shown to be representatives of separate paralogous lineages rather than simple genetic orthologs. The multiple subfunctionalization events that have occurred in this subfamily highlight the potential for gene duplication to lead to dissociation among genetic modules, thereby allowing an increase in morphological diversity.     

Impact of gene gains,losses and duplication modes on the origin and diversification of vertebrates

The study of the evolutionary origin of vertebrates has been linked to the study of genome duplications since Susumo Ohno suggested that the successful diversification of vertebrate innovations was facilitated by two rounds of whole-genome duplication (2R-WGD) in the stem vertebrate. Since then, studies on the functional evolution of many genes duplicated in the vertebrate lineage have provided the grounds to support experimentally this link. This article reviews cases of gene duplications derived either from the 2R-WGD or from local gene duplication events in vertebrates, analyzing their impact on the evolution of developmental innovations. We analyze how gene regulatory networks can be rewired by the activity of transposable elements after genome duplications, discuss how different mechanisms of duplication might affect the fate of duplicated genes, and how the loss of gene duplicates might influence the fate of surviving paralogs. We also discuss the evolutionary relationships between gene duplication and alternative splicing, in particular in the vertebrate lineage. Finally, we discuss the role that the 2R-WGD might have played in the evolution of vertebrate developmental gene networks, paying special attention to those related to vertebrate key features such as neural crest cells, placodes, and the complex tripartite brain. In this context, we argue that current evidences points that the 2R-WGD may not be linked to the origin of vertebrate innovations, but to their subsequent diversification in a broad variety of complex structures and functions that facilitated the successful transition from peaceful filter-feeding non-vertebrate ancestors to voracious vertebrate predators.     

Globin gene family evolution and functional diversification in annelids

Globins are the most common type of oxygen-binding protein in annelids. In this paper, we show that circulating intracellular globin (Alvinella pompejana and Glycera dibranchiata), noncirculating intracellular globin (Arenicola marina myoglobin) and extracellular globin from various annelids share a similar gene structure, with two conserved introns at canonical positions B12.2 and G7.0. Despite sequence divergence between intracellular and extracellular globins, these data strongly suggest that these three globin types are derived from a common ancestral globin-like gene and evolved by duplication events leading to diversification of globin types and derived functions. A phylogenetic analysis shows a distinct evolutionary history of annelid extracellular hemoglobins with respect to intracellular annelid hemoglobins and mollusc and arthropod extracellular hemoglobins. In addition, dehaloperoxidase (DHP) from the annelid, Amphitrite ornata, surprisingly exhibits close phylogenetic relationships to some annelid intracellular globins. We have characterized the gene structure of A. ornata DHP to confirm assumptions about its homology with globins. It appears that it has the same intron position as in globin genes, suggesting a common ancestry with globins. In A. ornata, DHP may be a derived globin with an unusual enzymatic function.     

Patterns of gene duplication and functional diversification during the evolution of the AP1/SQUA subfamily of plant MADS-box genes

Members of the AP1/SQUA subfamily of plant MADS-box genes play broad roles in the regulation of reproductive meristems, the specification of sepal and petal identities, and the development of leaves and fruits. It has been shown that AP1/SQUA-like genes are angiosperm-specific, and have experienced several major duplication events. However, the evolutionary history of this subfamily is still uncertain. Here, we report the isolation of 14 new AP1/SQUA-like genes from seven early-diverging eudicots and the identification of 11 previously uncharacterized ESTs and genomic sequences from public databases. Sequence comparisons of these and other published sequences reveal a conserved C-terminal region, the FUL motif, in addition to the known euAP1/paleoAP1 motif, in AP1/SQUA-like proteins. Phylogenetic analyses further suggest that there are three major lineages (euAP1, euFUL, and AGL79) in core eudicots, likely resulting from two close duplication events that predated the divergence of core eudicots. Among the three lineages, euFUL is structurally very similar to FUL-like genes from early-diverging eudicots and basal angiosperms, whereas euAP1 might have originally been generated through a 1-bp deletion in the exon 8 of an ancestral euFUL- or FUL-like gene. Because euFUL- and FUL-like genes usually have broad expression patterns, we speculate that AP1/SQUA-like genes initially had broad functions. Based on these observations, the evolutionary fates of duplicate genes and the contributions of the frameshift mutation and alternative splicing to functional diversity are discussed.