Search tips
Search criteria 


Logo of plosgenPLoS GeneticsSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
PLoS Genet. 2010 July; 6(7): e1001030.
Published online 2010 July 22. doi:  10.1371/journal.pgen.1001030
PMCID: PMC2908677

Tinkering Evolution of Post-Transcriptional RNA Regulons: Puf3p in Fungi as an Example

Jianzhi Zhang, Editor


Genome-wide studies of post-transcriptional mRNA regulation in model organisms indicate a “post-transcriptional RNA regulon” model, in which a set of functionally related genes is regulated by mRNA–binding RNAs or proteins. One well-studied post-transcriptional regulon by Puf3p functions in mitochondrial biogenesis in budding yeast. The evolution of the Puf3p regulon remains unclear because previous studies have shown functional divergence of Puf3p regulon targets among yeast, fruit fly, and humans. By analyzing evolutionary patterns of Puf3p and its targeted genes in forty-two sequenced fungi, we demonstrated that, although the Puf3p regulon is conserved among all of the studied fungi, the dedicated regulation of mitochondrial biogenesis by Puf3p emerged only in the Saccharomycotina clade. Moreover, the evolution of the Puf3p regulon was coupled with evolution of codon usage bias in down-regulating expression of genes that function in mitochondria in yeast species after genome duplication. Our results provide a scenario for how evolution like a tinker exploits pre-existing materials of a conserved post-transcriptional regulon to regulate gene expression for novel functional roles.

Author Summary

It is well known that the evolution of gene expression can account for significant phenotypic diversity among species. Gene expression is regulated at various levels. Many studies have demonstrated gene expression changes caused by mutations in transcription. Post-transcriptional regulation is also crucial for splicing, translation, localization, and degeneration of mRNAs in eukaryotes and, thus, is important in determining the abundance of gene expression. Changes in post-transcriptional regulons are also important in the evolution of gene expression. In this study we investigated the evolution of a particular post-transcriptional regulon, Puf3p, in fungal species. Our results illustrated an important evolutionary mode for evolution of post-transcriptional regulon, i.e., pre-existing materials of conserved post-transcriptional regulons can be recruited to regulate gene expression for novel functional roles.


Evolution of gene expression may account for significant phenotypic diversity among species [1][6]. Gene expression is regulated at various levels. Many studies have demonstrated gene expression changes caused by mutations in transcription [7], [8]. Post-transcriptional regulation is also crucial for splicing, translation, localization and degeneration of mRNAs in eukaryotes, and thus is important in determining the abundance of gene expression [9], [10]. It is likely that changes in post-transcriptional regulons are also important in the evolution of gene expression [11].

The control of post-transcriptional regulation is mediated by regulons such as mRNA-binding proteins (RBP) or RNAs (e.g., microRNAs) which usually bind to elements in the 3′ untranslated regions (UTR) and determine the fate of their targeted mRNAs. Evolution of microRNA post-transcriptional regulons has been well studied due to recent improvement in understanding their functions. It was shown that novel microRNAs can turn over rapidly during evolution [12], and for those that are highly conserved over long evolutionary distances, their targets can change dramatically even within populations [13]. Studies on microRNAs have revealed interesting information on evolution of this particular type of post-transcriptional regulon, whereas evolution of RBP regulons remains poorly understood. Furthermore, RBP regulons play major post-transcriptional roles in budding yeast because the species lost the microRNA regulatory machine [14].

One of the best-characterized RBP families is PUF (Pumilio and FBF, FBF represents for fem-3 binding factor), which is conserved in a wide variety of eukaryotes from yeast to humans [15][19]. The PUF post-transcriptional regulon regulates diverse gene sets in various model organisms. For example, in the budding yeast, Saccharomyces cerevisiae, genes most commonly targeted by Puf3p are in the mitochondria and play essential roles in mitochondrial biogenesis [17], [20]. In the fruit fly, Drosophila melanogaster, Pumilio (a PUF protein), which binds with the same element as Puf3p in the budding yeast, is necessary for early embryogenesis and development of primordial germ cells [21], . Genome-wide identification of the Pumilio targets in fruit flies uncovered genes involved particularly in nucleotide metabolism, transcriptional regulation and synthesis of membrane proteins [18]. In humans, two paralogous PUF proteins (Pum1p and Pum2p), which interact with the microRNA system in post-transcriptional regulation, share the same binding-element with yeast Puf3p and bind to mRNAs from genes that function in transcriptional regulation and cell proliferation [19], [23].

