Skip to main content
  • Research article
  • Open access
  • Published:

Complex fate of paralogs

Abstract

Background

Thanks to recent high coverage mass-spectrometry studies and reconstructed protein complexes, we are now in an unprecedented position to study the evolution of biological systems. Gene duplications, known to be a major source of innovation in evolution, can now be readily examined in the context of protein complexes.

Results

We observe that paralogs operating in the same complex fulfill different roles: mRNA dosage increase for more than a hundred cytosolic ribosomal proteins, mutually exclusive participation of at least 54 paralogs resulting in alternative forms of complexes, and 24 proteins contributing to bona fide structural growth. Inspection of paralogous proteins participating in two independent complexes shows that an ancient, pre-duplication protein functioned in both multi-protein assemblies and a gene duplication event allowed the respective copies to specialize and split their roles.

Conclusion

Variants with conditionally assembled, paralogous subunits likely have played a role in yeast's adaptation to anaerobic conditions. In a number of cases the gene duplication has given rise to one duplicate that is no longer part of a protein complex and shows an accelerated rate of evolution. Such genes could provide the raw material for the evolution of new functions.

Background

Gene duplication can be a major source of innovation in evolution [1], providing redundancy and additional genetic material to build upon and differentiate. In general, eukaryotic genomes contain a large fraction of gene duplicates, with paralogs stemming not only from single gene or segmental duplications, but, in the case of S. cerevisiae, also from a Whole-Genome Duplication event that occurred approximately 100 mln years ago (WGD; [2, 3]). Genomic instability and massive gene loss promptly followed WGD and purged most of the newly formed gene copies from the yeast genome, retaining approximately 10% of them [3]. Today, using multiple genomes of related fungal species with conserved synteny, we can unambiguously identify hundreds of gene pairs as WGD paralogs [4] in addition to normal small scale paralogs.

The identification of paralogs of WGD origin, in conjunction with the wealth of data on physical protein interactions and derived maps of protein complexes, puts us in an unprecedented position to test the fate of nascent duplicated genes and to potentially identify cases of duplication of whole complexes. Recently, it has been shown that, after gene duplication, protein interactions can be conserved [5, 6]. The data suggested that there exists a stepwise pathway of evolution for such functional modules [6], with duplications of homomeric interactions known to have a significant influence on the evolution of genes [5]. Moreover, it is known that gene duplicates can be found less often among the core components of protein complexes compared to sparse regions of protein interaction network [7]. For our study of the impact of gene duplication on protein complexes, we separated paralogs into two distinct, non-overlapping classes: genes that were duplicated at the WGD event, and non-WGD duplicates detectable by sequence similarity. Dubbed small scale duplications (SSD), these paralogs are the result of the most recent gene duplications, identified per event by employing a best bi-directional hit criterion on all yeast gene pairs (see Methods). From the analysis of the phylogenetic distribution and number of paralogs in related species, it appears that the time of duplication of SSD genes greatly predates the WGD event (see Methods). Both duplication types, WGD and SSD, cover together ~40% of yeast genes, providing a comprehensive overview of these evolutionary events. These two paralog types are already known to differ with respect to their expression pattern [8, 9] and synthetic lethality rate [10], by displaying different phenotypic effects when deleted [11] and occurrence across functional classes (e.g., stress responsive genes, [8]). Musso and colleagues [9] show that nearly half of WGD paralogs co-cluster in the same protein complex. Amoutzias and colleagues [12] indicate that whole genome duplication did not change the dimerization specificities of interacting homologs. Here, we show a much more detailed spectrum of evolutionary and functional fates of higher order protein complex subunits. This integrated overview, enables us to quantify the fates with respect to the duplication type and address questions related to protein specialization (subfunctionalization), as well as the emergence of novel functions related to complexes (neofunctionalization).

Our hypotheses were tested on various types of manually curated data: both complexes from MIPS consortium [13], and those annotated by SGD [14]. To avoid a possible bias introduced by manual curation, we also use computationally derived maps of complexes [15, 16], reconstruction of which was possible owing to recent mass-spectrometry studies [17, 18]. Integration of these datasets allowed us to systematically study the fates of all gene duplicates which are involved in protein complexes.

Results

The fates of duplicate genes in complexes

We carried out a systematic analysis of the fate of paralogs in protein complexes. From our first observations it became clear that the cytosolic ribosomal complex dominates the whole spectrum of gene duplications. In order to prevent this single protein complex to dominate our results, we analyze it separately (see Methods). The fates of other paralogs found within complexes fall into two other categories (Figure 1 and 2). Intra-complex paralogs (I) that are formed when both resulting genes remain within the same protein complex, whereas bi-complex paralogs (II) function within two separate complexes. The third class, which we define as overhangs (III), consists of subunits of complexes with a paralog possessing no association to a known protein complex whatsoever. SSD and WGD paralogs are equally divided over intra-complex and overhang classes, but differ with respect to the bi-complex class: many more SSD paralogs are present in two complexes compared to WGD paralogs (Figure 2b). We discuss this observation below.

Figure 1
figure 1

Complex fate of paralogs. a) Gene duplication and subsequent divergence, for cytosolic ribosomal proteins (cRP) followed by homogenizing gene conversion events. b) Impact of duplicated proteins on complexes. Intra-complex duplications include dosage increase, interacting homologs and module variants. Dosage increase requires many components of the complex to duplicate simultaneously (as in the case of cRP and the whole genome duplication). For interacting homologs, the two duplicated proteins become physically subunits of the complex (e.g., homomers turning into heterodimers after the duplication). In module variants only one of the two paralogs is present in the protein complex at a given time. Bi-complex paralogs operate in different protein complexes; two possible evolutionary routes are shown. Overhangs do not aggregate with other proteins in a non-transient manner, while their paralogs do.

