Evolution of the closely related, sex-related genes DM-W and DMRT1 in African clawed frogs (Xenopus)

DM-W is a dominant, female-specific, regulator of sex determination in the African clawed frog Xenopus laevis. This gene is derived from partial duplication of DMRT1, a male-related autosomal gene. We set out to better understand sex determination in Xenopus by studying this pair of genes. We found that DM-W evolved in Xenopus after divergence from the sister genus Silurana but before divergence of X. laevis and X. clivii, and that DM-W arose from partial duplication of DMRT1β, which is one of the two DMRT1 paralogs in the tetraploid ancestor of Xenopus. Using the rate ratio of nonsynonymous to synonymous substitutions per site and multilocus polymorphism data, we show that DM-W evolved non-neutrally. By cloning paralogs and using a pyrosequencing assay, we also demonstrate that DMRT1 underwent phylogenetically biased pseudogenization after polyploidization, and that expression of this gene is regulated by mechanisms that vary through development. One explanation for these observations is that the expression domain of DMRT1β was marginalized, which would explain why this paralog is dispensable in Xenopus polyploids and why DM-W has a narrow expression domain. These findings illustrate how evolution of the genetic control of stable phenotypes is facilitated by redundancy, degeneration, and compartmentalized regulation.     

Large-scale genome duplications and paralog divergence in fish

Based on fish genomic studies, we review mechanisms of divergence in duplicated genes (paralogs), resulted in small (“subfunctionalization”) or large (“neofunctionalization”) changes in paralogs. Gene divergence occurs due to several processes, such as non-synonymous substitutions, exon-intron structure rearrangement, and alterations in regulatory regions, which cause differential temporal or spatial expression of paralogous gene copies during ontogenesis.     

Questioning the Ubiquity of Neofunctionalization

Gene duplication provides much of the raw material from which functional diversity evolves. Two evolutionary mechanisms have been proposed that generate functional diversity: neofunctionalization, the de novo acquisition of function by one duplicate, and subfunctionalization, the partitioning of ancestral functions between gene duplicates. With protein interactions as a surrogate for protein functions, evidence of prodigious neofunctionalization and subfunctionalization has been identified in analyses of empirical protein interactions and evolutionary models of protein interactions. However, we have identified three phenomena that have contributed to neofunctionalization being erroneously identified as a significant factor in protein interaction network evolution. First, self-interacting proteins are underreported in interaction data due to biological artifacts and design limitations in the two most common high-throughput protein interaction assays. Second, evolutionary inferences have been drawn from paralog analysis without consideration for concurrent and subsequent duplication events. Third, the theoretical model of prodigious neofunctionalization is unable to reproduce empirical network clustering and relies on untenable parameter requirements. In light of these findings, we believe that protein interaction evolution is more persuasively characterized by subfunctionalization and self-interactions.     

Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis

Gene duplication plays an important role in the evolution of diversity and novel function and is especially prevalent in the nuclear genomes of flowering plants. Duplicate genes may be maintained through subfunctionalization and neofunctionalization at the level of expression or coding sequence. In order to test the hypothesis that duplicated regulatory genes will be differentially expressed in a specific manner indicative of regulatory subfunctionalization and/or neofunctionalization, we examined expression pattern shifts in duplicated regulatory genes in Arabidopsis. A two-way analysis of variance was performed on expression data for 280 phylogenetically identified paralogous pairs. Expression data were extracted from global expression profiles for wild-type root, stem, leaf, developing inflorescence, nearly mature flower buds, and seedpod. Gene, organ, and gene by organ interaction (G x O) effects were examined. Results indicate that 85% of the paralogous pairs exhibited a significant G x O effect indicative of regulatory subfunctionalization and/or neofunctionalization. A significant G x O effect was associated with complementary expression patterns in 45% of pairwise comparisons. No association was detected between a G x O effect and a relaxed evolutionary constraint as detected by the ratio of nonsynonymous to synonymous substitutions. Ancestral gene expression patterns inferred across a Type II MADS-box gene phylogeny suggest several cases of regulatory neofunctionalization and organ-specific nonfunctionalization. Complete linkage clustering of gene expression levels across organs suggests that regulatory modules for each organ are independent or ancestral genes had limited expression. We propose a new classification, regulatory hypofunctionalization, for an overall decrease in expression level in one member of a paralogous pair while still having a significant G x O effect. We conclude that expression divergence specifically indicative of subfunctionalization and/or neofunctionalization contributes to the maintenance of most if not all duplicated regulatory genes in Arabidopsis and hypothesize that this results in increasing expression diversity or specificity of regulatory genes after each round of duplication.     

