Genetic redundancy in vertebrates: polyploidy and persistence of genes encoding multidomain proteins

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

Genotype-Phenotype Relationships in the Context of Transcriptional Adaptation and Genetic Robustness

Genetic manipulations with a robust and predictable outcome are critical to investigate gene function, as well as for therapeutic genome engineering. For many years, knockdown approaches and reagents including RNA interference and antisense oligonucleotides dominated functional studies; however, with the advent of precise genome editing technologies, CRISPR-based knockout systems have become the state-of-the-art tools for such studies. These technologies have helped decipher the role of thousands of genes in development and disease. Their use has also revealed how limited our understanding of genotype-phenotype relationships is. The recent discovery that certain mutations can trigger the transcriptional modulation of other genes, a phenomenon called transcriptional adaptation, has provided an additional explanation for the contradicting phenotypes observed in knockdown versus knockout models and increased awareness about the use of each of these approaches. In this review, we first cover the strengths and limitations of different gene perturbation strategies. Then we highlight the diverse ways in which the genotype-phenotype relationship can be discordant between these different strategies. Finally, we review the genetic robustness mechanisms that can lead to such discrepancies, paying special attention to the recently discovered phenomenon of transcriptional adaptation. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

Paralog dependency indirectly affects the robustness of human cells

The protective redundancy of paralogous genes partly relies on the fact that they carry their functions independently. However, a significant fraction of paralogous proteins may form functionally dependent pairs, for instance, through heteromerization. As a consequence, one could expect these heteromeric paralogs to be less protective against deleterious mutations. To test this hypothesis, we examined the robustness landscape of gene loss-of-function by CRISPR-Cas9 in more than 450 human cell lines. This landscape shows regions of greater deleteriousness to gene inactivation as a function of key paralog properties. Heteromeric paralogs are more likely to occupy such regions owing to their high expression and large number of protein-protein interaction partners. Further investigation revealed that heteromers may also be under stricter dosage balance, which may also contribute to the higher deleteriousness upon gene inactivation. Finally, we suggest that physical dependency may contribute to the deleteriousness upon loss-of-function as revealed by the correlation between the strength of interactions between paralogs and their higher deleteriousness upon loss of function.

Selective Constraints on Coding Sequences of Nervous System Genes Are a Major Determinant of Duplicate Gene Retention in Vertebrates

The evolutionary history of vertebrates is marked by three ancient whole-genome duplications: two successive rounds in the ancestor of vertebrates, and a third one specific to teleost fishes. Biased loss of most duplicates enriched the genome for specific genes, such as slow evolving genes, but this selective retention process is not well understood. To understand what drives the long-term preservation of duplicate genes, we characterized duplicated genes in terms of their expression patterns. We used a new method of expression enrichment analysis, TopAnat, applied to in situ hybridization data from thousands of genes from zebrafish and mouse. We showed that the presence of expression in the nervous system is a good predictor of a higher rate of retention of duplicate genes after whole-genome duplication. Further analyses suggest that purifying selection against the toxic effects of misfolded or misinteracting proteins, which is particularly strong in nonrenewing neural tissues, likely constrains the evolution of coding sequences of nervous system genes, leading indirectly to the preservation of duplicate genes after whole-genome duplication. Whole-genome duplications thus greatly contributed to the expansion of the toolkit of genes available for the evolution of profound novelties of the nervous system at the base of the vertebrate radiation.

The Phenotypic Plasticity of Duplicated Genes in Saccharomyces cerevisiae and the Origin of Adaptations

Gene and genome duplication are the major sources of biological innovations in plants and animals. Functional and transcriptional divergence between the copies after gene duplication has been considered the main driver of innovations . However, here we show that increased phenotypic plasticity after duplication plays a more major role than thought before in the origin of adaptations. We perform an exhaustive analysis of the transcriptional alterations of duplicated genes in the unicellular eukaryote Saccharomyces cerevisiae when challenged with five different environmental stresses. Analysis of the transcriptomes of yeast shows that gene duplication increases the transcriptional response to environmental changes, with duplicated genes exhibiting signatures of adaptive transcriptional patterns in response to stress. The mechanism of duplication matters, with whole-genome duplicates being more transcriptionally altered than small-scale duplicates. The predominant transcriptional pattern follows the classic theory of evolution by gene duplication; with one gene copy remaining unaltered under stress, while its sister copy presents large transcriptional plasticity and a prominent role in adaptation. Moreover, we find additional transcriptional profiles that are suggestive of neo- and subfunctionalization of duplicate gene copies. These patterns are strongly correlated with the functional dependencies and sequence divergence profiles of gene copies. We show that, unlike singletons, duplicates respond more specifically to stress, supporting the role of natural selection in the transcriptional plasticity of duplicates. Our results reveal the underlying transcriptional complexity of duplicated genes and its role in the origin of adaptations.

Meiosis evolves: adaptation to external and internal environments

306 I. 306 II. 307 III. 312 IV. 317 V. 318 319 References 319 SUMMARY: Meiosis is essential for the fertility of most eukaryotes and its structures and progression are conserved across kingdoms. Yet many of its core proteins show evidence of rapid or adaptive evolution. What drives the evolution of meiosis proteins? How can constrained meiotic processes be modified in response to challenges without compromising their essential functions? In surveying the literature, we found evidence of two especially potent challenges to meiotic chromosome segregation that probably necessitate adaptive evolutionary responses: whole-genome duplication and abiotic environment, especially temperature. Evolutionary solutions to both kinds of challenge are likely to involve modification of homologous recombination and synapsis, probably via adjustments of core structural components important in meiosis I. Synthesizing these findings with broader patterns of meiosis gene evolution suggests that the structural components of meiosis coevolve as adaptive modules that may change in primary sequence and function while maintaining three-dimensional structures and protein interactions. The often sharp divergence of these genes among species probably reflects periodic modification of entire multiprotein complexes driven by genomic or environmental changes. We suggest that the pressures that cause meiosis to evolve to maintain fertility may cause pleiotropic alterations of global crossover rates. We highlight several important areas for future research.

