Yeast genome evolution in the post-genome era.

The Saccharomyces cerevisiae genome sequence, augmented by new data on gene expression and function, continues to yield new findings about eukaryote genome evolution. Analysis of the duplicate gene pairs formed by whole-genome duplication indicates that selection for increased levels of gene expression was a significant factor determining which genes were retained as duplicates and which were returned to a single-copy state, possibly in addition to selection for novel gene functions. Proteome comparisons between worm and yeast show that genes for core metabolic processes are shared among eukaryotes and unchanging in function, while comparisons between different yeast species identify 'orphan' genes as the most rapidly evolving fraction of the proteome. Natural hybridisation among yeast species is frequent, but its long-term evolutionary significance is unknown.     

Pervasive indels and their evolutionary dynamics after the fish-specific genome duplication

Insertions and deletions (indels) in protein-coding genes are important sources of genetic variation. Their role in creating new proteins may be especially important after gene duplication. However, little is known about how indels affect the divergence of duplicate genes. We here study thousands of duplicate genes in five fish (teleost) species with completely sequenced genomes. The ancestor of these species has been subject to a fish-specific genome duplication (FSGD) event that occurred approximately 350 Ma. We find that duplicate genes contain at least 25% more indels than single-copy genes. These indels accumulated preferentially in the first 40 my after the FSGD. A lack of widespread asymmetric indel accumulation indicates that both members of a duplicate gene pair typically experience relaxed selection. Strikingly, we observe a 30-80% excess of deletions over insertions that is consistent for indels of various lengths and across the five genomes. We also find that indels preferentially accumulate inside loop regions of protein secondary structure and in regions where amino acids are exposed to solvent. We show that duplicate genes with high indel density also show high DNA sequence divergence. Indel density, but not amino acid divergence, can explain a large proportion of the tertiary structure divergence between proteins encoded by duplicate genes. Our observations are consistent across all five fish species. Taken together, they suggest a general pattern of duplicate gene evolution in which indels are important driving forces of evolutionary change.     

Loss of LIN-35, the Caenorhabditis elegans ortholog of the tumor suppressor p105Rb, results in enhanced RNA interference

Background

Genome-wide RNA interference (RNAi) screening is a very powerful tool for analyzing gene function in vivo in Caenorhabditis elegans. The effectiveness of RNAi varies from gene to gene, however, and neuronally expressed genes are largely refractive to RNAi in wild-type worms.

Results

We found that C. elegans strains carrying mutations in lin-35, the worm ortholog of the tumor suppressor gene p105Rb, or a subset of the genetically related synMuv B family of chromatin-modifying genes, show increased strength and penetrance for many germline, embryonic, and post-embryonic RNAi phenotypes, including neuronal RNAi phenotypes. Mutations in these same genes also enhance somatic transgene silencing via an RNAi-dependent mechanism. Two genes, mes-4 and zfp-1, are required both for the vulval lineage defects resulting from mutations in synMuv B genes and for RNAi, suggesting a common mechanism for the function of synMuv B genes in vulval development and in regulating RNAi. Enhanced RNAi in the germline of lin-35 worms suggests that misexpression of germline genes in somatic cells cannot alone account for the enhanced RNAi observed in this strain.

Conclusion

A worm strain with a null mutation in lin-35 is more sensitive to RNAi than any other previously described single mutant strain, and so will prove very useful for future genome-wide RNAi screens, particularly for identifying genes with neuronal functions. As lin-35 is the worm ortholog of the mammalian tumor suppressor gene p105Rb, misregulation of RNAi may be important during human oncogenesis.     

Duplication of floral regulatory genes in the Lamiales

Duplication of some floral regulatory genes has occurred repeatedly in angiosperms, whereas others are thought to be single-copy in most lineages. We selected three genes that interact in a pathway regulating floral development conserved among higher tricolpates (LFY/FLO, UFO/FIM, and AP3/DEF) and screened for copy number among families of Lamiales that are closely related to the model species Antirrhinum majus. We show that two of three genes have duplicated at least twice in the Lamiales. Phylogenetic analyses of paralogs suggest that an ancient whole genome duplication shared among many families of Lamiales occurred after the ancestor of these families diverged from the lineage leading to Veronicaceae (including the single-copy species A. majus). Duplication is consistent with previous patterns among angiosperm lineages for AP3/DEF, but this is the first report of functional duplicate copies of LFY/FLO outside of tetraploid species. We propose Lamiales taxa will be good models for understanding mechanisms of duplicate gene preservation and how floral regulatory genes may contribute to morphological diversity.     

