Fungal Zuotin proteins evolved from MIDA1-like factors by lineage-specific loss of MYB domains

Structure and evolution of the 4-helix bundle domain of Zuotin, a J-domain protein co-chaperone of Hsp70

The J-domain protein Zuotin is a multi-domain eukaryotic Hsp70 co-chaperone. Though it is primarily ribosome-associated, positioned at the exit of the 60S subunit tunnel where it promotes folding of nascent polypeptide chains, Zuotin also has off-ribosome functions. Domains of Zuotin needed for 60S association and interaction with Hsp70 are conserved in eukaryotes. However, whether the 4-helix bundle (4HB) domain is conserved remains an open question. We undertook evolutionary and structural approaches to clarify this issue. We found that the 4HB segment of human Zuotin also forms a bundle of 4 helices. The positive charge of Helix I, which in Saccharomyces cerevisiae is responsible for interaction with the 40S subunit, is particularly conserved. However, the C-termini of fungal and human 4HBs are not similar. In fungi the C-terminal segment forms a plug that folds back into the bundle; in S. cerevisiae it plays an important role in bundle stability and, off the ribosome, in transcriptional activation. In human, C-terminal helix IV of the 4HB is extended, protruding from the bundle. This extension serves as a linker to the regulatory SANT domains, which are present in animals, plants and protists, but not fungi. Further analysis of Zuotin sequences revealed that the plug likely arose as a result of genomic rearrangement upon SANT domain loss early in the fungal lineage. In the lineage leading to S. cerevisiae, the 4HB was subjected to positive selection with the plug becoming increasingly hydrophobic. Eventually, these hydrophobic plug residues were coopted for a novel regulatory function—activation of a recently emerged transcription factor, Pdr1. Our data suggests that Zuotin evolved off-ribosome functions twice—once involving SANT domains, then later in fungi, after SANT domain loss, by coopting the hydrophobic plug. Zuotin serves as an example of complex intertwining of molecular chaperone function and cell regulation.

Protein quality control at the ribosome: focus on RAC, NAC and RQC

The biogenesis of new polypeptides by ribosomes and their subsequent correct folding and localization to the appropriate cellular compartments are essential key processes to maintain protein homoeostasis. These complex mechanisms are governed by a repertoire of protein biogenesis factors that directly bind to the ribosome and chaperone nascent polypeptide chains as soon as they emerge from the ribosomal tunnel exit. This nascent chain 'welcoming committee' regulates multiple co-translational processes including protein modifications, folding, targeting and degradation. Acting at the front of the protein production line, these ribosome-associated protein biogenesis factors lead the way in the cellular proteostasis network to ensure proteome integrity. In this article, I focus on three different systems in eukaryotes that are critical for the maintenance of protein homoeostasis by controlling the birth, life and death of nascent polypeptide chains.

The intriguing realm of protein biogenesis: Facing the green co-translational protein maturation networks

The ribosome is the cell's protein-making factory, a huge protein-RNA complex, that is essential to life. Determining the high-resolution structures of the stable "core" of this factory was among the major breakthroughs of the past decades, and was awarded the Nobel Prize in 2009. Now that the mysteries of the ribosome appear to be more traceable, detailed understanding of the mechanisms that regulate protein synthesis includes not only the well-known steps of initiation, elongation, and termination but also the less comprehended features of the co-translational events associated with the maturation of the nascent chains. The ribosome is a platform for co-translational events affecting the nascent polypeptide, including protein modifications, folding, targeting to various cellular compartments for integration into membrane or translocation, and proteolysis. These events are orchestrated by ribosome-associated protein biogenesis factors (RPBs), a group of a dozen or more factors that act as the "welcoming committee" for the nascent chain as it emerges from the ribosome. In plants these factors have evolved to fit the specificity of different cellular compartments: cytoplasm, mitochondria and chloroplast. This review focuses on the current state of knowledge of these factors and their interaction around the exit tunnel of dedicated ribosomes. Particular attention has been accorded to the plant system, highlighting the similarities and differences with other organisms.

Functional conservation and divergence of J-domain-containing ZUO1/ZRF orthologs throughout evolution