Tempo and mode of gene duplication in mammalian ribosomal protein evolution

Gene duplication has been widely recognized as a major driver of evolutionary change and organismal complexity through the generation of multi-gene families. Therefore, understanding the forces that govern the evolution of gene families through the retention or loss of duplicated genes is fundamentally important in our efforts to study genome evolution. Previous work from our lab has shown that ribosomal protein (RP) genes constitute one of the largest classes of conserved duplicated genes in mammals. This result was surprising due to the fact that ribosomal protein genes evolve slowly and transcript levels are very tightly regulated. In our present study, we identified and characterized all RP duplicates in eight mammalian genomes in order to investigate the tempo and mode of ribosomal protein family evolution. We show that a sizable number of duplicates are transcriptionally active and are very highly conserved. Furthermore, we conclude that existing gene duplication models do not readily account for the preservation of a very large number of intact retroduplicated ribosomal protein (RT-RP) genes observed in mammalian genomes. We suggest that selection against dominant-negative mutations may underlie the unexpected retention and conservation of duplicated RP genes, and may shape the fate of newly duplicated genes, regardless of duplication mechanism.

Gene dosage effects: nonlinearities, genetic interactions, and dosage compensation

High-throughput genomic analyses have shown that many mutations, including loss-of-function (LOF) mutations, are present in diseased as well as in healthy individuals. Gene dosage effects due to deletions, duplications, and LOF mutations provide avenues to explore oligo- and multigenic inheritance. Here, we focus on several mechanisms that mediate gene dosage effects and analyze biochemical interactions among multiple gene products that are sources of nonlinear relations connecting genotypes and phenotypes. We also explore potential mechanisms that compensate for gene dosage effects. Understanding these issues is critical to understanding why an individual bearing a few damaging mutations can be severely diseased, whereas others harboring tens of potentially deleterious mutations can appear quite healthy.

Genetic redundancy strengthens the circadian clock leading to a narrow entrainment range

Circadian clocks are internal timekeepers present in almost all organisms. Driven by a genetic network of highly conserved structure, they generate self-sustained oscillations that entrain to periodic external signals such as the 24 h light-dark cycle. Vertebrates possess multiple, functionally overlapping homologues of the core clock genes. Furthermore, vertebrate clocks entrain to a range of periods three times as narrow as that of other organisms. We asked whether genetic redundancies play a role in governing entrainment properties and analysed locomotor activity rhythms of genetically modified mice lacking one set of clock homologues. Exposing them to non-24 h light-dark cycles, we found that the mutant mice have a wider entrainment range than the wild types. Spectral analysis furthermore revealed nonlinear phenomena of periodically forced self-sustained oscillators for which the entrainment range relates inversely to oscillator amplitude. Using the forced oscillator model to explain the observed differences in entrainment range between mutant and wild-type mice, we sought to quantify the overall oscillator amplitude of their clocks from the activity rhythms and found that mutant mice have weaker circadian clocks than wild types. Our results suggest that genetic redundancy strengthens the circadian clock leading to a narrow entrainment range in vertebrates.

Prevalent role of gene features in determining evolutionary fates of whole-genome duplication duplicated genes in flowering plants

The evolution of genes and genomes after polyploidization has been the subject of extensive studies in evolutionary biology and plant sciences. While a significant number of duplicated genes are rapidly removed during a process called fractionation, which operates after the whole-genome duplication (WGD), another considerable number of genes are retained preferentially, leading to the phenomenon of biased gene retention. However, the evolutionary mechanisms underlying gene retention after WGD remain largely unknown. Through genome-wide analyses of sequence and functional data, we comprehensively investigated the relationships between gene features and the retention probability of duplicated genes after WGDs in six plant genomes, Arabidopsis (Arabidopsis thaliana), poplar (Populus trichocarpa), soybean (Glycine max), rice (Oryza sativa), sorghum (Sorghum bicolor), and maize (Zea mays). The results showed that multiple gene features were correlated with the probability of gene retention. Using a logistic regression model based on principal component analysis, we resolved evolutionary rate, structural complexity, and GC3 content as the three major contributors to gene retention. Cluster analysis of these features further classified retained genes into three distinct groups in terms of gene features and evolutionary behaviors. Type I genes are more prone to be selected by dosage balance; type II genes are possibly subject to subfunctionalization; and type III genes may serve as potential targets for neofunctionalization. This study highlights that gene features are able to act jointly as primary forces when determining the retention and evolution of WGD-derived duplicated genes in flowering plants. These findings thus may help to provide a resolution to the debate on different evolutionary models of gene fates after WGDs.

Convergent gene loss following gene and genome duplications creates single-copy families in flowering plants

The importance of gene gain through duplication has long been appreciated. In contrast, the importance of gene loss has only recently attracted attention. Indeed, studies in organisms ranging from plants to worms and humans suggest that duplication of some genes might be better tolerated than that of others. Here we have undertaken a large-scale study to investigate the existence of duplication-resistant genes in the sequenced genomes of 20 flowering plants. We demonstrate that there is a large set of genes that is convergently restored to single-copy status following multiple genome-wide and smaller scale duplication events. We rule out the possibility that such a pattern could be explained by random gene loss only and therefore propose that there is selection pressure to preserve such genes as singletons. This is further substantiated by the observation that angiosperm single-copy genes do not comprise a random fraction of the genome, but instead are often involved in essential housekeeping functions that are highly conserved across all eukaryotes. Furthermore, single-copy genes are generally expressed more highly and in more tissues than non-single-copy genes, and they exhibit higher sequence conservation. Finally, we propose different hypotheses to explain their resistance against duplication.

On the expansion of "dangerous" gene repertoires by whole-genome duplications in early vertebrates