Drosophila duplicate genes evolve new functions on the fly

Gene duplication is thought to play a key role in phenotypic innovation. While several processes have been hypothesized to drive the retention and functional evolution of duplicate genes, their genomic contributions have never been determined. We recently developed the first genome-wide method to classify these processes by comparing distances between expression profiles of duplicate genes and their ancestral single-copy orthologs. Application of our approach to spatial gene expression profiles in two Drosophila species revealed that a majority of young duplicate genes possess new functions, and that new functions are acquired rapidly—often within a few million years. Surprisingly, new functions tend to arise in younger copies of duplicate gene pairs. Moreover, we found that young duplicates are often specifically expressed in testes, whereas old duplicates are broadly expressed across several tissues, providing strong support for the hypothetical “out-of-testes” origin of new genes. In this Extra View, I discuss our findings in the context of theoretical predictions about gene duplication, with a particular emphasis on the importance of natural selection in the evolution of novel phenotypes.     

Gene Duplicability of Core Genes Is Highly Consistent across All Angiosperms

Gene duplication is an important mechanism for adding to genomic novelty. Hence, which genes undergo duplication and are preserved following duplication is an important question. It has been observed that gene duplicability, or the ability of genes to be retained following duplication, is a nonrandom process, with certain genes being more amenable to survive duplication events than others. Primarily, gene essentiality and the type of duplication (small-scale versus large-scale) have been shown in different species to influence the (long-term) survival of novel genes. However, an overarching view of “gene duplicability” is lacking, mainly due to the fact that previous studies usually focused on individual species and did not account for the influence of genomic context and the time of duplication. Here, we present a large-scale study in which we investigated duplicate retention for 9178 gene families shared between 37 flowering plant species, referred to as angiosperm core gene families. For most gene families, we observe a strikingly consistent pattern of gene duplicability across species, with gene families being either primarily single-copy or multicopy in all species. An intermediate class contains gene families that are often retained in duplicate for periods extending to tens of millions of years after whole-genome duplication, but ultimately appear to be largely restored to singleton status, suggesting that these genes may be dosage balance sensitive. The distinction between single-copy and multicopy gene families is reflected in their functional annotation, with single-copy genes being mainly involved in the maintenance of genome stability and organelle function and multicopy genes in signaling, transport, and metabolism. The intermediate class was overrepresented in regulatory genes, further suggesting that these represent putative dosage-balance-sensitive genes.     

Protein complexity, gene duplicability and gene dispensability in the yeast genome

Using functional genomic and protein structural data we studied the effects of protein complexity (here defined as the number of subunit types in a protein) on gene dispensability and gene duplicability. We found that in terms of gene duplicability the major distinction in protein complexity is between hetero-complexes, each of which includes at least two different types of subunits (polypeptides), and homo-complexes, which include monomers and complexes that consist of only subunits of one polypeptide type. However, gene dispensability decreases only gradually as the number of subunit types in a protein complex increases. These observations suggest that the dosage balance hypothesis can explain well gene duplicability of complex proteins, but cannot completely explain the difference in dispensabilities between hetero-complex subunits. It is likely that knocking out a gene coding for a hetero-complex subunit would disrupt the function of the whole complex, so that the deletion effect on fitness would increase with protein complexity. We also found that multi-domain polypeptide genes are less dispensable but more duplicable than single-domain polypeptide genes. Duplicate genes derived from the whole genome duplication event in yeast are more dispensable (except for ribosomal protein genes) than other duplicate genes. Further, we found that subunits of the same protein complex tend to have similar expression levels and similar effects of gene deletion on fitness. Finally, we estimated that in yeast the contribution of duplicate genes to genetic robustness against null mutation is approximately 9%, smaller than previously estimated. In yeast, protein complexity may serve as a better indicator of gene dispensability than do duplicate genes.     

Higher duplicability of less important genes in yeast genomes

Gene duplication plays an important role in evolution because it is the primary source of new genes. Many recent studies showed that gene duplicability varies considerably among genes. Several considerations led us to hypothesize that less important genes have higher rates of successful duplications, where gene importance is measured by the fitness reduction caused by the deletion of the gene. Here, we test this hypothesis by comparing the importance of two groups of singleton genes in the yeast Saccharomyces cerevisiae (Sce). Group S genes did not duplicate in four other yeast species examined, whereas group D experienced duplication in these species. Consistent with our hypothesis, we found group D genes to be less important than group S genes. Specifically, 17% of group D genes are essential in Sce, compared to 28% for group S. Furthermore, deleting a group D gene in Sce reduces the fitness by 24% on average, compared to 38% for group S. Our subsequent analysis showed that less important genes have more cis-regulatory motifs, which could lead to a higher chance of subfunctionalization of duplicate genes and result in an enhanced rate of gene retention. Less important genes may also have weaker dosage imbalance effects and cause fewer genetic perturbations when duplicated. Regardless of the cause, our observation indicates that the previous finding of a less severe fitness consequence of deleting a duplicate gene than deleting a singleton gene is at least in part due to the fact that duplicate genes are intrinsically less important than singleton genes and suggests that the contribution of duplicate genes to genetic robustness has been overestimated.     