Figure 2
figure 2

The roles of paralogs in protein complexes. a) Shaded areas mark a complex, dashed lines connect paralogs. I) Intra-complex paralogs: when both proteins participate in the same complex; ARG transcription complex includes an intra-complex duplication of genes encoding FUN80 and ARGR1 subunits. II) Bi-complex paralogs: two proteins are involved in different protein complexes; two small complexes are shown: zeta DNA polymerase complex (left) and delta DNA polymerase complex (right). Pair REV3/CDC2 are bi-complex paralogs. III) Overhangs: only one of the paralogs constitutes a subunit of a complex, while its homolog does not aggregate with other proteins in a non-transient manner; Vps4p ATPase transport complex. Here, CHM2 protein (a paralog of DID3) represents an overhang. b) Type of duplication and their contribution to protein complexes: left, whole genome duplication (cytoplasmic ribosomal proteins excluded), and right, small scale duplications. On the pie chart, fractions of all paralog pairs are denoted. Protein complex annotations after SGD consortium.

Intra-complex paralogs: retention is an important fate of paralogs within complexes

We observe a very strong preference for both duplicated proteins to function in the same module. Compared to a null model, where proteins are stochastically reshuffled between complexes, intra-complex paralogs are ~40-fold overrepresented (SGD modules, [14]). This preference is similar, and not statistically different for both duplication types (P = 0.97, chi-square test) and holds for other module definitions, including the computationally derived protein complexes from complex co-purification experiments (see additional file 1, Table S1). Paralog retention within the module is thus an important factor in shaping the map of protein complexes.

We thus recover the previously made observation that WGD and SSD paralogs are known to act within the ancestral protein complex after the duplication [7, 8]. Further analysis however revealed a wider spectrum of fates in which two paralogs can be involved in a single protein complex, as illustrated in the analysis of the essential yeast chromatin remodeling complex RSC. Owing to the availability of protein-protein interaction data [17, 18] we can distinguish between different modes of participation in a single complex. The first, a more "direct" bait-prey interaction mode, occurs when one protein was designated a bait and the other protein co-purified as a prey; this event is characterized by a high spoke value [17]. The second type of interaction, a prey-prey interaction, can be detected when two prey proteins were co-purified with the same bait in two independent purification experiments and does not provide the evidence for the two proteins to co-occur in the same protein complex at a given time. Hence, we were able to specify the following intra-complex fates:

Ia. mRNA dosage effect

Undoubtedly, the single most abundant fate of paralogs is the mRNA dosage increase of cytosolic ribosomal proteins (CRPs), as more than half of 200 intra-complex proteins are CRPs (see Table 1 and Methods for details). Their coordinate expression is primarily accomplished at the transcriptional level [19], signified by high expression correlation between the paralogs (our results show Spearman's rank correlation coefficient of 0.8). It is a markedly different type of evolutionary innovation from other classes, as paralogous CRPs, due to undergoing gene conversion [20], are highly similar in their protein sequence (see Methods).

Table 1 Paralogs in complexes.

Ib. Interacting homologs

This subclass consists of paralogs both present in a protein complex at the same time. Using protein-protein interaction data we identified 24 intra-complex paralogs with a bait-prey interaction type, signified by a high spoke value of least 5 (see Methods for details). This class is exemplified by the RSC3/RSC30 pair from RSC chromatin remodeling complex, known to form a stable heterodimer [21]. This kind of relationship between paralogs is likely to result from an ancestral homodimer, where a paralogous replacement of the dimer's components took place [5]. Strong positive co-expression (Spearman's correlation coefficient of 0.4), even though weaker than the tightly co-regulated CRPs, provides additional clues for simultaneous presence of both proteins in the functional module. Homomers undergoing this evolutionary route are probably the classic view on how two paralogs are involved in the same protein complex, as exemplified by the F1-ATPase alpha and beta subunits [22].

Ic. Module variants