The emergence and evolutionary expansion of gene families implicated in cancers and other severe genetic diseases is an evolutionary oddity from a natural selection perspective. Here, we show that gene families prone to deleterious mutations in the human genome have been preferentially expanded by the retention of "ohnolog" genes from two rounds of whole-genome duplication (WGD) dating back from the onset of jawed vertebrates. We further demonstrate that the retention of many ohnologs suspected to be dosage balanced is in fact indirectly mediated by their susceptibility to deleterious mutations. This enhanced retention of "dangerous" ohnologs, defined as prone to autosomal-dominant deleterious mutations, is shown to be a consequence of WGD-induced speciation and the ensuing purifying selection in post-WGD species. These findings highlight the importance of WGD-induced nonadaptive selection for the emergence of vertebrate complexity, while rationalizing, from an evolutionary perspective, the expansion of gene families frequently implicated in genetic disorders and cancers.

The odds of duplicate gene persistence after polyploidization

Background

Gene duplication is an important biological phenomenon associated with genomic redundancy, degeneration, specialization, innovation, and speciation. After duplication, both copies continue functioning when natural selection favors duplicated protein function or expression, or when mutations make them functionally distinct before one copy is silenced.

Results

Here we quantify the degree to which genetic parameters related to gene expression, molecular evolution, and gene structure in a diploid frog - Silurana tropicalis - influence the odds of functional persistence of orthologous duplicate genes in a closely related tetraploid species - Xenopus laevis. Using public databases and 454 pyrosequencing, we obtained genetic and expression data from S. tropicalis orthologs of 3,387 X. laevis paralogs and 4,746 X. laevis singletons - the most comprehensive dataset for African clawed frogs yet analyzed. Using logistic regression, we demonstrate that the most important predictors of the odds of duplicate gene persistence in the tetraploid species are the total gene expression level and evenness of expression across tissues and development in the diploid species. Slow protein evolution and information density (fewer exons, shorter introns) in the diploid are also positively correlated with duplicate gene persistence in the tetraploid.

Conclusions

Our findings suggest that a combination of factors contribute to duplicate gene persistence following whole genome duplication, but that the total expression level and evenness of expression across tissues and through development before duplication are most important. We speculate that these parameters are useful predictors of duplicate gene longevity after whole genome duplication in other taxa.

A generalized model of gene dosage and dominant negative effects in macromolecular complexes

Most genes and their corresponding products are supposed to be involved in genetic or biochemical interactions. A typical example is provided by macromolecular complexes, which may contain tens of proteins in defined stoichiometric proportions. Stoichiometric imbalances in such complexes can be a source of abnormal phenotypes. Comparable effects can also arise from negative dominance/transdominance, even though the underlying mechanisms are different. Here I propose a general yet simple biochemical model accounting for the effects of dosage changes and weak dominant/transdominant negative mutations in macromolecular complexes. The molecular alterations studied are predicted to lead to synergistic effects that can drive total multimer concentration and/or activity in a multiple heterozygote below a critical threshold required to ensure a normal phenotype, thus providing an explanation for the phenomenon of unlinked noncomplementation or nonallelic noncomplementation. The model also helps in understanding the basis of heterosis and the long-term consequences of gene dosage alterations and weak dominant/transdominant negative effects. Indeed, it can explain the observed extensive retention of paralogs in polyploids. Finally, because the effects of weak single-gene alterations can escape selection, they may accumulate in the population. This situation has important evolutionary consequences and may eventually lead to reproductive isolation and speciation.-Veitia, R. A. A generalized model of gene dosage and dominant negative effects in macromolecular complexes.

Genetics of athletic performance

Performance enhancing polymorphisms (PEPs) are examples of natural genetic variation that affect the outcome of athletic challenges. Elite athletes, and what separates them from the average competitor, have been the subjects of discussion and debate for decades. While training, diet, and mental fitness are all clearly important contributors to achieving athletic success, the fact that individuals reaching the pinnacle of their chosen sports often share both physical and physiological attributes suggests a role for genetics. That multiple members of a family often participate in highly competitive events, such as the Olympics, further supports this argument. In this review, we discuss what is known regarding the genes and gene families, including the mitochondrial genome, that are believed to play a role in human athletic performance. Where possible, we describe the physiological impact of the critical gene variants and consider predictions about other potentially important genes. Finally, we discuss the implications of these findings on the future for competitive athletics.

Dominant negative factors in health and disease

The idea of dominant mutations that interfere with the activity of a normal gene product has been known for more than 80 years-the famous Muller's antimorphs. However, only over half a century later, the mechanistic bases of dominant negative mutations (DNMs) were defined in a systematic way by Ira Herskowitz. Most analyses of DNMs consider only intralocus (interallelic) interactions. The typical textbook explanation invokes a defective subunit, which poisons a homo-dimer or a homo-oligomer. More complex cases exist and the quantitative dimension of this phenomenon will be explored here. The basic ideas underlying DN effects can be (and should be) extended to included epistatic (interloci) interactions. Indeed, poisoning heteromeric macromolecular complexes is per se a matter of 'transdominant' negative effects. In this context, non-allelic non-complementation is also considered. Given the importance of DNMs in human disease and in the study of gene function, understanding how they work is essential for understanding pathology and for the design of effective DN molecules that can also prove useful in therapeutics. Finally, the existence and potential relevance of an increasing number of physiological DN protein isoforms is briefly discussed.

Duplicated proteasome subunit genes in Drosophila and their roles in spermatogenesis

The proteasome is a large, multisubunit complex that acts as the cell's 'protein-degrading machine' in the ubiquitin-mediated proteolytic pathway for regulated protein turnover. Although proteasomes are usually thought of as being homogeneous structures, recent studies have revealed their more dynamic and heterogeneous nature. For example, in a number of plant and animal species, multiple isoforms of several proteasome subunits, encoded by paralogous genes, have been discovered, and in some cases, these alternative isoforms have been shown to be functionally distinct from their conventional counterparts. A particularly striking example of this phenomenon is seen in Drosophila melanogaster, where 12 of the 33 subunits that make up the 26S proteasome holoenzyme are represented in the genome by multiple paralogous genes. Remarkably, in every case, the 'extra' genes are expressed in a testis-specific manner. Here, we describe the extent and nature of these testis-specific gene duplications and discuss their functional significance, and speculate on why this situation might have evolved.