Parallel duplication and partial subfunctionalization of β-catenin/armadillo during insect evolution

β-Catenin is a multifunctional scaffolding protein with roles in Wnt signaling, cell adhesion, and centrosome separation. Here, we report on independent duplications of the insect β-Catenin ortholog armadillo (arm) in the red flour beetle Tribolium castaneum and the pea aphid Acyrthosiphon pisum. Detailed sequence analysis shows that in both species, one paralog lost critical residues of the α-Catenin binding domain, which is essential for cell adhesion, and accumulated a dramatically higher number of amino acid substitutions in the central Arm repeat domain. Residues associated with aspects of Wnt signaling, however, are conserved in both paralogs. Consistent with these molecular signatures, the effects of specific and combinatorial knockdown experiments in the Tribolium embryo indicate that the duplication resulted in redundant involvement in Wnt signaling of both β-Catenin paralogs but differential inheritance of the ancestral cell adhesion and centrosome separation functions. We conclude that the duplicated pea aphid and flour beetle β-catenin genes experienced partial subfunctionalization, which appears to be evolutionarily favored. Providing first evidence of genetic separability of the cell adhesion and centrosome separation functions, the duplicated Tribolium and Acyrthosiphon arm paralogs offer new inroads for context-specific analyses of β-Catenin. Our data also revealed the conservation of a C-terminally truncated Arm isoform in both singleton and duplicated homologs, suggesting an as yet unexplored role in Wnt signaling.     

Functional analysis of gene duplications in Saccharomyces cerevisiae

Gene duplication can occur on two scales: whole-genome duplications (WGD) and smaller-scale duplications (SSD) involving individual genes or genomic segments. Duplication may result in functionally redundant genes or diverge in function through neofunctionalization or subfunctionalization. The effect of duplication scale on functional evolution has not yet been explored, probably due to the lack of global knowledge of protein function and different times of duplication events. To address this question, we used integrated Bayesian analysis of diverse functional genomic data to accurately evaluate the extent of functional similarity and divergence between paralogs on a global scale. We found that paralogs resulting from the whole-genome duplication are more likely to share interaction partners and biological functions than smaller-scale duplicates, independent of sequence similarity. In addition, WGD paralogs show lower frequency of essential genes and higher synthetic lethality rate, but instead diverge more in expression pattern and upstream regulatory region. Thus, our analysis demonstrates that WGD paralogs generally have similar compensatory functions but diverging expression patterns, suggesting a potential of distinct evolutionary scenarios for paralogs that arose through different duplication mechanisms. Furthermore, by identifying these functional disparities between the two types of duplicates, we reconcile previous disputes on the relationship between sequence divergence and expression divergence or essentiality.     

Gene complexity and gene duplicability

Eukaryotic genes are on average more complex than prokaryotic genes in terms of expression regulation, protein length, and protein-domain structure [1-5]. Eukaryotes are also known to have a higher rate of gene duplication than prokaryotes do [6, 7]. Because gene duplication is the primary source of new genes [], the average gene complexity in a genome may have been increased by gene duplication if complex genes are preferentially duplicated. Here, we test this "gene complexity and gene duplicability" hypothesis with yeast genomic data. We show that, on average, duplicate genes from either whole-genome or individual-gene duplication have longer protein sequences, more functional domains, and more cis-regulatory motifs than singleton genes. This phenomenon is not a by-product of previously known mechanisms, such as protein function [10-13], evolutionary rate [14, 15], dosage [11], and dosage balance [16], that influence gene duplicability. Rather, it appears to have resulted from the sub-neo-functionalization process in duplicate-gene evolution [11]. Under this process, complex genes are more likely to be retained after duplication because they are prone to subfunctionalization, and gene complexity is regained via subsequent neofunctionalization. Thus, gene duplication increases both gene number and gene complexity, two important factors in the origin of genomic and organismal complexity.     