This perhaps somewhat less explicitly recognized category embraces paralogs with a seemingly intrinsic contradiction: operating within the same "complex", yet never present together with only a prey-prey evidence for their interaction. Such mutually exclusive presence implies existence of different variants of the same complex. To assign proteins to the module variants class, we select intra-complex paralogs with no evidence of direct interaction. That includes paralogs never purified together, or with a negative spoke value (see Methods for details). Our analysis yields 54 intra-complex paralogs that belong to this class. Lower co-expression of these genes, likely resulting from the functional role undertaken by paralogs, confirms that these subunits are alternatively present in a module, thus not required to operate simultaneously (average co-expression Spearman's correlation 0.2, statistically different from other classes, one sided Wilcoxon ranked sum test, P < 0.02). More divergent expression also suggests a mechanism of control of complex activity by conditional assembly (analogous to just-in-time assembly for cell cycle complexes, [23]).

The literature provides many tantalizing clues to the conditions in which these alternative assemblies (module variants) may operate. For example, the WGD paralogs TPK1/TPK3 have a negatively correlated co-expression (Spearman's correlation coefficient -0.3), although they are both part of the cAMP-dependent protein kinase complex. Glucose-induced hyperaccumulation of cAMP was observed in exponential-phase glucose-grown cells of the TPK1-deficient but not the TPK3-deficient strain [24]. Moreover, investigation of mitochondrial respiration by in vivo 31P nuclear magnetic resonance spectroscopy showed the tpk1- and not tpk3-mutant, to be defective in glucose repression [24]. Another clue hinting at specialization for carbon source utilization comes from the WGD paralogs COX5a/COX5b. The two subunits are encoded by divergent sequences, but are functionally interchangeable forms of yeast cytochrome c oxidase subunit V [25]. COX5a/COX5b paralogs are oppositely regulated (Spearman's expression correlation -0.3) and are known to be expressed in a mutually exclusive manner under aerobic (COX5a) and anaerobic conditions (COX5b) [26]. Taking all module variants together, we observe their enrichment among WGD paralogs (P < 0.02) and also module variants are 2.5 less likely to be essential when considering viability of single-gene knockouts (P < 0.01, Table 2). This apparent redundancy of module variants in rich medium does not exclude their possible contribution to cell's survival in other growth media. To test this hypothesis we analyzed the data on growth rates of yeast deletion mutants in nine fermentable and non-fermentable substrates [27]. 12 module variants and none of interacting homologs show differential and partly complementary pattern of growth rates in various carbon sources (Table 3). Additionally, phenotype of single-gene deletions of module variants does not correlate between paralogs (average correlation coefficient ~0.1) when tested in hundreds chemical and environmental stress conditions [28]. Together, all this evidence leads us to suggest that the WGD event might have facilitated the evolution of anaerobic fermentation in S. cerevisiae via introduction of many specialized module variants.

Table 2 Type of intra-complex paralogs and viability of single-gene knockouts in rich medium.
Table 3 Deletion mutants of module variants exhibit differential growth patterns when cultured on various carbon sources.

Bi-complex paralogs: proteins functioning in different complexes

As opposed to intra-complex paralogs, where both proteins function in the same module, bi-complex paralogs each participate in distinct ones. Depending on the map of protein complexes, 44 or more genes fall into this category (see Table 1 and additional file 1, Table S2). We confirmed the lack of interaction between this type of paralogs with protein-complex purification data (only two out of 31 pairs were ever purified together, significantly less than intra-complex paralogs, Fisher exact test P < 4e-5, odds ratio 15).

Interestingly, for bi-complex paralogs, a significant difference between WGD and SSD duplicates can be seen. The majority of them are SSD duplicates (see additional file 1, Table S2). This strong bias, with SSD constituting more than 80% of the bi-complex class, contrasts with handful of WGD paralogs split between different complexes. We propose this to be an effect associated with the age of duplication. The lion's share of SSD paralogs not only predate the WGD event but are older than the divergence with S. pombe. While none of the eight post-S. pombe SSD duplications is bi-complex, three duplications are intra-complex (see additional file 1, Table S5), a hint that not the type of duplication (SSD versus WGD), but its age has a greater influence on the paralog's fate. Over extended evolutionary time since the ancient duplication of majority of SSD paralogs, many specialized (subfunctionalization), join or even established a new complex (neofunctionalization), ultimately leading to the bi-complex relationship. The conservative nature of interaction evolution after gene duplication is confirmed by the underrepresentation of bi-complex paralogs, compared to a null model where proteins are free to change their complex following the duplication (Table 1 and additional file 1, Table S2).

What evolutionary route would lead to the emergence of bi-complex paralogs? Two possible scenarios are shown in Figure 3a. Unfortunately, we do not have a map of protein complexes prior to gene duplications and we have to rely on the indirect evidence of ancestral state. Manual inspection of the five cases of WGD bi-complex paralogs indicates that their complexes overlap (e.g. still posses shared subunits) and the bi-complex paralogs tend to show some genetic redundancy. For example, the DPB4 subunit shared between epsilon DNA Polymerase and CHRAC complexes (see Figure 3b), is indication that before the WGD event the ancestral DPB3/DLS1 subunit functioned, similarly, in the two complexes (see additional file 1, "Tracing evolutionary past of paralogs"). Based on these clues we conclude that the most common evolutionary scenario for bi-complex paralogs is as follows: a single ancestral protein operated in two independent complexes and a duplication event allowed the respective copies to specialize and split their roles. This preferred route would be an indication of subfunctionalization, rather than of neofunctionalization.

Figure 3
figure 3

Evolution of bi-complex paralogs. a) Bi-complex paralogs (right) are potentially cases of sub-functionalization (upper route; bi-complex protein becomes bi-complex paralogous pair) neo-functionalization (single protein in which one paralogous copy is added to another complex). b) A possible evolutionary scenario for a gene duplication event involving epsilon DNA Polymerase and CHRAC complexes. Left, proposed pre-duplication relationships between complexes (ancestral DPB3/DLS1 gene is pictured here before its duplication). Right, the map of modules after a DLS1-DPB3 protein duplication (current S. cerevisiae protein complexes) suggests subfunctionalization.

Examples of whole-complex duplications

A dramatic result of numerous bi-complex duplications could be a duplication of a whole complex. Particularly, a whole-genome duplication, by simultaneous duplication of all subunits, provides the necessary material to make a "carbon copy" of complexes. Instead, we observed a statistical under-representation of bi-complex paralogs. Nevertheless, single events of whole-complex duplication can be identified. For example, the farnesyltransferase and geranylgeranyltransferase complexes both consists of two subunits: RAM1-RAM2 and BET2-BET4 respectively, originating as ancient, pre-WGD duplications of RAM1/BET2 and RAM2/BET4 ancestral genes. Although it is not know which was the original complex (or whether the ancestral complex served both functions) this suggests that two stepwise small-scale duplications occurred, ultimately copying the whole complex that went on to evolve distinct functions (see Figure 4a).

Figure 4
figure 4

Examples of a whole-complex duplication. a) Farnesyltransferase complex (left) prenylates a-factor mating pheromone and Ras proteins, is required for their membrane localization [48]. Geranylgeranyltransferase complex (right) catalyzes the prenylation of Rab proteins, is required for membrane attachment and their subsequent biological activity [49], and is indispensable for vesicular transport between the endoplasmic reticulum and the Golgi. b) Endomplasmic reticulum membrane translocons Sec61 (left) and Ssh1 (right).