KEPE--a motif frequently superimposed on sumoylation sites in metazoan chromatin proteins and transcription factors

Motivation: We noted that the sumoylation site in C/EBP homologues is conserved beyond the canonical consensus sequence for sumoylation. Therefore, we investigated whether this pattern might define a more general protein motif.

Results: We undertook a survey of the human proteome using a regular expression based on the C/EBP motif. This revealed significant enrichment of the motif using different Gene Ontology terms (e.g. ‘transcription’) that pertain to the nucleus. When considering requirements for the motif to be functional (evolutionary conservation, structural accessibility of the motif and proper cell localization of the protein), more than 130 human proteins were retrieved from the UniProt/Swiss-Prot database. These candidates were particularly enriched in transcription factors, including FOS, JUN, Hif-1α, MLL2 and members of the KLF, MAF and NFATC families; chromatin modifiers like CHD-8, HDAC4 and DNA Top1; and the transcriptional regulatory kinases HIPK1 and HIPK2. The KEPEmotif appears to be restricted to the metazoan lineage and has three length variants—short, medium and long—which do not appear to interchange.

Contact:toby.gibson@embl.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Preferential subfunctionalization of slow-evolving genes after allopolyploidization in Xenopus laevis

As paleopolyploid genomes evolve, the expression profiles of retained gene pairs are expected to diverge. To examine this divergence process on a large scale in a vertebrate system, we compare Xenopus laevis, which has retained approximately 40% of loci in duplicate after a recent whole-genome duplication (WGD), with its unduplicated relative Silurana (Xenopus) tropicalis. This comparison of ingroup pairs to an outgroup allows the direction of change in expression profiles to be inferred for a set of 1,300 X. laevis pairs, relative to their single orthologs in S. tropicalis, across 11 tissues. We identify 68 pairs in which X. laevis is inferred to have undergone a significant reduction of expression in at least two tissues since WGD. Of these pairs, one-third show evidence of subfunctionalization, with decreases in the expression levels of different gene copies in two different tissues. Surprisingly, we find that genes with slow rates of evolution are particularly prone to subfunctionalization, even when the tendency for highly expressed genes to evolve slowly is controlled for. We interpret this result to be an effect of allopolyploidization. We then compare the outcomes of this WGD with an independent one that happened in the teleost fish lineage. We find that if a gene pair was retained in duplicate in X. laevis, the orthologous pair is more likely to have been retained in duplicate in zebrafish, suggesting that similar factors, among them subfunctionalization, determined which gene pairs survived in duplicate after the two WGDs.

Duplicate gene evolution and expression in the wake of vertebrate allopolyploidization

Background

The mechanism by which duplicate genes originate – whether by duplication of a whole genome or of a genomic segment – influences their genetic fates. To study events that trigger duplicate gene persistence after whole genome duplication in vertebrates, we have analyzed molecular evolution and expression of hundreds of persistent duplicate gene pairs in allopolyploid clawed frogs (Xenopus and Silurana). We collected comparative data that allowed us to tease apart the molecular events that occurred soon after duplication from those that occurred later on. We also quantified expression profile divergence of hundreds of paralogs during development and in different tissues.

Results

Our analyses indicate that persistent duplicates generated by allopolyploidization are subjected to strong purifying selection soon after duplication. The level of purifying selection is relaxed compared to a singleton ortholog, but not significantly variable over a period spanning about 40 million years. Despite persistent functional constraints, however, analysis of paralogous expression profiles indicates that quantitative aspects of their expression diverged substantially during this period.

Conclusion

These results offer clues into how vertebrate transcriptomes are sculpted in the wake of whole genome duplication (WGD), such as those that occurred in our early ancestors. That functional constraints were relaxed relative to a singleton ortholog but not significantly different in the early compared to the later stage of duplicate gene evolution suggests that the timescale for a return to pre-duplication levels is drawn out over tens of millions of years – beyond the age of these tetraploid species. Quantitative expression divergence can occur soon after WGD and with a magnitude that is not correlated with the rate of protein sequence divergence. On a coarse scale, quantitative expression divergence appears to be more prevalent than spatial and temporal expression divergence, and also faster or more frequent than other processes that operate at the protein level, such as some types of neofunctionalization.

The evolutionary dynamics of the Saccharomyces cerevisiae protein interaction network after duplication

Gene duplication is an important mechanism in the evolution of protein interaction networks. Duplications are followed by the gain and loss of interactions, rewiring the network at some unknown rate. Because rewiring is likely to change the distribution of network motifs within the duplicated interaction set, it should be possible to study network rewiring by tracking the evolution of these motifs. We have developed a mathematical framework that, together with duplication data from comparative genomic and proteomic studies, allows us to infer the connectivity of the preduplication network and the changes in connectivity over time. We focused on the whole-genome duplication (WGD) event in Saccharomyces cerevisiae. The model allowed us to predict the frequency of intergene interaction before WGD and the post duplication probabilities of interaction gain and loss. We find that the predicted frequency of self-interactions in the preduplication network is significantly higher than that observed in today's network. This could suggest a structural difference between the modern and ancestral networks, preferential addition or retention of interactions between ohnologs, or selective pressure to preserve duplicates of self-interacting proteins.

Comparative genomics of transcription factors and chromatin proteins in parasitic protists and other eukaryotes