Heat shock protein 40s (Hsp40s), also known as J-proteins, are conserved in prokaryotes and eukaryotes. The Zuotin/Zuotin-related factor (ZUO1/ZRF) family belongs to a novel Hsp40 clade exclusively found in eukaryotes. Zuotin/Zuotin-related factor proteins are characterized by a large N terminal ZUO1 domain originally identified in the yeast ZUO1 protein. The ZUO1 domain is characterized by a highly conserved J-domain, together with an atypical UBD domain first identified in the human ZRF1 protein. Furthermore, ZUO1/ZRF protein families in animals and plants harbor a pair of C terminal SANT domains, suggesting the divergence of their functions with those in fungi. Zuotin/Zuotin-related factor proteins retain the ancestral function as an Hsp70co-chaperone implicated in protein folding and renaturation after stress; these proteins also perform diverse neofunctions in the cytoplasm and transcriptional and/or epigenetic regulatory functions in the nucleus. Therefore, these proteins are involved in translational fidelity control, ribosomal biogenesis, asymmetric cell division, cell cycle, apoptosis, differentiation, and tumorigenesis. The results of sequence and domain organization analysis of proteins from diverse organisms provided valuable insights into the evolutionary conservation and diversity of ZUO1/ZRF protein family. Further, phylogenetic analysis provides a platform for future functional investigation on the ZUO1/ZRF protein family, particularly in higher plants.

Experimental evidence validating the computational inference of functional associations from gene fusion events: a critical survey

More than a decade ago, a number of methods were proposed for the inference of protein interactions, using whole-genome information from gene clusters, gene fusions and phylogenetic profiles. This structural and evolutionary view of entire genomes has provided a valuable approach for the functional characterization of proteins, especially those without sequence similarity to proteins of known function. Furthermore, this view has raised the real possibility to detect functional associations of genes and their corresponding proteins for any entire genome sequence. Yet, despite these exciting developments, there have been relatively few cases of real use of these methods outside the computational biology field, as reflected from citation analysis. These methods have the potential to be used in high-throughput experimental settings in functional genomics and proteomics to validate results with very high accuracy and good coverage. In this critical survey, we provide a comprehensive overview of 30 most prominent examples of single pairwise protein interaction cases in small-scale studies, where protein interactions have either been detected by gene fusion or yielded additional, corroborating evidence from biochemical observations. Our conclusion is that with the derivation of a validated gold-standard corpus and better data integration with big experiments, gene fusion detection can truly become a valuable tool for large-scale experimental biology.

Role for the molecular chaperones Zuo1 and Ssz1 in quorum sensing via activation of the transcription factor Pdr1

Zuo1 functions as a J-protein cochaperone of its partner Hsp70. In addition, the C terminus of Zuo1 and the N terminus of Ssz1, with which Zuo1 forms a heterodimer, can independently activate the Saccharomyces cerevisiae transcription factor pleiotropic drug resistance 1 (Pdr1). Here we report that activation of Pdr1 by Zuo1 or Ssz1 causes premature growth arrest of cells during the diauxic shift, as they adapt to the changing environmental conditions. Conversely, cells lacking Zuo1 or Ssz1 overgrow, arresting at a higher cell density, an effect overcome by activation of Pdr1. Cells lacking the genes encoding plasma membrane transporters Pdr5 and Snq2, two targets of Pdr1, also overgrow at the diauxic shift. Adding conditioned medium harvested from cultures of wild-type cells attenuated the overgrowth of both zuo1Δssz1Δ and pdr5Δsnq2Δ cells, suggesting the extracellular presence of molecules that signal growth arrest. In addition, our yeast two-hybrid analysis revealed an interaction between Pdr1 and both Zuo1 and Ssz1. Together, our results support a model in which (i) membrane transporters, encoded by Pdr1 target genes act to promote cell-cell communication by exporting quorum sensing molecules, in addition to playing a role in pleiotropic drug resistance; and (ii) molecular chaperones function at promoters to regulate this intercellular communication through their activation of the transcription factor Pdr1.

Evolutionary genomics implies a specific function of Ant4 in mammalian and anole lizard male germ cells