Another case of whole-complex duplication involves a three-protein Sec61 complex (also referred to as a translocon, Figure 4b). This essential complex forms a channel in the ER membrane and mediates translocation of secretory and membrane proteins into the ER and also retrograde transport of misfolded proteins to the cytoplasm for degradation [29, 30]. The complex has duplicated in the course of evolution to form an Ssh1 translocon complex [31]. The Ssh1 complex, a result of small scale duplications, also functions in co-translational import to the endoplasmic reticulum (an essential paralogous subunit Sec61p plays a post-translational role as well), and is required for normal growth rates.

Overhangs – lone paralogs

The final class of paralogs are overhangs, proteins without an association to a functional module, but with a paralog known to be involved in a protein complex (Figure 1). For SGD protein complexes, we found 58 such proteins, with no significant difference in contributions of WGD and SSD duplication types for most of the protein complex maps (see additional file 1, Table S3). Validation with TAP protein complex purification data shows virtually no association of overhangs with their paralog's module (average interaction spoke value for overhangs is 0.06 compared to 2.6 for their "in-module" partners). Additionally, compared to their paralogs, less functional data about overhangs is available. Perhaps predictably, 11/58 overhangs genes are unnamed genes (i.e. not described in a scientific publication), compared to all of their paralogs being named (Fisher's exact test, P < 0.01). Naming roughly reflects the state of our knowledge about the gene, and we further observe absence of annotation in Molecular Function (P < 0.02) and Biological Process (P < 0.01) classes of GO.

To further validate the role overhangs play in the cellular processes we counted the essential genes (inviable null mutants) among them. Even after excluding unnamed genes from this analysis, we have only 4/48 essential overhangs, compared to 17/58 of their in-module paralogs (Fisher's exact test P < 0.01, odds ratio 4.5). This corroborates with the hypothesis of Hart et al. [15] that essentiality is a product of the protein complex, rather than the individual protein. We conclude that overhangs play a much less important role in cell biology, at least in the rich medium conditions in which most of the functional studies are performed.

We observed that overhangs are less constrained by evolution on the sequence level. For WGD overhangs, we compared amino acid identity levels of paralogs against their Kluyveromyces waltii ortholog (there is a single ortholog in K. waltii, as this species diverged before the WGD event). The amino acid sequence of overhangs diverges significantly faster compared to their in-module paralogs (34% vs 40% global amino acid identity, one sided paired Wilcoxon signed rank test P < 0.02). We therefore conclude that being a part of the protein complex imposes certain constraints on divergence, and the process of orphaning coincides with an increased rate of sequence evolution.

A higher rate of protein sequence evolution and almost complete loss of interactions with an ancestral protein complex are manifestations of rapid functional divergence. The orphaned proteins are involved in different cellular processes: e.g., an overhang SSD1 (suppressor of SIT4 deletion, YDR293C), interacts with a TOR pathway, and functions in sustaining cell wall integrity [32], while its paralog DIS3 is a catalytic component of exosome [33], also involved in mitotic control [34]. We measured the degree of function divergence of overhangs and their paralogs. Using semantic similarity based on Gene Ontologies (see Methods), genes were assigned values between 0 (for different function) and 1 (highly similar or identical function). We observe a rapid divergence of functionality for overhangs (additional file 1, Figure S1). This analysis hints to the overhangs as one of nature's methods to gene neofunctionalization.

Discussion and conclusion

For the paralogs participating in different complexes (bi-complex paralogs), we see a quantitative difference between duplicates of different age, with only a minority of bi-complex paralogs stemming from WGD. We attribute the higher representation of SSD paralogs to the time of the duplication. The mixture of functional data and the knowledge of their evolutionary history enabled us to reconstruct the evolutionary past of WGD paralogs. As bi-complex paralogs might have potentially undergone either neo- or subfunctionalization (see Figure 3a), we suggest, based on the examination of the association between complexes, that bi-complex paralogs could be examples of function specialization in the protein interaction network.

As observed in [8] there is no overrepresentation of whole modules being duplicated at the WGD event. A massive duplication is a unique opportunity for an organism to replicate components of its cellular machinery (e.g., protein complexes) and let it subsequently evolve independently, with each complex following its own evolutionary path. And even though it appears that gene pairs [35] and transcriptional network show features of partitioning into heavily intra-connected, but sparsely inter-connected clusters [36] at the protein complex level we did not observe large-scale duplications. Is it maybe that the ancestor of S. cerevisiae around 100 mln years ago had a chance to duplicate complexes as a whole, but missed the unique opportunity? Certainly the case of cytoplasmic ribosomes is an example of the ancestral yeast cell taking advantage of WGD event and doubling the subunit count in this protein complex. In fact, the completeness of the duplication of the cytoplasmic ribosome (both the small and large subunit) allowed the cell to maintain required equimolar concentrations of CRPs [37, 38] while doubling the gene repertoire, a goal not attainable by stepwise module growth and multiple small duplications.

Many parallel functional modules has been identified across all life kingdoms, including helicases and heat shock proteins in yeast [39]. Nevertheless, the absence of complete complex duplication at the WGD event (with a notable exception of CRPs) indicates that a stepwise duplication of modules [6], rather than whole-complex duplications, is a major mode of protein complex evolution in eukaryotes. Subunit-by-subunit module expansions amend the cellular machinery with the introduction of module variants. Subsequent duplications may give rise to bi-complex paralogs, which can be seen as intermediate phases on the evolutionary path leading to whole-complex duplication (see Figure 5). The process may be accompanied by attachment of additional subunits to one of the complexes or a differential loss of existing ones.

Figure 5
figure 5

A possible evolutionary scenario for incremental duplication of whole complexes. a) Ancestral three-subunit protein complex before duplications. b) Intra-complex duplication, gives rise to paralogous proteins, mutually exclusive as a part of the complex. Two module variants arise, with either A' or A" as one of subunits. c) Another intra-complex duplication makes possible assembly of up to four distinct module variants. d) Paralogs co-evolve to work together, specializing to accept only one of the interacting paralogs. This gives rise to more accentuated module variants, which now can be interpreted as separate complexes with a shared subunit. e) Eventually, all subunits may get duplicated, creating two independent protein complexes.