Predicting the Proportion of Essential Genes in Mouse Duplicates Based on Biased Mouse Knockout Genes

In the yeast or nematode, the proportion of essential genes in duplicates is lower than in singletons (single-copy genes), due to the functional redundancy. One may expect that it should be the same in the mouse genome. However, based on the publicly available mouse knockout data, it was observed that the proportion of essential genes in duplicates is similar to that in singletons. The most straightforward interpretation, as claimed in a recent study, is that duplicate genes may have a negligible role in the mouse genetic robustness. Here we show that in the current mouse knockout dataset, recently duplicated genes have been highly underrepresented, leading to an overestimation of the proportion of essential genes in duplicates. After estimating the duplication time of mouse duplication events, we have developed a simple bias-correcting procedure and shown that the bias-corrected proportion of essential genes in mouse duplicates is significantly lower than that in singletons.     

Uncover Genetic Interactions in Caenorhabditis elegans by RNA Interference

RNA-mediated interference (RNAi) has emerged recently as one of the most powerful functional genomics tools. RNAi has been particularly effective in the nematode worm C. elegans where RNAi has been used to analyse the loss-of-function phenotypes of almost all predicted genes. In this review, we illustrate how RNAi has been used to analyse gene function in C. elegans as well as pointing to some future directions for using RNAi to examine genetic interactions in a systematic manner.     

The TRiC/CCT Chaperone Is Implicated in Alzheimer's Disease Based on Patient GWAS and an RNAi Screen in Aβ-Expressing Caenorhabditis elegans

The human Aβ peptide causes progressive paralysis when expressed in the muscles of the nematode worm, C. elegans. We have exploited this model of Aβ toxicity by carrying out an RNAi screen to identify genes whose reduced expression modifies the severity of this locomotor phenotype. Our initial finding was that none of the human orthologues of these worm genes is identical with the genome-wide significant GWAS genes reported to date (the “white zone”); moreover there was no identity between worm screen hits and the longer list of GWAS genes which included those with borderline levels of significance (the “grey zone”). This indicates that Aβ toxicity should not be considered as equivalent to sporadic AD. To increase the sensitivity of our analysis, we then considered the physical interactors (+1 interactome) of the products of the genes in both the worm and the white+grey zone lists. When we consider these worm and GWAS gene lists we find that 4 of the 60 worm genes have a +1 interactome overlap that is larger than expected by chance. Two of these genes form a chaperonin complex, the third is closely associated with this complex and the fourth gene codes for actin, the major substrate of the same chaperonin.     

Predominant gain of promoter TATA box after gene duplication associated with stress responses

TATA box, the core promoter element, exists in a broad range of eukaryotes, and the expression of TATA-containing genes usually responds to various environmental stresses. Hence, the evolution of TATA-box in duplicate genes may provide some clues for the interrelationship among environmental stress, expression differentiation, and duplicate gene preservation. In the present study, we observed that the TATA box is significantly overrepresented in duplicate genes compared with singletons in human, worm, Arabidopsis, and yeast genomes. We then conducted an extensive functional genomic analysis to investigate the evolution of TATA box along over 700 yeast gene family phylogenies. After reconstructing the ancestral TATA-box states (presence or absence), we found that significantly higher numbers of TATA box gain events than loss events had occurred after yeast gene duplications-the overall gain-loss ratio is about 3-4 to 1. Interestingly, these TATA-gain duplicate genes on average have experienced greater expression divergence from the ancestral expression states than their most closely related TATA-less duplicate partners, but only under environmental stress conditions (asymmetric evolution); indeed, under normal physiological conditions, they have similar expression divergence (symmetric evolution). Moreover, we showed that TATA-gain duplicates are enriched in stress-associated functional categories but that is not the case for TATA-ancestral duplicates (those inherited from their ancestors prior to duplication). Together, we conclude that after the gene duplication, gain of the TATA box in duplicate promoters may have played an important role in yeast duplicate preservation by accelerating expression divergence that may facilitate the adaptive evolution of the organism in response to environmental changes.     