Comparative genomics of parasitic protists and their free-living relatives are profoundly impacting our understanding of the regulatory systems involved in transcription and chromatin dynamics. While some parts of these systems are highly conserved, other parts are rapidly evolving, thereby providing the molecular basis for the variety in the regulatory adaptations of eukaryotes. The gross number of specific transcription factors and chromatin proteins are positively correlated with proteome size in eukaryotes. However, the individual types of specific transcription factors show an enormous variety across different eukaryotic lineages. The dominant families of specific transcription factors even differ between sister lineages, and have been shaped by gene loss and lineage-specific expansions. Recognition of this principle has helped in identifying the hitherto unknown, major specific transcription factors of several parasites, such as apicomplexans, Entamoeba histolytica, Trichomonas vaginalis, Phytophthora and ciliates. Comparative analysis of predicted chromatin proteins from protists allows reconstruction of the early evolutionary history of histone and DNA modification, nucleosome assembly and chromatin-remodeling systems. Many key catalytic, peptide-binding and DNA-binding domains in these systems ultimately had bacterial precursors, but were put together into distinctive regulatory complexes that are unique to the eukaryotes. In the case of histone methylases, histone demethylases and SWI2/SNF2 ATPases, proliferation of paralogous families followed by acquisition of novel domain architectures, seem to have played a major role in producing a diverse set of enzymes that create and respond to an epigenetic code of modified histones. The diversification of histone acetylases and DNA methylases appears to have proceeded via repeated emergence of new versions, most probably via transfers from bacteria to different eukaryotic lineages, again resulting in lineage-specific diversity in epigenetic signals. Even though the key histone modifications are universal to eukaryotes, domain architectures of proteins binding post-translationally modified-histones vary considerably across eukaryotes. This indicates that the histone code might be "interpreted" differently from model organisms in parasitic protists and their relatives. The complexity of domain architectures of chromatin proteins appears to have increased during eukaryotic evolution. Thus, Trichomonas, Giardia, Naegleria and kinetoplastids have relatively simple domain architectures, whereas apicomplexans and oomycetes have more complex architectures. RNA-dependent post-transcriptional silencing systems, which interact with chromatin-level regulatory systems, show considerable variability across parasitic protists, with complete loss in many apicomplexans and partial loss in Trichomonas vaginalis. This evolutionary synthesis offers a robust scaffold for future investigation of transcription and chromatin structure in parasitic protists.

Molecular evolution and genome divergence at RPB2 gene of the St and H genome in Elymus species

Molecular evolution of the second largest subunit of low copy nuclear RNA polymerase II (RPB2) in allotetrploid StH genomic species of Elymus is characterized here. Our study first reported a 39-bp MITE stowaway element insertion in the genic region of RPB2 gene for all tetraploid Elymus St genome and diploid Pseudoroegneria spicata and P. stipifolia St genome. The sequences on 3'-end are highly conserved, with AGTA in all sequences but H10339 (E. fibrosis), in which the AGTA was replaced with AGAA. All 12 Stowaway-containing sequences encompassed a 9 bp conserved TIRs (GAGGGAGTA). Interestingly, the 5'-end sequence of GGTA which was changed to AGTA or deleted resulted in Stowaway excision in the H genome of Elymus sepcies, in which Stowaway excision did not leave footprint. Another two large insertions in all St genome sequences are also transposable-like elements detected in the genic region of RPB2 gene. Our results indicated that these three transposable element indels have occurred prior to polyploidization, and shaped the homoeologous RPB2 loci in St and H genome of Eymus species. Nucleotide diversity analysis suggested that the RPB2 sequence may evolve faster in the polyploid species than in the diploids. Higher level of polymorphism and genome-specific amplicons generated by this gene indicated that RPB2 is an excellent tool for investigating the phylogeny and evolutionary dynamics of speciation, and the mode of polyploidy formation in Elymus species.

Functionally redundant SHI family genes regulate Arabidopsis gynoecium development in a dose-dependent manner

Gene duplication events, and the subsequent functional divergence of duplicates, are believed to be important evolutionary agents, driving morphological diversification. We have studied the structural and functional diversification of members of a plant-specific gene family in Arabidopsis thaliana by analysing mutant phenotypes, expression patterns and phylogeny. The SHI gene family comprises ten members that encode proteins with a RING finger-like zinc finger motif. We show that, despite being highly divergent in sequence, except in two conserved regions, many of the SHI-related genes are partially redundant in function and synergistically promote gynoecium, stamen and leaf development in Arabidopsis. Gynoecia of the loss-of-function sty1-1 mutant display subtle morphological defects, and, although mutations in the related STY2, SHI, SRS3, SRS4, SRS5, SRS7 and LRP1 genes have no apparent effect on gynoecium development, the sty1-1 mutant phenotype is gradually enhanced in double, triple, quadruple and quintuple mutant combinations, suggesting a remarkably extensive functional conservation within the family, which appears to be based on dosage dependency and protection against dominant negative mutations. In multiple mutant lines, all marginal tissues in the apical part of the gynoecium are dramatically reduced or missing, and our data indicate that SHI family members may promote formation of these tissues downstream of the transcriptional co-repressor LEUNIG (LUG).

Multiple mechanisms promote the retained expression of gene duplicates in the tetraploid frog Xenopus laevis

