Paralogs in Polyploids: One for All and All for One?

Plants are remarkable with regard to their capability to go beyond diploidy and exist as stable polyploids. It is estimated that ∼50 to 70% of flowering plants are polyploids (Masterson, 1994; Otto and Whitton, 2000). Theory predicts that gene duplicates degenerate and disappear, but the extensive retention of paralogs in polyploids or paleopolyploids runs counter to this prediction. What is the fate of gene duplicates in polyploids? Multiple evolutionary fates for paralogs have been explored from different perspectives (Lynch and Conery, 2000; Wolfe, 2001). Here, this question is revisited. Specifically, several scenarios will be explored based on the notion that point mutations or deletions of the paralogs can be deleterious if their products are involved in protein–protein or protein–DNA complexes. This perspective is not exclusive but rather complementary to the ones developed in previous studies. The extent of paralog survival in plants is so important that protein diversity in this kingdom has been hypothesized to be generated primarily through gene duplication rather than by alternative splicing, as has been proposed for vertebrates (Hodgkin, 2000; Yu et al., 2002). However, one cannot overlook that much of vertebrate diversity comes from the fishes, major lineages of which appear to be the result of genome duplication events. While still controversial, it would not be unreasonable to say that data are suggestive of a role for polyploidy in the majority of vertebrate diversity (Makalowski, 2001). Thus, while the focus of this review is the fate of paralogs in plants, the ideas developed are applicable to other organisms. Force et al. (1999) proposed that random mutation and drift can alter the expression patterns of paralogous genes and segregate their regulation so that both (together) cover the full range of functions of the ancestral gene. Accordingly, both copies would be needed. For instance, each paralogous gene might have a specific domain of expression. This process is known as subfunctionalization. Once subfunctional mutations arise and become fixed, selection will maintain both paralogous copies. Acquisition of a new function by one of the paralogs, neofunctionalization (a rare event at least in small populations), can also explain evolutionary retention (Ohno, 1970; Force et al., 1999). Another mechanism that can lead to the retention of paralogs, not exclusive of subfunctionalization, is a requirement for the preservation of stoichiometry within complexes or pathways. Gene products that participate together in a macromolecular complex or pathway are thought to require a dosage balance for proper function (Veitia, 2002, 2004). The importance of gene dosage balance is obvious in organisms that alternate between haploid and diploid phases, especially haplodiplobiontics, such as plants, which undergo development in both phases. During the haploid gametophyte phase, deletion of genes required for viability cannot be tolerated, and the balance among interacting components is achieved with one gene copy of each locus. In the diploid state, a heterozygous deletion of a gene can be deleterious to diploid development if it leads to a stoichiometric imbalance among components of a complex (Veitia, 2002, 2003a). General theories of dominance have stated that one of the advantages of diploidy is that one allele can mask deleterious mutations in the other allele at a particular locus (Orr and Otto, 1994). At dosage sensitive genes, however, one active allele in a diploid cannot mask a complete lack of activity of the other. Accordingly, a recent genomic survey has shown that deletion of one allele of 21% of the yeast genes analyzed leads to defective growth (i.e., these genes are haplo-insufficient; Kondrashov and Koonin, 2004). In humans, the effects of dosage sensitivity might be more dramatic because haplo-insufficient mutations responsible for disease seem to outnumber the recessive ones (Kondrashov and Koonin, 2004). Further evidence for the importance of dosage balance has been found in yeast, which is able to undergo vegetative growth in both haploid and diploid conditions. Dosage-sensitive genes (whose underexpression or overexpression induces a growth defect) are more likely to encode proteins involved in macromolecular complexes than genes with low dosage sensitivity (Papp et al., 2003). Birchler et al. (2001) also presented evidence suggesting that the stoichiometric relationships among the components of regulatory complexes affect target gene expression. The contribution of polyploids to evolution, either with multisomic inheritance or as diploidized forms, presumably would be minor if most duplicated copies were to be silenced (nonfunctionalized). Walsh (1995) estimated that ∼99% of duplicate genes would become pseudogenes under certain conditions of population size and mutation rates. The analysis of ancient polyploids shows that a much larger fraction of functional paralogs is actually retained. In the paleopolyploid yeast Saccharomyces cerevisiae, ∼12% of paralogs have retained functionality for 100 million years (Kellis et al., 2004). The fraction of retained paralogs rises to 72% in maize, through 11 million years of evolution (Ahn and Tanksley, 1993). However, their functionality is yet to be assessed. As discussed below, Arabidopsis thaliana also provides many examples of functional paralogs retained over several cycles of polyploidization. Gene duplicates generated by genome duplication are likely to be retained if they undergo selection for absolute gene product dosage (Otto and Whitton, 2000). This simple idea entails more complexity than might be suspected at first sight. Thus, before proceeding with the discussion, we must consider that on the one hand, in plant ploidy series the mRNA level of most genes is roughly proportional to ploidy (Guo et al., 1996). On the other hand, when similar cell types are compared, cell volume tends to be proportional to DNA content (Cavalier-Smith, 1982; Gregory, 2001). For instance, pollen grain volume and stomatal cell size correlate with ploidy (Bennett, 1972; Masterson, 1994). Even subtle changes in genomic DNA content (i.e., by adding or substracting small amounts of DNA) have an impact on cell volume (Gregory, 2001). Thus, when correlations between gene expression, cell volume, and ploidy are considered together, an increase in gene dosage after autopolyploidization is not expected to be paralleled by an increase in the concentrations of gene products. An increase in the absolute gene product concentration can only be accomplished by a reduction in cell volume if expression rates of the paralogous alleles do not change and if there is some degree of saturation in protein degradation. A reduction in cell volume could be achieved by the deletion of DNA. A number of other characteristics of the three-dimensional cell will be altered with an increase in ploidy. For example, the surface-to-volume ratio of the plasma membrane to the cytosol, or of the nuclear envelope to the nucleoplasm and cytoplasm, might be altered. Because these effects have been less well studied, the discussion will focus on gene products in suspension in the cytosol or the nucleoplasm. Absolute increases in the quantity of functional gene products are straightforward when the proteins act as monomers poorly connected within the cellular interaction network, for which increased dosage can be achieved through segmental DNA duplication. However, for dosage-sensitive genes, increases in the expression level of one component can lead to an alteration in the amount of functional complex. In such cases, a coordinated increase of dosage of the interacting gene products must be achieved, possibly through a reduction in cellular volume after genome duplication. The survival of one paralog will depend on the retention of paralogs encoding the other stoichiometric interactors (see below). In yeast, homozygous deletions of many genes are associated with little or no apparent phenotypic effects. At least one-quarter of those gene deletions that have no growth phenotype are expected to be compensated for by duplicate genes. Gu et al. (2003) found a higher probability of phenotypic compensation for a gene with two recognizable copies than for singletons, and higher sequence similarities between duplicates correlated with higher frequencies of compensation. Interestingly, when the more highly expressed paralog was deleted, there was a higher probability of a severe fitness effect. This observation points to the existence of absolute or relative dosage requirements for many duplicated genes in yeast. Dominant Negative Effects and Retention of Paralogs. (A) Punnett square showing all possible combinations of dimers produced from three normal A alleles and one mutant allele A′ that renders the dimer inactive. The presence of only one mutant allele out of four leads to 56% of normal dimer AA. (B) Passage through a reduced ploidy phase is even worse because only 25% of dimers will be active. This fraction corresponds, in genetic terms, to the proportion of homozygotes AA, which are favored by selection. Removal of the deleterious A′ allele from the population can be extremely effective during this phase if fitness depends tightly on AA. Similar results are obtained for the trimer A-B-A. Dosage Effects and Retention of Paralogs. (A) In the context of the homodimer AA, consider that the synthesis of A results from the translation of the mRNA(A), present at a stable concentration, then A undergoes reversible dimerization or is degraded with a rate proportional to the concentration of monomers. According to this scheme, the concentration of AA in the steady state (ss) is proportional to the square of the ratio of specific rates of synthesis of A (Sa) and its degradation (Da; derivation in Veitia, 2004). Deletion of one allele out of four leads to 56% of normal AA. The graph shows the effects of deleting one, two, or three alleles expressed as the percentage of active AA with respect to the normal tetraploid system (the dotted line represents 50% of concentration of dimer). In a diploid phase, obliteration of one copy of A leads to 25% AA. Results are different if, for instance, AA becomes the degradable entity, as in a fast irreversible dimerization. When AA is degraded with a rate proportional to its concentration, a change in monomer input will have a proportional impact on the concentration of the dimer. (B) What happens to the trimer A-B-A in the steady state? Under similar assumptions as above, the yield of ABA depends again on the square of Sa. The black bars in the graph show the steady state amounts of ABA for different doses of the wild-type A allele (for the maximum/tetraploid dose of B). The gray bars show how deleting one allele of B increases dosage sensitivity of A. Again, the results are slightly different if the dimers or trimer are allowed to be degraded, and less dramatic effects are obtained when ABA is the only degradable entity. There is increasing evidence for nonequivalent contributions of the two parents to gene expression in allopolyploids and hybrids. A classical example is nucleolar dominance, in which one parental set of rRNA genes is silenced in an interspecific hybrid or allopolyploid (Pikaard, 2000). For instance, in A. suecica, an allotetraploid of A. thaliana and Cardaminopsis arenosa, there is silencing of the rRNA genes contributed by the A. thaliana subgenome (Chen et al., 1998). A similar phenomenon has been observed for many genes in other polyploids. For example, a recent analysis by Adams et al. (2003) has shown that expression levels of paralogous genes in synthetic allopolyploids simulating the ancestor of the most widely cultivated cotton, Gossypium hirsutum, were not additive. That is, the parental subgenomes did not contribute equally to the transcriptome of the polyploid. Indeed, expression of each paralog varied among organs and was suggestive of differential developmental regulation. Interestingly, the expression patterns in the synthetic polyploids can recapitulate those found in the naturally occurring cotton (Adams et al., 2003). This suggests that certain patterns of expression appear just after a polyploidization event and persist over time. Classical epigenetic mechanisms, such as altered DNA methylation or chromatin structure, have been proposed to explain this phenomenon (Kellogg, 2003). Heritable changes, such as DNA methylation, deletions, and chromosomal rearrangements, have been documented in other cases of allopolyploidy (Chen et al., 1998; Ozkan et al., 2001; Pires et al., 2004). Alternatively and not exclusively, sudden changes in transcription can be attributable to a regulatory mismatch between the effectors and the target genes contributed by the merging genomes. This mismatch can arise from the divergence of the relevant regulatory mechanisms in the parents and/or from global effects of the genome on gene regulation. Potential Impact of Allopolyploidization or Hybridization on Transcription. (A) Effect of an allopolyploidization event involving the two parental cells P1 and P2 on transcription from homologous promoters responding to the transactivator A. The black triangles represent the binding sites for A. (B) In the simplest case, the parental cells have similar volumes (Vol) and similar DNA contents but might contain different amounts of A and promoters with different thresholds (AThr). After the merger, let us suppose that the volume doubles and that the regulation of the expression of A is unaltered. It is easy to see that, in this particular case, the expected concentration of A is intermediate. For instance, this concentration can be suitable to trigger expression from the promoter coming from P1 but not from P2. (C) Transcriptional responses from promoters with three or four binding sites allowing cooperative interactions with the transactivator A. The solid curves are redrawn from Figure 5 of Veitia (2003b), with the same parameters. Decreasing the number of similar binding sites (n) shifts the position of the threshold to higher concentrations of A. In the case of a newly formed allopolyploid, if the resulting concentration of A is higher than the thresholds of both types of promoters (case A1), transcription is expected from both. If, on the contrary, the concentration of A is not high enough (i.e., A2), the red promoter will be on, while the green one will be off. Of course, cases of incomplete silencing are expected depending on the actual A concentration with respect to AThr. The thin black sigmoid represents a promoter with four binding sites that produces a response similar to that with three binding sites. It can be shown that there is an infinity of ways to obtain such a result, simply by changing the strength of the interaction between A and its binding sites, as long as cooperativity is maintained. This can explain the outcome of compensatory mutations. At another extreme, the transcriptional responses of the parents can be very similar, while the underlying system's promoter-activator can be different. Indeed, similar responses are obtained with mutations decreasing (or increasing) the number of binding sites compensated by increasing (or decreasing) the affinity of A for DNA (Figures 3B and 3C). Accordingly, an allopolyploid or hybrid can carry a mixture of A molecules that have different affinities for DNA in a volume that is different from that of the parents. Several outcomes are possible, but as a rule the promoter with a higher number of binding sites is more likely to be triggered by preferential interaction with A molecules having the highest affinity for DNA. In these simple cases involving only one activator, the promoter having the lowest AThr, whatever its parental origin, is more likely to remain active after the genomes merge. Although this does not involve any mark on the chromatin, it provides a mechanism allowing transmission of regulatory patterns. In the case of nucleolar dominance, chromatin modifications seem to be involved in epigenetic transmission (Chen et al., 1998). In this case, the factor A might attract a histone acetylase to promote transcription of the dominant genotype, or it might somehow induce repression of the silent genotype. The results of Adams et al. (2003) show that the relative expression of many paralogous gene pairs is dependent on tissue, organ, and/or developmental stage. This complexity could easily stem from the interaction of many loci, inhibitors, and activators that have an impact on the same promoter (Birchler et al., 2001). Indeed, this must also be the case for nucleolar dominance because a model relying on a promoter responding to a single factor cannot account for all the findings described by Chen et al. (1998), but this analysis is beyond the scope of this essay. Sudden alterations in gene expression patterns in a newly formed allopolyploid are not likely to be adaptive or optimal. However, the evolutionary persistence of a similar pattern of expression in a polyploid, as noted by Adams et al. (2003), is suggestive of some underlying selective advantage or constraint. The emergence of a new pattern of expression provides a further explanation for the retention of the paralogs, even if the initial event was a regulatory mismatch. Genes involved in the same pathway or macromolecular complex are often coregulated, responding basically to the same regulators. If a regulatory mismatch arises after it is expected that genes from each paralogous pathway or complex will undergo a change in their regulation interacting genes from one genome might have each of these paralogous or complexes will to a copy of each of their dosage Thus, an alteration in the mechanism of can a simple explanation for both the coordinated divergence of duplicated or complexes and paralog This would be to one paralogous system from the of selective between duplicated can be if they are expressed in different or at different developmental (i.e., This process could be thought of as an of and contribute to the of polyploids to A. a has at least three of polyploidy et al., 2003). and have patterns of functional divergence of paralogs formed during the most recent polyploidy found of retained paralogs encoding proteins involved in and paralogs of genes have been Another genomic analysis in the same plant has shown that the same of genes have a higher probability of retained from one of genome duplication to another and 2004). According to and more than 50% of the paralogous pairs in A. thaliana different patterns of expression. Similar results have been obtained for gene duplicates in and and Wolfe, suggesting that this might be a more in the context of A. genes expression were which are likely to from the paralogs in the same and Wolfe, 2004). The results in terms, with those of et al. found more paralogs than expressed ones in A. this also found that of duplicated gene pairs expression between promoters of duplicated genes, which suggests that the process of transcriptional neofunctionalization and is to a fraction of or that we are to expression and subtle sequence and through a are complex and similar can be obtained with different combinations at different In such compensatory mutations in the responses while allowing the existence of high levels of in the As a phenotypic and can be obtained when selection The components of such can in levels or patterns of expression, if they to the deleterious effects of dosage the of this process have been For instance, et al. have how the number of interactions a protein has within a genetic of gene expression in yeast and found that gene expression was correlated with the number of protein–protein the expression levels of interactors were correlated and the extent of in expression among genes encoding interactors was than that of random pairs of genes. This might result from the of balance between interactors within complexes or among complexes of the same It is that selection might lead to the retention of specific of at interacting involved in two paralogous so as to the pathways. That is, before selection the paralogous to similar in basically the same After the will to different in different ways and This can be with a process of not of a single gene but of a the simple case of components poorly connected within the cellular interaction network, when selection on one of a pathway or complex involved in a with other with its is This explain retention (or of or or complexes through of polyploidization. in the fate of one paralog must be considered from the perspective of the selective advantage (or by the to the and especially for their and on many have the and of one to and were also of this for with a of on several that to the