Most vertebrates have three paralogous genes with identical intron-exon structures and a high degree of sequence identity that encode mitochondrial adenine nucleotide translocase (Ant) proteins, Ant1 (Slc25a4), Ant2 (Slc25a5) and Ant3 (Slc25a6). Recently, we and others identified a fourth mammalian Ant paralog, Ant4 (Slc25a31), with a distinct intron-exon structure and a lower degree of sequence identity. Ant4 was expressed selectively in testis and sperm in adult mammals and was indeed essential for mouse spermatogenesis, but it was absent in birds, fish and frogs. Since Ant2 is X-linked in mammalian genomes, we hypothesized that the autosomal Ant4 gene may compensate for the loss of Ant2 gene expression during male meiosis in mammals. Here we report that the Ant4 ortholog is conserved in green anole lizard (Anolis carolinensis) and demonstrate that it is expressed in the anole testis. Further, a degenerate DNA fragment of putative Ant4 gene was identified in syntenic regions of avian genomes, indicating that Ant4 was present in the common amniote ancestor. Phylogenetic analyses suggest an even more ancient origin of the Ant4 gene. Although anole lizards are presumed male (XY) heterogametic, like mammals, copy numbers of the Ant2 as well as its neighboring gene were similar between male and female anole genomes, indicating that the anole Ant2 gene is either autosomal or located in the pseudoautosomal region of the sex chromosomes, in contrast to the case to mammals. These results imply the conservation of Ant4 is not likely simply driven by the sex chromosomal localization of the Ant2 gene and its subsequent inactivation during male meiosis. Taken together with the fact that Ant4 protein has a uniquely conserved structure when compared to other somatic Ant1, 2 and 3, there may be a specific advantage for mammals and lizards to express Ant4 in their male germ cells.

Evolutionary and comparative analysis of MYB and bHLH plant transcription factors

The expansion of gene families encoding regulatory proteins is typically associated with the increase in complexity characteristic of multi-cellular organisms. The MYB and basic helix-loop-helix (bHLH) families provide excellent examples of how gene duplication and divergence within particular groups of transcription factors are associated with, if not driven by, the morphological and metabolic diversity that characterize the higher plants. These gene families expanded dramatically in higher plants; for example, there are approximately 339 and 162 MYB and bHLH genes, respectively, in Arabidopsis, and approximately 230 and 111, respectively, in rice. In contrast, the Chlamydomonas genome has only 38 MYB genes and eight bHLH genes. In this review, we compare the MYB and bHLH gene families from structural, evolutionary and functional perspectives. The knowledge acquired on the role of many of these factors in Arabidopsis provides an excellent reference to explore sequence-function relationships in crops and other plants. The physical interaction and regulatory synergy between particular sub-classes of MYB and bHLH factors is perhaps one of the best examples of combinatorial plant gene regulation. However, members of the MYB and bHLH families also interact with a number of other regulatory proteins, forming complexes that either activate or repress the expression of sets of target genes that are increasingly being identified through a diversity of high-throughput genomic approaches. The next few years are likely to witness an increasing understanding of the extent to which conserved transcription factors participate at similar positions in gene regulatory networks across plant species.

The chaperone network connected to human ribosome-associated complex

Mammalian ribosome-associated complex (mRAC), consisting of the J-domain protein MPP11 and the atypical Hsp70 homolog (70-homolog) Hsp70L1, can partly complement the function of RAC, which is the homologous complex from Saccharomyces cerevisiae. RAC is the J-domain partner exclusively of the 70-homolog Ssb, which directly and independently of RAC binds to the ribosome. We here show that growth defects due to mRAC depletion in HeLa cells resemble those of yeast strains lacking RAC. Functional conservation, however, did not extend to the 70-homolog partner of mRAC. None of the major human 70-homologs was able to complement the growth defects of yeast strains lacking Ssb or was bound to ribosomes in an Ssb-like manner. Instead, our data suggest that mRAC was a specific partner of human Hsp70 but not of its close homolog Hsc70. On a mechanistic level, ATP binding, but not ATP hydrolysis, by Hsp70L1 affected mRAC's function as a J-domain partner of Hsp70. The combined data indicate that, while functionally conserved, yeast and mammalian cells have evolved distinct solutions to ensure that Hsp70-type chaperones can efficiently assist the biogenesis of newly synthesized polypeptide chains.

Phylogenomic analyses reveal the evolutionary origin of the inhibin alpha-subunit, a unique TGFbeta superfamily antagonist