Density peaks of paralog pairs in human and mouse genomes

Paralog gene trees, which reflect the increase of genomic complexity in the evolution, can be complicated and ambiguous. A simpler complementary approach is analysis of density distribution of paralog pairs. It can reveal general features of genome evolution, which may be hidden in the forest of gene trees. It is known that distribution of human paralog pairs along the axis of protein divergence between pair members forms two main peaks. Here I show that there are three main peaks in the mouse genome. Thus, the multimodality of paralog pair distribution seems to be a fundamental feature of mammalian genomes. Despite the great diversity of domains presented in small amounts or in multidomain architectures with a few predominant domains, both in human and mouse the first peak consists mostly of gene pairs with zinc finger domains or olfactory receptor domain. In the mouse the olfactory receptor predominates, which stipulates the three-peak distribution (since in the olfactory receptors the second peak is closer to the first peak than in other genes). The mammalian-wide zinc finger orthologs are biased towards the second peak. Thus, the marsupial orthologs are nearly absent in the first peak of human and mouse. The gene pairs in the first peak show a lower ratio of nonsynonymous to synonymous substitutions, which suggests that their evolution is more constrained. The plausible explanation is that they are in subfunctionalization state (partition of initial function of ancestral gene), whereas the second peak contains gene pairs that are already in neofunctionalization state (acquiring of novel functions). These data suggest that the adaptive radiation of mammals was accompanied by a burst of duplication of zinc finger genes, which are located in the first (most recent) peak of paralog pairs.     

In silico evidence for functional specialization after genome duplication in yeast

A fairly recent whole-genome duplication (WGD) event in yeast enables the effects of gene duplication and subsequent functional divergence to be characterized. We examined 15 ohnolog pairs (i.e. paralogs from a WGD) out of c . 500 Saccharomyces cerevisiae ohnolog pairs that have persisted over an estimated 100 million years of evolution. These 15 pairs were chosen for their high levels of asymmetry, i.e. within the pair, one ohnolog had evolved much faster than the other. Sequence comparisons of the 15 pairs revealed that the faster evolving duplicated genes typically appear to have experienced partially – but not fully – relaxed negative selection as evidenced by an average nonsynonymous/synonymous substitution ratio (d N /d S _avg=0.44) that is higher than the slow-evolving genes' ratio (d N /d S _avg=0.14) but still <1. Increased number of insertions and deletions in the fast-evolving genes also indicated loosened structural constraints. Sequence and structural comparisons indicated that a subset of these pairs had significant differences in their catalytically important residues and active or cofactor-binding sites. A literature survey revealed that several of the fast-evolving genes have gained a specialized function. Our results indicate that subfunctionalization and even neofunctionalization has occurred along with degenerative evolution, in which unneeded functions were destroyed by mutations.     

Gene duplication within the Green Lineage: the case of TEL genes