Previous studies have reported that the binding site of Puf3p is highly conserved in sensu stricto yeasts [24][26]. Taking advantage of a large number of genomic sequences, in this study we investigated the evolution of the Puf3p post-transcriptional regulon in fungi. Our results show continuous steps of functional innovation in the Puf3p regulon despite its ultra conservation in these fungal species. First, the regulation of mitochondrial biogenesis by the Puf3p regulon originated in the Saccharomycotina subdivision; second, the Puf3p regulon was coupled with codon usage bias to modulate expression of genes that function in mitochondria in yeasts after whole genome duplication (WGD). Our work and reports from other labs show that mitochondria underwent significant functional changes during the origin of an efficient aerobic fermentation system in the yeast species that went through WGD [27][30]. This current report provides evidence suggesting that the Puf3p post-transcriptional regulon was involved in the evolution of this novel life history in yeasts.


The binding element of the Puf3p regulon is conserved in fungi

The RNA-binding domain of Puf3p, called the PUF homology domain (PUF-HD), consists of eight repeated peptide motifs [31][33]. In order to study the evolution of the RNA-binding domain, orthologous genes of PUF3 were identified from the forty-two sequenced fungal species (Table S1). Domain alignments in SMART [34], [35] indicate that almost all Puf3p orthologs contain the eight repeated motifs even though some repeats are not highly conserved (Figure 1A). Puf3p in Lodderomyces elongisporus and Rhizopus oryzae lost one repeat which may have resulted from insufficient genome sequencing or assembly because both orthologs are located at the end of the assembled contigs. This result demonstrates that the binding domain of Puf3p is conserved among all studied fungi.

Figure 1
Conservation of the Puf3p regulon in fungi.

The evolutionary trajectory of the Puf3p-binding element in its target genes was further investigated. Because the 8nt-core motif of P3E is conserved from yeast to human, we used the 8nt-core P3E profiling in the budding yeast as a reference to identify all possible puf3p targeted genes that contain at least one P3E at their 3′ downstream sequences (Figure 1B) [19]. Because GC content is very low in the P3E, the genomic GC content would inevitably affect the frequency of P3E in each species. Indeed, as shown in Figure 1C, at the genome level, the number of genes with P3E is negatively correlated with the genomic GC content among the studied species.

In order to exclude the impact of genomic GC content on our results, we generated 10,000 random sequences for each species based on the average GC content of 1,000 bp downstream sequences of all the annotated genes in this species. Using the occurrence of P3E in these random sequences as background in each species, we estimated the probability of having the observed P3E frequency in the 3′ downstream sequences of all the annotated genes in the same species. The probability was calculated in each sliding window of 50 bp in the 3′ downstream sequences of each species. As shown in Figure 1D, all studied species exhibited significant enrichment of P3E in the first several sliding windows in the 3′ downstream regions. As the 3′ UTR in yeast is usually shorter than 250 bp [36], our results indicate that the enrichment of P3E in the studied fungi results from P3E conservation in the 3′ UTR sequences. When we used the GC contents in the 250 bp regions after the stop codon or in each 50 bp sliding window to calculate the background P3E motif distribution, similar enrichment of P3E motif in the 3′UTR regions are still observed in most fungi species (Figure S1).

The Puf3p target genes are significantly enriched for mitochondrial functions only in the Saccharomycotina subdivision

