Research article | Open | Published:
Identification of novel conserved peptide uORF homology groups in Arabidopsis and rice reveals ancient eukaryotic origin of select groups and preferential association with transcription factor-encoding genes
BMC Biologyvolume 5, Article number: 32 (2007)
Upstream open reading frames (uORFs) can mediate translational control over the largest, or major ORF (mORF) in response to starvation, polyamine concentrations, and sucrose concentrations. One plant uORF with conserved peptide sequences has been shown to exert this control in an amino acid sequence-dependent manner but generally it is not clear what kinds of genes are regulated, or how extensively this mechanism is invoked in a given genome.
By comparing full-length cDNA sequences from Arabidopsis and rice we identified 26 distinct homology groups of conserved peptide uORFs, only three of which have been reported previously. Pairwise K a /K s analysis showed that purifying selection had acted on nearly all conserved peptide uORFs and their associated mORFs. Functions of predicted mORF proteins could be inferred for 16 homology groups and many of these proteins appear to have a regulatory function, including 6 transcription factors, 5 signal transduction factors, 3 developmental signal molecules, a homolog of translation initiation factor eIF5, and a RING finger protein. Transcription factors are clearly overrepresented in this data set when compared to the frequency calculated for the entire genome (p = 1.2 × 10-7). Duplicate gene pairs arising from a whole genome duplication (ohnologs) with a conserved uORF are much more likely to have been retained in Arabidopsis (Arabidopsis thaliana) than are ohnologs of other genes (39% vs 14% of ancestral genes, p = 5 × 10-3). Two uORF groups were found in animals, indicating an ancient origin of these putative regulatory elements.
Conservation of uORF amino acid sequence, association with homologous mORFs over long evolutionary time periods, preferential retention after whole genome duplications, and preferential association with mORFs coding for transcription factors suggest that the conserved peptide uORFs identified in this study are strong candidates for translational controllers of regulatory genes.
Upstream open reading frames (uORFs) are small open reading frames found in the 5' UTR of a mature mRNA, and can mediate translational regulation of the largest, or major, ORF (mORF). Regulation by uORFs has been studied in several individual transcripts demonstrating the importance of uORFs in such processes as polyamine production , amino acid production [2, 3], and sucrose response , but the biological effect of uORFs in the vaste majority of transcripts of the genome is still unclear. Upstream start codons (uAUGs) occur in 20–30% of yeast, mammalian, and plant transcript 5' UTRs [5–7] therefore potentially thousands of genes are regulated in this manner.
The majority of characterized uORFs appear to act in an amino acid sequence-independent manner, regulating mORF translation by the uORF start codon nucleotide context, by the uORF length, or by the distance between the uORF stop codon and the mORF start codon, rather than by uORF-encoded peptides [8–11]. Some uORFs, however, do rely on peptide sequences to mediate translational regulation of the associated mORF, but few examples have been identified and characterized to date. In fungi and animals, a few genes have been shown to contain uORFs whose amino acid sequences are similar between two or more species [12–17], but only two cases, CPA1  and SAMDC1/AdoMetDC1 , have demonstrated uORF sequence-dependent regulation. In plants two groups of genes, S-Adenosylmethionine decarboxylases (AdoMetDCs; EC 188.8.131.52) and group S basic region leucine zipper (bZIP) transcription factors, have been shown to contain uORFs with similar amino acids between monocots and dicots [19, 20]. In the former group, mORF translational regulation is dependent on the sequence of the uORF peptide [1, 4] and overexpression of the mORF in either group results in stunted or lethal phenotypes, suggesting that these genes play a critical role in growth and/or development. Indeed, AdoMetDC is required for polyamine synthesis, molecules that are implicated in essential plant functions such as cell division, embryogenesis, leaf, root, and flower development, and stress responses [21, 22].
In general, it has been difficult to carry out genome-wide surveys of conserved peptide uORFs due to poor annotation of 5' UTRs. The availability of expressed sequence tags (ESTs) has improved exon and intron annotation of the genomic sequence, but they are relatively short and often do not predict the entire mRNA molecule, even when several ESTs overlap the same genomic region and can be assembled to predict one transcript. As there are very few introns in yeast transcripts, prediction of uORF conservation has been attempted in S. cerevisiae by analyzing genomic sequence upstream of predicted mORF start sites , but it is still not clear whether these uORFs are truly conserved (i.e., are under negative selection pressures), or are simply undergoing evolutionary drift. With the sequencing of the Aspergillus nidulans genome, comparison to A. fumigatus and A. oryzae has identified 38 uORFs with putatively conserved start and stop codon positions relative to the mORF, 14 of which are conserved in one of Neurospora crassa, Fusarium graminaerum, or Magnaporthe grisea , but the authors did not comment on whether the uORF amino acid sequences are also conserved.
With the emergence of large plant full-length cDNA sequence collections [24–26], it is now possible to adopt a comparative genomics approach to determine the prevalence of conserved amino acid uORFs in the genome and the persistence of these elements throughout eukaryotic evolution. Because rice and Arabidopsis shared a common ancestor 140–200 million years ago (Mya) [27–29], sequence similarity retained over this amount of time provides good candidates for truly conserved peptide uORF sequences. In this study we have used Oryza sativa (rice) and Arabidopsis thaliana (Arabidopsis) full-length cDNA sequence collections to estimate the incidence of conserved peptide uORFs in the rice and Arabidopsis genomes, to determine the prevalence of uORFs within regulatory genes, and to compare evolutionary rates for uORFs versus mORFs. By examining more distantly related sequences, we posit an ancient origin for select uORFs and we provide evidence for one mechanism by which uORFs can arise within genes.
Identification of conserved peptide uORFs by comparison of rice and Arabidopsis transcripts
To identify conserved peptide uORFs, we developed "uORF-Finder", a Perl program that first compares the mORF amino acid sequence of each cDNA from one collection with the mORF sequences of another species' collection to identify putative mORF homologs, and then compares the uORFs in the 5' UTRs of the two paired sequences to identify uORFs with conserved amino acid sequences (see Methods). Comparison by uORF-Finder of a corrected set of 34000 full-length cDNA sequences from Arabidopsis with a similar set from rice resulted in the identification of conserved peptide uORFs in 44 Arabidopsis genes and 36 rice genes, which together comprise 19 homology groups based on uORF amino acid similarity (Tables 1, 2, 3; Figures 1, 2, 3, 4, 5). All three of the homology groups that had been previously reported were identified by uORF-Finder [1, 4]. The other 16 conserved uORFs have not been reported previously. Homologs of these 19 conserved uORF groups also exist in other angiosperm species (Figures 1, 2, 3, 4, 5).
Comparison of Arabidopsis homologs detects additional conserved uORFs
Conserved uORFs that are not sufficiently well conserved to be detected in a rice-Arabidopsis comparison could conceivably be detected in ohnologs, homologous genes arising by whole-genome duplication (WGD) , and in paralogs, homologous genes arising from segmental duplication or tandem duplication. Modification of uORF-Finder allowed comparison of each full-length cDNA to all other cDNAs in the same collection (see Methods), and identified seven additional conserved uORF homology groups (Tables 1 and 4; Figures 6, 7, 8). Six of these pairs are ohnologs, created by the most recent WGD (24–40 Mya) in an ancestor of Arabidopsis [31–33]. The seventh pair is not found in syntenic regions and is most likely a paralogous pair. It appears to have arisen at about the same time as the recent WGD event because its synonymous substitution frequency (K s value) of 0.7 is similar to the median K s of recent duplicate pairs (0.8) and is within their K s range (0.4–1.6) . The corresponding rice genes in four of the seven homology groups possess uORFs, but lack sufficient uORF sequence similarity to have been detected in the Arabidopsis-rice comparison (Figures 6, 7, 8).
Purifying selection maintains uORF amino acid sequences
Pairwise K a /K s tests for selection on amino acid sequences were applied to each uORF homology group and their associated mORFs to determine whether uORF amino acid sequences are under selective constraints similar to their associated mORFs. Both an approximate method (Yn00) and a maximum likelihood method (codeml) were used to calculate mean pairwise K a /K s ratios for each group. A K a /K s ratio less than 1 implies that negative, or purifying, selection has acted on the sequence, a ratio equal to 1 suggests drift, and a ratio greater than 1 indicates that positive selection has acted on an amino acid sequence. It is also true that conservation at the nucleotide level, not the amino acid level, can drive the K a /K s ratio to one. Analysis of all 26 homology groups showed that generally both uORFs and mORFs have been under mild to strong purifying selection since the divergence of each gene pair (Table 5) and these low K a /K s ratios suggest that the conservation is at the amino acid level, not simply at the nucleotide level.
One possible explanation for low K a /K s ratios in the putative uORFs invokes an incomplete splicing of the full-length cDNAs for which the uORF and mORF are normally fused. To address this possibility, all Genbank Arabidopsis ESTs were screened for evidence of uORF-mORF translational fusions. No ORFs were found to run continuously between the uORF and mORF, with one exception. A fusion product (Genbank accession no. DR353698) was identified between the N-terminal and central region of the uORF and the central and C-terminal region of the mORF found at locus At5g03190 (group 17). Classification of this putative uORF is shown in Table 1 for two reasons. Firstly, the four uORF C-terminal amino acids that are excluded in the fusion EST are perfectly conserved in monocot and dicot members, and the position of their stop codon is perfectly conserved, therefore it is difficult to explain this conservation if the uORF is not translated. Secondly, the N-terminal portion of the mORF that is removed in the fusion EST is similar between three Arabidopsis loci of the same homology group, with the start codon position also being conserved in these three members. It is likely, therefore, that the fusion EST represents an alternatively spliced form of this transcript, but further characterization of this locus will be needed to support this conclusion. Most of the homology groups show uORFs with conserved amino acid residues at the C-terminus and an identical positioning of the uORF stop codon (Figures 1, 2, 3, 4, 5, 6, 7, 8). This would suggest that the full-length cDNAs are fully spliced and are not erroneously predicting uORF sequences due to incomplete splicing.
Conserved features of uORF sequences
The lengths of uORFs vary to differing degrees within and among homology groups, but in amino acid sequence alignments nearly all groups exhibit considerable conservation of the position of the N-terminus and/or the C-terminus, i.e., length variation is usually due to a variable region in the middle or at one end of the uORF (Table 6; Figures 1, 2, 3, 4, 5, 6, 7, 8). The amino acid sequences of some uORFs possess potentially interesting features. Notably, some uORF groups possess regions rich in serine, threonine, and/or tyrosine, and others possess regions rich in lysine and/or arginine. Two homology groups are particularly noteworthy: Group 8 uORFs specify peptides with a coiled coil-helix, coiled coil-helix (CHCH) domain (Pfam accession number PF06747; Figure 9), and group 13 uORFs encode peptides that are extremely serine/arginine-rich (Figure 10). Both of these unusual peptides will be discussed in further detail below.
Most genes with conserved uORFs appear to have regulatory functions
A total of 31% of mORFs encoded by conserved peptide uORF loci in Arabidopsis were predicted to be a transcription factor, as determined by GO molecular function terms (Tables 2 and 4), whereas only 5.9% of all Arabidopsis loci are predicted to encode transcription factors . Thus, genes predicted to encode transcription factors are significantly overrepresented (p = 1.2 × 10-7) among conserved peptide uORF loci. In each case, GO terms were validated by manual annotation of protein functions using domain predictions from NCBI Conserved Domain and InterProScan Database searches [35, 36]. A variety of different types of transcription factors, including bZIP, Ethylene Response Factor/Apetala 2-like (ERF/AP2-like), basic helix-loop-helix (bHLH), and homeobox proteins, are represented among conserved peptide uORF loci with no demonstrable bias. No other GO terms were found to be significantly over- or under-represented in the uORF data set.
Biological functions could be inferred for 16 of the 26 uORF homology groups (Table 1). Six groups encode transcription factor homologs and so are presumably involved in transcriptional control (1, 2, 14, 15, 18, and 21). Five groups are likely to be involved in signal transduction, including four protein kinases and a putative calmodulin-binding protein involved in auxin response (groups 10, 16, 19, 23, 25). Two groups are involved in the metabolism of small molecules that regulate plant development: polyamines (group 3)  and trehalose (group 11) . One group (13) encodes the key enzyme in the biosynthesis of phosphocholine, which is an intermediate in biosynthesis of phosphatidylcholine and phosphatidic acid; phosphocholine levels influence levels of phosphatidic acid, an important physiological and developmental signal molecule [38–40]. Group 7 putatively encodes translation initiation factor eIF5, which influences start codon selection, and Group 26 encodes a RING finger protein, suggesting a role in targeted protein turnover by ubiquitination. Of the remaining 10 groups, 8 encode predicted proteins of unknown function, 1 encodes an ankyrin-repeat protein, and 1 encodes an amine oxidase. Thus, all but two families of conserved uORF genes whose functions are known or can be inferred potentially play a regulatory role in the biology of plants.
Genes with conserved uORFs were preferentially retained after whole genome duplication
Since the most recent WGD event in the Arabidopsis lineage, only 14% of the original gene pairs present in the ancestral tetraploid have been retained as a duplicate pair in the extant Arabidopsis genome, i.e., for the remaining 86% of ancestral gene pairs, one member has been lost . Among 31 ancestral gene pairs that possessed conserved uORFs at the time immediately following the genome duplication, 12 (39%) pairs have been retained in the present Arabidopsis genome (Table 2), which is significantly higher than the genome-wide average (p = 0.0005). The conserved uORF was retained in both copies of each of the twelve retained duplicate pairs. Retention of these 12 uORFs in both paralogs suggests that they act in cis, consistent with the expectation that uORFs typically control translation of downstream mORFs on the same RNA molecule .
The overrepresentation of transcription factors among conserved uORF loci could be due, in part, to preferential retention of transcription factor recent duplicates (22.7% retention of transcription factor duplicates vs 14.4% retention genome-wide) , but this alone does not account for the high frequency of predicted transcription factors among the uORF loci. When duplicate history bias is removed by calculating GO term frequencies of the pre-genome-duplication set of loci, transcription factors are still overrepresented (11/31 loci, or 35%).
Conserved angiosperm uORF peptide sequences in primitive plants and other eukaryotes
To determine whether any of the 19 uORF homology groups conserved between rice and Arabidopsis might also be present in other eukaryotes, we searched for uORF sequences in all Genbank eukaryotic ESTs. Amino acid sequences similar to four homology groups (3, 8, 13, and 15) were detected in non-angiosperms. Group 15 was found only as distantly as a fern (Adiantum); group 3 was found as far from angiosperms as the green algae (Ulva); group 13 was found in an animal (Xenopus tropicalis); and group 8 uORF sequence was found in primitive plants, animals, fungi, and a slime mold (Figures 9 and 10). Another algal sequence (Chlamydomonas) from the Genbank non-redundant database was identified belonging to group 3 (Genbank: AJ841703). The group 13 uORF homolog found in a X. tropicalis EST was also found in a genomic contig sequence  in which the uORF homolog is flanked by genes that are more similar to animal sequences than to any known plant sequences. Thus, this group 13 uORF homolog most likely exists in the Xenopus genome rather than being an EST library contaminant.
Sequences similar to group 8 Arabidopsis and rice uORFs were found in most eukaryotes, but transcript sequence following the uORF varied among the different lineages. All land plant uORFs were associated with macrophage inhibitory cytokine-1-like (Mic1-like) mORF sequences while the mORFs downstream of the group 8 uORF homologs in nematodes and arthropods code for an unknown protein and a putative mannosyl transferase, respectively (Figure 11). Available EST sequences for each of the group 8 uORF homologs in mammals, fungi, algae, and slime mold end shortly after the conserved peptide uORF, suggesting that in these eukaryotes the uORF homolog is not associated with a mORF and is simply a short ORF. This is further supported by more than 10 human ESTs that end at the same position and include a polyA sequence. In the sea squirt lineage a putative mORF is present in the EST sequences, but a full-length cDNA sequence will be needed to further investigate this possibility.
Although there is variability in the sequences found downstream of group 8 uORFs, three features of these uORF homologs are relatively well conserved: the length of the predicted uORF, the relative positions of four cysteine codons, and the positions of two introns (Figure 9). The length of the uORF peptide ranges from 51 amino acids in Haemonchus (nematode), to 74 amino acids in humans, and length is even more highly conserved within each of the land plant, arthropod, nematode, fungal, and vertebrate lineages (59–62, 65–69, 51–68, 54–66, and 69–74 amino acids, respectively). Four cysteine residues consistently align in all eukaryotes, with nine amino acids separating the first and second cysteine residues, as well as the third and fourth cysteine residues, whereas 11–15 residues separate the second and third cysteines. Two intron positions are perfectly conserved among the land plants, vertebrates, and at least one member of the fungal lineage. The first intron lies between the third and fourth amino acids following the first conserved cysteine position, and the second intron lies between the fourth and fifth amino acids following the fourth conserved cysteine position (Figure 11). The first and/or second intron positions are present in Dictyostelium, algae, and some fungi, but are absent in nematodes, arthropods, and sea squirts.
The four cysteines are part of a putative coiled coil-helix, coiled coil-helix (CHCH) domain (Pfam accession number PF06747), also found in three small yeast proteins, Cox17p, Cox19p, and Mrp10p. Cox17p and Cox19p are required for assembly of functional cytochrome oxidase and Mrp10p is homologous to a nuclear-encoded mitochondrial ribosomal protein. A hypothetical human gene, CHCH domain 7 (CHCHD7), is also similar to the group 8 uORF, as determined by BLAST similarity searches.
Phylogenetic relationships among group 8-like ORFs
Fungal, animal, and plant representatives of each CHCH-containing ORF were identified using a BLAST search, and their evolutionary relationships were inferred using a Bayesian phylogenetic analysis (Figure 12; Additional file 1). Animal Mrp10p-like (Genbank: BC075310, DR155443 and BX935835), Debaryomyces group 8-like (Genbank: NC_006045), and Dictyostelium Cox19p-like (Genbank: XM_631387) sequences were more divergent than other sequences, causing long branch attraction . Thus, these sequences were removed from the analysis to prevent tree topology distortion. Five distinct clades were observed, which we refer to as Cox17p-like, Cox19p-like, Mrp10p-like, CHCHD7-like, and uORF group 8-like (Figure 12). All clades but one (Mrp10p-like) contain representatives from fungi, animals, and plants and are strongly supported, showing branch order probabilities greater than 0.8, which suggests that these sequences emerged in a common eukaryotic ancestor and have since diverged in the three lineages. Mrp10p-like sequences do not strongly group independently of other branches (P = 0.57), which could be due to highly divergent amino acid sequence represented by relatively long branches. The tree shows that the group 8-like proteins are a distinct clade from other CHCH domain proteins (P = 1.0), and that CHCHD7-like proteins are more closely related to group 8-like members than to other CHCH-containing proteins (P = 0.94). The tree topology also indicates that Cox17p-like and Cox19p-like genes are more closely related to each other than to other CHCH proteins (P = 0.97).
A separate phylogenetic analysis of the 46 group 8-like sequences shows that most cluster into five taxonomic groups (plants and green algae, arthropods, nematodes, vertebrates, and fungi) with strong branch support (0.85–1.00) in all but the fungal lineage (0.58; Figure 13). Sea squirt sequences group with one of two Branchiostoma sequences with weak branch support (0.53). Dictyostelium, sea urchin (Strongylocentrotus), and one further Branchiostoma sequence do not group with any of these with weak support (0.53). Sea squirt, Branchiostoma, and sea urchin sequences should be more similar to other deuterastomes (includes the vertebrate lineage) than other organisms, but the short group 8-like sequence alignment could prevent resolution of correct evolutionary relationships of some groups (Additional file 2). Despite weakly supported branches, there is strong support for independent clustering of the arthropods, nematodes, vertebrates and plants, as expected.
Although two Branchiostoma group 8-like sequences (Brafl1 and 2) suggest that there has been a duplication event within this lineage, there is no evidence for maintenance of ancient group 8-like gene duplications occurring within the plant, vertebrate, nematode, arthropod, or fungal lineages. In Arabidopsis both the recent and ancient duplicates from two WGD events have been lost from the genome. Only the Mesostigma genome contains two group 8-like transcripts. Their short branch lengths indicate that this duplication occurred relatively recently and it is possible that insufficient time has passed for loss of the second copy.
Discussion and conclusion
Comparative analysis by uORF-Finder of 5' UTRs in full-length cDNAs from two distantly related plant species, rice and Arabidopsis, identified conserved peptide uORFs in 58 Arabidopsis loci that comprised 26 uORF homology groups and in 36 rice loci that comprised 19 homology groups, increasing the number of known conserved uORF homology groups from two to 26 and providing useful, new information for investigations of regulatory biology. Because full-length cDNAs derived from both Arabidopsis and rice only represent a fraction of all nuclear genes, not all conserved uORFs are expected to be detected by this approach. Extrapolation to the whole Arabidopsis genome suggests that it possesses approximately 61 to 102 genes with conserved peptide uORFs that are also conserved in the rice genome (see Methods for calculation). An additional 24 conserved peptide uORF genes are predicted among Arabidopsis loci with retained duplicates from the most recent WGD event. In all, there are likely to be approximately 99–140 genes, or 0.38–0.53% of all protein-coding genes, with conserved peptide uORFs in the Arabidopsis genome. Because short conserved uORFs (<20 amino acids) would not have been detected by uORF-Finder, this is a conservative estimate.
To find additional conserved uORFs, more extensive collections of full-length cDNA sequences will need to be developed and/or 5' UTRs predicted from genomic sequence will be required. As full-length cDNA sequence resources become available for other plant species, such as maize  and poplar , it should be possible to identify additional conserved uORFs that might be specific to taxonomic groups, such as monocotyledons or dicotyledons. Similarly, analysis of ancient tetraploidy events in species such as poplar and maize might be able to identify uORFs conserved between retained duplicates.
Conserved uORF genes are regulatory genes
Based on the study of a few hundred genes, it has been suggested that uORFs are usually associated with mORFs that encode proteins that regulate cell growth [41, 46], but a genome-wide study of upstream AUGs (uAUGs) found no correlation of uAUG-containing transcripts with any particular gene ontology (GO) molecular function term in mammalian transcripts . These observations did not differentiate between sequence-dependent and sequence-independent uORFs. Our analysis shows that genes encoding transcription factors are overrepresented among genes predicted to encode conserved peptide uORFs, representing almost one third of the 58 Arabidopsis loci as compared to 6% of all genes. Moreover, nearly all genes whose function can be reasonably inferred appear to play some regulatory role in the biology of plants.
Do conserved peptide uORFs mediate feedback translational regulation by small regulatory molecules?
Certain eukaryotic conserved peptide uORFs are known to control translation of a downstream mORF in response to a metabolic product such as arginine or polyamines [4, 14, 47]. In the case of the fungal arginine-regulated carbamoyl-phosphate synthase subunit, a uORF codes for the arginine attenuator peptide that responds to increased arginine concentrations by causing ribosomes to stall near the 3' end of the uORF, interfering with ribosome scanning and translation of the downstream mORF . A similar mechanism has been elucidated for the regulation of AdoMetDC in which the uORF peptide interferes with the termination of uORF translation in a polyamine-dependent manner [48, 49]. In plants, sucrose is a signaling molecule that controls not only the transcription of many genes, but also translation of a class of bZIP transcription factors via their conserved uORF, suggesting the possibility of sucrose interaction with a uORF-encoded peptide to regulate translation downstream .
Our analysis identified not only these previously known examples of genes involved in pathways exhibiting small molecule feedback in a uORF sequence-dependent manner, but several additional genes that might also act via this mechanism. One is the conserved group 13 uORFs, which are present in genes that encode phosphoethanolamine N-methyltransferase (PEAMT/NMT), the key enzyme in phosphocholine (PCho) biosynthesis. Recently, NMT1 has been shown to contain a uORF that differentially affects translation of the mORF in response to exogenously added choline . This effect is observed when the uORF start codon is abolished but it remains to be determined whether the response to choline is uORF sequence-specific. Intriguingly, the group 13 uORF peptide is rich in arginine and serine (40–48% in Arabidopsis and rice genes; Table 6). A variety of arginine-rich peptides 15–20 amino acids long with 5 or more arginines bind to specific RNA sequences . The predicted group 13 uORF peptide has 5–7 arginines in a 16–17 amino acid region, well within this range, suggesting the possibility that it might bind to a specific RNA sequence, perhaps in PEAMT/NMT transcripts. The fact that the group 13 uORF peptide was also found in Xenopus suggests that its regulatory role is widespread in eukaryotes.
Another example is homology group 11, whose mORFs are predicted to encode trehalose-6-phosphate phosphatase (TPPase); trehalose-6-phosphate is postulated to regulate sugar metabolism in plants . In summary, sucrose, polyamines, phosphatidic acid, and trehalose-6-phosphate are possible regulators of translation of downstream mORFs through interaction with conserved uORFs. Also interesting in this light are group 19, which specifies an auxin-induced calmodulin-binding homolog, and group 15, which encodes a bHLH transcription factor that is believed to be subject to translational control through its conserved uORF by spermine synthase . Spermine is a polyamine signal molecule necessary for normal plant growth and defense responses.
As mentioned, six conserved uORF families specify transcription factors, one of which is regulated by the small signaling molecule sucrose. In plants, transcription factors often act quantitatively to control target gene expression proportionate to transcription factor concentration . Therefore, it is interesting to consider the possibility that translational control of transcription factor protein levels could be mediated by interaction of a conserved uORF peptide with a metabolite. This might be an effective means for quantitatively modulating the levels of expression of a pathway or network of downstream genes, for instance, in response to changing physiological or environmental conditions. This logic can equally be applied to other key control proteins and their uORFs.
How is translational control mediated by conserved peptide uORFs?
If conserved uORF peptides can regulate mORF levels in response to small molecules, they are clearly analogous to RNA sensors and riboswitches that sense small molecules and regulate transcript translation accordingly [28, 54]. It is interesting to think of conserved peptide uORFs too as sensors of cellular, physiological, or developmental conditions. Although the role of conserved uORFs as 'sensors' of cellular metabolites has been clearly established in the cases of polyamine, sucrose, and arginine concentration, it is still not clear how uORF peptides gauge cellular conditions. uORF peptides could affect mORF translation by interacting directly with the ribosomal complex, by associating with other proteins that influence the translational machinery, and/or by stabilizing or destabilizing RNA secondary structures in the 5' UTR that impede or promote mORF translation. Given the variety of uORF peptides represented in the 26 homology groups, each of these possibilities could occur one or more times.
It is perhaps interesting to note also that the uORFs of 9 homology groups are rich in serine, threonine, and/or tyrosine. These amino acids are potential targets for phosphorylation that conceivably could promote or inhibit ribosome stalling or initiation at downstream mORFs. As mentioned above, lysine/arginine-rich motifs could function in RNA binding .
Effect of nonsense-mediated decay on uORF transcripts
Because uORFs create a premature termination codon (PTC), the nonsense-mediated decay (NMD) system might target uORF transcripts for degradation. Yoine et al  carried out a microarray analysis of plants mutant in the UPF1 ortholog, which is required for NMD.
Among 75 genes that Yoine et al identified that accumulate transcripts at more than twice the level in the upf1 mutant as in wild type Arabidopsis, we found representatives of seven uORF homology groups (1, 7, 10, 12, 13, 15, and 17), suggesting that these uORF transcripts are susceptible to nonsense-mediated decay. The uORFs in these groups might work in a manner analogous to the uORF arginine attenuator protein (AAP) in the fungal CPA1 transcript. The CPA1 transcript exclusively exhibits increased levels of degradation via NMD when the AAP inhibits translation termination in response to high levels of arginine, ultimately decreasing translation using a two-pronged approach . Similarly, the above-identified plant uORFs could intensify translational inhibition of their associated mORFs by both blocking the ribosome physically and inducing the NMD pathway.
Evolutionary emergence of uORFs and a 'transcriptional fusion' model
Very little is known about how uORFs arise. In the extant rice and Arabidopsis genomes, sequences homologous to uORFs identified by uORF-Finder were observed only in 5' UTRs and never as part of another mORF, within 3' UTRs, within introns, or in non-transcribed regions. Possible origins of 5' UTR ORFs include (a) fragmentation of mORF sequences, (b) creation of an AUG or alternate start codon by random mutation within the 5' UTR and subsequent selection for the peptide sequence, and (c) relocation of other ORF sequences within the genome to the 5' UTR or upstream region of a given gene and subsequent transcriptional fusion of the two ORFs.
Transcriptional fusions occur in an estimated 2% of adjacently transcribed mRNAs in the human genome . The evolutionary history of uORF homology group 8 suggests a stable transcriptional fusion model leading to uORF emergence in plants, arthropods and nematodes. Group 8 uORFs are associated with three independent mORFs in the land plant, arthropod and nematode lineages, while the vertebrate, slime mold, algal, and fungal small ORFs that are orthologous to group 8 uORFs do not seem to be associated with mORFs. Given the phylogenetic relationships among these species , the most parsimonious explanation for the evolutionary origin of group 8 uORFs is that they originated as a small ORF transcribed independently of a mORF. Subsequently, this small ORF gene was displaced via genome rearrangements or transposition events to regions upstream of three independent large ORFs resulting in transcriptional fusions of the two previously independent transcripts. The uORFs and mORFs in the plant, nematode, and arthropod lineages have remained associated within the same transcript for 300–500 My, therefore these transcriptional fusion events seem to be stable and perhaps biologically advantageous. Evidence for other uORF emergence models, such as mORF fragmentation or de novo creation, will require further analysis of closely related organisms.
Potential dual role for uORF proteins
uORFs can regulate specific mORF protein expression in trans when the cis uORF is intact [59, 60] but it is still unclear whether uORF proteins can play additional roles in the cell. Small proteins, similar in length to uORFs, play a role in plant development and could also be involved in plant defense [61, 62]. Potentially, uORFs could affect such processes independently of their role as a translational regulator. Homology group 8 uORFs are largely conserved in length, sequence, and intron position across most eukaryotes, but in fungi, algae, slime mold, and vertebrates, the associated mORF seems to be absent. The absence of the mORF and strong conservation of the uORF amino acid sequence over one billion years in these eukaryotes indicates that, in plants, this protein could act as both a regulator of mORF expression and as a trans acting factor in the cell.
Group 13 uORFs contain peptides similar to RS motifs found in SR proteins. SR proteins are a family of proteins required for alternative and constitutive pre-mRNA splicing [63, 64]. A subset of these proteins, shuttling SR proteins, have not only been implicated in splicing but have also been shown to stimulate translation of a reporter gene when fused to the same transcript , analogous to a uORF-mORF associated pair. It is possible then, that group 13 uORF proteins could also play a dual role, as a translational regulator and trans factor.
Similarly, some uORFs in mammalian genomes might adopt these dual roles and further characterization of conserved mammalian uORFs  could resolve a dual role model.
K a /K s analyses suggest that conserved peptide uORFs are under mild to strong negative selection and might therefore be useful for resolving orthology and paralogy of specific gene pairs. For example, phylogenetic studies have sometimes failed to identify all members within a uORF homology group when only considering the mORF sequence (e.g. homology group 2). Although the bHLH transcription factor domain occurs in the mORF of all three group 2 members, none were identified in the original studies, and only two of the three members have been included in the latest description of Arabidopsis bHLH family members [67–69].
Further characterization of conserved peptide uORFs and their functional mechanisms might also provide useful tools for creating inducible or repressible expression vectors in plants. AdoMetDC1, bZIP11, and PEAMT/NMT1 protein levels are regulated by conserved uORFs in a metabolite-dependent manner (polyamine, sucrose, and choline, respectively) and other conserved uORFs might also regulate mORF translation in response to cellular compounds, such as TPPases. If this is the case, further functional characterization of conserved peptide uORFs could provide the tools necessary to build constructs that are quickly inducible or repressible at the translational level under various conditions.
Identifying conserved uORFs in rice and Arabidopsis
Corrected RIKEN and Genoscope Arabidopsis thaliana ecotype Columbia and NIAS, FAIS and RIKEN Oryza sativa spp. japonica cv Nipponbare full-length cDNA collections were used for all analysis . A cDNA's major ORF (mORF) was defined as the longest ORF starting with an AUG, the sequence upstream of this AUG was designated the 5' UTR, and upstream ORFs (uORFs) were any ORFs found in the 5' UTR starting with an AUG. All ORFs were identified using getorf . Arabidopsis mORFs were aligned to rice cDNAs using tBLASTn with an E-value cutoff = 1e-5 [72, 73] to find putative homologs. Rice cDNAs with hits below this threshold were paired with their respective Arabidopsis transcript, 5' UTR sequences extracted from both, uORFs determined using getorf, and all combinations of rice and Arabidopsis uORF peptide pairs aligned using needle . The reciprocal analysis was also performed, starting with rice full-length cDNA sequences and comparing them to Arabidopsis transcript sequences. All uORFs greater than 100 amino acids were excluded from this analysis.
All pairs with scores >50 were kept and examined manually against existing Arabidopsis transcript annotations (TAIR and TIGR) and existing ESTs to determine whether aligned peptides fall within a probable 5' UTR. To validate the putative uORFs, the first 100 amino acids of the Arabidopsis mORF were aligned to Genbank plant ESTs using tBLASTn (E-value = 1e-10, limit: Viridiplantae [orgn] NOT Arabidopsis [orgn], complexity filter off), and all retrieved plant uORF sequences were aligned to rice and Arabidopsis uORFs using ClustalW , manually adjusted, and visualized using Jalview  (Figures 1, 2, 3, 4, 5, 6, 7, 8). There were two exceptions to this procedure. Because the uORFs in group 10 are 400–600 bp upstream of the mORF AUG, only the first 25 mORF amino acids were used to search Genbank plant ESTs (first 25 amino acids are very highly conserved). Secondly, high identity was limited to the 3' end of mORFs in group 17, therefore the Arabidopsis transcript's terminal 50 amino acids were aligned to Genbank non-EST plant sequences. Support for a conserved uORF was found in the Medicago truncatula and Lotus corniculatus genomic sequences.
To test whether uORFs appear upstream of non-homologous genes, Arabidopsis uORF sequences were aligned to the entire Arabidopsis genome (version 5)  using tBLASTn (E-value = 10). Predicted conserved uORFs were found to lie upstream of the annotated gene instead of in the annotated 5' UTR in approximately 10% of Arabidopsis and 25% of rice genes (Tables 2, 3, 4). The discrepancies with the accepted annotations, found at TAIR  and TIGR , respectively, demonstrate the benefit of using full-length cDNA sequences for this analysis.
To determine whether sequences similar to these conserved uORFs reside elsewhere in the rice and Arabidopsis genomes, uORF amino acid sequences were aligned with sequences translated from the genome sequence using tBLASTn . Sequences similar to these uORFs were found within 5' UTRs of homologous mORF loci, and were absent from non-homologous transcripts, intronic regions, and intergenic regions with only one exception, Arabidopsis NMT3 (AGI locus identifier At1g73600). The annotated mORF for NMT3  is not covered by any available full-length cDNA and has no EST support at its 5' end. Thus, we annotated NMT3 by comparison with its paralog, NMT1 (At3g18000) . NMT3 possesses sequences similar to the NMT1 uORF, as well as sequences similar to the NMT1 mORF, but the TAIR annotation fuses these into a single ORF. However, NMT3 possesses potential splice sites that would produce transcripts with uORF and mORF sequences similar to those in NMT1. The NMT3 uORF predicted by one alternative splice model is the same length as, and is 72% identical to, the NMT1 uORF amino acid sequence (Group 13 in Figure 4).
The TAIR website was used to assign locus numbers for each Arabidopsis transcript and the TIGR website for rice locus numbers. The Arabidopsis locus numbers were then used to search for retained duplicates from the recent and ancient whole genome duplications as defined on the Arabidopsis Paralogon website .
Calculating K a /K s
For homology groups 1–19, K a /K s values for homologous rice and Arabidopsis mORFs and uORFs were determined using pairwise_kaks.PLS (version 1.7) . Both the approximate method (option-kaks yn00) and the maximum likelihood method (-kaks codeml) were used. Any K a /K s values resulting from a K a or K s value >10 was excluded from the analysis, as these values result in inaccurate predictions of K a /K s [80, 81]. The K a /K s values for homology groups 20–26 were determined with the same approach using Arabidopsis sequences only.
GO molecular function terms
GO molecular function terms  were retrieved from TAIR Locus History pages . GO terms for all Arabidopsis loci were downloaded from the TAIR website and used to compare genome-wide GO molecular function term frequencies to those found in the conserved uORF-containing loci. Statistically significant differences were detected using the Exact Binomial test as described in the R program package . This analysis was also carried out by GeneMerge, a program that incorporates a Bonferroni corrected P-value .
Identification of Arabidopsis ohnologs and paralogs with conserved uORF
Conserved uORFs were found in Arabidopsis duplicates in much the same way as conserved uORFs were found between rice and Arabidopsis. uORFs and mORFs were defined in the same way, and mORF sequences were aligned to the entire Arabidopsis full-length cDNA collection using BLASTp (E-value cutoff = 1e-5) to detect transcripts deriving from a duplicated locus. mORFs aligning with ≥ 99% identity were discarded, and uORFs of all remaining pairs were aligned using needle and validated as above.
Generation of phylogenetic trees
Sequences similar to Cox17p, Cox19p, Mrp10, CHCHD7, and uORF homology group 8 (as determined by tBLASTn and analyzed for conservation of the CHCH motif) were aligned using Muscle , trimmed of non-informative sites, and analyzed using Mr. Bayes v. 3.0  (rates = gamma, aamodel = mixed, ngen = 2000000). Phylogenetic trees were visualized using PHYLIP's DRAWTREE program v. 3.65 .
Sequences similar to uORF homology group 8 were aligned, edited, and analyzed in the same manner with one exception, ngen = 3000000.
Estimate of conserved peptide uORF prevalence
Number of Arabidopsis-rice loci
There is an average of 2.23 full-length cDNAs per uORF locus identified (excluding loci identified by BLAST alignment), which suggests that 15200 Arabidopsis genes are represented in the cDNA collections (34000 cDNAs/2.23 cDNAs per locus), representing approximately 60% of all Arabidopsis genes (assuming 26000 genes) . In addition, Kikuchi et al  report that the 28000 rice full-length cDNA sequences represent 20000 transcription units (TUs) and that 64% of these (12800) have a homolog in Arabidopsis. Assuming that 60–100% of these homologs are represented in the Arabidopsis cDNA collections, the estimated number of Arabidopsis homologs screened for uORF conservation is 7800–13000. Only 80% of Arabidopsis genes also have a homolog in rice (~21000) , therefore the uORF-Finder program has identified 37–62% of all conserved upstream ORFs (7800/21000 to 13000/21000) when comparing rice and Arabidopsis full-length cDNAs. Therefore, there should be 61–102 loci that contain conserved uORFs: 38 loci found by uORF-Finder, 6 additional loci found by aligning known uORF sequences with the Arabidopsis genome using BLAST, and 17–58 presently unidentified loci. Using both uORF-Finder and BLAST algorithms we estimate that between 43% and 72% of conserved peptide uORFs between monocots and dicots have been identified.
Number of Arabidopsis-Arabidopsis loci
A total of 60% of Arabidopsis genes are represented in the full-length cDNA collections used for this study. Therefore, the probability of selecting two loci that have conserved peptide uORFs from the pool of known sequences is 0.6*0.6 = 0.36. This translates to a total of 38 loci that have conserved uORFs using an Arabidopsis-Arabidopsis comparison (14 identified (36%), and 24 unidentified).
We therefore predict that there are between 99 and 140 loci in the Arabidopsis genome that contain conserved peptide uORFs, 41–58% of which have been identified.
Abbreviated species names and Genbank accession number, cDNA clone number, or genome identifier
Group 1: Arath1 (CNS0ABWH); Arath2 (CNS09Y87); Arath3 (CNS0A364); Arath4 (CNS0A728); Arath5 (RAFL11-10-D10); Orysa1 (AK070887); Orysa2 (AK065180); Orysa3 (AK064903); Orysa4 (AK109929); Orysa5 (LOC_Os12g37410).
Group2: Arath1 (At2g31280); Arath2 (At1g06150); Arath3 (RAFL04-15-e03); Lacsa (BQ869454); Lyces (AW621910); Medtr, (BF643643); Orysa (AK074015.1).
Group3: Arath1 (CNS0A7A6); Arath2 (RAFL04-16-A04); Arath3 (RAFL09-22-L13); Cycru (CB092297); Orysa1 (AK072162); Orysa2 (AK100397); Orysa3 (AK070259); Selmo (DN838497); Torru (CN201012); Ulvli (AJ892634).
Group4: Arath1 (RAFL09-11-P17); Arath2 (RAFL09-63-H05); Arath3 (RAFL06-76-P19); Brana (CD823274); Goshi (AI730427); Gosra (CO113165); Medtr (AW689516); Orysa (AK060830); Poptt (BU896557); Prupe (BU045695).
Group 5: Arath1 (RAFL05-05-C03); Arath2 (CNS0A9PN); Gosra (CO130855); Hevbr (CB376393); Lacse (BU011020); Orysa (AK103103); Phaac (BU791117); Triae (BJ233459).
Group 6: Arath1 (RAFL05-17-I08); Arath2 (CNS0A6ZP); Aspof (CV291431); Glyma (BM143067); Gosar (BG442153); Orysa (AK064902); Pontr (CD576165); Triae (CK161649); Vitvi (CB980452).
Group 7: Arath (RAFL09-25-N17); Brana (CD836460); Mescr (BM301482); Nicbe (CK290710); Orysa1 (AK067685); Orysa2 (LOC_Os06g48350); Triae (CV066319).
Group 8: Arath (RAFL07-08-P17); Chlre (BE121764); Mesvi1 (DN255332); Mesvi2 (DN261354); Orysa (AK072620); Phypa (BJ174896); Popca (CX178804); Popeu (AJ776458); Sachc (CF573523); Triae (CA499582).
Group 9: Allce (CF443194); (Arath1 (RAFL07-09-G06); Arath2 (RAFL09-23-F23); Arath3 (At1g64140); Gosra (CO081490); Orysa1 (AK101398); Orysa2 (AK105763); Orysa3 (AK068099); Orysa4 (AK099577).
Group 10: Arath1 (RAFL07-11-O11); Arath2 (RAFL09-17-I10); Brana (CN732239); Orysa1 (AK069526); Orysa2 (AK100056); Poptt (BI131713); Sorbi (CN139168); Theha (BE758596).
Group 11: Arath1 (RAFL07-14-D12); Arath2 (CNS0A404); Glyma (CA783255); Jugre (CV197923); Medtr (AW691064); Orysa1 (AK103391); Orysa2 (AK069361); Soltu (BQ113418).
Group 12: Arath1 (RAFL07-18-F03); Arath2 (CNS0AB39); Brana (CD812479); Citse (CN185367); Jugre (CV196770); Orysa (AK060405); Popde (CK319714); Triae (BQ752938); Zeama (CD433782).
Group 13: Arath1 (RAFL08-10-M03); Arath2 (At1g48600.2); Arath3 (At1g73600); Cycru (CB093136); Gosra (CO080661); Iponi (BJ562806); Linus (CA483285); Medtr (AW587372); Orysa1 (LOC_Os05g47540); Orysa2 (AK102037); Phypa (BJ204269); Xenla (CA792398); Xentr (CX412233); Zeama (AY103779).
Group14: Allce (CF450799); Arath (RAFL09-10-M04); Medtr (AW267817); Nicbe (CK295530); Orysa (AK101569); Soltu (CK258175); Zeama (CO519993). Group 15: Adica (BP914226); Arath1 (CNS0ADY7); Arath2 (RAFL08-17-G21); Arath3 (RAFL04-17-N21); Arath4 (RAFL16-69-M04); Citpa (DN959636); Gosra (CO125506); Maldo (CV082382); Medtr (CX528608); Orysa1 (AK102703); Orysa2 (AK101749); Orysa3 (AK071582); Orysa4 (AK065674); Sacof (CA154823); Vitvi (CB001711); Welma (DT579937).
Group 15: Adica (BP914226); Arath1 (CNS0ADY7); Arath2 (RAFL08-17-G21); Arath3 (RAFL04-17-N21); Arath4 (RAFL16-69-M04); Citpa (DN959636); Gosra (CO125506); Maldo (CV082382); Medtr (CX528608); Orysa1 (AK102703); Orysa2 (AK101749); Orysa3 (AK071582); Orysa4 (AK065674); Sacof (CA154823); Vitvi (CB001711); Welma (DT579937).
Group 16: Arath (CNS0A4RC); Medtr (AW693231); Orysa1 (AK071885); Orysa2 (AK067447).
Group 17: Arath1 (RAFL09-25-E19); Arath2 (At5g03190); Arath3 (RAFL19-67-G09); Arath4 (At5g01710); Gosra (CO108440); Lyces (AW738430); Medtr1 (BQ149694); Medtr2 (AC144517); Orysa1 (AK69088); Orysa2 (AK070250); Sacof (CA191644).
Group 18: Arath (RAFL08-18-B11); Gosra (CO115325); Nicbe (CK286574); Orysa (AK061433).
Group 19: Arath (CNS09ZXM); Eupes (DV113097); Helan (AJ541596); Medtr (BI309364); Orysa (AK068270); Triae (CD927685); Vitvi (CB918939); Zeama (DV166198).
Group 20: Arath1 (RAFL04-17-G13); Arath2 (CNS0A8YX); Brana (CD835762); Brugy (BP941533); Gosar (BF274209); Maldo (CN940921); Medtr (BE316669); Styhu (L36823).
Group 21: Allce (CF450138); Arath1 (RAFL07-08-G04); Arath2 (RAFL21-49-G19); Betvu (BQ594525); Brana (CD835573); Erate (DN481483); Escca (CD481239); Eupti (BP958766); Gosra (CO074819); Glyma (BU761432); Horvu (AV834976); Lacse (BQ998418); Maldo (CV881926); Medtr (CA991201); Orysa (AK100575); Popca (CX182168).
Group 22: Arath1 (RAFL07-11-D20); Arath2 (RAFL11-03-J07); Brana (CD836422); Horvu (CA023398); Orysa (CK041713); Sacof (CA242575); Triae (BJ247925); Zeama (CO458204).
Group 23: Arath1 (RAFL07-11-L03); Arath2 (RAFL09-07-L11); Citsi (CV720092); Glyma (BI892512).
Group 24: Arath1 (RAFL07-14-J09); Arath2 (CNS0A44P); Brana (CD828343); Glyma (BI471587); Horvu (BQ471053); Orysa (AK119634); Sacof (CA118382); Sorbi (CB928687); Triae (CA483985); Zeama (CO520078).
Group 25: Arath1 (RAFL09-94-P19); Arath2 (CNS0A6N0); Brana (CD835519); Citsi (CN191447); Escca (CD481312); Glyma (BE805986); Phaco (CA913939); Soltu (DN940765); Vitsh (CV098492).
Group 26: Arath1 (CNS0A7NI); Arath2 (CNS0A1F5); Citja (CO912573); Pethy (CV298852); Poptd (CN521002); Prupe (BU045483).
Arath (RAFL07-08-P17); Caeel (U10402); Ciosa (BW577210); Danre (CO350578); Dicdi (AU072562); Drome (AI297387); Homsa (BU541024); Mesvi (DN255332); Neucr (BX284746); Orysa, (AK072620); Phypa (BJ174896); Strpu (CX079489); Ustma (CF644197).
Arath1 (RAFL08-10-M03); Arath2 (At1g48600.2); Arath3 (At1g73600); Cycru (CB093136); Gosra (CO080661); Iponi (BJ562806); Linus (CA483285); Medtr (AW587372); Orysa1 (LOC_Os05g47540); Orysa2 (AK102037); Phypa (BJ204269); Xenla (CA792398); Xentr (CX412233); Zeama (AY103779).
Acypi (CV847404); Ajeca (CV605785); Anoga1 (BX617953), Anoga2 (XM_552406); Apime (NW_622706); Arath1 (BP562704), Arath2 (AY065264), Arath3 (RAFL07-08-P17), Arath4 (NM_179521), Arath5 (NM_112400); Ascsu (BM964977); Bosta (CO877216); Brafl1 (BW786058), Brafl2 (BW840607); Caeel (U10402); Chlre1 (BE121764), Chlre2 (AF280543); Cicli (CV156944); Ciosa (BW577210); Cocpo (CO006101); Cryne (XM_572394); Danre (CO350578); Debha (NC_006045); Dicdi1 (AU072562), Dicdi2 (XM_631387); Drome1 (AI297387), Drome2 (AY102691); Drops (DR121964), Erate (DN481021); Galga1 (BX935835), Galga2 (CR407540); Gibze (BI750032); Glomo (BX557417); Glyso (BG045953); Haeco (CA956938); Hetgl (CB299856); Homsa1 (DR155443), Homsa2 (CR607136), Homsa3 (BU541024), Homsa4 (AY957566), Homsa5 (NM_005694); Hordvu (BF628344); Locmi1 (CO854527), Locmi2 (CO825844); Mesvi1 (DN255332), Mesvi2 (DN261354); Molte (CJ368011); Musmu1 (BC030366), Musmu2 (AK010111); Neucr (BX284746); Oncmy (BX081024); Oryla (BJ737531); Orysa1 (XM_482456), Orysa2 (AK072620), Orysa3 (AK120143), Orysa4 (XM_468245); Parbr (CA581923); Phypa1 (BJ966696), Phypa2 (BJ174896); Popca (CX178804); Popeu (AJ776458); Sacce1 (NC_001136), Sacce2 (AY692601), Sacce3 (NC_001144), Sacce4 (NC_001144); Sachc (CF573523); Schma (CD081475); Schpo1 (NM_001019463), Schpo2 (NM_001022867), Schpo3 (NM_001022571); Sorbi (CD423660); Strpu (CX079489); Strra (BI323578); Tetni1 (CR709012), Tetni2 (CNS0G27U); Triae (CA499582); Trisp (BQ693345); Ustma1 (CF644197), Ustma2 (XM_754796); Xenla1 (BI477811), Xenla2 (BC084847); Xentr1 (BC075310), Xentr2 (CN119217); Yarli (XM_500713).
Acypi (CV847404); Ajeca (CV605785); Anoga (BX617953); Apime (NW_622706); Arath (RAFL07-08-P17); Ascsu (BM964977); Bosta (CO877216); Brafl1 (BW840607), Brafl2 (BW786058); Caeel (U10402); Chlre (BE121764); Cicli (CV156944); Ciosa (BW577210); Danre (CO350578); Debha (NC_006045); Dicdi (AU072562); Drome (AI297387); Drops (DR121964); Galga (CR407540); Gibze (BI750032); Glomo (BX557417); Haeco (CA956938); Hetgl (CB299856); Homsa (BU541024); Locmi (CO825844); Mesvi1 (DN255332), Mesvi2 (DN261354); Molte (CJ368011); Musmu (AK010111); Neucr (BX284746); Oncmy (BX081024); Oryla (BJ737531); Orysa (AK072620); Phypa (BJ174896); Popca (CX178804); Popeu (AJ776458); Sachc (CF573523); Schma (CD081475); Sorbi (CD423660); Strpu (CX079489); Strra (BI323578); Tetni (CR709012); Triae (CA499582); Trisp1 (BQ693345), Trisp2 (BQ692350); Ustma (CF644197); Xenla (BI477811); Xentr (CN119217).
Homo sapiens Horvu, Hordeum vulgare
Populus trichocarpa ×Populus deltoides
Populus tremula ×Populus tremuloides
Saccharum hybrid cultivar
Hanfrey C, Franceschetti M, Mayer MJ, Illingworth C, Michael AJ: Abrogation of upstream open reading frame-mediated translational control of a plant S-adenosylmethionine decarboxylase results in polyamine disruption and growth perturbations. J Biol Chem. 2002, 277: 44131-44139. 10.1074/jbc.M206161200.
Hinnebusch AG: Translational regulation of yeast GCN4. J Biol Chem. 1997, 272: 21661-21664. 10.1074/jbc.272.35.21661.
Werner M, Feller A, Messenguy F, Pierard A: The leader peptide of yeast gene CPA1 is essential for the translational repression of its expression. Cell. 1987, 49: 805-813. 10.1016/0092-8674(87)90618-0.
Wiese A, Elzinga N, Wobbes B, Smeekens S: A conserved upstream open reading frame mediates sucrose-induced repression of translation. Plant Cell. 2004, 16: 1717-1729. 10.1105/tpc.019349.
Galagan JE, Calvo SE, Cuomo C, Ma LJ, Wortman JR, Batzoglou S, Lee SI, Basturkmen M, Spevak CC, Clutterbuck J, et al: Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature. 2005, 438: 1105-1115. 10.1038/nature04341.
Churbanov A, Rogozin IB, Babenko VN, Ali H, Koonin EV: Evolutionary conservation suggests a regulatory function of AUG triplets in 5'-UTRs of eukaryotic genes. Nucleic Acids Res. 2005, 33: 5512-5520. 10.1093/nar/gki847.
Kawaguchi R, Bailey-Serres J: mRNA sequence features that contribute to translational regulation in Arabidopsis. Nucleic Acids Res. 2005, 33: 955-965. 10.1093/nar/gki240.
Futterer J, Hohn T: Role of an upstream open reading frame in the translation of polycistronic mRNAs in plant cells. Nucleic Acids Res. 1992, 20: 3851-3857. 10.1093/nar/20.15.3851.
Kozak M: Effects of intercistronic length on the efficiency of reinitiation by eucaryotic ribosomes. Mol Cell Biol. 1987, 7: 3438-3445.
Kozak M: Pushing the limits of the scanning mechanism for initiation of translation. Gene. 2002, 299: 1-34. 10.1016/S0378-1119(02)01056-9.
Luukkonen BG, Tan W, Schwartz S: Efficiency of reinitiation of translation on human immunodeficiency virus type 1 mRNAs is determined by the length of the upstream open reading frame and by intercistronic distance. J Virol. 1995, 69: 4086-4094.
Gopfert U, Kullmann M, Hengst L: Cell cycle-dependent translation of p27 involves a responsive element in its 5'-UTR that overlaps with a uORF. Hum Mol Genet. 2003, 12: 1767-1779. 10.1093/hmg/ddg177.
Gray TA, Saitoh S, Nicholls RD: An imprinted, mammalian bicistronic transcript encodes two independent proteins. Proc Natl Acad Sci USA. 1999, 96: 5616-5621. 10.1073/pnas.96.10.5616.
Fang P, Wang Z, Sachs MS: Evolutionarily conserved features of the arginine attenuator peptide provide the necessary requirements for its function in translational regulation. J Biol Chem. 2000, 275: 26710-26719. 10.1074/jbc.275.12.8945.
Jin X, Turcott E, Englehardt S, Mize GJ, Morris DR: The two upstream open reading frames of oncogene mdm2 have different translational regulatory properties. J Biol Chem. 2003, 278: 25716-25721. 10.1074/jbc.M300316200.
Lee J, Park EH, Couture G, Harvey I, Garneau P, Pelletier J: An upstream open reading frame impedes translation of the huntingtin gene. Nucleic Acids Res. 2002, 30: 5110-5119. 10.1093/nar/gkf664.
Lincoln AJ, Monczak Y, Williams SC, Johnson PF: Inhibition of CCAAT/enhancer-binding protein alpha and beta translation by upstream open reading frames. J Biol Chem. 1998, 273: 9552-9560. 10.1074/jbc.273.16.9552.
Hill JR, Morris DR: Cell-specific translational regulation of S-adenosylmethionine decarboxylase mRNA. J Biol Chem. 1993, 268: 726-731.
Lee MM, Lee SH, Park KY: Characterization and expression of two members of the S-adenosylmethionine decarboxylase gene family in carnation flower. Plant Mol Biol. 1997, 34: 371-382. 10.1023/A:1005811229988.
Martinez-Garcia JF, Moyano E, Alcocer MJ, Martin C: Two bZIP proteins from Antirrhinum flowers preferentially bind a hybrid C-box/G-box motif and help to define a new sub-family of bZIP transcription factors. Plant J. 1998, 13: 489-505. 10.1046/j.1365-313X.1998.00050.x.
Evans PT, Malmberg RL: Do polyamines have roles in plant development?. Annu Rev Plant Phys. 1989, 40: 235-269.
Walden R, Cordeiro A, Tiburcio AF: Polyamines: small molecules triggering pathways in plant growth and development. Plant Physiol. 1997, 113: 1009-1013. 10.1104/pp.113.4.1009.
Zhang Z, Dietrich FS: Identification and characterization of upstream open reading frames (uORF) in the 5' untranslated regions (UTR) of genes in Saccharomyces cerevisiae. Curr Genet. 2005, 48: 77-87. 10.1007/s00294-005-0001-x.
Castelli V, Aury JM, Jaillon O, Wincker P, Clepet C, Menard M, Cruaud C, Quetier F, Scarpelli C, Schachter V, et al: Whole genome sequence comparisons and "full-length" cDNA sequences: a combined approach to evaluate and improve Arabidopsis genome annotation. Genome Res. 2004, 14: 406-413. 10.1101/gr.1515604.
Kikuchi S, Satoh K, Nagata T, Kawagashira N, Doi K, Kishimoto N, Yazaki J, Ishikawa M, Yamada H, Ooka H, et al: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301: 376-379. 10.1126/science.1081288.
Seki M, Narusaka M, Kamiya A, Ishida J, Satou M, Sakurai T, Nakajima M, Enju A, Akiyama K, Oono Y, et al: Functional annotation of a full-length Arabidopsis cDNA collection. Science. 2002, 296: 141-145. 10.1126/science.1071006.
Wolfe KH, Gouy M, Yang YW, Sharp PM, Li WH: Date of the monocot-dicot divergence estimated from chloroplast DNA sequence data. Proc Natl Acad Sci USA. 1989, 86: 6201-6205. 10.1073/pnas.86.16.6201.
Sanderson MJ: A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol Biol Evol. 1997, 14: 1218-1231.
Chaw SM, Chang CC, Chen HL, Li WH: Dating the monocot-dicot divergence and the origin of core eudicots using whole chloroplast genomes. J Mol Evol. 2004, 58: 424-441. 10.1007/s00239-003-2564-9.
Wolfe K: Robustness – it's not where you think it is. Nat Genet. 2000, 25: 3-4. 10.1038/75560.
Schranz ME, Mitchell-Olds T: Independent ancient polyploidy events in the sister families Brassicaceae and Cleomaceae. Plant Cell. 2006, 18: 1152-1165. 10.1105/tpc.106.041111.
Blanc G, Wolfe KH: Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell. 2004, 16: 1679-1691. 10.1105/tpc.021410.
Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 2003, 13: 137-144. 10.1101/gr.751803.
Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, et al: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290: 2105-2110. 10.1126/science.290.5499.2105.
Marchler-Bauer A, Anderson JB, Cherukuri PF, DeWeese-Scott C, Geer LY, Gwadz M, He S, Hurwitz DI, Jackson JD, Ke Z, et al: CDD: a Conserved Domain Database for protein classification. Nucleic Acids Res. 2005, 33: D192-196. 10.1093/nar/gki069.
Zdobnov EM, Apweiler R: InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848. 10.1093/bioinformatics/17.9.847.
Eastmond PJ, Graham IA: Trehalose metabolism: a regulatory role for trehalose-6-phosphate?. Curr Opin Plant Biol. 2003, 6: 231-235. 10.1016/S1369-5266(03)00037-2.
Cruz-Ramirez A, Lopez-Bucio J, Ramirez-Pimentel G, Zurita-Silva A, Sanchez-Calderon L, Ramirez-Chavez E, Gonzalez-Ortega E, Herrera-Estrella L: The xipotl mutant of Arabidopsis reveals a critical role for phospholipid metabolism in root system development and epidermal cell integrity. Plant Cell. 2004, 16: 2020-2034. 10.1105/tpc.103.018648.
Mou Z, Wang X, Fu Z, Dai Y, Han C, Ouyang J, Bao F, Hu Y, Li J: Silencing of phosphoethanolamine N-methyltransferase results in temperature-sensitive male sterility and salt hypersensitivity in Arabidopsis. Plant Cell. 2002, 14: 2031-2043. 10.1105/tpc.001701.
Wang X: Regulatory functions of phospholipase D and phosphatidic acid in plant growth, development, and stress responses. Plant Physiol. 2005, 139: 566-573. 10.1104/pp.105.068809.
Geballe AP, Sachs MS: Translational control by upstream open reading frames. Translational control of gene expression. Edited by: Sonenberg N, Hershey JWB, Mathews MB. 2000, Cold Spring Harbor, New York: Cold Spring Harbor Laboratory Press, 595-614.
UCSC X. tropicalis BLAT search. [http://genome.ucsc.edu/cgi-bin/hgBlat]
Felsenstein J: Cases in which parsimony or compatibility methods will be positively misleading. Syst Zool. 1978, 27: 401-410. 10.2307/2412923.
The Maize full length cDNA Project. [http://www.maizecdna.org]
RIKEN Poplar full-length cDNA clones. [http://www.brc.riken.jp/lab/epd/Eng/catalog/poplar.shtml]
Kozak M: An analysis of 5'-noncoding sequences from 699 vertebrate messenger RNAs. Nucleic Acids Res. 1987, 15: 8125-8148. 10.1093/nar/15.20.8125.
Hanfrey C, Elliott KA, Franceschetti M, Mayer MJ, Illingworth C, Michael AJ: A dual upstream open reading frame-based autoregulatory circuit controlling polyamine-responsive translation. J Biol Chem. 2005, 280: 39229-39237. 10.1074/jbc.M509340200.
Law GL, Raney A, Heusner C, Morris DR: Polyamine regulation of ribosome pausing at the upstream open reading frame of S-adenosylmethionine decarboxylase. J Biol Chem. 2001, 276: 38036-38043.
Raney A, Law GL, Mize GJ, Morris DR: Regulated translation termination at the upstream open reading frame in S-adenosylmethionine decarboxylase mRNA. J Biol Chem. 2002, 277: 5988-5994. 10.1074/jbc.M108375200.
Tabuchi T, Okada T, Azuma T, Nanmori T, Yasuda T: Posttranscriptional regulation by the upstream open reading frame of the phosphoethanolamine N-methyltransferase gene. Biosci Biotechnol Biochem. 2006, 70: 2330-2334. 10.1271/bbb.60309.
Bayer TS, Booth LN, Knudsen SM, Ellington AD: Arginine-rich motifs present multiple interfaces for specific binding by RNA. RNA. 2005, 11: 1848-1857. 10.1261/rna.2167605.
Imai A, Hanzawa Y, Komura M, Yamamoto KT, Komeda Y, Takahashi T: The dwarf phenotype of the Arabidopsis acl5 mutant is suppressed by a mutation in an upstream ORF of a bHLH gene. Development. 2006, 133: 3575-3585. 10.1242/dev.02535.
Hollick JB, Patterson GI, Asmundsson IM, Chandler VL: Paramutation alters regulatory control of the maize pl locus. Genetics. 2000, 154: 1827-1838.
Lai EC: RNA sensors and riboswitches: self-regulating messages. Curr Biol. 2003, 13: R285-291. 10.1016/S0960-9822(03)00203-3.
Yoine M, Ohto MA, Onai K, Mita S, Nakamura K: The lba1 mutation of UPF1 RNA helicase involved in nonsense-mediated mRNA decay causes pleiotropic phenotypic changes and altered sugar signalling in Arabidopsis. Plant J. 2006, 47: 49-62. 10.1111/j.1365-313X.2006.02771.x.
Gaba A, Jacobson A, Sachs MS: Ribosome occupancy of the yeast CPA1 upstream open reading frame termination codon modulates nonsense-mediated mRNA decay. Mol Cell. 2005, 20: 449-460. 10.1016/j.molcel.2005.09.019.
Akiva P, Toporik A, Edelheit S, Peretz Y, Diber A, Shemesh R, Novik A, Sorek R: Transcription-mediated gene fusion in the human genome. Genome Res. 2006, 16: 30-36. 10.1101/gr.4137606.
Douzery EJ, Snell EA, Bapteste E, Delsuc F, Philippe H: The timing of eukaryotic evolution: does a relaxed molecular clock reconcile proteins and fossils?. Proc Natl Acad Sci USA. 2004, 101: 15386-15391. 10.1073/pnas.0403984101.
Parola AL, Kobilka BK: The peptide product of a 5' leader cistron in the beta 2 adrenergic receptor mRNA inhibits receptor synthesis. J Biol Chem. 1994, 269: 4497-4505.
Pendleton LC, Goodwin BL, Solomonson LP, Eichler DC: Regulation of endothelial argininosuccinate synthase expression and NO production by an upstream open reading frame. J Biol Chem. 2005, 280: 24252-24260. 10.1074/jbc.M500106200.
Marton ML, Cordts S, Broadhvest J, Dresselhaus T: Micropylar pollen tube guidance by egg apparatus 1 of maize. Science. 2005, 307: 573-576. 10.1126/science.1104954.
Wen JQ, Lease KA, Walker JC: DVL, a novel class of small polypeptides: overexpression alters Arabidopsis development. Plant J. 2004, 37: 668-677. 10.1111/j.1365-313X.2003.01994.x.
Graveley BR: Sorting out the complexity of SR protein functions. RNA. 2000, 6: 1197-1211. 10.1017/S1355838200000960.
Fu XD: The superfamily of arginine/serine-rich splicing factors. RNA. 1995, 1: 663-680.
Sanford JR, Gray NK, Beckmann K, Caceres JF: A novel role for shuttling SR proteins in mRNA translation. Genes Dev. 2004, 18: 755-768. 10.1101/gad.286404.
Crowe ML, Wang XQ, Rothnagel JA: Evidence for conservation and selection of upstream open reading frames suggests probable encoding of bioactive peptides. BMC Genomics. 2006, 7: 16-10.1186/1471-2164-7-16.
Toledo-Ortiz G, Huq E, Quail PH: The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003, 15: 1749-1770. 10.1105/tpc.013839.
Bailey PC, Martin C, Toledo-Ortiz G, Quail PH, Huq E, Heim MA, Jakoby M, Werber M, Weisshaar B: Update on the basic helix-loop-helix transcription factor gene family in Arabidopsis thaliana. Plant Cell. 2003, 15: 2497-2502. 10.1105/tpc.151140.
Heim MA, Jakoby M, Werber M, Martin C, Weisshaar B, Bailey PC: The basic helix-loop-helix transcription factor family in plants: a genome-wide study of protein structure and functional diversity. Mol Biol Evol. 2003, 20: 735-747. 10.1093/molbev/msg088.
Hayden CA, Wheeler TJ, Jorgensen RA: Evaluating and improving cDNA sequence quality with cQC. Bioinformatics. 2005, 21: 4414-4415. 10.1093/bioinformatics/bti709.
EMBOSS-European Molecular Biology Open Software Suite. [http://emboss.sourceforge.net]
Gish W, States DJ: Identification of protein coding regions by database similarity search. Nat Genet. 1993, 3: 266-272. 10.1038/ng0393-266.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Clamp M, Cuff J, Searle SM, Barton GJ: The Jalview Java Alignment Editor. Bioinformatics. 2004, 20: 426-427. 10.1093/bioinformatics/btg430.
The Arabidopsis Information Resource. [http://www.arabidopsis.org]
TIGR Rice Genome Annotation Project – web BLASTserver. [http://tigrblast.tigr.org/euk-blast/index.cgi?project=osa1]
TAIR locus At1g73600. [http://www.arabidopsis.org/servlets/TairObject?id=29540&type=locus]
Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15: 568-573.
Anisimova M, Bielawski JP, Yang Z: Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution. Mol Biol Evol. 2001, 18: 1585-1592.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
The R Project for statistical computing. [http://www.r-project.org]
Castillo-Davis CI, Hartl DL: GeneMerge – post-genomic analysis, data mining, and hypothesis testing. Bioinformatics. 2003, 19: 891-892. 10.1093/bioinformatics/btg114.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Felsenstein J: PHYLIP-Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.
AGI: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815. 10.1038/35048692.
Rook F, Gerrits N, Kortstee A, van Kampen M, Borrias M, Weisbeek P, Smeekens S: Sucrose-specific signalling represses translation of the Arabidopsis ATB2 bZIP transcription factor gene. Plant J. 1998, 15: 253-263. 10.1046/j.1365-313X.1998.00205.x.
Satoh-Nagasawa N, Nagasawa N, Malcomber S, Sakai H, Jackson D: A trehalose metabolic enzyme controls inflorescence architecture in maize. Nature. 2006, 441: 227-230. 10.1038/nature04725.
Verhagen BW, Glazebrook J, Zhu T, Chang HS, van Loon LC, Pieterse CM: The transcriptome of rhizobacteria-induced systemic resistance in Arabidopsis. Mol Plant Microbe Interact. 2004, 17: 895-908. 10.1094/MPMI.2004.17.8.895.
Henriksson E, Olsson AS, Johannesson H, Johansson H, Hanson J, Engstrom P, Soderman E: Homeodomain leucine zipper class I genes in Arabidopsis. Expression patterns and phylogenetic relationships. Plant Physiol. 2005, 139: 509-518. 10.1104/pp.105.063461.
Aoyama T, Dong CH, Wu Y, Carabelli M, Sessa G, Ruberti I, Morelli G, Chua NH: Ectopic expression of the Arabidopsis transcriptional activator Athb-1 alters leaf cell fate in tobacco. Plant Cell. 1995, 7: 1773-1785. 10.1105/tpc.7.11.1773.
Bharti K, Von Koskull-Doring P, Bharti S, Kumar P, Tintschl-Korbitzer A, Treuter E, Nover L: Tomato heat stress transcription factor HsfB1 represents a novel type of general transcription coactivator with a histone-like motif interacting with the plant CREB binding protein ortholog HAC1. Plant Cell. 2004, 16: 1521-1535. 10.1105/tpc.019927.
Nover L, Scharf KD, Gagliardi D, Vergne P, Czarnecka-Verner E, Gurley WB: The Hsf world: classification and properties of plant heat stress transcription factors. Cell Stress Chaperones. 1996, 1: 215-223. 10.1379/1466-1268(1996)001<0215:THWCAP>2.3.CO;2.
Yang T, Poovaiah BW: Molecular and biochemical evidence for the involvement of calcium/calmodulin in auxin action. J Biol Chem. 2000, 275: 3137-3143. 10.1074/jbc.275.5.3137.
Nakano T, Suzuki K, Fujimura T, Shinshi H: Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006, 140: 411-432. 10.1104/pp.105.073783.
Gutterson N, Reuber TL: Regulation of disease resistance pathways by AP2/ERF transcription factors. Curr Opin Plant Biol. 2004, 7: 465-471. 10.1016/j.pbi.2004.04.007.
PlantsP kinase classification. [http://plantsp.genomics.purdue.edu/html/families.html]
Sakuma Y, Liu Q, Dubouzet JG, Abe H, Shinozaki K, Yamaguchi-Shinozaki K: DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration- and cold-inducible gene expression. Biochem Biophys Res Commun. 2002, 290: 998-1009. 10.1006/bbrc.2001.6299.
We thank Dr Yadegari for suggesting this project. We would also like to recognize N Merchant and S Miller at the Biotechnology Computing Facility, as well as T Wheeler in the Department of Computer Science for invaluable programming suggestions. This research was supported by a University of Arizona NSF IGERT Genomics Initiative fellowship (DGE-0114420) to CAH and the NSF Plant Genome Program Grant No. DBI-0421679 (RAJ).
Both CAH and RAJ designed and implemented the analyses for the present study. CAH drafted the manuscript and RAJ provided critical comments. Both authors have read and approved the final manuscript.