Transforming growth factor-beta (TGFβ) homologues form a diverse superfamily that arose early in animal evolution and control cellular function through membrane-spanning, conserved serine-threonine kinases (RII and RI receptors). Activin and inhibin are related dimers within the TGFβ superfamily that share a common β-subunit. The evolution of the inhibin α-subunit created the only antagonist within the TGFβ superfamily and the only member known to act as an endocrine hormone. This hormone introduced a new level of complexity and control to vertebrate reproductive function. The novel functions of the inhibin α-subunit appear to reflect specific insertion-deletion changes within the inhibin β-subunit that occurred during evolution. Using phylogenomic analysis, we correlated specific insertions with the acquisition of distinct functions that underlie the phenotypic complexity of vertebrate reproductive processes. This phylogenomic approach presents a new way of understanding the structure-function relationships between inhibin, activin, and the larger TGFβ superfamily.

Functional analysis of the Volvox carteri asymmetric division protein GlsA

The Zuotin-family J protein chaperone GlsA is essential for the asymmetric divisions that establish germ and somatic cell initials during embryogenesis in the green alga Volvox carteri, but it is not known on what cellular process GlsA acts to carry out this function. Most GlsA protein is nuclear, and GlsA possesses two SANT domains, suggesting that GlsA may function as a transcriptional regulator. On the other hand, close homologs from yeast and mice are ribosome-associated factors that regulate translation fidelity, implying GlsA might also regulate translation. Here we set out to gain additional evidence regarding the function of GlsA, specifically with respect to its possible involvement in transcription and translation. We found that like zuotin mutants, glsA mutants are ultrasensitive to both cold and to the ribosome-binding aminoglycoside antibiotic paromomycin, so some fraction of GlsA is likely to be ribosome associated. We also found that GlsA co-immunoprecipitates with histones and that this interaction is dependent on the presence of intact SANT domains. Through rescue experiments using transgenes that encode GlsA variants, we determined that the growth and asymmetric division defects of the glsA mutant are separable-a GlsA variant that rescued the growth defects did not completely rescue the asymmetric division phenotype. Considered in total, our results suggest that GlsA acts both at the level of translation and transcription, but the function that is essential for tolerance to paromomycin and cold is not sufficient for asymmetric cell division.

Arabidopsis thaliana J-class heat shock proteins: cellular stress sensors

Plants are sessile organisms that have evolved a variety of mechanisms to maintain their cellular homeostasis under stressful environmental conditions. Survival of plants under abiotic stress conditions requires specialized group of heat shock protein machinery, belonging to Hsp70:J-protein family. These heat shock proteins are most ubiquitous types of chaperone machineries involved in diverse cellular processes including protein folding, translocation across cell membranes, and protein degradation. They play a crucial role in maintaining the protein homeostasis by reestablishing functional native conformations under environmental stress conditions, thus providing protection to the cell. J-proteins are co-chaperones of Hsp70 machine, which play a critical role by stimulating Hsp70s ATPase activity, thereby stabilizing its interaction with client proteins. Using genome-wide analysis of Arabidopsis thaliana, here we have outlined identification and systematic classification of J-protein co-chaperones which are key regulators of Hsp70s function. In comparison with Saccharomyces cerevisiae model system, a comprehensive domain structural organization, cellular localization, and functional diversity of A. thaliana J-proteins have also been summarized.

GRASSIUS: a platform for comparative regulatory genomics across the grasses

Transcription factors (TFs) are major players in gene regulatory networks and interactions between TFs and their target genes furnish spatiotemporal patterns of gene expression. Establishing the architecture of regulatory networks requires gathering information on TFs, their targets in the genome, and the corresponding binding sites. We have developed GRASSIUS (Grass Regulatory Information Services) as a knowledge-based Web resource that integrates information on TFs and gene promoters across the grasses. In its initial implementation, GRASSIUS consists of two separate, yet linked, databases. GrassTFDB holds information on TFs from maize (Zea mays), sorghum (Sorghum bicolor), sugarcane (Saccharum spp.), and rice (Oryza sativa). TFs are classified into families and phylogenetic relationships begin to uncover orthologous relationships among the participating species. This database also provides a centralized clearinghouse for TF synonyms in the grasses. GrassTFDB is linked to the grass TFome collection, which provides clones in recombination-based vectors corresponding to full-length open reading frames for a growing number of grass TFs. GrassPROMDB contains promoter and cis-regulatory element information for those grass species and genes for which enough data are available. The integration of GrassTFDB and GrassPROMDB will be accomplished through GrassRegNet as a first step in representing the architecture of grass regulatory networks. GRASSIUS can be accessed from www.grassius.org.

Functional characterization of the atypical Hsp70 subunit of yeast ribosome-associated complex