In eukaryotes a single protein events (loss, duplication or gain) dominate the evolution of functional modules. Even though here we do not quantify the prokaryote/eukaryote difference, scientific literature indicates that multiple copies of a protein complex can be found in bacteria. In the case of Complex I submodules, homologs of some of the recruited proteins already performed a function together previous to their involvement in the new pathway, and were duplicated in parallel of shortly after each other. This type of modular evolution in prokaryotes includes a duplication of, sometimes sizable, complexes: we know that a formate hydrogenlyase complex (FHL) of E. coli is in close evolutionary relation to Complex I [40]. Additionally, a duplication-prone FHL complex can be found in two copies in E. coli (FHL-1 and FHL-2), differing by only three subunits [41]. This observations lead to the hypothesis that appearance of copies of protein complexes in prokaryotes may be associated with the operon structure. The whole module encoded by an operon could duplicate by means of a single, segmental duplication. Alternatively, related complexes could evolve independently in different bacterial species and then be brought together by the horizontal gene transfer of the whole operon. Either way separate, independently functioning copy of a module could, for example, become recruited as a submodule of a bigger protein complex [42]. Interestingly, both E. coli FHL complexes are encoded by two operons, Hyf and Hyc [41].

RSC is an ATP-dependent chromatin remodeling complex of S. cerevisiae, essential for mitotic growth [43], that plays a role in expression regulation by activating and repressing the transcription [21]. This single complex exemplifies almost all aforementioned fates awaiting duplicated proteins (see Figure 6 and Table 4).

Figure 6
figure 6

The RSC protein complex (right) together with the related SWI/SNF complex (left) exhibit intricate homologous relationships between complexes' components. SWI/SNF and RSC complexes share 3 common proteins (ARP7, ARP9, RTT102). Additional three genes from each complex are paralogs (dashed lines), suggesting series of duplications (this includes a WGD pair RSC6-SNF12). For these bi-complex paralogs, it is very likely that the pre-duplication/ancestral forms of these genes were participating in both complexes, and the duplication event allowed paralogs to assume more specific roles in either complex. Zinc-cluster domain containing RSC3 and RSC30 are paralogs found in RSC complexes in equal proportions [43] and are known to interact physically, forming a stable heteromeric complex [21]. Despite operating as a heterodimer (possibly within the RSC complex as well), the two genes have different functions, with only RSC3 being an essential gene. The intra-complex paralog pair originated at the WGD event. At a given time, both proteins can be found in the complex, hence increase the subunit count by one above the initial complex, given that the ancestral interaction was probably a homomeric one. RSC1 and RSC2 (WGD paralogs), associate with the complex in a mutually exclusive manner [47] (confirmed by a low spoke value of 2.1 compared to average 5.0 for other subunits). Strains deficient in either RSC1 or RSC2 do not display significant growth deficiency, while synthetic lethality of paralogs suggests a certain degree of function redundancy [47]. A similar, but non-identical function of the two paralogs is revealed by differences in phenotypes of yeast strains with RSC1 (growth deficiency in hydroxyurea) or RSC2 deletions (temperature sensitive, [47]). Both complexes are required for certain functions such as proper expression of mid-late sporulation specific genes (e.g. spore wall formation), yet distinctly required, and specialized, for other operations [48]. The two module variants can be seen as a remarkable example of subfunctionalization.

Table 4 Different fates of WGD paralogs involved in RSC complex.

Overhangs do not, unlike their paralogs, participate in a protein complex. Direct interaction data confirm that overhangs do not seem to be associated with their paralog's protein complex. In our opinion features such as lower fraction of essential proteins or faster sequence evolution, make overhangs likely to be cases of neofunctionalization, initially working under relaxed evolutionary constraints. We hypothesize that overhangs released from the control of the ancestral protein complex, which are not purged from the genome (such as aforementioned SSD1/DIS3 pair), may form "seeds" for emerging complexes. This, accompanied with draft of additional subunits may to form novel complexes and ultimately become more embedded in the core cellular machinery.

Methods

WGD, SSD gene sets and the assessment of the age of duplications

WGD paralogs, genes that duplicated at the time of whole-genome duplication were taken from [4]. SSD paralogs (Small Scale Duplications) represent the most recent, non-WGD gene duplications. The SSD list consists of best bi-directional hits, i.e., gene pairs (A, B), such that their alignment score is higher than alignments of A against any other gene in the genome, and higher than alignment of B with any other gene in the genome (self-alignments excluding). 87% of WGD genes pass the criterion of best bi-directional similarity and were excluded from the SSD dataset.

WGD and SSD types of paralogs both stem from the most recent duplication of a given gene. To determine whether SSD duplications preceded or followed the WGD event, it is enough to assess the phylogenetic distribution of paralogs in multiple fungi species. Using orthology data from [44] we established that SSD paralogs were present in two copies before the ancestor of yeast underwent the WGD event. More specifically, the analysis of the fungal gene trees [8] shows that among SSD paralogs which participate in complexes, out of 84 pairs only a single gene pair duplicated after the WGD event (RSC30/YHR054C) and only for eight SSD paralog pairs gene trees imply duplication after the divergence with S. pombe (see additional file 1, Table S5 and Methods).

Ribosomal paralogs