Previous reports showed that Puf3p plays an essential role in mitochondrial biogenesis in S. cerevisiae [17], [37], [38]. This observation prompted us to investigate whether the functional profile of the Puf3p regulon is also conserved among fungi species. Accordingly, we identified all of orthologous genes between each studied fungal species and budding yeast. Genes in each species are classified into categories based on the sub-cellular localization of their orthologs in the budding yeast [39]. We then estimated the enrichment of genes with P3E in each localization category. As shown in Figure 2A, we discovered that all of the studied species in the Saccharomycotina subdivision show significant enrichment of P3E in genes that function in mitochondria. Other clades of the studied fungi did not have this pattern. Indeed, close to 50% of genes that function in mitochondria have P3E in the Saccharomycotina subdivision, which is significantly higher than that of species in other clades (Figure 2B, P -value = 6×10−22). A similar pattern was observed when a slightly different P3E motif profile [17] was used to calculate the motif frequency (Figure S2).

Figure 2
The evolution of the puf3p regulon and mitochondrial biogenesis.

Because regulation of the mitochondrial translational machine is essential for mitochondrial biogenesis [40], we further investigated conservation of P3E among the orthologous genes with this particular function in the Saccharomycotina subdivision. Our results showed that ~80% of genes with highly conserved P3E were involved in mitochondrial translation, whereas this number is only ~4% for genes with little P3E conservation, indicating that the Puf3p regulation of genes that are involved in the mitochondrial translational machine is highly conserved in these Saccharomycotina species (Figure 2C).

Relaxation of codon usages in Puf3p-regulated mitochondrial genes

Gene expression is regulated at multiple levels. Biased usage of preferred codons can result in enhanced accuracy and speed of protein synthesis in highly expressed genes [41], [42]. Previous studies reported that codon usage bias in mitochondrial genes is relaxed, possibly due to a relaxed function of the organelle with the origin of an efficient aerobic fermentation system in the fungal lineage with WGD [27], [28]. We predicted that Puf3p-regulated mitochondrial genes, due to their importance in mitochondrial biogenesis and functions, would experience more relaxation of codon usage bias than other mitochondrial genes in the post-WGD yeast species. In order to test this, we calculated the average codon bias adaptation index (CAI) for the mitochondrial genes, with and without P3E, for each species. As shown in Figure 3A, the mitochondrial genes with P3E in the post-WGD species show significantly smaller CAI than those genes in the fungal species that diverged from the common ancestor before the WGD event (student t-test, P = 3×10−8), whereas mitochondrial genes without P3E did not show such a pattern (Figure 3B, student t-test, P = 0.1).

Figure 3
The codon adaptation index (CAI) of mitochondrial genes in species in the clade Saccharomycotina.

Puf3p regulon in mitochondrial gene regulation