Gene duplication provides a window of opportunity for biological variants to persist under the protection of a co-expressed copy with similar or redundant function. Duplication catalyzes innovation (neofunctionalization), subfunction degeneration (subfunctionalization), and genetic buffering (redundancy), and the genetic survival of each paralog is triggered by mechanisms that add, compromise, or do not alter protein function. We tested the applicability of three types of mechanisms for promoting the retained expression of duplicated genes in 290 expressed paralogs of the tetraploid clawed frog, Xenopus laevis. Tests were based on explicit expectations concerning the ka/ks ratio, and the number and location of nonsynonymous substitutions after duplication. Functional constraints on the majority of paralogs are not significantly different from a singleton ortholog. However, we recover strong support that some of them have an asymmetric rate of nonsynonymous substitution: 6% match predictions of the neofunctionalization hypothesis in that (1) each paralog accumulated nonsynonymous substitutions at a significantly different rate and (2) the one that evolves faster has a higher ka/ks ratio than the other paralog and than a singleton ortholog. Fewer paralogs (3%) exhibit a complementary pattern of substitution at the protein level that is predicted by enhancement or degradation of different functional domains, and the remaining 13% have a higher average ka/ks ratio in both paralogs that is consistent with altered functional constraints, diversifying selection, or activity-reducing mutations after duplication. We estimate that these paralogs have been retained since they originated by genome duplication between 21 and 41 million years ago. Multiple mechanisms operate to promote the retained expression of duplicates in the same genome, in genes in the same functional class, over the same period of time following duplication, and sometimes in the same pair of paralogs. None of these paralogs are superfluous; degradation or enhancement of different protein subfunctions and neofunctionalization are plausible hypotheses for the retained expression of some of them. Evolution of most X. laevis paralogs, however, is consistent with retained expression via mechanisms that do not radically alter functional constraints, such as selection to preserve post-duplication stoichiometry or temporal, quantitative, or spatial subfunctionalization.

Gene duplication plays a fundamental role in biological innovation but it is not clear how both copies of a duplicated gene manage to circumvent degradation by mutation if neither is unique. This study explores genetic mechanisms that could make each copy of a duplicate gene different, and therefore distinguishable and potentially preserved by natural selection. It is based on DNA sequences of the protein-coding region of 290 expressed duplicated genes in a frog, Xenopus laevis, that underwent complete duplication of its entire genome. Results provide evidence for multiple mechanisms acting within the same genome, within the same functional classes of genes, within the same period of time following duplication, and even on the same set of duplicated genes. Each copy of a duplicate gene may be subject to distinct evolutionary constraints, and this could be associated with degradation or enhancement of function. Functional constraints of most of these duplicates, however, are not substantially different from a single copy gene; their persistence in the first dozens of millions of years after duplication may more frequently be explained by mechanisms acting on their expression rather than their function.

Many genes in fish have species-specific asymmetric rates of molecular evolution

BACKGROUND

Gene and genome duplication events increase the amount of genetic material that might then contribute to an increase in the genomic and phenotypic complexity of organisms during evolution. Thus, it has been argued that there is a relationship between gene copy number and morphological complexity and/or species diversity. This hypothesis implies that duplicated genes have subdivided or evolved novel functions compared to their pre-duplication proto-orthologs. Such a functional divergence might be caused by an increase in evolutionary rates in one ortholog, by changes in expression, regulatory evolution, insertion of repetitive elements, or due to positive Darwinian selection in one copy. We studied a set of 2466 genes that were present in Danio rerio, Takifugu rubripes, Tetraodon nigroviridis and Oryzias latipes to test (i) for forces of positive Darwinian selection; (ii) how frequently duplicated genes are retained, and (iii) whether novel gene functions might have evolved.

RESULTS

25% (610) of all investigated genes show significantly smaller or higher genetic distances in the genomes of particular fish species compared to their human ortholog than their orthologs in other fish according to relative rate tests. We identified 49 new paralogous pairs of duplicated genes in fish, in which one of the paralogs is under positive Darwinian selection and shows a significantly higher rate of molecular evolution in one of the four fish species, whereas the other copy apparently did not undergo adaptive changes since it retained the original rate of evolution. Among the genes under positive Darwinian selection, we found a surprisingly high number of ATP binding proteins and transcription factors.

CONCLUSION

The significant rate difference suggests that the function of these rate-changed genes might be essential for the respective fish species. We demonstrate that the measurement of positive selection is a powerful tool to identify divergence rates of duplicated genes and that this method has the capacity to identify potentially interesting candidates for adaptive gene evolution.

Evolution of myelin proteolipid proteins: Gene duplication in teleosts and expression pattern divergence

The coevolution of neurons and their supporting glia to the highly specialized axon-myelin unit included the recruitment of proteolipids as neuronal glycoproteins (DMbeta, DMgamma) or myelin proteins (DMalpha/PLP/DM20). Consistent with a genome duplication at the root of teleosts, we identified three proteolipid pairs in zebrafish, termed DMalpha1 and DMalpha2, DMbeta1 and DMbeta2, DMgamma1 and DMgamma2. The paralogous amino acid sequences diverged remarkably after gene duplication, indicating functional specialization. Each proteolipid has adopted a distinct spatio-temporal expression pattern in neural progenitors, neurons, and in glia. DMalpha2, the closest homolog to mammalian PLP/DM20, is coexpressed with P0 in oligodendrocytes and upregulated after optic nerve lesion. DMgamma2 is expressed in multipotential stem cells, and the other four proteolipids are confined to subsets of CNS neurons. Comparing protein sequences and gene structures from birds, teleosts, one urochordate species, and four invertebrates, we have reconstructed major steps in the evolution of proteolipids.

Gene complexity and gene duplicability

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

Linear motifs: evolutionary interaction switches

Linear motifs are short sequence patterns associated with a particular function. They differ fundamentally from longer, globular protein domains in terms of their binding affinities, evolution and in how they are found experimentally or computationally. In this Minireview, we discuss various aspects of these critically important functional regions.

Natural selection during functional divergence to LMP7 and proteasome subunit X (PSMB5) following gene duplication

The LMP7 and PSMB5 genes were created through an ancient gene duplication event of their ancestral locus. These proteins contain an active site of proteolysis, and LMP7 replaces PSMB5 as a component of the 20S proteasome after stimulation of cells by interferon-gamma. Replacement of PSMB5 by LMP7 changes the profile of the products of 20S proteasome processing, predisposing digested peptides for transport to and display by the immune system. The purpose of this study is to investigate evolutionary forces influencing functional divergence between LMP7 and PSMB5 following duplication. Levels of synonymous and nonsynonymous substitution rates are estimated to infer differences in levels of natural selection. Estimates of substitution rates indicate that natural selection elevated rates of nonsynonymous substitution in LMP7 following gene duplication, whereas PSMB5 experienced an increase in substitution rate that was not likely due to diversifying natural selection following duplication. Following initial divergence, nearly neutral mutations have dominated gene evolution in both lineages. The LMP7 gene locus provides a rare example of a protein with specialized function arising from duplication and divergence of a housekeeping protein by way of natural selection.