From the initial analysis of our dataset of paralogs, it is apparent that genes involved in translation followed a distinct evolutionary route. Indeed, it is known from the literature that ribosomal protein sequences are highly constrained and many of the ribosomal protein pairs show exceptionally high levels of identity, likely subject to periodic gene conversion [20]. Interestingly, almost all cytoplasmic ribosomal proteins (CRPs), in a stark contrast to mitochondrial RP, were retained in duplicate after the WGD event. What could be the raison d'être this massive duplication? The set of CRP paralogs has many distinguishing properties, such as (a) a very similar amino acid sequence to paralogs, (b) a high mRNA expression correlation between paralogs (Spearman's correlation coefficient 0.8, see Methods for details), (c) the whole functional class, with few exceptions, was duplicated at WGD. These features imply a low level of functional differentiation and possibly an mRNA dosage increase as an explanation for the retention of both duplicates in CRPs, although new provocative evidence suggests more functional divergence than expected [45]. Nevertheless, such "Whole Ribosome Duplication" may signify the role the WGD event played in the evolution of anaerobic fermentation in yeast (compare with mitochondrial RP, additional file 1, Table S4).

Expression data were taken from the Gene Expression Omnibus (GEO) database [46] of the National Center for Biotechnology Information (NCBI), downloaded on 21 December 2006. Only multi-array datasets were considered, resulting in 357 microarray samples from 12 experiments, subsequently normalized (see additional file 1, Methods). We used Spearman's rank correlation coefficients to calculate the degree of co-expression between all gene pairs.

