Modeling the evolution of a classic genetic switch

Background  

The regulatory network underlying the yeast galactose-use pathway has emerged as a model system for the study of regulatory network evolution. Evidence has recently been provided for adaptive evolution in this network following a whole genome duplication event. An ancestral gene encoding a bi-functional galactokinase and co-inducer protein molecule has become subfunctionalized as paralogous genes (GAL1 and GAL3) in Saccharomyces cerevisiae, with most fitness gains being attributable to changes in cis-regulatory elements. However, the quantitative functional implications of the evolutionary changes in this regulatory network remain unexplored.     

Global Expression Analysis of the Yeast Lachancea (Saccharomyces) kluyveri Reveals New URC Genes Involved in Pyrimidine Catabolism

Pyrimidines are important nucleic acid precursors which are constantly synthesized, degraded, and rebuilt in the cell. Four degradation pathways, two of which are found in eukaryotes, have been described. One of them, the URC pathway, has been initially discovered in our laboratory in the yeast Lachancea kluyveri. Here, we present the global changes in gene expression in L. kluyveri in response to different nitrogen sources, including uracil, uridine, dihydrouracil, and ammonia. The expression pattern of the known URC genes, URC1-6, helped to identify nine putative novel URC genes with a similar expression pattern. The microarray analysis provided evidence that both the URC and PYD genes are under nitrogen catabolite repression in L. kluyveri and are induced by uracil or dihydrouracil, respectively. We determined the function of URC8, which was found to catalyze the reduction of malonate semialdehyde to 3-hydroxypropionate, the final degradation product of the pathway. The other eight genes studied were all putative permeases. Our analysis of double deletion strains showed that the L. kluyveri Fui1p protein transported uridine, just like its homolog in Saccharomyces cerevisiae, but we demonstrated that is was not the only uridine transporter in L. kluyveri. We also showed that the L. kluyveri homologs of DUR3 and FUR4 do not have the same function that they have in S. cerevisiae, where they transport urea and uracil, respectively. In L. kluyveri, both of these deletion strains grew normally on uracil and urea.     

Pervasive and persistent redundancy among duplicated genes in yeast

The loss of functional redundancy is the key process in the evolution of duplicated genes. Here we systematically assess the extent of functional redundancy among a large set of duplicated genes in Saccharomyces cerevisiae. We quantify growth rate in rich medium for a large number of S. cerevisiae strains that carry single and double deletions of duplicated and singleton genes. We demonstrate that duplicated genes can maintain substantial redundancy for extensive periods of time following duplication (~100 million years). We find high levels of redundancy among genes duplicated both via the whole genome duplication and via smaller scale duplications. Further, we see no evidence that two duplicated genes together contribute to fitness in rich medium substantially beyond that of their ancestral progenitor gene. We argue that duplicate genes do not often evolve to behave like singleton genes even after very long periods of time.     

Additions, Losses, and Rearrangements on the Evolutionary Route from a Reconstructed Ancestor to the Modern Saccharomyces cerevisiae Genome

Comparative genomics can be used to infer the history of genomic rearrangements that occurred during the evolution of a species. We used the principle of parsimony, applied to aligned synteny blocks from 11 yeast species, to infer the gene content and gene order that existed in the genome of an extinct ancestral yeast about 100 Mya, immediately before it underwent whole-genome duplication (WGD). The reconstructed ancestral genome contains 4,703 ordered loci on eight chromosomes. The reconstruction is complete except for the subtelomeric regions. We then inferred the series of rearrangement steps that led from this ancestor to the current Saccharomyces cerevisiae genome; relative to the ancestral genome we observe 73 inversions, 66 reciprocal translocations, and five translocations involving telomeres. Some fragile chromosomal sites were reused as evolutionary breakpoints multiple times. We identified 124 genes that have been gained by S. cerevisiae in the time since the WGD, including one that is derived from a hAT family transposon, and 88 ancestral loci at which S. cerevisiae did not retain either of the gene copies that were formed by WGD. Sites of gene gain and evolutionary breakpoints both tend to be associated with tRNA genes and, to a lesser extent, with origins of replication. Many of the gained genes in S. cerevisiae have functions associated with ethanol production, growth in hypoxic environments, or the uptake of alternative nutrient sources.     

Alternative Splicing and Subfunctionalization Generates Functional Diversity in Fungal Proteomes

Alternative splicing is commonly used by the Metazoa to generate more than one protein from a gene. However, such diversification of the proteome by alternative splicing is much rarer in fungi. We describe here an ancient fungal alternative splicing event in which these two proteins are generated from a single alternatively spliced ancestral SKI7/HBS1 gene retained in many species in both the Ascomycota and Basidiomycota. While the ability to express two proteins from a single SKI7/HBS1 gene is conserved in many fungi, the exact mechanism by which they achieve this varies. The alternative splicing was lost in Saccharomyces cerevisiae following the whole-genome duplication event as these two genes subfunctionalized into the present functionally distinct HBS1 and SKI7 genes. When expressed in yeast, the single gene from Lachancea kluyveri generates two functionally distinct proteins. Expression of one of these proteins complements hbs1, but not ski7 mutations, while the other protein complements ski7, but not hbs1. This is the first known case of subfunctionalization by loss of alternative splicing in yeast. By coincidence, the ancestral alternatively spliced gene was also duplicated in Schizosaccharomyces pombe with subsequent subfunctionalization and loss of splicing. Similar subfunctionalization by loss of alternative splicing in fungi also explains the presence of two PTC7 genes in the budding yeast Tetrapisispora blattae, suggesting that this is a common mechanism to preserve duplicate alternatively spliced genes.     