From 2R to 3R: evidence for a fish-specific genome duplication (FSGD)

An important mechanism for the evolution of phenotypic complexity, diversity and innovation, and the origin of novel gene functions is the duplication of genes and entire genomes. Recent phylogenomic studies suggest that, during the evolution of vertebrates, the entire genome was duplicated in two rounds (2R) of duplication. Later, approximately 350 mya, in the stem lineage of ray-finned (actinopterygian) fishes, but not in that of the land vertebrates, a third genome duplication occurred-the fish-specific genome duplication (FSGD or 3R), leading, at least initially, to up to eight copies of the ancestral deuterostome genome. Therefore, the sarcopterygian (lobe-finned fishes and tetrapods) genome possessed originally only half as many genes compared to the derived fishes, just like the most-basal and species-poor lineages of extant fishes that diverged from the fish stem lineage before the 3R duplication. Most duplicated genes were secondarily lost, yet some evolved new functions. The genomic complexity of the teleosts might be the reason for their evolutionary success and astounding biological diversity.

Computational approaches to unveiling ancient genome duplications

Recent analyses of complete genome sequences have revealed that many genomes have been duplicated in their evolutionary past. Such events have been associated with important biological transitions, major leaps in evolution and adaptive radiations of species. Here, we consider recently developed computational methods to detect such ancient large-scale gene duplication events. Several new approaches have been used to show that large-scale gene duplications are more common than previously thought.

Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution

To study the evolutionary effects of polyploidy on plant gene functions, we analyzed functional genomics data for a large number of duplicated gene pairs formed by ancient polyploidy events in Arabidopsis thaliana. Genes retained in duplicate are not distributed evenly among Gene Ontology or Munich Information Center for Protein Sequences functional categories, which indicates a nonrandom process of gene loss. Genes involved in signal transduction and transcription have been preferentially retained, and those involved in DNA repair have been preferentially lost. Although the two members of each gene pair must originally have had identical transcription profiles, less than half of the pairs formed by the most recent polyploidy event still retain significantly correlated profiles. We identified several cases where groups of duplicated gene pairs have diverged in concert, forming two parallel networks, each containing one member of each gene pair. In these cases, the expression of each gene is strongly correlated with the other nonhomologous genes in its network but poorly correlated with its paralog in the other network. We also find that the rate of protein sequence evolution has been significantly asymmetric in >20% of duplicate pairs. Together, these results suggest that functional diversification of the surviving duplicated genes is a major feature of the long-term evolution of polyploids.

Preferential duplication of conserved proteins in eukaryotic genomes

A central goal in genome biology is to understand the origin and maintenance of genic diversity. Over evolutionary time, each gene's contribution to the genic content of an organism depends not only on its probability of long-term survival, but also on its propensity to generate duplicates that are themselves capable of long-term survival. In this study we investigate which types of genes are likely to generate functional and persistent duplicates. We demonstrate that genes that have generated duplicates in the C. elegans and S. cerevisiae genomes were 25%-50% more constrained prior to duplication than the genes that failed to leave duplicates. We further show that conserved genes have been consistently prolific in generating duplicates for hundreds of millions of years in these two species. These findings reveal one way in which gene duplication shapes the content of eukaryotic genomes. Our finding that the set of duplicate genes is biased has important implications for genome-scale studies.

Testican-1: A Differentially Expressed Proteoglycan with Protease Inhibiting Activities

Testican-1 is a highly conserved, multidomain proteoglycan that is most prominently expressed in the thalamus of the brain, and is upregulated in activated astroglial cells of the cerebrum. Several functions of this gene product have now been demonstrated in vitro including membrane-type matrix metalloproteinase inhibition, cathepsin L inhibition, and low-affinity calcium binding. The purified gene product has been shown to inhibit cell attachment and neurite extensions in culture. Functions of testican in vivo have yet to be demonstrated in knockout mice or other models. Testican has been shown to carry substantial amounts of chondroitin sulfate as well as other oligosaccharides, but the biological significance of these embellishments is not yet known.

Asymmetric sequence divergence of duplicate genes

Much like humans, gene duplicates may be created equal, but they do not stay that way for long. For four completely sequenced genomes we show that 20%-30% of duplicate gene pairs show asymmetric evolution in the amino acid sequence of their protein products. That is, one of the duplicates evolves much faster than the other. The greater this asymmetry, the greater the ratio Ka/Ks of amino acid substitutions (Ka) to silent substitutions (Ks) in a gene pair. This indicates that most asymmetric divergence may be caused by relaxed selective constraints on one of the duplicates. However, we also find some candidate duplicates where positive (directional) selection of beneficial mutations (Ka/Ks > 1) may play a role in asymmetric divergence. Our analysis rests on a codon-based model of molecular evolution that allows a test for asymmetric divergence in Ka. The method is also more sensitive in detecting positive selection (Ka/Ks > 1) than models relying only on pairwise gene comparisons.

Evolution of glutamine synthetase in vertebrates: multiple glutamine synthetase genes expressed in rainbow trout (Oncorhynchus mykiss)