K. waltii orthologs were downloaded from Yeast Genome Order Browser (http://wolfe.gen.tcd.ie/browser[4]). Protein identity levels were calculated from alignments available in http://compbio.mit.edu.

Multiple module definitions were used to avoid bias of a certain protein complex annotation and make sure that results obtained hold independent of various protein complex maps used. MIPS data on protein complexes were downloaded from The MIPS Comprehensive Yeast Genome Database (CYGD, http://mips.gsf.de/[14]). The SGD GO complexes (in total 233 complexes, 1705 proteins) were generated by using the SGD GO component annotations (as of 9 May 2007) and then keeping only those components that have a GO description containing one of the following strings: complex, subunit, ribosome, proteasome, nucleosome, repairosome, degradosome, apoptosome, replisome, holoenzyme, snRNP. Only the lowest possible annotation level was maintained. Associations that where obtained from large scale experiments were removed.

References

  1. Ohno S: Evolution by Gene Duplication. 1970, London: Allen & Unwin

    Book  Google Scholar 

  2. Wolfe KH, Shields DC: Molecular evidence for an ancient duplication of the entire yeast genome. Nature. 1997, 387: 708-13. 10.1038/42711.

    Article  CAS  PubMed  Google Scholar 

  3. Kellis M, Patterson N, Endrizzi M, Birren B, Lander ES: Sequencing and comparison of yeast species to identify genes and regulatory elements. Nature. 2003, 423: 241-54. 10.1038/nature01644.

    Article  CAS  PubMed  Google Scholar 

  4. Byrne KP, Wolfe KH: The Yeast Gene Order Browser: Combining curated homology and syntenic context reveals gene fate in polyploid species. Genome Res. 2005, 15: 1456-146110. 10.1101/gr.3672305.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Pereira-Leal JB, Levy ED, Kamp C, Teichmann SA: Evolution of protein complexes by duplication of homomeric interactions. Genome Biol. 2007, 8: R51-10.1186/gb-2007-8-4-r51.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Pereira-Leal JB, Teichmann SA: Novel specificities emerge by stepwise duplication of functional modules. Genome Res. 2005, 15: 552-9. 10.1101/gr.3102105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Li L, Huang Y, Xia X, Sun Z: Preferential duplication in the sparse part of yeast protein interaction network. Mol Biol Evol. 2006, 23: 2467-73msl121. 10.1093/molbev/msl121.

    Article  CAS  PubMed  Google Scholar 

  8. Wapinski I, Pfeffer A, Friedman N, Regev A: Natural history and evolutionary principles of gene duplication in fungi. Nature. 2007, 449: 54-61. 10.1038/nature06107.

    Article  CAS  PubMed  Google Scholar 

  9. Musso G, Zhang Z, Emili A: Retention of protein complex membership by ancient duplicated gene products in budding yeast. Trends Genet. 2007, 23:

    Google Scholar 

  10. Guan Y, Dunham MJ, Troyanskaya OG: Functional Analysis of Gene Duplications in Saccharomyces cerevisiae. Genetics. 2007, 175: 933-94310. 10.1534/genetics.106.064329.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Hakes L, Pinney J, Lovell S, Oliver S, Robertson D: All duplicates are not equal: the difference between small-scale and genome duplication. Genome Biology. 2007, 8: R20910-10.1186/gb-2007-8-10-r209.

    Article  Google Scholar 

  12. Amoutzias GD, Veron AS, Weiner J, Robinson-Rechavi M, Bornberg-Bauer E, Oliver SG, Robertson DL: One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity. Mol Biol Evol. 2007, 24: 827-35msl211. 10.1093/molbev/msl211.

    Article  CAS  PubMed  Google Scholar 

  13. Mewes HW, Amid C, Arnold R, Frishman D, Güldener U, Mannhaupt G, Münsterkötter M, Pagel P, Strack N, Stümpflen V, Warfsmann J, Ruepp A: MIPS: analysis and annotation of proteins from whole genomes. Nucleic Acids Res. 2004, 32: D41-4. 10.1093/nar/gkh092.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Hong EL, Balakrishnan R, Dong Q, Christie KR, Park J, Binkley G, Costanzo MC, Dwight SS, Engel SR, Fisk DG, Hirschman JE, Hitz BC, Krieger CJ, Livstone MS, Miyasato SR, Nash RS, Oughtred R, Skrzypek MS, Weng S, Wong ED, Zhu KK, Dolinski K, Botstein D, Cherry JM: Gene Ontology annotations at SGD: new data sources and annotation methods. Nucleic Acids Res. 2008, 36: D577-81gkm909. 10.1093/nar/gkm909.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Hart GT, Lee I, Marcotte E: A high-accuracy consensus map of yeast protein complexes reveals modular nature of gene essentiality. BMC Bioinformatics. 2007, 8: 236-10.1186/1471-2105-8-236.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Collins SR, Kemmeren P, Zhao X, Greenblatt JF, Spencer F, Holstege FCP, Weissman JS, Krogan NJ: Toward a comprehensive atlas of the physical interactome of Saccharomyces cerevisiae. Mol Cell Proteomics. 2007, 6: 439-50.

    Article  CAS  PubMed  Google Scholar 

  17. Gavin A, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, Rau C, Jensen LJ, Bastuck S, Dümpelfeld B, Edelmann A, Heurtier M, Hoffman V, Hoefert C, Klein K, Hudak M, Michon A, Schelder M, Schirle M, Remor M, Rudi T, Hooper S, Bauer A, Bouwmeester T, Casari G, Drewes G, Neubauer G, Rick JM, Kuster B, Bork P, et al: Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006, 440: 631-6. 10.1038/nature04532.

    Article  CAS  PubMed  Google Scholar 

  18. Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, Li J, Pu S, Datta N, Tikuisis AP, Punna T, Peregrín-Alvarez JM, Shales M, Zhang X, Davey M, Robinson MD, Paccanaro A, Bray JE, Sheung A, Beattie B, Richards DP, Canadien V, Lalev A, Mena F, Wong P, Starostine A, Canete MM, Vlasblom J, Wu S, Orsi C, et al: Global landscape of protein complexes in the yeast Saccharomyces cerevisiae. Nature. 2006, 440: 637-64310. 10.1038/nature04670.

    Article  CAS  PubMed  Google Scholar 

  19. Planta RJ, Mager WH: The list of cytoplasmic ribosomal proteins of Saccharomyces cerevisiae. Yeast. 1998, 14: 471-7. 10.1002/(SICI)1097-0061(19980330)14:5<471::AID-YEA241>3.0.CO;2-U.

    Article  CAS  PubMed  Google Scholar 

  20. Kellis M, Birren BW, Lander ES: Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae. Nature. 2004, 428: 617-24. 10.1038/nature02424.

    Article  CAS  PubMed  Google Scholar 

  21. Angus-Hill ML, Schlichter A, Roberts D, Erdjument-Bromage H, Tempst P, Cairns BR: A Rsc3/Rsc30 Zinc Cluster Dimer Reveals Novel Roles for the Chromatin Remodeler RSC in Gene Expression and Cell Cycle Control. Molecular Cell. 2001, 7: 741-751. 10.1016/S1097-2765(01)00219-2.

    Article  CAS  PubMed  Google Scholar 

  22. Iwabe N, Kuma K, Hasegawa M, Osawa S, Miyata T: Evolutionary relationship of archaebacteria, eubacteria, and eukaryotes inferred from phylogenetic trees of duplicated genes. Proc Natl Acad Sci USA. 1989, 86: 9355-9. 10.1073/pnas.86.23.9355.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. de Lichtenberg U, Jensen LJ, Brunak S, Bork P: Dynamic Complex Formation During the Yeast Cell Cycle. Science. 2005, 307: 724-72710. 10.1126/science.1105103.

    Article  CAS  PubMed  Google Scholar 

  24. Mbonyi K, van Aelst L, Arguelles JC, Jans AW, Thevelein JM: Glucose-induced hyperaccumulation of cyclic AMP and defective glucose repression in yeast strains with reduced activity of cyclic AMP-dependent protein kinase. Mol Cell Biol. 1990, 10: 4518-4523.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Hodge MR, Kim G, Singh K, Cumsky MG: Inverse regulation of the yeast COX5 genes by oxygen and heme. Mol Cell Biol. 1989, 9: 1958-1964.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Burke PV, Raitt DC, Allen LA, Kellogg EA, Poyton RO: Effects of Oxygen Concentration on the Expression of Cytochrome c and Cytochrome c Oxidase Genes in Yeast. J Biol Chem. 1997, 272: 14705-1471210. 10.1074/jbc.272.23.14705.

    Article  CAS  PubMed  Google Scholar 

  27. Steinmetz LM, Scharfe C, Deutschbauer AM, Mokranjac D, Herman ZS, Jones T, Chu AM, Giaever G, Prokisch H, Oefner PJ, Davis RW: Systematic screen for human disease genes in yeast. Nat Genet. 2002, 31: 400-4.

    CAS  PubMed  Google Scholar 

  28. Hillenmeyer ME, Fung E, Wildenhain J, Pierce SE, Hoon S, Lee W, Proctor M, St Onge RP, Tyers M, Koller D, Altman RB, Davis RW, Nislow C, Giaever G: The chemical genomic portrait of yeast: uncovering a phenotype for all genes. Science. 2008, 320: 362-5. 10.1126/science.1150021.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Romisch K: Surfing the Sec61 channel: bidirectional protein translocation across the ER membrane. J Cell Sci. 1999, 112: 4185-4191.

    CAS  PubMed  Google Scholar 

  30. Sommer T, Wolf DH: Endoplasmic reticulum degradation: reverse protein flow of no return. FASEB J. 1997, 11: 1227-33.

    CAS  PubMed  Google Scholar 

  31. Robb A, Brown JD: Protein transport: two translocons are better than one. Mol Cell. 2001, 8: 484-6. 10.1016/S1097-2765(01)00339-2.

    Article  CAS  PubMed  Google Scholar 

  32. Kaeberlein M, Guarente L: Saccharomyces cerevisiae MPT5 and SSD1 function in parallel pathways to promote cell wall integrity. Genetics. 2002, 160: 83-95.

    PubMed Central  CAS  PubMed  Google Scholar 

  33. Mitchell P, Petfalski E, Shevchenko A, Mann M, Tollervey D: The exosome: a conserved eukaryotic RNA processing complex containing multiple 3'-->5' exoribonucleases. Cell. 1997, 91: 457-66. 10.1016/S0092-8674(00)80432-8.

    Article  CAS  PubMed  Google Scholar 

  34. Noguchi E, Hayashi N, Azuma Y, Seki T, Nakamura M, Nakashima N, Yanagida M, He X, Mueller U, Sazer S, Nishimoto T: Dis3, implicated in mitotic control, binds directly to Ran and enhances the GEF activity of RCC1. EMBO J. 1996, 15: 5595-605.

    PubMed Central  CAS  PubMed  Google Scholar 

  35. van Noort V, Snel B, Huynen MA: Predicting gene function by conserved co-expression. Trends in Genetics. 2003, 19: 238-242. 10.1016/S0168-9525(03)00056-8.

    Article  CAS  PubMed  Google Scholar 

  36. Conant GC, Wolfe KH: Functional partitioning of yeast co-expression networks after genome duplication. PLoS Biol. 2006, 4: e109-10.1371/journal.pbio.0040109.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Warner JR: The economics of ribosome biosynthesis in yeast. Trends Biochem Sci. 1999, 24: 437-40. 10.1016/S0968-0004(99)01460-7.

    Article  CAS  PubMed  Google Scholar 

  38. Planta RJ: Regulation of ribosome synthesis in yeast. Yeast. 1997, 13: 1505-18. 10.1002/(SICI)1097-0061(199712)13:16<1505::AID-YEA229>3.0.CO;2-I.

    Article  CAS  PubMed  Google Scholar 

  39. Li H, Pellegrini M, Eisenberg D: Detection of parallel functional modules by comparative analysis of genome sequences. Nat Biotechnol. 2005, 23: 253-60nbt1065. 10.1038/nbt1065.

    Article  CAS  PubMed  Google Scholar 

  40. Weiss H, Friedrich T, Hofhaus G, Preis D: The respiratory-chain NADH dehydrogenase (complex I) of mitochondria. Eur J Biochem. 1991, 197: 563-76. 10.1111/j.1432-1033.1991.tb15945.x.

    Article  CAS  PubMed  Google Scholar 

  41. Finel M: Organization and evolution of structural elements within complex I. Biochim Biophys Acta. 1998, 1364: 112-219593850. 10.1016/S0005-2728(98)00022-X.

    Article  CAS  PubMed  Google Scholar 

  42. Huynen MA, Gabaldón T, Snel B: Variation and evolution of biomolecular systems: searching for functional relevance. FEBS Lett. 2005, 579 (8): 1839-1845. 10.1016/j.febslet.2005.02.004.

    Article  CAS  PubMed  Google Scholar 

  43. Cairns BR, Lorch Y, Li Y, Zhang M, Lacomis L, Erdjument-Bromage H, Tempst P, Du J, Laurent B, Kornberg RD: RSC, an essential, abundant chromatin-remodeling complex. Cell. 1996, 87: 1249-60. 10.1016/S0092-8674(00)81820-6.

    Article  CAS  PubMed  Google Scholar 

  44. Scannell DR, Byrne KP, Gordon JL, Wong S, Wolfe KH: Multiple rounds of speciation associated with reciprocal gene loss in polyploid yeasts. Nature. 2006, 440: 341-5. 10.1038/nature04562.

    Article  CAS  PubMed  Google Scholar 

  45. Komili S, Farny NG, Roth FP, Silver PA: Functional Specificity among Ribosomal Proteins Regulates Gene Expression. Cell. 2007, 131: 557-571. 10.1016/j.cell.2007.08.037.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-10. 10.1093/nar/30.1.207.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Collins SR, Miller KM, Maas NL, Roguev A, Fillingham J, Chu CS, Schuldiner M, Gebbia M, Recht J, Shales M, Ding H, Xu H, Han J, Ingvarsdottir K, Cheng B, Andrews B, Boone C, Berger SL, Hieter P, Zhang Z, Brown GW, Ingles CJ, Emili A, Allis CD, Toczyski DP, Weissman JS, Greenblatt JF, Krogan NJ: Functional dissection of protein complexes involved in yeast chromosome biology using a genetic interaction map. Nature. 2007, 446: 806-10. 10.1038/nature05649.

    Article  CAS  PubMed  Google Scholar 

  48. He B, Chen P, Chen SY, Vancura KL, Michaelis S, Powers S: RAM2, an essential gene of yeast, and RAM1 encode the two polypeptide components of the farnesyltransferase that prenylates a-factor and Ras proteins. Proc Natl Acad Sci USA. 1991, 88: 11373-7. 10.1073/pnas.88.24.11373.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Witter DJ, Poulter CD: Yeast geranylgeranyltransferase type-II: steady state kinetic studies of the recombinant enzyme. Biochemistry. 1996, 35: 10454-63. 10.1021/bi960500y.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Ken Wolfe and Gavin Conant for insightful comments. Authors are grateful to Like Fokkens and Jos Boekhorst for discussions, Martin Oti for co-expression dataset, Joanna Parmley for carefully reading the manuscript and Patrick Kemmeren for sharing protein complex data. We would also like to thank anonymous reviewers for their valuable comments. This work was supported by the Netherlands Genomics Initiative (Horizon programme).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Berend Snel.

Additional information

Authors' contributions

RS and BS designed the study. RS performed the analysis. MH contributed analytical methods. RS and BS wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Szklarczyk, R., Huynen, M.A. & Snel, B. Complex fate of paralogs. BMC Evol Biol 8, 337 (2008). https://doi.org/10.1186/1471-2148-8-337

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-8-337

Keywords