Effects of cis and trans regulatory variations on the expression divergence of heat shock response genes between yeast strains

Phenotypic variation among individuals in a population can be due to DNA sequence variation in protein coding regions or in regulatory elements. Recently, many studies have indicated that mutations in regulatory elements may be the major cause of phenotypic evolution. However, the mechanisms for evolutionary changes in gene expression are still not well understood. Here, we studied the relative roles of cis and trans regulatory changes in Saccharomyces cerevisiae cells to cope with heat stress. It has been found that the expression level of ~ 300 genes was induced at least two fold and that of ~ 500 genes was repressed at least two fold in response to heat shock. From the former set of genes, we randomly selected 65 genes that showed polymorphism(s) between the BY and RM strains for pyrosequencing analysis to explore the relative contributions of cis and trans regulatory variations to the expression divergence between BY and RM. Our data indicated that the expression divergence between BY and RM was mainly due to trans regulatory variations under either the normal condition or the heat stress condition. However, the relative contribution of trans regulatory variation was decreased from 76.9% to 61.5% after the heat shock stress. These results indicated that the cis regulatory variation may play an important role in the adaption to heat stress. In our data, 43.1% (28 genes) of the 65 genes showed the same trend of cis or trans variation effect after the heat shock stress, 35.4% (23 genes) showed an increased cis variation effect and 21.5% (14 genes) showed an increased trans variation effect after the heat shock stress. Thus, our data give insights into the relative roles of cis and trans variations in response to heat shock in yeast.     

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.     

Paralogous ALT1 and ALT2 Retention and Diversification Have Generated Catalytically Active and Inactive Aminotransferases in Saccharomyces cerevisiae

Background

Gene duplication and the subsequent divergence of paralogous pairs play a central role in the evolution of novel gene functions. S. cerevisiae possesses two paralogous genes (ALT1/ALT2) which presumably encode alanine aminotransferases. It has been previously shown that Alt1 encodes an alanine aminotransferase, involved in alanine metabolism; however the physiological role of Alt2 is not known. Here we investigate whether ALT2 encodes an active alanine aminotransferase.

Principal Findings

Our results show that although ALT1 and ALT2 encode 65% identical proteins, only Alt1 displays alanine aminotransferase activity; in contrast ALT2 encodes a catalytically inert protein. ALT1 and ALT2 expression is modulated by Nrg1 and by the intracellular alanine pool. ALT1 is alanine-induced showing a regulatory profile of a gene encoding an enzyme involved in amino acid catabolism, in agreement with the fact that Alt1 is the sole pathway for alanine catabolism present in S. cerevisiae. Conversely, ALT2 expression is alanine-repressed, indicating a role in alanine biosynthesis, although the encoded-protein has no alanine aminotransferase enzymatic activity. In the ancestral-like yeast L. kluyveri, the alanine aminotransferase activity was higher in the presence of alanine than in the presence of ammonium, suggesting that as for ALT1, LkALT1 expression could be alanine-induced. ALT2 retention poses the questions of whether the encoded protein plays a particular function, and if this function was present in the ancestral gene. It could be hypotesized that ALT2 diverged after duplication, through neo-functionalization or that ALT2 function was present in the ancestral gene, with a yet undiscovered function.

Conclusions

ALT1 and ALT2 divergence has resulted in delegation of alanine aminotransferase activity to Alt1. These genes display opposed regulatory profiles: ALT1 is alanine-induced, while ALT2 is alanine repressed. Both genes are negatively regulated by the Nrg1 repressor. Presented results indicate that alanine could act as ALT2 Nrg1-co-repressor.     

A role for gene duplication and natural variation of gene expression in the evolution of metabolism

Background

Most eukaryotic genomes have undergone whole genome duplications during their evolutionary history. Recent studies have shown that the function of these duplicated genes can diverge from the ancestral gene via neo- or sub-functionalization within single genotypes. An additional possibility is that gene duplicates may also undergo partitioning of function among different genotypes of a species leading to genetic differentiation. Finally, the ability of gene duplicates to diverge may be limited by their biological function.

Methodology/Principal Findings

To test these hypotheses, I estimated the impact of gene duplication and metabolic function upon intraspecific gene expression variation of segmental and tandem duplicated genes within Arabidopsis thaliana. In all instances, the younger tandem duplicated genes showed higher intraspecific gene expression variation than the average Arabidopsis gene. Surprisingly, the older segmental duplicates also showed evidence of elevated intraspecific gene expression variation albeit typically lower than for the tandem duplicates. The specific biological function of the gene as defined by metabolic pathway also modulated the level of intraspecific gene expression variation. The major energy metabolism and biosynthetic pathways showed decreased variation, suggesting that they are constrained in their ability to accumulate gene expression variation. In contrast, a major herbivory defense pathway showed significantly elevated intraspecific variation suggesting that it may be under pressure to maintain and/or generate diversity in response to fluctuating insect herbivory pressures.

Conclusion

These data show that intraspecific variation in gene expression is facilitated by an interaction of gene duplication and biological activity. Further, this plays a role in controlling diversity of plant metabolism.