Glutamine synthetase (GSase) is a key enzyme in nitrogen metabolism and encoded by a single gene in mammals. Using PCR cloning techniques, including RT-PCR from total RNA and PCR from a cDNA library, we find evidence of four expressed GSase mRNAs for the tetraploid rainbow trout. For two of these mRNAs (Onmy-GS01, -GS02) we characterize the full-length coding regions, and for two others (Onmy-GS03, -GS04), we describe partial sequences. Northern analysis of Onmy-GS01, -GS02, -GS03 and -GS04 indicates that (1) Onmy-GS02 is expressed at higher levels relative to the other transcripts in most adult tissues, with the exception of brain and gill, where Onmy-GS01 is at the highest level, and (2) the tissue with the highest level of expression of all four transcripts is the brain, with decreasing levels in the intestine, liver, red muscle, gill/kidney, white muscle and heart. Clearly, rainbow trout possess multiple GSase genes with differing levels of tissue expression, implying manifold potential routes of regulation for this octameric enzyme. Our data also indicate that caution should be taken when interpreting mRNA expression data of a single gene, unless multiple genes have been ruled out. Consistent with a southern blot, phylogenetic and intron sequence analyses imply that the trout genes are encoded by at least four separate loci, belonging to two distinct evolutionary branches. Our data on rainbow trout, together with those from two full-length zebrafish Danio rerio GSase genes compiled from GenBank ESTs, support the idea that fish GSases are polyphyletic and that gene duplications have occurred at multiple points and in independent lineages throughout the evolution of bony fishes.

Role of duplicate genes in genetic robustness against null mutations

Deleting a gene in an organism often has little phenotypic effect, owing to two mechanisms of compensation. The first is the existence of duplicate genes: that is, the loss of function in one copy can be compensated by the other copy or copies. The second mechanism of compensation stems from alternative metabolic pathways, regulatory networks, and so on. The relative importance of the two mechanisms has not been investigated except for a limited study, which suggested that the role of duplicate genes in compensation is negligible. The availability of fitness data for a nearly complete set of single-gene-deletion mutants of the Saccharomyces cerevisiae genome has enabled us to carry out a genome-wide evaluation of the role of duplicate genes in genetic robustness against null mutations. Here we show that there is a significantly higher probability of functional compensation for a duplicate gene than for a singleton, a high correlation between the frequency of compensation and the sequence similarity of two duplicates, and a higher probability of a severe fitness effect when the duplicate copy that is more highly expressed is deleted. We estimate that in S. cerevisiae at least a quarter of those gene deletions that have no phenotype are compensated by duplicate genes.

Detecting the undetectable: uncovering duplicated segments in Arabidopsis by comparison with rice

Genome analysis shows that large-scale gene duplications have occurred in fungi, animals and plants, creating genomic regions that show similarity in gene content and order. However, the high frequency of gene loss reduces colinearity resulting in duplicated regions that, in the extreme, no longer share homologous genes. Here, we show that by comparison with an appropriate second genome, such paralogous regions can still be identified.

The hidden duplication past of Arabidopsis thaliana

Analysis of the genome sequence of Arabidopsis thaliana shows that this genome, like that of many other eukaryotic organisms, has undergone large-scale gene duplications or even duplications of the entire genome. However, the high frequency of gene loss after duplication events reduces colinearity and therefore the chance of finding duplicated regions that, at the extreme, no longer share homologous genes. In this study we show that heavily degenerated block duplications that can no longer be recognized by directly comparing two segments because of differential gene loss, can still be detected through indirect comparison with other segments. When these so-called hidden duplications in Arabidopsis are taken into account, many homologous genomic regions can be found in five to eight copies. This finding strongly implies that Arabidopsis has undergone three, but probably no more, rounds of genome duplications. Therefore, adding such hidden blocks to the duplication landscape of Arabidopsis sheds light on the number of polyploidy events that this model plant genome has undergone in its evolutionary past.

Asymmetric Functional Divergence of Duplicate Genes in Yeast

Most duplicate genes are eliminated from a genome shortly after duplication, but those that remain are an important source of biochemical diversity. Here, I present evidence from genome-scale protein-protein interaction data, microarray expression data, and large-scale gene knockout data that this diversification is often asymmetrical: one duplicate usually shows significantly more molecular or genetic interactions than the other. I propose a model that can explain this divergence pattern if asymmetrically diverging duplicate gene pairs show increased robustness to deleterious mutations.

The Hox Paradox: More complex(es) than imagined

An understanding of the origin of different body plans requires knowledge of how the genes and genetic pathways that control embryonic development have evolved. The Hox genes provide an appealing starting point for such studies because they play a well-understood causal role in the regionalization of the body plan of all bilaterally symmetric animals. Vertebrate evolution has been characterized by gene, and possibly genome, duplication events, which are believed to have provided raw genetic material for selection to act upon. It has recently been established that the Hox gene organization of ray-finned fishes, such as the zebrafish, differs dramatically from that of their lobe-finned relatives, a group that includes humans and all the other widely used vertebrate model systems. This unusual Hox gene organization of zebrafish is the result of a duplication event within the ray-finned fish lineage. Thus, teleosts, such as zebrafish, have more Hox genes arrayed over more clusters (or "complexes") than do tetrapod vertebrates. Here, I review our understanding of Hox cluster architecture in different vertebrates and consider the implications of gene duplication for Hox gene regulation and function and the evolution of different body plans.

Gene family evolution and homology: genomics meets phylogenetics

With the advent of high-throughput DNA sequencing and whole-genome analysis, it has become clear that the coding portions of the genome are organized hierarchically in gene families and superfamilies. Because the hierarchy of genes, like that of living organisms, reflects an ancient and continuing process of gene duplication and divergence, many of the conceptual and analytical tools used in phylogenetic systematics can and should be used in comparative genomics. Phylogenetic principles and techniques for assessing homology, inferring relationships among genes, and reconstructing evolutionary events provide a powerful way to interpret the ever increasing body of sequence data. In this review, we outline the application of phylogenetic approaches to comparative genomics, beginning with the inference of phylogeny and the assessment of gene orthology and paralogy. We also show how the phylogenetic approach makes possible novel kinds of comparative analysis, including detection of domain shuffling and lateral gene transfer, reconstruction of the evolutionary diversification of gene families, tracing of evolutionary change in protein function at the amino acid level, and prediction of structure-function relationships. A marriage of the principles of phylogenetic systematics with the copious data generated by genomics promises unprecedented insights into the nature of biological organization and the historical processes that created it.