Eukaryotic ribosomes carry a stable chaperone complex termed ribosome-associated complex consisting of the J-domain protein Zuo1 and the Hsp70 Ssz1. Zuo1 and Ssz1 together with the Hsp70 homolog Ssb1/2 form a functional triad involved in translation and early polypeptide folding processes. Strains lacking one of these components display slow growth, cold sensitivity, and defects in translational fidelity. Ssz1 diverges from canonical Hsp70s insofar that neither the ability to hydrolyze ATP nor binding to peptide substrates is essential in vivo. The exact role within the chaperone triad and whether or not Ssz1 can hydrolyze ATP has remained unclear. We now find that Ssz1 is not an ATPase in vitro, and even its ability to bind ATP is dispensable in vivo. Furthermore, Ssz1 function was independent of ribosome-associated complex formation, indicating that Ssz1 is not merely a structural scaffold for Zuo1. Finally, Ssz1 function in vivo was inactivated when both nucleotide binding and Zuo1 interaction via the C-terminal domain were disrupted in the same mutant. The two domains of this protein thus cooperate in a way that allows for severe interference in either but not in both of them.

The evolutionary relationship of the domain architectures in the RhoGEF-containing proteins

Domain insertions and deletions lead to variations in the domain architectures of the proteins from their common ancestor. In this work, we investigated four groups of the RhoGEF-containing proteins from different organisms with domain architectures RhoGEF-PH-SH3, SH3-RhoGEF-PH, RhoGEF-PH, and SH3-RhoGEF defined in the Pfam database. The phylogenetic trees were constructed using each individual domain and/or the combinations of all the domains. The phylogenetic analysis suggests that RhoGEF-PH-SH3 and SH3-RhoGEF-PH might have evolved from RhoGEF-PH through the insertion of SH3 independently, while SH3-RhoGEF of proteins in fruit fly might have evolved from SH3-RhoGEF-PH by the degeneration of PH domain.

The chaperones MPP11 and Hsp70L1 form the mammalian ribosome-associated complex

Soluble Hsp70 homologs cotranslationally interact with nascent polypeptides in all kingdoms of life. In addition, fungi possess a specialized Hsp70 system attached to ribosomes, which in Saccharomyces cerevisiae consists of the Hsp70 homologs Ssb1/2p, Ssz1p, and the Hsp40 homolog zuotin. Ssz1p and zuotin are assembled into a unique heterodimeric complex termed ribosome-associated complex. So far, no such specialized chaperones have been identified on ribosomes of higher eukaryotes. However, a family of proteins characterized by an N-terminal zuotin-homology domain fused to a C-terminal two-repeat Myb domain is present in animals and plants. Members of this family, like human MPP11 and mouse MIDA1, have been implicated in the regulation of cell growth. Specific targets of MPP11/MIDA1, however, have remained elusive. Here, we report that MPP11 is localized to the cytosol and associates with ribosomes. Purification of MPP11 revealed that it forms a stable complex with Hsp70L1, a distantly related homolog of Ssz1p. Complementation experiments indicate that mammalian ribosome-associated complex is functional in yeast. We conclude that despite a low degree of homology on the amino acid level cooperation of ribosome-associated chaperones with the translational apparatus is well conserved in eukaryotic cells.

The SANT2 domain of the murine tumor cell DnaJ-like protein 1 human homologue interacts with alpha1-antichymotrypsin and kinetically interferes with its serpin inhibitory activity

The murine tumor cell DnaJ-like protein 1 or MTJ1/ERdj1 is a membrane J-domain protein enriched in microsomal and nuclear fractions. We previously showed that its lumenal J-domain stimulates the ATPase activity of the molecular chaperone BiP/GRP78 (Chevalier, M., Rhee, H., Elguindi, E. C., and Blond, S. Y. (2000) J. Biol. Chem. 275, 19620-19627). MTJ1/ERdj1 also contains a large carboxyl-terminal cytosolic extension composed of two tryptophan-mediated repeats or SANT domains for which the function(s) is unknown. Here we describe the cloning of the human homologue HTJ1 and its interaction with alpha(1)-antichymotrypsin (ACT), a member of the serine proteinase inhibitor (serpin) family. The interaction was initially identified in a two-hybrid screening and further confirmed in vitro by dot blots, native electrophoresis, and fluorescence studies. The second SANT domain of HTJ1 (SANT2) was found to be sufficient for binding to ACT, both in yeast and in vitro. Single tryptophan-alanine substitutions at two strictly conserved residues significantly (Trp-497) or totally (Trp-520) abolished the interaction with ACT. SANT2 binds to human ACT with an intrinsic affinity equal to 0.5 nm. Preincubation of ACT with nearly stoichiometric concentrations of SANT2 wild-type but not SANT2: W520A results in an apparent loss of ACT inhibitory activity toward chymotrypsin. Kinetic analysis indicates that the formation of the covalent inhibitory complex ACT-chymotrypsin is significantly delayed in the presence of SANT2 with no change on the catalytic efficiency of the enzyme. This work demonstrates for the first time that the SANT2 domain of MTJ1/HTJ1/ERdj1 mediates stable and high affinity protein-protein interactions.