Recent years have witnessed a breathtaking increase in the availability of genome sequence data, providing evidence of the highly duplicate nature of eukaryotic genomes. Plants are exceptional among eukaryotic organisms in that duplicate loci compose a large fraction of their genomes, partly because of the frequent occurrence of polyploidy (or whole-genome duplication) events. Tandem gene duplication and transposition have also contributed to the large number of duplicated genes in plant genomes. Evolutionary analyses allowed the dynamics of duplicate gene evolution to be studied and several models were proposed. It seems that, over time, many duplicated genes were lost and some of those that were retained gained new functions and/or expression patterns (neofunctionalization) or subdivided their functions and/or expression patterns between them (subfunctionalization). Recent studies have provided examples of genes that originated by duplication with successive diversification within plants. In this review, we focused on the TEL (TERMINAL EAR1-like) genes to illustrate such mechanisms. Emerged from the mei2 gene family, these TEL genes are likely to be land plant-specific. Phylogenetic analyses revealed one or two TEL copies per diploid genome. TEL gene degeneration and loss in several Angiosperm species such as in poplar and maize seem to have occurred. In Arabidopsis thaliana, whose genome experienced at least three polyploidy events followed by massive gene loss and genomic reorganization, two TEL genes were retained and two new shorter TEL-like (MCT) genes emerged. Molecular and expression analyses suggest for these genes sub- and neofunctionalization events, but confirmation will come from their functional characterization.     

Gene duplication in the evolution of sexual dimorphism

Males and females share most of the same genes, so selection in one sex will typically produce a correlated response in the other sex. Yet, the sexes have evolved to differ in a multitude of behavioral, morphological, and physiological traits. How did this sexual dimorphism evolve despite the presence of a common underlying genome? We investigated the potential role of gene duplication in the evolution of sexual dimorphism. Because duplication events provide extra genetic material, the sexes each might use this redundancy to facilitate sex‐specific gene expression, permitting the evolution of dimorphism. We investigated this hypothesis at the genome‐wide level in Drosophila melanogaster, using the presence of sex‐biased expression as a proxy for the sex‐specific specialization of gene function. We expected that if sexually antagonistic selection is a potent force acting upon individual genes, duplication will result in paralog families whose members differ in sex‐biased expression. Gene members of the same duplicate family can have different expression patterns in males versus females. In particular, duplicate pairs containing a male‐biased gene are found more frequently than expected, in agreement with previous studies. Furthermore, when the singleton ortholog is unbiased, duplication appears to allow one of the paralog copies to acquire male‐biased expression. Conversely, female‐biased expression is not common among duplicates; fewer duplicate genes are expressed in the female‐soma and ovaries than in the male‐soma and testes. Expression divergence exists more in older than in younger duplicates pairs, but expression divergence does not correlate with protein sequence divergence. Finally, genomic proximity may have an effect on whether paralogs differ in sex‐biased expression. We conclude that the data are consistent with a role of gene duplication in fostering male‐biased, but not female‐biased, gene expression, thereby aiding the evolution of sexual dimorphism.     

Evolution of the vertebrate twist family and synfunctionalization: a mechanism for differential gene loss through merging of expression domains

Twist genes are essential for embryonic development and are conserved from jellyfish to human. To study the vertebrate twist family and its evolution, the entire complement of twist genes was obtained for 9 representative species. Phylogenetic analysis showed that a single protochordate twist gene was duplicated at least twice before the teleost-tetrapod split to give rise to 3 ancestral genes, which were further duplicated or deleted, resulting in fluctuating number of twist paralogs in different vertebrate lineages. To find whether changes in gene copy number were associated with changes in gene function, embryonic expression patterns of twist orthologs were evaluated against the number of twist paralogs in different species. The results showed evidence for both neo- and subfunctionalization, and, in addition, for loss of an ancestral regulatory gene. For example, in Xenopus, twist2 was lost, but the twist1 paralog acquired, and therefore preserved, twist2 function. A general model is proposed to explain the data. In this process, termed synfunctionalization, one paralog acquires the expression domain(s) of another. The merging may lead to function shuffle. Alternatively, it may leave one paralog redundant and thus subject to deletion--while its function is retained by the surviving paralog(s). Synfunctionalization is a mechanism that, together with neo- and subfunctionalization, may work to establish equilibrium in the number of genes that regulate developmental processes; it may regulate the complexity of regulatory regions as well as gene copy number and therefore may play a role in evolution of gene function and the structure of genome.     

