- © 2011 American Society of Plant Biologists
Abstract
The abundance of detected ancient polyploids in extant genomes raises questions regarding evolution after whole-genome duplication (WGD). For instance, what rules govern the preservation or loss of the duplicated genes created by WGD? We explore this question by contrasting two possible preservation forces: selection on relative and absolute gene dosages. Constraints on the relative dosages of central network genes represent an important force for maintaining duplicates (the dosage balance hypothesis). However, preservation may also result from selection on the absolute abundance of certain gene products. The metabolic network of the model plant Arabidopsis thaliana is a powerful system for comparing these hypotheses. We analyzed the surviving WGD-produced duplicate genes in this network, finding evidence that the surviving duplicates from the most recent WGD (WGD-α) are clustered in the network, as predicted by the dosage balance hypothesis. A flux balance analysis suggests an association between the survival of duplicates from a more ancient WGD (WGD-β) and reactions with high metabolic flux. We argue for an interplay of relative and absolute dosage constraints, such that the relative constraints imposed by the recent WGD are still being resolved by evolution, while they have been essentially fully resolved for the ancient event.
INTRODUCTION
Duplication of genetic material is an engine of evolutionary innovation (Taylor and Raes, 2004). In particular, Ohno (1970) argued for an important role for whole-genome duplications (WGD) in the shaping of the genetic and morphological complexity of multicellular eukaryotes. While such associations are of considerable interest, it is quite difficult to associate changes in rates of diversification or complexity with particular WGD events (Sémon and Wolfe, 2007). This difficulty is compounded by the fact that most duplicate genes produced by WGD are quickly lost (Scannell et al., 2007). As a result, understanding the interplay of duplicate gene loss and adaptation after WGD is an important open problem in evolutionary biology.
The model plant Arabidopsis thaliana (thale cress) is an attractive system for studying this question because, despite its small size (157 Mb), the Arabidopsis genome includes the remnants of at least three ancient WGDs (Vision et al., 2000; Simillion et al., 2002; Blanc et al., 2003; Bowers et al., 2003; Maere et al., 2005; Barker et al., 2009). These WGDs (Bowers et al., 2003), termed α (the most recent event, with 3947 retained duplicate genes in the genome), β (an intermediate event in time with 2765 retained duplicates), and γ (the most ancient event, 771 retained duplicates), have been validated by both gene synteny and comparative duplicate age (Ks) analyses (Maere et al., 2005). While the γ event is quite widely distributed, having occurred near the origin of all eudicots (Soltis et al., 2009), the two more recent events (α and β) are known to be restricted to the order Brassicales, as papaya (Carica papaya), a member of this order, possesses neither (Soltis et al., 2009). Chronologically, α and β are estimated to have occurred roughly 47 to 65 and 65 to 115 million years ago, respectively (Beilstein et al., 2010). Here, we explore the patterns of duplicate preservation after WGD to gain a better appreciation of the evolutionary role of such duplications. That their influence is important is suggested by several observations. For instance, recent studies have demonstrated that the two recent WGD events gave rise to novel pathways for the synthesis of indole and Met-derived glucosinolates (see Kliebenstein, 2008; Schranz et al., 2011, and references therein) and may have increased species diversification rates (Beilstein et al., 2010). Likewise, it has been argued that repeated WGD events have contributed to the great diversity of angiosperm lineages (De Bodt et al., 2005; Soltis et al., 2008; Fawcett et al., 2009).
While the loss of duplicated genes after polyploidy may be rapid (Scannell et al., 2007), it is also distinctly nonrandom. Similar types of genes (kinases, transcription factors, and ribosomal proteins, among others) have been retained in duplicate after independent WGD events in ancestors of modern Arabidopsis, rice (Oryza sativa), paramecium (Paramecium tetraurelia), bakers’ yeast (Saccharomyces cerevisiae), and maize (Zea mays) (Seoighe and Wolfe, 1999; Blanc and Wolfe, 2004; Seoighe and Gehring, 2004; Tian et al., 2005; Aury et al., 2006; Conant and Wolfe, 2008; Schnable et al., 2009; Wu and Qi, 2010). In addition, other genetic features have been shown to associate with increased frequency of duplicate survival after WGD: they include greater numbers of protein interactions (Guan et al., 2007; Hakes et al., 2007), higher levels of gene expression (Seoighe and Wolfe, 1999), and greater numbers of phosphorylation sites (Amoutzias et al., 2010). Collectively, these trends suggest a pattern of anticorrelation in the functional classes of duplicate genes created by the two major types of gene duplications: those from small-scale duplications (SSDs) and those created by WGD (Maere et al., 2005; Dopman and Hartl, 2007; Wapinski et al., 2007). In contrast with WGD duplicates, SSD-produced duplicate genes tend to have fewer protein interactions and not to have core regulatory functions (Maere et al., 2005; Freeling, 2009). Given this strong differentiation between the two duplication mechanisms, there is a need for theories explaining the rules by which duplicate genes survive, especially following WGD.
In fact, the study of duplicate genes is a flourishing field in molecular evolution, and numerous mechanisms of preservation have been proposed (Innan and Kondrashov, 2010). Two of the most important are neofunctionalization (Hughes, 1994), whereby one of the two duplicate copies acquires a new, selectively beneficial function, and subfunctionalization, whereby the functions of ancestrally multifunctioned genes are subdivided between duplicated copies by a variety of mechanisms (Force et al., 1999; Stoltzfus, 1999; Des Marais and Rausher, 2008; Innan and Kondrashov, 2010). As powerful as these explanations are for exploring SSD events, they are less useful for understanding WGD, both because we lack the detailed historical and functional knowledge to employ them at the genome scale and because they do not seem to account for the kinds of functional biases observed in WGD-produced duplicates (Freeling, 2009).
Instead, at least two (likely complementary) hypotheses have been proposed to explain which genes survive in duplicate after WGD. The first is the dosage balance hypothesis (DBH) also known as the gene balance hypothesis (Papp et al., 2003; Freeling and Thomas, 2006; Birchler and Veitia, 2007). It states that, in eukaryotes, there is selection operating to disfavor duplications of central network genes due to the imbalance in network stoichiometry that results (Maere et al., 2005; Conant and Wolfe, 2008; Edger and Pires, 2009). This situation is reversed in the context of the duplication of the entire genome, since in that case it would be the loss of a second copy of a particular gene that would introduce imbalances relative to the remaining, duplicated, genes. Genes such as ribosomal proteins, kinases, and transcription factors that have pervasive and extensive functional interactions with other biomolecules are precisely the classes of genes that the DBH predicts to be maintained in duplicate post-WGD, and empirical evidence indicates they are dosage sensitive (Birchler et al., 2001). Likewise, the observation that genes coding for proteins of high interaction degree are overretained is in agreement with the DBH. We refer to these retention patterns as relative dosage constraints to distinguish them from cases of absolute dosage constraint, the second proposed explanation for duplicate gene overretention post-WGD.
An absolute dosage constraint implies that an absolute increase in the concentration of a gene product is beneficial. Such selection has been observed for duplicate genes from a variety of organisms (Kondrashov and Kondrashov, 2006), including plants (van Hoof et al., 2001; Widholm et al., 2001). At first blush, it might appear that this type of selection would be agnostic as to the type of duplication employed. However, as Kacser and Burns (1981) demonstrated in their classic article, flux through a metabolic pathway is unlikely to be limited by the availability of a single enzyme, meaning that no single gene duplication is likely to dramatically increase flux. We have therefore argued that when a dosage increase is favored not for a single reaction but an entire pathway (such as glycolysis in yeast), WGD is a more rapid adaptation than is awaiting the necessary sequence of single gene duplications (Conant and Wolfe, 2007). Our observation is in keeping with a number of analyses that have shown an association between the yeast WGD and the propensity to ferment glucose even in the presence of oxygen (Blank et al., 2005; Piškur et al., 2006; Chen et al., 2008). In an elegant computational evolution experiment, van Hoek and Hogeweg (2009) found that overretention of certain enzymes was a recurrent feature of post-WGD evolution under (simulated) selection for rapid growth in a defined environment. This result suggests an interesting answer to Gould’s famous question of the results of replaying the tape of life (Gould, 1989): at least in certain environments, both phenotypic and genotype patterns can recur.
Here, we examine the above two hypotheses of relative dosage and absolute dosage in the context of the patterns of duplicated metabolic gene retention after two of the nested WGD events (α and β) from Arabidopsis. In addition, we consider the question of how repeatable post-WGD genome evolution actually is. Maere et al. (2005) have already shown that categories of duplicate genes surviving after the α and β differ in their functional distributions, and our results suggest that such differences are seen even within metabolism. On the other hand, Seoighe and Gehring (2004) found that, over the genome as a whole, genes duplicated in one WGD were also likely to survive in duplicate after a second. One can imagine differing hypotheses in this regard. van Hoek and Hogeweg’s results could imply a situation where repeated polyploidy gives rise to increasingly large but similarly constituted gene families. A similar pattern might be observed if the primary selection acting on WGD is internal, as implied by the DBH. In both cases, one argues that the selective environment post-WGD is relatively constant, meaning that one expects the same retention patterns. Differences would be expected if instead each polyploidy is a unique adaption to a local environment.
Using the primary metabolic network from Arabidopsis (de Oliveira Dal’Molin et al., 2010), we explored whether the WGD-produced α and β duplicates tend to cluster together in the metabolic network and whether genes associated with high flux metabolic reactions were more likely to remain duplicated. Our work follows others in finding a general, but not universal, trend for metabolic genes to survive in duplicate post-WGD at higher frequency than genes in the genome at large (Maere et al., 2005; Kliebenstein, 2008). We argue that natural selection immediately following WGD operates to maintain relative dosage balance in the metabolic network (i.e., clustering of duplicates in the network for the WGD-α duplicates), but, for the older WGD-β duplicates, selection instead operates on absolute gene product abundance. Thus, our results suggest that all polyploids may not be created equal.
RESULTS
Network Construction and Validation
We employed a published Arabidopsis primary metabolic network (de Oliveira Dal’Molin et al., 2010) to explore the patterns of post-WGD evolution in duplicated metabolic genes. A metabolic network is an abstract representation of the relationships between the biochemical reactions an organism is capable of (which are generally enzyme catalyzed) and the metabolites used or produced by those reactions (Wagner and Fell, 2001). In this case, we represent each metabolic reaction as a node in the network (a circle in Figure 1). Two nodes (reactions) are connected by an edge if they share a metabolite (lines in Figure 1). For Arabidopsis, two groups have produced models describing the metabolites and enzymes of primary metabolism (Poolman et al., 2009; de Oliveira Dal’Molin et al., 2010). The resulting networks include such pathways as photosynthesis, core anabolic and catabolic processes, and bulk biomass synthesis. Generally absent are secondary metabolic processes such as glucosinolate synthesis. The resulting network from de Oliveira Dal'Molin and colleagues consists of 1709 actual metabolites and 1550 metabolic reactions. The full network has 161,309 edges and the following network statistics: diameter, 11; average shortest path, 2.50; and density, 0.067 (Sabidussi, 1966; Coleman and More, 1983). There are 1398 genes annotated in the network, and 1168 reactions are associated with at least one gene. The metabolic network includes 334 reactions annotated with redundant genes (the redundant reactions in Supplemental Data Set 1 online): these are distinct reactions that share the same genes. Many of these genes are close paralogs that have undergone postduplication sub- or neofunctionalization (Table 1; Conant and Wolfe, 2008). For the purposes of this analysis, we thus employed a nonredundant 1217-reaction primary metabolic network created by merging nodes with the same genes and compartmental annotations but different metabolites (see Methods). This filtered network has 102,591 edges and the following network statistics: diameter, 8; average shortest path, 2.47; and density, 0.069.
Gene Numbers and Duplication Status in the Arabidopsis Primary Metabolic Network.
(A) and (B) For each pair of reactions (nodes), an edge is created if the two reactions share a metabolite. The number of associated enzyme genes (A) or SSD events (B) is indicated for each reaction by a color gradient. The observed mean value was set to the gradient mid color to highlight the actual distribution, shown at the top of the color gradient.
(C) The WGD-α and WGD-β duplication events are shown separately. Nodes with at least one such event are reported (the detailed distribution is also illustrated in the associated pie charts). White nodes are nodes without associated genes (31.5% of the reactions).
Selected List of Sub- and Neofunctionalized Duplicates
Duplication Mapping
Arabidopsis duplicate genes were identified and attributed either to the (recent) SSDs (245 duplicates) of Rizzon et al. (2006) or to WGD, as identified by Bowers et al. (2003) and modified by Thomas et al. (2006). Duplications attributed to genome duplication were further subdivided as surviving from the early β duplication (156 duplicates) and/or from the subsequent α event (420 duplicates). We then individually mapped these distinct types of duplicate genes onto the reaction-centered metabolic network (Figure 1; see Supplemental Data Set 2 online).
Propensity for Duplicate Preservation after WGD Varies by Functional Role
We first evaluated the frequency of duplicate genes (of both SSD and WGD origins) in the metabolic network compared with the genome at large. There are more surviving duplicate genes from the metabolic network than would be expected: duplicates produced by SSD are observed 29% more often than expected (compared with the duplication level observed in the overall genome), while the corresponding figures for the α and β events are 30 and 12%, respectively (Fisher's exact test, P < 0.001). We were further curious whether duplication propensity was associated with the cellular compartment in which an enzyme acts, especially given that these compartments are specialized for differing types of pathways. Furthermore, two of the compartments (the mitochondria and the chloroplast) possess their own genomes whose ploidy level was unaltered by nuclear WGD events (Figure 2). Except for a deficit of duplicates annotated as transporters (P < 0.003), the proportion of genes retained following WGD does not differ significantly among the compartments. By contrast, the frequency of SSD was significantly different among the compartments, with the highest frequency of events seen in the nuclear-encoded fraction of the enzymes active in the mitochondria and chloroplast. No SSDs of transporters were detected. Similarly, Table 2 summarizes the frequency of duplicate retention after the WGD-α and WGD-β events for genes of differing functions, inferred using the Gene Ontology framework (Carbon et al., 2009). Retention rates for genes from the metabolic network are high relative to the genome at large but are lower than that for transcription factors, a class of genes known to be highly retained post-WGD (Maere et al., 2005).
Duplication Frequencies by Cellular Compartment.
For each compartment, the proportion of reactions with duplicate genes surviving from SSD and WGD events is shown. The diameter and shade of gray of each pie is proportional to the frequency of duplicates (inner numbers). The proportion of surviving WGD-produced duplicates is not significantly different between compartments save for the transporters (P < 0.003). The cellular compartments fall into three distinct groups/categories with respect to the frequency of SSD events (bottom brackets): [Mitochondria, Chloroplast], [Cytosol], and [Peroxisome] (P < 0.001).
Retained Functions after WGD
Duplicated Genes Form Clusters in the Metabolic Network
Figure 1 suggests that duplicated metabolic genes are clustered nonrandomly in the network. We used a novel cluster detection algorithm to assess this apparent pattern (see Methods; Bekaert and Conant, 2011). With this algorithm, we analyzed five distinct sets of duplicated genes: (1) SSD duplicates, (2 and 3) genomic loci with duplicates derived from the α and β events (WGD-α and WGD-β, respectively), (4) genomic loci with duplicates surviving from both events (WGD-α∩β), and (5) a pooled data set of loci with duplicates from either WGD event (WGD-α∪β). To analyze each set, we first removed from the network all nodes not possessing genes from that set. We then calculated the number of connected components (Watts and Strogatz, 1998) among the remaining nd nodes with appropriate duplicates. We compared this number to the number of components obtained by retaining nd nodes selected at random from the network (see Methods; Figure 3). For the sets WGD-α (Figure 3B) and WGD-α∪β, we found significantly fewer and larger clusters than expected (Figure 3D; see Supplemental Table 1 online). These two duplicate gene sets also have higher in and out degrees than would be expected (i.e., their products are more likely to be reactants in other reactions, while their reactants are produced by more other reactions than those of the average reaction). While the WGD-β set (Figure 3C) does not show a significant bias in clustering or node degree compared with the remainder of the network, this may be due to the small number (84 nodes total) of known β duplicates: the pooled set of α and β duplicates (WGD-α∪β) show strong biases. We see no tendency for the set of genes produced by SSD to cluster in the network (Figures 1B and 3D). Note that we see less clustering than expected among the reactions catalyzed by only a single gene (Figure 3D). This result is expected because our analysis is essentially zero sum: an excess of clustering among some reactions will be balanced by a deficit among other classes. This fact should be clear when one recognizes that our randomization approach compares clustering in a subset of the network to the average clustering in the whole network (given by the random subsets created for statistical analysis).
WGD-α Duplicate Genes Cluster in the Metabolic Network.
(A) Overview of the cluster detection method. The number and size of the connected components (shaded clusters in the figure) for the real network are calculated after nonduplicated nodes (white) are removed. These clusters are then compared with those seen in randomized networks with the same number of duplicated nodes (see Methods).
(B) The clusters observed for the 247 reactions with retained WGD-α duplicates (dark nodes). The inset shows a detail of the chloroplast glycolysis. EC 2.7.7.27, glucose-1-phosphate adenylyltransferase; EC 2.7.1.11, ATP:d-fructose-6-phosphate 1-phosphotransferase; EC 2.7.2.3, ATP:3-phospho-d-glycerate 1-phosphotransferase; EC 6.3.1.2, glutamate-ammonia ligase; EC 2.7.1.40, ATP:pyruvate 2-O-phosphotransferase.
(C) The clusters observed for the 84 reactions with retained WGD-β duplicates (dark nodes).
(D) Results of the clustering tests for the five data sets described in Results. Significant P values (one-sided tests, α < 0.05) are shown along with the direction of the test (higher or lower than the expectation). Thus, the WGD-α duplicates show larger clusters than seen in the randomized networks (P < 0.05). Comparisons indicated with a dash were not significantly different from randomized networks.
Enzymes Duplicated in Ancient Polyploidy Events Are Associated with High Flux
A natural question is how the flux an enzyme carries influences its chances of surviving in duplicate after WGD. To explore this question, we used flux balance analysis (Orth et al., 2010) to estimate the flux through every reaction in the Arabidopsis primary metabolic network under both photosynthetic and nonphotosynthetic conditions (see Supplemental Figure 1 online). The published network does not allow the synthesis of either DNA/RNA nucleotides (ATP, dATP, etc.) or the amino acids Arg and His. However, we were able to model the biosynthesis of these molecules by adding a transport reaction that imports succinate in the mitochondria in exchange for fumarate (Catoni et al., 2003). This transporter is encoded by the known Arabidopsis gene At5g01340 (Catoni et al., 2003). Using these flux values (covering 40% of all reactions; see Supplemental Figure 1 online) and the above five sets of duplicate genes (SSD, WGD-α, WGD-β, WGD-α∩β, and WGD-α∪β), we built logistic regression models to evaluate whether reaction flux is predictive of membership in a duplication set. Among these sets, high metabolic flux is strongly predictive of membership in WGD-β (P ≈ 0.002; Figure 4); no other set shows a significant association between flux and duplication status.
High Flux Reactions Are Enriched for Surviving WGD-β Duplicate Genes.
The association between metabolic flux (x axis) and the likelihood of a surviving WGD-produced duplication (y axis) as predicted by logistic regression (see Methods). Thus, high reaction flux is a predictor of the presence of a WGD-β, but not a WGD-α, duplicate. *Significant P values (α < 0.05).
DISCUSSION
The degree to which genome duplication is associated with functional or morphological innovation is a long-standing evolutionary question (Kliebenstein, 2009). Ohno (1970) was a strong proponent of such a relationship, but direct links between the two types of event are difficult to establish. Nonetheless, a number of authors have speculated that repeated polyploidy events, especially in plants, might have driven both speciation and increases in morphological complexity (Freeling and Thomas, 2006; Scannell et al., 2006; Edger and Pires, 2009; Wood et al., 2009). One of the few cases where the details of such links are available is in S. cerevisiae, where it has been argued that a WGD helped give rise to yeast’s propensity to ferment glucose even in the presence of oxygen (Blank et al., 2005; Conant and Wolfe, 2007; Merico et al., 2007; Chen et al., 2008). Importantly, this adaption appears to have arisen through selection for high absolute dosage among metabolic genes, even though relative dosage constraints appear to have been active in other parts of the genome (Seoighe and Wolfe, 1999). Here, we tested the roles of relative dosage constraints and absolute dosage in Arabidopsis in influencing the fates of two classes of gene duplications: ancient polyploidy events (WGD-β and -α) and tandem duplications, with the prediction that they would differ (Edger and Pires, 2009; Freeling, 2009).
Is it possible that, as in yeast, absolute dosage selection was active at the time of the WGD-β in Arabidopsis? Figure 4 suggests that high flux reactions have been overretained duplicated after this event, and there should be a strong association between high flux and the need for multiple copies of the associated enzymes. However, no such association is detected for the more recent α duplication (Figure 4), suggesting, at a minimum, that selection on absolute dosage is not the dominant force in all polyploidy events. And of course, flux high should not be read as a synonym for gene essentiality or importance: some high flux reactions can be bypassed by alternative metabolic routes (Wagner, 2005), while low flux reactions may be essential. Alternatively, we suggest that both absolute and relative dosage constraints have acted during the post-WGD evolution of Arabidopsis. It is intuitive that, under the DBH, clusters of interacting metabolic genes would tend to either all remain duplicated or all return to single copy. Such relative dosage constraints are exactly the pattern we observe after the α event (e.g., the clustering in Figure 3). We suggest that the pattern of flux distribution among WGD-β duplicates can also be partly understood in terms of the DBH, despite the fact that little clustering is seen among these duplicates.
Because polyploidy imparts an immediate and relatively predictable increase in cell volume (Galitski et al., 1999; Leitch and Bennett, 2004; Veitia, 2005), it is reasonable to argue that, on average, mRNA and protein concentrations change relatively little after genome doubling. Thereafter though, the rapid and nonrandom loss of duplicated genes after WGD (Scannell et al., 2007) should alter the dosages of the new genome’s genes relative to each other (Conant and Wolfe, 2007). We argue that selection on relative dosage balance is responsible for the immediate preservation of duplicates (consistent with the clustering results discussed above). One caveat to this argument is that it is at least possible that the lack of clustering among the single-enzyme reactions (Figure 3C) is the true evolutionary signal from our analysis and not the clustering among the duplicates (see Results). From this perspective, rather than selection to maintain duplicates with many interactions after WGD, there is actually selection to return certain enzymes to the single copy state (Paterson et al., 2006; Rong et al., 2010). We are skeptical of this argument, however, because it is unclear why such selection would be in force.
While the DBH is an attractive unifying theory, its predictive power is not expected to be absolute (Veitia et al., 2008). One reason is that much of the cellular regulatory apparatus that acts to buffer gene expression noise can also partly buffer ploidy changes (Stelling et al., 2004; Raser and O’Shea, 2005). It seems plausible therefore, that, over the longer term, changes in gene expression and regulation (Wray et al., 2003; Prud’homme et al., 2007; Wagner and Lynch, 2008) will release relative dosage constraints (Carroll, 2000; Lynch and Wagner, 2008; Veitia et al., 2008; Birchler and Veitia, 2010). However, at least for metabolic genes, because relative expression levels remain tuned to the current genome, there may be a set of high flux enzymes originally preserved by selection on relative dosage but for which absolute dosage selection is now acting. Absolute dosage selection will oppose the loss of a duplicate copy in cases where reduced gene expression interferes with the high enzyme copy number needed to maintain flux. Again, note that the history of WGD means that these dosage constraints are acting at the pathway and not the single gene level, meaning that we do not necessarily expect further duplication-based dosage amplifications. Examples of selection on gene copy number to achieve high enzyme activity are known from a number of species (Brown et al., 1998; Guillemaud et al., 1999; Price et al., 2004; Perry et al., 2007). Note that the argument is not that all genes behave in this dosage-dependent manner; merely that it is a general trend. This hypothesis of two types of dosage selection is attractive because it explains the difference in flux bias seen between the α and β duplicates: the α event is not yet fully resolved, meaning that enzymes other than those of high flux still survive due to selection on relative dosage balance, while the β duplicates are being acted on primarily by (late-arriving) selection on absolute dosage.
Once relative dosage constraints are lifted, duplicate genes have a number of alternate fates beyond absolute dosage selection. These fates include gene loss as well as subfunctionalization and neofunctionalization. We previously illustrated an example of post-WGD neofunctionalization (Kliebenstein, 2008; Schranz et al., 2011): the ability to synthesize two major classes of glucosinolates, which gives members of the order Brassicales their unique sharp taste, arose through the α and β WGDs. These compounds are synthesized from novel substrates by retained α or β duplicates, and the novelty is clear from the retained ancestral function in the other duplicate. Here, we identified several other cases of putative neo- and subfunctionalization post-WGD (Table 1). In these cases, functional genomics, gene expression and knockout, and protein localization data indicate that functional divergence has been layered on top of conserved enzyme structures and reaction mechanisms. In this view, the dosage balance hypothesis provides one absolutely key ingredient for functional change: time. By preserving duplicates, it allows a sufficient period for the necessary function-altering mutations to appear, meaning that the relaxation of relative dosage constraints represents the beginning, not the end, of the evolutionary impact of WGD.
An important caveat is that our study did not examine glucosinolates or any other potentially ecologically important secondary metabolic processes. Genome-scale metabolic models allow us to explore the potential for very interesting links between metabolism, genome duplication, and evolution. However, these models do not include most secondary metabolites. One might think this omission would tend to artificially lower certain fluxes in the network and hence mislead us when looking for links between metabolism and evolution. However, even the relatively abundant glucosinolates represent <2% of biomass for adult Arabidopsis plants (Brown et al., 2003). This proportion is an order of magnitude less than the fraction of protein by mass (Poolman et al., 2009). Since the highest flux reactions considered here are not even biomass synthesis reactions but those for energy production, it is unlikely that including the synthesis of secondary metabolites would greatly alter the numeric flux values considered here. Indeed, even the artificial doubling of flux through amino acid synthesis pathways has no discernable effect on the results presented in Figure 4 (see Supplemental Table 2 online). The more serious concern is that omitting such evolutionary labile processes gives a false sense of stability as to the role of WGD in metabolic evolution (Kliebenstein, 2008). It will be very interesting in future studies to collect this missing metabolic data to explore the potential for functional changes in secondary metabolism enzymes duplicated at WGD.
Here, we tried to integrate our understanding of the roles of relative dosage and absolute dosage in driving evolution after genome duplication. It is still an open question whether the WGD events in the ancestors of Arabidopsis conferred some immediate benefit that led to their fixation. What is clearer is that these events profoundly affected the subsequent evolution of this lineage. One effect that should not be underestimated is the role of WGD in relaxing the epistatic constraints all genomes acquire over their history. Subfunctionalization is one of the best illustrations of this principle: the divergence between the genes in Table 1 occurred post-WGD, and it is intriguing that such changes could take place even in this relatively limited window of divergence. Whether such relaxation could be a route to the increased phylogenetic and phenotypic diversity thought to characterize polyploid lineages is an exciting open question.
METHODS
Data Collection
The complete Arabidopsis thaliana primary metabolic network v1.0 was obtained from de Oliveira Dal’Molin et al. (2010). It includes a list of metabolites and their respective cellular compartments as well as the biochemical equation for each reaction. In cases where the enzyme catalyzing a reaction is known, the encoding gene is also noted. Each reaction is defined as a node in our reaction-centric network. Edges between these nodes are defined by shared metabolites between the reactions (see Supplemental Figure 2A online). The network is directed: for irreversible reactions if the product of one reaction is a reactant in the second, we define an edge. Reversible reactions are treated similarly, except that both directions of the reaction are allowed and handled independently. A specific compartment was created for the transporters (which we defined as reactions having metabolites in two compartments). The networks were visualized with Gephi v0.7 alpha (Bastian et al., 2009) using the Force-based algorithm ForceAtlas. Our goal was to use this network to assign gene duplication events. To do so, we must account for the fact that several of the reactions are associated with the same set of genes (e.g., At1g04710, At2g33150, and At5g48880 all encode enzymes responsible two different reactions: the conversion of 2-methylacetoacetyl-CoA into acetyl-CoA and the conversion of 3-α 7-α dihydroxy-5-β 24-oxocholestanoyl-CoA into chenodeoxyglycocholoyl-CoA) (see Supplemental Figure 2B online). To account for these ambiguous annotations, we merged 334 such reactions (see Supplemental Figure 2B online) in order to not overestimate patterns of gene duplication and duplicate clustering. The result is the simplified metabolic network (with only 1217 nodes rather than 1550) that we use for mapping the gene duplications (see Supplemental Figure 2C online).
Identifying Enzyme Copy Numbers and Duplicate Gene Origins
The mean, median, and maximum number of genes associated with a reaction were 2.56, 1, and 84, respectively. For visualization purposes (Figure 1A), we use two color gradients in order to distinguish nodes with fewer (blue to yellow shades) or more (yellow to red shades) genes than the mean. To identify duplicate genes produced by SSD, we employed the data from Rizzon et al. (2006) and mapped the presence of SSDs (4043 events under a low stringency criteria with up to 10 spacer genes; Rizzon et al., 2006) onto the set of reactions (mean, median, and maximum number of duplicate genes were 4.30, 1, and 47, respectively). We plotted nodes with more or fewer SSD genes than the mean in Figure 1B. Duplicate genes produced by the α and β events were taken from the list given by Bowers et al. (2003) as modified by Thomas et al. (2006). We then mapped the presence of single copy genes or WGD-α or WGD-β duplicated genes onto the metabolic network (Figure 1C).
Retained Functions
We investigated the functional biases in the set of genes retained in duplicate following WGD. We collected the highest level of the Gene Ontology molecular function ontology (Ashburner et al., 2000). Genes might have more than one distinct function and therefore might be annotated with more than one Gene Ontology term. As a result, the total number of functional classifications is greater than the total number of genes. We calculated the significance of the enrichment with two-tailed Fisher’s exact tests under the null hypothesis that there is no bias in the proportion of duplicate genes in the specified category.
Clustering Tests
We were interested in whether reactions with the same history of duplications cluster in the metabolic network. We thus calculated the number and maximal size of the network components containing nodes with each tested event (Figure 2C). To assess whether these components were bigger than would be expected, we used network randomization. We begin by copying the original network and reassigning the gene number or duplication status at random (Figure 3A). We then computed connected components for the random networks. We performed 10,000 permutations and used the distribution of component sizes to determine whether the clusters in the real network were larger than expected. The procedure was implemented in C++ using the Boost Libraries (http://www.boost.org/).
Flux Balance Analysis
We used the Systems Biology Research Tool v2.0.0 (Wright and Wagner, 2008) to perform flux balance analysis on the Arabidopsis primary metabolic network. Because the original network lacked the ability to synthesize nucleotides, DNA and RNA were not included in the biomass reaction (see Results). Thus, the nucleic acid requirements used here were calculated based on the weight of the each nucleotide and the composition of the Arabidopsis genome. Fluxes computed with this tool were also verified with our own implementation of flux balance analysis. Using the network and the stoichiometry provided (de Oliveira Dal’Molin et al., 2010), we estimated the optimal biomass production rate under photosynthetic conditions (photon import allowed and sugar imports forbidden) and under nonphotosynthetic conditions (photon import forbidden but sugar imports allowed). In each case, we also made every possible reaction knockout whereby the reaction flux is constrained to 0, and the remainder of the network is reoptimized (see Supplemental Figure 1 online). For each such knockout analysis, all fluxes were first normalized by the value of the computed biomass flux for that analysis. For each reaction across all analyses, we then selected the maximum flux (i.e., across all possible conditions). We compared this maximal flux to the duplication status of that reaction node.
Logistic Regression
A logistic regression model (Sokal and Rohlf, 1995) was used to evaluate the relationship between flux and gene duplications. In this framework, reaction flux is used as a predictor of a binary duplication status variable (i.e., to infer whether a given reaction, drawn at random, possesses a duplicate given its flux value). We can use a likelihood ratio test to ask whether adding flux information significantly improves our ability to make such a prediction as opposed to simply using the overall duplication frequency as our predictor. The analysis was implemented in the package mLogit v0.18 in R.
Supplemental Data
The following materials are available in the online version of this article.
Supplemental Figure 1. Metabolic Flux Maps.
Supplemental Figure 2. Initial Metabolic Network and Processed Network after Merging Redundant Reactions.
Supplemental Table 1. Clustering Test Results.
Supplemental Table 2. Amino Acids Drain the Secondary Metabolism.
Supplemental Data Set 1. List of All the Case of Reactions Encoded by the Same Set of Genes.
Supplemental Data Set 2. Reactions’ Duplication Status and Flux.
Acknowledgments
We thank T. Arias, J. Birchler, C. Hudson, D. Mayfield, and J. Schultz for helpful discussions. This work was supported by a Research Board grant from the University of Missouri (M.B.), the Reproductive Biology Group of the Food for the 21st Century program at the University of Missouri (G.C.C.), and by the U.S. National Science Foundation (DBI 0501712 and DBI 0638536; P.P.E. and J.C.P.).
Footnotes
The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Michaël Bekaert (bekaertm{at}missouri.edu).
↵[W] Online version contains Web-only data.
- Received November 15, 2010.
- Revised April 7, 2011.
- Accepted April 15, 2011.
- Published May 3, 2011.
References