Molecular cloning of the estrogen and progesterone receptors of the American alligator

Steroid hormones perform many essential roles in vertebrates during embryonic development, reproduction, growth, water balance, and responses to stress. The estrogens are essential for normal reproductive activity in female and male vertebrates and appear to have direct actions during sex determination in some vertebrates. To begin to understand the molecular mechanisms of estrogen action in alligators, we have isolated cDNAs encoding the estrogen receptors (ER) from the ovary. Degenerate PCR primers specific to ER were designed and used to amplify alligator ovary RNA. Two different DNA fragments (ERalpha and ERbeta) were obtained and the full-length alligator ERalpha cDNA was obtained using 5' and 3' RACE. The inferred amino acid sequence of alligator ERalpha (aERalpha) was very similar to the chicken ERalpha (91% identity), although phylogenetic analyses suggested profound differences in the rate of sequence evolution for vertebrate ER sequences. We also isolated partial DNA fragments encoding ERbeta and the progesterone receptor (PR) of the alligator, both of which show strong sequence similarities to avian ERbeta and PR. We examined the expression levels of these three steroid receptors (ERalpha, ERbeta, and PR) in the ovary of juvenile alligators and observed detectable levels of all three receptors. Quantitative RT-PCR showed that gonadal ERalpha transcript levels in juvenile alligators decreased after E2 treatment whereas ERbeta and PR transcripts were not changed. These results provide tools that will allow future studies examining the regulation and ontogenic expression of steroid receptors in alligators and expand our knowledge of vertebrate steroid receptor evolution.

The role of GlsA in the evolution of asymmetric cell division in the green alga Volvox carteri

Volvox carteri, a green alga in the order Volvocales, contains two completely differentiated cell types, small motile somatic cells and large reproductive cells called gonidia, that are set apart from each other during embryogenesis by a series of visibly asymmetric cell divisions. Mutational analysis has revealed a class of genes (gonidialess, gls) that are required specifically for asymmetric divisions in V. carteri, but that are dispensable for symmetric divisions. Previously we cloned one of these genes, glsA, and showed that it encodes a chaperone-like protein (GlsA) that has close orthologs in a diverse set of eukaryotes, ranging from fungi to vertebrates and higher plants. In the present study we set out to explore the role of glsA in the evolution of asymmetric division in the volvocine algae by cloning and characterizing a glsA ortholog from one of the simplest members of the group, Chlamydomonas reinhardtii, which does not undergo asymmetric divisions. This ortholog (which we have named gar1, for glsA related) is predicted to encode a protein that is 70% identical to GlsA overall, and that is most closely related to GlsA in the same domains that are most highly conserved between GlsA and its other known orthologs. We report that a gar1 transgene fully complements the glsA mutation in V. carteri, a result that suggests that asymmetric division probably arose through the modification of a gene whose product interacts with GlsA, but not through a modification of glsA itself.

Genome and protein evolution in eukaryotes

The past year has seen the completion of the genome sequence of the flowering plant Arabidopsis thaliana and the initial sequence reports of the human genome. The availability of completely sequenced eukaryotic genomes from disparate phylogenetic lineages has opened the door to comparative analyses and a better understanding of the evolutionary processes shaping genomes. Complex many-to-many relationships between genes from different species appear to be the norm, suggesting that transfer of detailed functional annotation will not be straightforward. In addition to expansion and contraction of gene families, new genes evolve from recombination of pre-existing domains, although some domain families do appear to have evolved recently and to be specific to restricted phylogenetic lineages. The overall picture is of a huge diversity of gene content within eukaryotic genomes, reflecting different functional demands in different species.