Evolution of Olfactory Receptors Tuned to Mustard Oils in Herbivorous Drosophilidae

The diversity of herbivorous insects is attributed to their propensity to specialize on toxic plants. In an evolutionary twist, toxins betray the identity of their bearers when herbivores coopt them as cues for host-plant finding, but the evolutionary mechanisms underlying this phenomenon are poorly understood. We focused on Scaptomyza flava, an herbivorous drosophilid specialized on isothiocyanate (ITC)-producing (Brassicales) plants, and identified Or67b paralogs that were triplicated as mustard-specific herbivory evolved. Using in vivo heterologous systems for the expression of olfactory receptors, we found that S. flava Or67bs, but not the homologs from microbe-feeding relatives, responded selectively to ITCs, each paralog detecting different ITC subsets. Consistent with this, S. flava was attracted to ITCs, as was Drosophila melanogaster expressing S. flava Or67b3 in the homologous Or67b olfactory circuit. ITCs were likely coopted as olfactory attractants through gene duplication and functional specialization (neofunctionalization and subfunctionalization) in S. flava, a recently derived herbivore.     

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, “marking” one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.     

Duplicated genes evolve slower than singletons despite the initial rate increase

Background

Gene duplication is an important mechanism that can lead to the emergence of new functions during evolution. The impact of duplication on the mode of gene evolution has been the subject of several theoretical and empirical comparative-genomic studies. It has been shown that, shortly after the duplication, genes seem to experience a considerable relaxation of purifying selection.

Results

Here we demonstrate two opposite effects of gene duplication on evolutionary rates. Sequence comparisons between paralogs show that, in accord with previous observations, a substantial acceleration in the evolution of paralogs occurs after duplication, presumably due to relaxation of purifying selection. The effect of gene duplication on evolutionary rate was also assessed by sequence comparison between orthologs that have paralogs (duplicates) and those that do not (singletons). It is shown that, in eukaryotes, duplicates, on average, evolve significantly slower than singletons. Eukaryotic ortholog evolutionary rates for duplicates are also negatively correlated with the number of paralogs per gene and the strength of selection between paralogs. A tally of annotated gene functions shows that duplicates tend to be enriched for proteins with known functions, particularly those involved in signaling and related cellular processes; by contrast, singletons include an over-abundance of poorly characterized proteins.

Conclusions

These results suggest that whether or not a gene duplicate is retained by selection depends critically on the pre-existing functional utility of the protein encoded by the ancestral singleton. Duplicates of genes of a higher biological import, which are subject to strong functional constraints on the sequence, are retained relatively more often. Thus, the evolutionary trajectory of duplicated genes appears to be determined by two opposing trends, namely, the post-duplication rate acceleration and the generally slow evolutionary rate owing to the high level of functional constraints.

     

Preferential retention of the slowly evolving gene in pairs of duplicates in angiosperm genomes

Gene duplication provides raw material for functional innovation, but gene duplicability varies considerably. Previous studies have found widespread asymmetrical sequence evolution between paralogs. However, it remains unknown whether the rate of evolution among paralogs affects their propensity of being retained after another round of whole-genome duplication (WGD). In this study, we investigated gene groups that have experienced two successive WGDs to determine which of two older duplicates with different evolutionary rates was more likely to retain both younger duplicates. To uncouple the measurement of evolutionary rates from any assignment of duplicate or singleton status, we measured the evolutionary rates of singleton genes in out-lineages but classified these singleton genes according to whether they are retained or not in a crown group of species. We found that genes that retained younger duplicates in the crown group of genomes were more constrained prior to the younger duplication event than those that failed to leave duplicates. In addition, we also found that the retained clades have more genes in out-lineages. Subsequent analyses showed that genes in the retained clades were expressed more broadly and highly than genes in the singleton clades. We concluded that the set of repeatedly retained genes after two WGDs is biased toward slowly evolving genes in angiosperms, suggesting that the potential of genes for both functional conservation and divergence likely affects their propensity of being retained after WGD in angiosperms.