A forward genetic screen reveals essential and non-essential RNAi factors in Paramecium tetraurelia

In most eukaryotes, small RNA-mediated gene silencing pathways form complex interacting networks. In the ciliate Paramecium tetraurelia, at least two RNA interference (RNAi) mechanisms coexist, involving distinct but overlapping sets of protein factors and producing different types of short interfering RNAs (siRNAs). One is specifically triggered by high-copy transgenes, and the other by feeding cells with double-stranded RNA (dsRNA)-producing bacteria. In this study, we designed a forward genetic screen for mutants deficient in dsRNA-induced silencing, and a powerful method to identify the relevant mutations by whole-genome sequencing. We present a set of 47 mutant alleles for five genes, revealing two previously unknown RNAi factors: a novel Paramecium-specific protein (Pds1) and a Cid1-like nucleotidyl transferase. Analyses of allelic diversity distinguish non-essential and essential genes and suggest that the screen is saturated for non-essential, single-copy genes. We show that non-essential genes are specifically involved in dsRNA-induced RNAi while essential ones are also involved in transgene-induced RNAi. One of the latter, the RNA-dependent RNA polymerase RDR2, is further shown to be required for all known types of siRNAs, as well as for sexual reproduction. These results open the way for the dissection of the genetic complexity, interconnection, mechanisms and natural functions of RNAi pathways in P. tetraurelia.     

Distributed robustness versus redundancy as causes of mutational robustness

A biological system is robust to mutations if it continues to function after genetic changes in its parts. Such robustness is pervasive on different levels of biological organization, from macromolecules to genetic networks and whole organisms. I here ask which of two possible causes of such robustness are more important on a genome-wide scale, for systems whose parts are genes, such as metabolic and genetic networks. The first of the two causes is redundancy of a system's parts: A gene may be dispensable if the genome contains redundant, back-up copies of the gene. The second cause, distributed robustness, is more poorly understood. It emerges from the distributed nature of many biological systems, where many (and different) parts contribute to system functions. I will here discuss evidence suggesting that distributed robustness is equally or more important for mutational robustness than gene redundancy. This evidence comes from the functional divergence of redundant genes, as well as from large-scale gene deletion studies. I also ask whether one can quantify the extent to which redundancy or distributed robustness contribute to mutational robustness.     

Selection for distinct gene expression properties favours the evolution of mutational robustness in gene regulatory networks

Mutational robustness is a genotype's tendency to keep a phenotypic trait with little and few changes in the face of mutations. Mutational robustness is both ubiquitous and evolutionarily important as it affects in different ways the probability that new phenotypic variation arises. Understanding the origins of robustness is specially relevant for systems of development that are phylogenetically widespread and that construct phenotypic traits with a strong impact on fitness. Gene regulatory networks are examples of this class of systems. They comprise sets of genes that, through cross‐regulation, build the gene activity patterns that define cellular responses, different tissues or distinct cell types. Several empirical observations, such as a greater robustness of wild‐type phenotypes, suggest that stabilizing selection underlies the evolution of mutational robustness. However, the role of selection in the evolution of robustness is still under debate. Computer simulations of the dynamics and evolution of gene regulatory networks have shown that selection for any gene activity pattern that is steady and self‐sustaining is sufficient to promote the evolution of mutational robustness. Here, I generalize this scenario using a computational model to show that selection for different aspects of a gene activity phenotype increases mutational robustness. Mutational robustness evolves even when selection favours properties that conflict with the stationarity of a gene activity pattern. The results that I present support an important role for stabilizing selection in the evolution of robustness in gene regulatory networks.     

Validation of RNAi silencing specificity using synthetic genes: salicylic acid-binding protein 2 is required for innate immunity in plants

RNA interference (RNAi) is widely used to specifically silence the expression of any gene to study its function and to identify and validate therapeutic targets. Despite the popularity of this technology, recent studies have shown that RNAi may also silence non-targeted genes. Here we demonstrate the utility of a quick, efficient and robust approach to directly validate the specificity of RNAi as an alternative to indirect validation of RNAi through gene expression profiling. Our approach involves reversing (complementing) the RNAi-induced phenotype by introducing a synthetic version of the target gene that is designed to escape silencing. This synthetic gene complementation approach can also be used for mutational analysis of the target gene, or to provide a functional version of a defective protein after silencing the defective gene by RNAi. Using this approach we demonstrate that the loss of systemic acquired resistance, a form of innate immunity in plants, is indeed due to the silencing of salicylic acid-binding protein 2 rather than to off-target effects.