Due to the importance of Puf3p regulation in mitochondrial gene degradation, we further investigated expression of its target mitochondrial genes under the fermentative condition. Using gene-expression profiling measured by microarray data [43], we discovered that significantly more mitochondrial genes with P3E were down-regulated in the fermentative medium (YPD) than those without P3E (Figure 4A, Fisher's exact test, P = 1.2×10−4). Furthermore, we found that mitochondrial genes with P3E tend to be co-regulated because the average correlation coefficients of gene expression among mitochondrial genes with P3E in different conditions is significantly higher than that of genes without P3E (Figure 4B, student t test, P = 0). Therefore Puf3p regulon plays an important role in regulating mitochondrial genes in different conditions.

Figure 4
Regulation of mitochondrial genes by Puf3p regulon.


Function of the Puf3p regulon under various environmental conditions

PUF protein was first characterized in Drosophila as an mRNA-binding factor that recruits other proteins to inhibit the translation of the bound mRNA [21]. Subsequently, many studies in yeast revealed that the PUF family regulates specific mRNA degeneration by their RNA-binding domains [17], [31], [44], [45]. It was shown that the function of targeting mRNA for degeneration by Puf3p is much more efficient in 2% glucose (YPD, fermentative) than in 3% ethanol (YPE, non-fermentative) [17], [46]. Furthermore, it was shown that Puf3p is crucial for mitochondrial biogenesis and motility under non-fermentative conditions in budding yeast [43]. Saint-Georges and his colleagues reported that Puf3p can transfer its target mRNAs to the peripheral mitochondria in the non-fermentative growth medium [37]. The expression of PUF3 gene is significantly higher in yeast growing in YPE than in YPD (Figure 5A, P-value<0.05). We speculate that this is true because the positive regulation of mitochondrial biogenesis might not be as important for Puf3p in fermentative conditions as that in respiratory conditions: First, based on gene deletion data, the mitochondrial genes with P3E are significantly more important (having more severe growth defects after gene deletion) than those genes without P3E (P = 4×10−6) under non-fermentative conditions, but these two gene groups do not show obvious difference in deletion phenotype under fermentative conditions (Figure 5B). Second, severe growth defect after PUF3 gene deletion was observed in YPE, but not in YPD (Figure 5C). Therefore in non-fermentative condition Puf3p regulates both mitochondrial biogenesis and mRNA degradation, but in fermentative condition, it might only regulate mRNA degradation, albeit more efficiently in this condition. The expression difference of PUF3 in two growth conditions can also be explained by the fact that mitochondrial biogenesis in non-fermentative conditions is extremely important for yeast because the organism relies on respiration, and therefore mitochondria, to generate cellular energy in these conditions. In contrast, the function of mRNA degradation might not be as essential to the organism under fermentative conditions.

Figure 5
Function of the Puf3p regulon under different conditions.

Although loss of a functional PUF3 gene shows a negligible effect on organism growth in YPD (Figure 5C), our results in this study indicate that Puf3p regulation of mitochondrial gene degradation might be important for yeast fermentative growth during evolution. After whole genome duplication, the post-WGD species (including budding yeast) evolved efficient fermentative ability [47]. Mitochondrial function became relaxed in these species [27], [28]. Most post-WGD yeast species can live even without a functional mitochondrial genome [47]. A large number of mitochondrial genes are down-regulated in yeast fermentative growth [48]. Degradation of mRNA by Puf3p may accelerate this gene-expression regulation process during environmental switches. Interestingly, our results showed that mitochondrial genes having P3E had significantly relaxed codon usage bias in the post-WGD species, which is not true for other mitochondrial genes.

Evolution of the Puf3p post-transcriptional regulon and reuse of extant modules for novel functional roles

Understanding the origin of genetic novelties is a challenging issue in evolutionary biology. One of the prominent models proposed by Francois Jacob for evolution of genetic novelties in gene regulatory network is tinkering evolution, in which evolution reorganizes pre-existing networks to produce novelties [6], [49][52]. Our results provide an interesting paradigm for the evolution of post-transcriptional regulons that is consistent with this model. The Puf3p post-transcriptional regulon might have changed significantly at least twice in fungal evolution. First, although it is conserved in all the studied fungal species, the dedicated function of the Puf3p post-transcriptional regulon in mitochondrial biogenesis independently evolved in the Saccharomycotina subdivision. Second, although the regulation of mitochondrial mRNAs by Puf3p is conserved across the Saccharomycotina subdivision, the Puf3p target genes evolved a reduced codon usage bias in the post-WGD species, which might be consistent with the functional relaxation of mitochondrial genes in the post-WGD species due to the emergence of their fermentative life-style during evolution.

It was shown that although some microRNAs are highly conserved, their target networks can change dramatically during evolution [13]. Our results, together with those from previous studies, indicate that post-transcriptional regulons by RBP and microRNAs might share similar evolutionary patterns: i.e., the interaction mechanisms between the regulators and their target genes are conserved, whereas the target network is plastic during evolution. As post-transcriptional regulation plays an important role in regulating gene expression, this evolutionary scenario involving post-transcriptional regulons could lead to significant gene-expression divergence among species.

Materials and Methods

Sequence data

Sequences for the forty-two sequenced fungal species were downloaded from the Fungal Comparative Genomics database [53] and the National Center for Biotechnology Information (

Orthologous gene definition and motif detection

Using the InParanoid software package [54], orthologs between budding yeast and other fungal species were identified. The PUF3 orthologs in other species of fungi were identified manually based on the best alignment. Eight repeated motifs of PUF protein were detected by the SMART sequence analysis (

Because the Puf3p binding motif in its target genes (P3E) is conserved between yeast and humans, based on the profile of P3E in budding yeast, we used a Perl script to detect the target locus of Puf3p by fixing all seven invariable sites and allowing flexibility in the fifth site (Figure 1B). For each species, we scanned 1,000 bp of DNA sequence downstream of all annotated genes to discern the occurrence of the P3E motif. The percentage of motif occurrences in each 50-bp window among all genes in each species was calculated. Multiple occurrences of the motif in the same sliding window were regarded as independent events. To see whether the motif occurrences in a species are different from random expectation, we calculated GC content of the 1,000bp downstream sequences for all genes in each species. Based on the observed GC content, 10,000 random sequences of 1,000bp were generated by a perl script and the occurrences of PUF3 motif were calculated. The significance of PUF3 motif occurrences in each sliding window was calculated by comparing its motif frequency against the frequencies of P3E in these random sequences by Fisher's exact test (more focused tests were also conducted in Figure S1). The Bonferroni correction was used to correct for multiple comparisons. Genes with P3E in their 250bp 3′ downstream regions were defined as the target of Puf3p.

Functional enrichment analysis

The sub-localization information for genes in budding yeast was downloaded from the Saccharomyces Genome Database ( [39]. Based on the identified orthologous genes between each studied fungal species and budding yeast, genes in each species are classified into different localization categories based on the sub-cellular localization of their orthologs in the budding yeast. The hypergeometric test was used to test the enrichment of genes with P3E in each localization category. The Bonferroni correction was used to correct for multiple comparisons.

Codon usage bias calculation

According to the InParanoid results, we identified all the orthologous clusters that contained the known cytoplasmic ribosomal protein genes in S. cerevisiae, regardless of gene copy number in each species. The ribosomal protein genes in each species were used as references to calculate the codon adaptation index (CAI) value for each individual gene in the same species by CodonW ( [55]. In order to compare codon usage bias between different species, CAI values in each species were standardized so that the mean and standard deviation were 0 and 1, respectively.

Microarray data analysis

Gene expression in fermentative conditions

Microarray data were downloaded from ( The expression profiles of mitochondrial genes during cell growth in fermentative (YPD) media were used to do the analysis. We counted the number of the up-regulated and down-regulated genes (defined as their expression level was higher or lower by 1.5 folds than the reference pool, respectively). Fisher's exact test was used to assess whether the number of up or down-regulated genes were significantly different between mitochondrial genes with P3E and without P3E.

Gene co-regulation analysis

The collected microarray data under different experimental conditions were downloaded from ( [29]. The dataset compiles microarray gene expression under multiple conditions and contains 1,011 data points for each gene. Correlation coefficients of gene expression among the studied mitochondrial genes were calculated in R.

Quantitative real-time PCR

Cells were grown in the YPD and YPE media until optical densities at 600 nm reaches 1. Total RNA was extracted using the Trizol protocol [56] and cDNA was synthesized using an Invitrogen kit (Cat. No.18080-051). Using the ACT1 gene as reference, the expressional levels of PUF3 in fermentative and non-fermentative conditions were measured by quantitative real-time PCR.

Gene deletion analysis and growth measurement

The rate of growth for deletion mutants were downloaded from [57]. Two-tails student t test was used to compare the average fitness contribution of mitochondrial genes with and without P3E in fermentative and non-fermentative growth conditions.

The deletion of the PUF3 gene was conducted in the BY4741 strain (Mat a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0) by homologous recombination and ura- was used as the selection marker (the primer sequences for gene deletion are available upon request). The mutant and the wild type strains were grown overnight in YPD (2% glucose) and YPE (3% ethanol) media. Cells were then transferred into fresh media and grown until the optical density (600 nm) reached 0.2. Then 4ul of growth media were dotted onto the YPD and YPE plates with ten-fold dilutions. YPD and YPE plates were incubated at 30°C for 48 and 72h, respectively.

Supporting Information

Figure S1

The occurrences of the P3E at the downstream of all annotated genes in each species. (A) The randomly generated sequences based on the GC contents of 250 bp regions after the stop codon were used as background to calculate motif enrichment; (B) the randomly generated sequences based on the GC contents in each sliding window were used as background to calculate motif enrichment. Red denotes significant enrichment of the observed P3E motif in the real sequences in comparison to the random sequences.

(0.21 MB TIF)

Figure S2

The occurrence of P3E for genes functioning in mitochondria in yeast species. The PUF3 motif profile from Gerber et al, (2004) [17], which has two extra nucleotides at the 5′ of the motif in the paper, was used. As shown in the figure, species in Saccharomycotina subdivision have significantly higher percentages of mitochondrial genes that have this motif than the other species. The hypergeometric test was used to test the enrichment of mitochondrial genes with P3E in each species. The P value is smaller than 1×10−8 for each species in the Saccharomycotina subdivision, but larger than 0.2 in each of the other species.

(0.74 MB TIF)

Table S1

The forty-two sequenced fungal species used in this study.

(0.09 MB DOC)


We thank Jian Lu, Lin Xu, and Hong Chen for discussions and Haijing Yu, Jun Sun, Roberto Flores, and Allan Eaglesham for reading the manuscript.


The authors have declared that no competing interests exist.

This work was supported by startup funds from Cornell University and NSF DEB-0949556 awarded to ZG. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Lynch VJ, Wagner GP. Resurrecting the role of transcription factor change in developmental evolution. Evolution. 2008;62:2131–2154. [PubMed]
2. Khaitovich P, Enard W, Lachmann M, Paabo S. Evolution of primate gene expression. Nat Rev Genet. 2006;7:693–702. [PubMed]
3. Sung HM, Wang TY, Wang D, Huang YS, Wu JP, et al. Roles of Trans and Cis Variation in Yeast Intra-species Evolution of Gene Expression. Mol Biol Evol 2009 [PMC free article] [PubMed]
4. Tuch BB, Li H, Johnson AD. Evolution of eukaryotic transcription circuits. Science. 2008;319:1797–1799. [PubMed]
5. King MC, Wilson AC. Evolution at two levels in humans and chimpanzees. Science. 1975;188:107–116. [PubMed]
6. Wray GA, Hahn MW, Abouheif E, Balhoff JP, Pizer M, et al. The evolution of transcriptional regulation in eukaryotes. Mol Biol Evol. 2003;20:1377–1419. [PubMed]
7. Wagner GP, Lynch VJ. The gene regulatory logic of transcription factor evolution. Trends Ecol Evol. 2008;23:377–385. [PubMed]
8. Wray GA. The evolutionary significance of cis-regulatory mutations. Nat Rev Genet. 2007;8:206–216. [PubMed]
9. Mata J, Marguerat S, Bahler J. Post-transcriptional control of gene expression: a genome-wide perspective. Trends Biochem Sci. 2005;30:506–514. [PubMed]
10. Moore MJ. From birth to death: the complex lives of eukaryotic mRNAs. Science. 2005;309:1514–1518. [PubMed]
11. Keene JD. RNA regulons: coordination of post-transcriptional events. Nat Rev Genet. 2007;8:533–543. [PubMed]
12. Lu J, Shen Y, Wu Q, Kumar S, He B, et al. The birth and death of microRNA genes in Drosophila. Nat Genet. 2008;40:351–355. [PubMed]
13. Chen K, Rajewsky N. Natural selection on human microRNA binding sites inferred from SNP data. Nat Genet. 2006;38:1452–1456. [PubMed]
14. Drinnenberg IA, Weinberg DE, Xie KT, Mower JP, Wolfe KH, et al. RNAi in budding yeast. Science. 2009;326:544–550. [PMC free article] [PubMed]
15. Wickens M, Bernstein DS, Kimble J, Parker R. A PUF family portrait: 3′UTR regulation as a way of life. Trends Genet. 2002;18:150–157. [PubMed]
16. Spassov DS, Jurecic R. The PUF family of RNA-binding proteins: does evolutionarily conserved structure equal conserved function? IUBMB Life. 2003;55:359–366. [PubMed]
17. Gerber AP, Herschlag D, Brown PO. Extensive association of functionally and cytotopically related mRNAs with Puf family RNA-binding proteins in yeast. PLoS Biol. 2004;2:e79. doi: 10.1371/journal.pbio.0020079. [PMC free article] [PubMed]
18. Gerber AP, Luschnig S, Krasnow MA, Brown PO, Herschlag D. Genome-wide identification of mRNAs associated with the translational regulator PUMILIO in Drosophila melanogaster. Proc Natl Acad Sci U S A. 2006;103:4487–4492. [PubMed]
19. Galgano A, Forrer M, Jaskiewicz L, Kanitz A, Zavolan M, et al. Comparative analysis of mRNA targets for human PUF-family proteins suggests extensive interaction with the miRNA regulatory system. PLoS One. 2008;3:e3164. doi: 10.1371/journal.pone.0003164. [PMC free article] [PubMed]
20. Tu BP, Kudlicki A, Rowicka M, McKnight SL. Logic of the yeast metabolic cycle: temporal compartmentalization of cellular processes. Science. 2005;310:1152–1158. [PubMed]
21. Murata Y, Wharton RP. Binding of pumilio to maternal hunchback mRNA is required for posterior patterning in Drosophila embryos. Cell. 1995;80:747–756. [PubMed]
22. Lin H, Spradling AC. A novel group of pumilio mutations affects the asymmetric division of germline stem cells in the Drosophila ovary. Development. 1997;124:2463–2476. [PubMed]
23. Morris AR, Mukherjee N, Keene JD. Ribonomic analysis of human Pum1 reveals cis-trans conservation across species despite evolution of diverse mRNA target sets. Mol Cell Biol. 2008;28:4093–4103. [PMC free article] [PubMed]
24. Jacobs Anderson JS, Parker R. Computational identification of cis-acting elements affecting post-transcriptional control of gene expression in Saccharomyces cerevisiae. Nucleic Acids Res. 2000;28:1604–1617. [PMC free article] [PubMed]
25. Shalgi R, Lapidot M, Shamir R, Pilpel Y. A catalog of stability-associated sequence elements in 3′ UTRs of yeast mRNAs. Genome Biol. 2005;6:R86. [PMC free article] [PubMed]
26. Kellis M, Birren BW, Lander ES. Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae. Nature. 2004;428:617–624. [PubMed]
27. Jiang H, Guan W, Pinney D, Wang W, Gu Z. Relaxation of yeast mitochondrial functions after whole-genome duplication. Genome Res. 2008;18:1466–1471. [PubMed]
28. Man O, Pilpel Y. Differential translation efficiency of orthologous genes is involved in phenotypic divergence of yeast species. Nat Genet. 2007;39:415–421. [PubMed]
29. Ihmels J, Bergmann S, Gerami-Nejad M, Yanai I, McClellan M, et al. Rewiring of the yeast transcriptional network through the evolution of motif usage. Science. 2005;309:938–940. [PubMed]
30. Field Y, Fondufe-Mittendorf Y, Moore IK, Mieczkowski P, Kaplan N, et al. Gene expression divergence in yeast is coupled to evolution of DNA-encoded nucleosome organization. Nat Genet. 2009;41:438–445. [PMC free article] [PubMed]
31. Olivas W, Parker R. The Puf3 protein is a transcript-specific regulator of mRNA degradation in yeast. Embo J. 2000;19:6602–6611. [PubMed]
32. Wang X, McLachlan J, Zamore PD, Hall TM. Modular recognition of RNA by a human pumilio-homology domain. Cell. 2002;110:501–512. [PubMed]
33. Miller MT, Higgin JJ, Hall TM. Basis of altered RNA-binding specificity by PUF proteins revealed by crystal structures of yeast Puf4p. Nat Struct Mol Biol. 2008;15:397–402. [PMC free article] [PubMed]
34. Schultz J, Milpetz F, Bork P, Ponting CP. SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A. 1998;95:5857–5864. [PubMed]
35. Letunic I, Doerks T, Bork P. SMART 6: recent updates and new developments. Nucleic Acids Res. 2009;37:D229–232. [PMC free article] [PubMed]
36. Mignone F, Gissi C, Liuni S, Pesole G. Untranslated regions of mRNAs. Genome Biol. 2002;3:REVIEWS0004. [PMC free article] [PubMed]
37. Saint-Georges Y, Garcia M, Delaveau T, Jourdren L, Le Crom S, et al. Yeast mitochondrial biogenesis: a role for the PUF RNA-binding protein Puf3p in mRNA localization. PLoS One. 2008;3:e2293. doi: 10.1371/journal.pone.0002293. [PMC free article] [PubMed]
38. Garcia-Rodriguez LJ, Gay AC, Pon LA. Puf3p, a Pumilio family RNA binding protein, localizes to mitochondria and regulates mitochondrial biogenesis and motility in budding yeast. J Cell Biol. 2007;176:197–207. [PMC free article] [PubMed]
39. Huh WK, Falvo JV, Gerke LC, Carroll AS, Howson RW, et al. Global analysis of protein localization in budding yeast. Nature. 2003;425:686–691. [PubMed]
40. Jacobs HT, Turnbull DM. Nuclear genes and mitochondrial translation: a new class of genetic disease. Trends Genet. 2005;21:312–314. [PubMed]
41. Sharp PM, Li WH. The codon Adaptation Index–a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987;15:1281–1295. [PMC free article] [PubMed]
42. Akashi H. Translational selection and yeast proteome evolution. Genetics. 2003;164:1291–1303. [PubMed]
43. Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, et al. Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000;11:4241–4257. [PMC free article] [PubMed]
44. Tadauchi T, Matsumoto K, Herskowitz I, Irie K. Post-transcriptional regulation through the HO 3′-UTR by Mpt5, a yeast homolog of Pumilio and FBF. Embo J. 2001;20:552–561. [PubMed]
45. Jackson JS, Jr, Houshmandi SS, Lopez Leban F, Olivas WM. Recruitment of the Puf3 protein to its mRNA target for regulation of mRNA decay in yeast. Rna. 2004;10:1625–1636. [PubMed]
46. Foat BC, Houshmandi SS, Olivas WM, Bussemaker HJ. Profiling condition-specific, genome-wide regulation of mRNA stability in yeast. Proc Natl Acad Sci U S A. 2005;102:17675–17680. [PubMed]
47. Merico A, Sulo P, Piskur J, Compagno C. Fermentative lifestyle in yeasts belonging to the Saccharomyces complex. Febs J. 2007;274:976–989. [PubMed]
48. Santangelo GM. Glucose signaling in Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 2006;70:253–282. [PMC free article] [PubMed]
49. Jacob F. Evolution and tinkering. Science. 1977;196:1161–1166. [PubMed]
50. Carroll SB. Evolution at two levels: on genes and form. PLoS Biol. 2005;3:e245. doi: 10.1371/journal.pbio.0030245. [PMC free article] [PubMed]
51. Jacob F. Complexity and tinkering. Ann N Y Acad Sci. 2001;929:71–73. [PubMed]
52. Long M. Evolution of novel genes. Curr Opin Genet Dev. 2001;11:673–680. [PubMed]
53. Fitzpatrick DA, Logue ME, Stajich JE, Butler G. A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis. BMC Evol Biol. 2006;6:99. [PMC free article] [PubMed]
54. Remm M, Storm CE, Sonnhammer EL. Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001;314:1041–1052. [PubMed]
55. Sharp PM, Matassi G. Codon usage and genome evolution. Curr Opin Genet Dev. 1994;4:851–860. [PubMed]
56. Chomczynski P, Mackey K. Short technical reports. Modification of the TRI reagent procedure for isolation of RNA from polysaccharide- and proteoglycan-rich sources. Biotechniques. 1995;19:942–945. [PubMed]
57. Steinmetz LM, Scharfe C, Deutschbauer AM, Mokranjac D, Herman ZS, et al. Systematic screen for human disease genes in yeast. Nat Genet. 2002;31:400–404. [PubMed]
58. Bailey TL, Elkan C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994;2:28–36. [PubMed]
59. Workman CT, Yin Y, Corcoran DL, Ideker T, Stormo GD, et al. enoLOGOS: a versatile web tool for energy normalized sequence logos. Nucleic Acids Res. 2005;33:W389–392. [PMC free article] [PubMed]

Articles from PLoS Genetics are provided here courtesy of Public Library of Science