<title>Background</title> <p>RNA-omics, or more simply R-nomics, is the large-scale study of RNA structure and function [##REF##15980545##1##]. One of the major challenges faced by R-nomics is to understand the regulatory mechanisms of complex signals found in the untranslated regions (UTRs) of messenger RNAs. In particular, the control signals found in the 5'-UTR of some eukaryotic mRNAs play a crucial role in translational control that can result in rapid changes to the proteome [##REF##15520293##2##]. These post-transcriptionally regulated mRNAs frequently encode important regulatory proteins (e.g., proto-oncogenes, growth factors, and transcription factors) [##REF##11897027##3##] that need to be strongly or precisely regulated for normal cellular activity. In other cases, control signals in the 5'-UTR provide continuous regulation of essential mRNAs by providing an alternative route for translation when cap-dependent translation is compromised (e.g., under stress conditions) [##REF##11050335##4##].</p> <p>Translational control signals are often found in long 5'-UTRs (&gt;100 nt) [##REF##3313277##5##] where they can contain either a single control signal [##REF##10829027##6##] or multiple control signals that function independently [##REF##11244117##7##] or in a coordinated fashion [##REF##12757712##8##, ####REF##11139406##9##, ##REF##12730202##10####12730202##10##]. One important translational control signal found in both prokaryotes and eukaryotes is the upstream open reading frame (uORF), a small open reading frame located upstream of the main coding region [##REF##8801438##11##].</p> <p>Two types of functional upstream open reading frames have been described that have a demonstrated activity either <italic>in-vitro </italic>or <italic>in-vivo</italic>: a) uORFs encoding bioactive peptides [##REF##17663791##12##, ####REF##16438715##13##, ##REF##18237443##14##, ##REF##15777708##15####15777708##15##] that either act on translation or have biological roles other than reducing the translation of the main ORF, and therefore can be described as sequence-dependent, and b) sequence-independent uORFs. A sequence-dependent uORF encodes a small peptide, and some of these uORF-encoded peptides have been shown to directly affect translation via either ribosomal stalling during translation of the uORF or termination of translation by inhibiting the peptidyl transferase activity of the ribosome and thus peptide bond formation [##REF##11707416##16##,##REF##8636015##17##]. For sequence-independent uORFs, the uORF-encoded peptide is not important for translational control, but other factors like uORF recognition, length, stop codon environment, and the downstream intercistronic sequence (length and structure) can affect reinitiation efficiency at the downstream ORF [##REF##12890013##18##,##REF##12117416##19##]. Sequence-independent uORFs can also indirectly affect translation by allowing ribosomes to bypass inhibitory stem structures [##REF##10938098##20##] or activate dormant internal ribosome entry sites (IRES) [##REF##12757712##8##] via conformational changes induced by the translation of the uORF. These distinct mechanisms of translational control have been proven to be important through <italic>in-vitro </italic>genetic (mutational analyses) and biochemical (toe-printing) assays [##REF##11707416##16##].</p> <p>There are two known pathways where uORFs can influence mRNA stability. Studies in yeast have indicated that both sequence-dependent and sequence-independent uORFs can cause mRNA destabilisation by the nonsense-mediated mRNA decay pathway [##REF##10892745##21##]. Mutations in the mRNA 5'-UTR that insert an uORF trigger the nonsense-mediated decay pathway and lead to decapping of the mRNA. Alternatively, mRNA destabilisation can occur via the termination dependent decay pathway [##REF##10357825##22##]. In this pathway, the 40S ribosomal units are released from the mRNA due to features such as stop codon environment (e.g., GC rich) or short intercistronic sequence containing a secondary structure. Release of the 40S ribosomal units prevent reinitiation of translation downstream of the uORF, and the mRNA becomes susceptible to decay. The mechanisms underlying both uORF nonsense-mediated decay and post-termination mediated decay remain unclear.</p> <p>Identifying uORFs involved in regulation of gene expression remains a challenge [##REF##11707416##16##,##REF##17913625##23##,##REF##17045969##24##]. Recently it has been estimated that it would take 20 man-months to find a single functional uORF by random selection and testing of yeast mRNAs [##UREF##0##25##]. To overcome this problem Selpi <italic>et al</italic>. [##UREF##0##25##] used an artificial intelligence approach called inductive logic programming to identify likely functional uORFs. The approach used rules based on background knowledge of uORFs in yeast mRNAs and as such may not be applicable to other organisms such as plants.</p> <p>Another approach for identifying sequence-independent uORFs was recently described [##REF##18358843##26##]. Kochetov <italic>et al</italic>. [##REF##18358843##26##] selected human mRNAs with specific sequence organisation (i.e., uORF overlapping the main ORF) that could facilitate reinitiation at downstream start codons. If the downstream start codons were nested in-frame with the main ORF then potentially N-terminally truncated variants of the main protein could be produced via reinitiation. Kochetov <italic>et al</italic>. [##REF##18358843##26##] reported that 297 out of 754 mRNAs (39% of the sub-sample) contained this specific sequence organisation with an average intercistronic spacer of 66 ± 77 nt, which provides sufficient space for reinitiation. This novel approach highlights another way in which uORFs can be functional via the generation of novel protein isoforms.</p> <p>The number of characterised uORFs in plants is apparently less than 100 (0.3%) based on a PubMed search, and about four have been identified in cereals. They include the uORFs of the S-adenosylmethionine decarboxylase gene (AdoMetDC) in both monocots and dicots [##REF##11139406##9##,##REF##15821146##27##,##REF##16982115##28##], rice myb7 gene [##REF##11855732##29##]; transcription factorssuch as maize Opaque-2 [##REF##8439744##30##], maize <italic>R </italic>[##REF##9761799##31##], and maize <italic>Lc </italic>[##REF##11244117##7##]. Also, uORFs have been found in dicot plant genes that include <italic>AtB2/AtbZIP11 </italic>[##REF##15208401##32##], <italic>ABI3 </italic>[##REF##15159632##33##], and <italic>CpbZIP2 </italic>[##REF##16897481##34##]; and auxin responsive factor genes <italic>ETT </italic>and <italic>MP </italic>[##REF##16227452##35##]. These characterised uORFs (&lt;0.3%) in plants are much lower than the estimated number of genes that contain uORFs, which can vary from 11% [##REF##11164040##36##] to 60% [##REF##17663791##12##].</p> <p>One strategy for identifying functional uORFs in plants is to use a comparative approach [##REF##17663791##12##,##REF##16438715##13##,##REF##17286865##37##]. There are extensive assembled EST datasets for five important cereal crops and Arabidopsis. The cereals include rice (<italic>Oryza sativa </italic>L.), wheat (<italic>Triticum aestivum</italic>), barley (<italic>Hordeum vulgare</italic>), maize (<italic>Zea mays</italic>), and sorghum (<italic>Sorghum bicolor</italic>). Rice is the best characterised of these cereals with a sequenced genome [##REF##11935017##38##] and a cDNA database containing 32,000 clones that were enriched for 5' full-length sequences [##REF##12869764##39##]. Cereals such as wheat are unlikely to be fully sequenced in the near future because of their large genome size. Wheat has a hexaploid genome of 16000 Mb [##REF##12830387##40##] that is 37 times that of rice (430 Mb), and 5.5 times (2900 Mb) the size of the human genome [##REF##11181995##41##]. A comparative approach is likely to identify sequence dependent uORFs [##REF##11707416##16##] where the encoded peptide of the uORF is involved in regulation of gene expression.</p> <p>In this study, we used comparative R-nomics to identify conserved uORF motifs in cereals and Arabidopsis. We constructed a bioinformatics pipeline called uORFSCAN that performs a comparative analysis on the important agronomic crops rice, wheat, barley, maize, and sorghum; and the well studied dicot plant Arabidopsis. To account for the variable quality of assembled EST data, we have used orthologous sequence clustering, iterative sequence analysis, and manual curation. Our comparative method is easily transferable to uORF identification in other species.</p>
<title>Methods</title> <title>Data material</title> <p>KOME (Knowledge-based Oryza Molecular biological Encyclopedia) full-length rice cDNA sequences were obtained from <ext-link ext-link-type="ftp" xlink:href="ftp://cdna01.dna.affrc.go.jp/pub/data/CURRENT/INE_FULL_SEQUENCE_DB.zip"/>. This file is dated Tuesday, 24 January 2006, and contains 32,127 full-length cDNA clones (originally 28,469). The TIGR plant gene indices database <ext-link ext-link-type="ftp" xlink:href="ftp://ftp.tigr.org/pub/data/tgi/"/> was used to obtain tentative contigs (TCs) from wheat (release 10.0, Jan 05, 580155 ESTs, 44954 TCs), barley (release 9.0, Sept 04, 370546 ESTs, 23176 TCs), maize (release 17.0, Nov 06, 695811 ESTs, 56687 TCs), and sorghum (release 8.0, Nov 05, 187282 ESTs, 20029 TCs). Data cleaning was performed on the TIGR dataset to select for sequences that are designated as tentative contigs (identifiers prefixed with "TC"), thereby excluding all singletons. All data files were imported and managed using Microsoft Access 2003. We also re-ran the analysis using the TIGR Plant Transcript Assemblies (last updated on October 17<sup>th</sup>, 2006) for wheat (840871 ESTs), barley (456410 ESTs), maize (1084701 ESTs), and sorghum (203575 ESTs) on the uORFSCAN pipeline, but did not find any additional conserved upstream open reading frames (uORFs).</p> <title>Orthologue searches</title> <p>The reciprocal best hit method (rbh) was adapted to account for alternative splice forms that are present in the KOME dataset that would otherwise give many false negatives. The problem with alternative splice forms is that they will never have the highest score in the reverse BLAST because the presence of a longer alternative splice form will always be listed higher on the hit list due to the way BLAST [##REF##2231712##42##] ranks hits (according to score and e-value). To account for alternative splice forms, we examined not only the top hit but also similar hits (percent identity to top hit: Δ -5%, similar length to top hit: +/- 20%) for symmetry with the top hit in the forward blast. If there is symmetry between the forward and reverse blasts then we considered the reciprocal pair to be orthologous. General parameters for similarity searches were: tblastx program, expect threshold value at 1.0e-50, scoring parameters, BLOSUM62 matrix; gap costs (existence, 11; extension 1), and filter and masking, off. Only sequence alignments with at least 70% sequence coverage were considered further. Similarity searches were performed at The South Australian Partnership for Advanced Computing (SAPAC) <ext-link ext-link-type="uri" xlink:href="http://www.sapac.edu.au/"/>.</p> <title>Verification of main ORF</title> <p>The rice cDNA sequences containing conserved uORFs were used in a blastn search against NCBI Non-redundant database to identify uORFs predicted from ribosomal RNA genes, chloroplastic genes, and mitochondrial genes. These genes do not represent coding genes derived from the nuclear genome, and therefore have been removed from this study. Also, the main open reading frames, predicted by uORFSCAN were used to search (blastn) the coding sequence (CDS) annotations from TIGR rice pseudomolecules database <ext-link ext-link-type="uri" xlink:href="http://www.tigr.org/tdb/e2k1/osa1/data_download.shtml"/>. Alignments not starting from the beginning of the CDS were regarded as suspicious. As additional verification, the rice main open reading frame predictions were also compared with protein data from The UniProt Knowledgebase (UniProtKB) <ext-link ext-link-type="uri" xlink:href="http://www.ebi.uniprot.org/database/download.shtml"/>. Translations of the rice cDNA sequences in the same frame as the predicted main open reading frame, starting from the 5'-untranslated region to the end of the main open reading frame, were used to search (blastp) against UniProtKB. Aligments not beginning from the start of the protein sequence were discarded if they also did not have TIGR CDS support.</p> <title>Statistical analysis of codon usage</title> <p>The <italic>p-values </italic>were calculated according to the following formulas:</p> <p>The probability to observe the number of times each codon was present in the uORFs (n<sub>obs</sub>) that was less than or equal to the expected (n<sub>av</sub>) by chance alone is:</p> <p></p> <p>The probability to observe the number of times each codon was present in the uORFs (n<sub>obs</sub>) that was greater than or equal to the expected (n<sub>av</sub>) by chance alone is:</p> <p></p> <p>Where,</p> <p>n<sub>obs </sub>= The observed number of times a codon was present in the uORFs.</p> <p>n<sub>av </sub>= The average number of times a codon was present in the uORFs based on the frequency of this codon in the mORF and the sample size (the observed number of codons for the set of codons for an amino in the uORFs).</p>
<title>Results</title> <title>The uORFSCAN pipeline for discovering uORFs</title> <p>The uORFSCAN pipeline used rice full-length cDNAs [##REF##12869764##39##] and wheat, barley, maize, and sorghum assembled EST data for comparative analysis (Figure ##FIG##0##1##). In the first step of the pipeline, we identified rice genes that had orthologues in wheat, barley, maize, and sorghum. The use of orthologous sequences allowed us to more accurately predict the main coding region and define the 5'-UTR that is necessary to identify conserved uORFs.</p> <p>The reciprocal best hit (rbh) method was used to find true orthologues by a process of eliminating paralogues [##REF##9790834##43##,##REF##9381173##44##]. The principle of rbh is that a pair of sequences are orthologues if they are each others best hit. We modified the rbh method to find orthologues of rice genes in wheat, barley, maize, and sorghum such that it allowed us to keep alternative splice forms while at the same time eliminating paralogues. Alternative splice forms of a gene were distinguished by changes in gene length while still maintaining high sequence identity. We found that the modified reciprocal best hit method eliminated 70–75% of paralogue sequences. For example, in the one directional BLAST against the barley assembled EST database 19,655 sequences were identified, however this number was reduced to 5,115 (26%) sequences when the reciprocal best hit method was used (Figure ##FIG##0##1##, Step 1).</p> <p>Only 1723 of the rice genes had conserved orthologues in the other four cereals (wheat, barley, maize, and sorghum), most likely because none of the assembled EST datasets contained the entire transcriptome. To account for missing or erroneous sequences, we grouped the orthologues into three datasets for 5'-UTR analysis (Figure ##FIG##0##1##, Step 2). The datasets included rice genes that had orthologues in four other cereals (5 out of 5 dataset), in three other cereals (4 out of 5 dataset), and in two other cereals (3 out of 5 dataset).</p> <p>In Figure ##FIG##0##1## (Step 3), we developed a program called uORFSCAN (see Additional file ##SUPPL##8##9##)  to find conserved uORFs. uORFSCAN takes as input a FASTA file containing the rice cDNA sequence and its orthologues, and identifies for each of these sequences all the possible open reading frames (ORFs). In the first iteration, the longest conserved ORF was designated as the main coding region. However, the longest ORF is not always the main coding region when there are other ORFs of similar length. Therefore, a comparative approach was used to identify the main coding region (Figure ##FIG##0##1##, Step 3.1). This involved finding the longest ORF that was present in all orthologous sequences, and then iteratively reducing the number of orthologous sequences, one at a time, to determine if a longer conserved set of ORFs could be found, and finally terminating when there was no improvement. The longest ORF in at least three out of five cereals was considered the main coding region. In Figure ##FIG##0##1## (Step 3.2), uORFSCAN attempts to align rice uORFs with similar length orthologous uORFs (+/- 5%) at the protein level using ClustalW ([##REF##7984417##45##], see Additional file ##SUPPL##7##8##) . Finally, uORFSCAN analysed each alignment file to determine the average conservation of the uORFs, and grouped the alignments based on the number of conserved orthologous uORFs found. For example, using the 4 out of 5 dataset generated the 4 out of 4 and the 3 out of 4 datasets (Figure ##FIG##0##1##). We report uORFs from orthologous genes that shared sequence similarity because of our interest in finding functional uORFs.</p> <p>The final step (Figure ##FIG##0##1##, Step 4) was manual curation to verify the predicted rice main coding region of each gene by comparing it with the genome annotation and other protein data. This was necessary, as uORFSCAN is expected to be sensitive to inaccurate (e.g., frame-shifts) and/or incomplete sequence data. For example, rice full length sequences can be incomplete because of failure of the 5' capping method [##REF##12869764##39##]. If the coding region is truncated, this can result in an internal methionine selected as the start codon and therefore the derived 5'-UTR is actually coding sequence, which is often highly conserved and can lead to false positive predictions.</p> <title>Conserved upstream open reading frames appear to be rare</title> <p>The uORFSCAN pipeline identified nine cDNAs containing uORFs that were conserved in the five out of five (5/5) dataset containing orthologous sequences in all 5 cereals (Table ##TAB##0##1##). Three of these cDNAs encoded multiple uORFs, one of the cDNAs being AdoMetDC, which has previously been reported to contain two uORFs [##REF##12205086##46##]. We manually curated all nine cDNAs and showed that they were all reliable based on our validation criteria (Table ##TAB##1##2##), which included the removal of the uORFs predicted from ribosomal rRNA genes (data not shown). The cDNAs included the multiple uORFs in S-adenosylmethionine decarboxylase cDNA [##REF##12205086##46##], alkaline phytoceramidase cDNA, calcineurin B-like (CBL)-interacting protein kinase cDNA; and a single conserved uORF in a cDNA encoding an oxidoreductase protein, ribosomal protein S6 kinase, trehalose-6-phosphate phosphatase, ubiquitin-fold protein, F9L1.29 protein, and an ankyrin-3 protein.</p> <p>To account for variable quality in assembled EST data, we also looked for cases where the uORFs (4/5, 3/5, 4/4, 3/4, and 3/3 result set) were conserved in only four out of five (4/5) orthologues and three out of five orthologues (see Additional files ##SUPPL##0##1##, ##SUPPL##1##2##, ##SUPPL##2##3##, ##SUPPL##3##4##, ##SUPPL##4##5##; Figure ##FIG##0##1##, Step 3.2). In brief, the 4/5 result set contains 16 rice genes with a total of 20 conserved uORFs in orthologous cereal genes, the 3/5 result set contains 44 rice genes with a total of 79 conserved uORFs in orthologous cereal genes, the 4/4 result set contains 16 rice genes with a total of 23 conserved uORFs in orthologous genes, the 3/4 result set contains 113 rice genes with a total of 129 conserved uORFs in orthologous genes, and finally the 3/5 result set contains 65 genes with a total of 93 conserved uORFs in orthologous genes.</p> <p>In order to identify sequence dependent uORFs, we extended our search for cereal uORFs that might also be conserved in the dicot Arabidopsis by using the rice cDNAs that contained conserved uORFs in four other cereals (5/5 result set) and the Arabidopsis Tair 7 cDNA dataset (see Methods). The uORFSCAN pipeline identified 13 rice cDNAs containing uORFs that were conserved in Arabidopsis (Table ##TAB##2##3##). Four of these cDNAs encoded multiple uORFs. Of the 13 cDNAs with uORFs, only 11 were verified as reliable based on manual curation (Table ##TAB##3##4##) that removed the uORFs predicted from a cDNA encoding a helicase. Manual curation of the helicase cDNA revealed that the genome and protein annotation for the coding region extended further upstream than predicted by uORFSCAN, highlighting the limitations of using assembled EST data where frame-shift errors was the likely reason for the false positive prediction. The reliable predictions included the multiple uORFs found in a cDNA encoding ww domain containing protein, trehalose-6-phosphate phosphatase, GAMYB-binding protein, and ankyrin-3. The latter three cDNAs contained a combination of uORFs that were conserved between the cereals (rice and at least two other cereals) and Arabidopsis, and uORFs conserved between rice and Arabidopsis (Table ##TAB##2##3##). uORFSCAN also identified seven rice cDNAs containing a single uORF that were conserved in Arabidopsis and in almost all cases (except cDNA encoding an auxilin-like protein) the cereals as well (Table ##TAB##2##3##). They included the uORFs found in a cDNA encoding phosphatase 2a protein, homeodomain containing protein, S-Adenosylmethionine decarboxylase, auxilin-like protein, CBL-interacting protein kinase, protein kinase ATN1, and a hypothetical protein.</p> <title>Position and occupation of uORFs in 5'-UTRs</title> <p>Studies have shown that the position of an uORF within its 5'-UTR, which determines the pre-orf and intercistronic distances, can have profound effects on its function [##REF##12890013##18##,##REF##10357825##22##]. When we examined the positions of cereal uORFs within their 5'-UTRs we found that there was no positional preference with the exception that they were not positioned too closely to the start of their individual 5'-UTR and coding region (Figure ##FIG##1##2##). For example, all of the uORFs conserved in five orthologous cereals (5/5 result set) and in Arabidopsis were at least positioned 20 nucleotides from the start of their 5'-UTR, which is thought to be the minimum number of nucleotides required for a functional uORF [##REF##12890013##18##]. We noticed that the intercistronic distances for these uORFs were generally shorter than the pre-orf distance. We also found seven uORFs, which included the functional small AdoMetDC uORF, that occupied greater than 20% of their individual 5'-UTR.</p> <title>The length distribution of uORFs</title> <p>Since earlier reports showed that plant uORFs can vary in length from 6 to 156 nucleotides [##REF##11244117##7##,##REF##11139406##9##,##REF##11855732##29##, ####REF##8439744##30##, ##REF##9761799##31####9761799##31##,##REF##12205086##46##], we examined the length distribution of the cereal uORFs. There are two peaks in the distribution that were found between 1 to 40 nucleotides, and 81 to 120 nucleotides (Figure ##FIG##2##3##). The uORFs found in the first peak are tiny with 9 (out of 14) uORFs having a length of nine nucleotides. Some of these tiny uORFs could be artefactual as a result of point mutations that insert an in-frame start and/or stop codon in the 5'-UTR. If these artefactual uORFs were removed then the uORF length distribution would move towards a normal distribution. Seventy six percent of the uORFs in the length distribution are shorter than 100 nucleotides, and 48% are shorter than 40 nucleotides. The shortest conserved uORF found in four independent cDNAs was nine nucleotides, even though the cut-off length used by uORFSCAN to identify uORFs was six nucleotides (a start and a stop codon). One of the nine nucleotide uORFs was the 5' tiny uORF found in the S-adenosylmethionine decarboxylase cDNA [##REF##11139406##9##], and three new uORFs, two found in a cDNA encoding alkaline phytoceramidase, and one in a cDNA encoding oxidoreductase, (Table ##TAB##0##1##). Two long conserved uORFs (&gt;181 nucleotides) were found in cDNAs encoding protein kinases that included one uORF found in a cDNA encoding a CBL-interacting protein kinase and another uORF found in a cDNA encoding a ribosomal protein S6 kinase.</p> <title>Sequence conservation in uORFs</title> <p>The level of amino acid sequence conservation in cereal uORFs was generally high, as expected, based on our approach of reporting similar length orthologous uORFs that shared sequence similarity. For example, in the 5 out of 5 dataset the median value is 50% sequence similarity. When the two main datasets were included (uORFs conserved in all five cereals and uORFs conserved between rice and Arabidopsis), the median value is 36% sequence similarity. The uORFs conserved between the cereals (rice and at least two others) and Arabidopsis (median value of 36% sequence similarity) generally had a higher amino acid sequence similarity than those uORFs conserved between rice and Arabidopsis (median value of 28% sequence similarity). Given that the uORFs from orthologous genes were selected to be within a given length interval for alignment purposes, the high amino acid sequence similarity may suggest that these uORFs have a functional role (e.g., ribosomal stalling) that is mediated by the encoded uORF peptide.</p> <title>Start codon context and codon usage of uORFs</title> <p>The presence of uORFs does not mean that they will be translated. The sequence context of some plant uORFs has been shown to be sub-optimal for efficient initiation [##REF##9426620##47##,##REF##14990743##48##]. We therefore examined the sequence context of our cereal uORF AUG codons to see if there was any sequence conservation that may aid in their ribosome initiation. We found that there were no informative positions in the uORF consensus sequence context (see Additional file ##SUPPL##5##6## – Tran_FigureS1.eps) based on the observed number of positions that showed sequence conservation was not greater than expected by chance alone. However when the context of the AUGs demarcating the conserved uORFs were compared with the context of the AUG at the main ORF it was evident that the main ORF generally had a better sequence context denoted by a purine in the -3 position and a guanine in the +4 position (Table ##TAB##4##5##).</p> <p>Recent work showed that ribosome stalling could occur at rare codons [##REF##15694341##49##, ####REF##12799445##50##, ##REF##16703415##51##, ##REF##17017124##52##, ##REF##17683808##53####17683808##53##]. We therefore examined the codon usage of the identified uORFs to determine if they contained an increased number of rare codons. We showed that the frequencies of some codons had a <italic>p-Value </italic>less than &lt;<italic>0.05 </italic>in the rice uORF codon usage compared to the rice main coding region based on a significant deviation of observed from expected numbers of uORF codons (Equations 1 and 2); however, the number of codons that had these <italic>p-Values </italic>were not greater than expected by chance (see Additional file ##SUPPL##6##7## – Tran_FigureS2.eps).</p>
<title>Discussion and conclusion</title> <title>Conserved uORFs are rare</title> <p>This study provides a method to identify conserved uORFs from large assembled EST datasets. We developed a pipeline that used a modified reciprocal best hit method to identify putative orthologous sequences that were then analysed by a comparative R-nomics program called uORFSCAN to find conserved uORFs. We showed that this pipeline was successful in identifying 29 rice uORFs that are conserved at the amino acid level (median value of 36% sequence similarity) in wheat, barley, maize, sorghum, and in some cases (33%) Arabidopsis.</p> <p>The number of conserved uORFs that share sequence similarity in the transcriptome of cereals appears to be low. This is consistent with reports of conserved uORFs in distantly related plants (i.e., rice and Arabidopsis) [##REF##17663791##12##] and in <italic>Drosophila melanogaster </italic>[##REF##18237443##14##]. One explanation is that genes controlled at the level of translation by uORFs have low levels of transcription [##REF##15821146##27##] and therefore are under- represented in cDNA and assembled EST databases. Another explanation for the low numbers of conserved cereal uORFs is that the uORFs have evolved in both length and sequence such that they no longer share sequence similarity across minor taxonomic groups (i.e., within the cereals) (Table ##TAB##1##2## and ##TAB##3##4##). Furthermore, if the codon usage of cereal uORFs rather than the uORF-encoded peptide were a major controlling mechanism then amino acid sequence may not be conserved.</p> <title>Cereal uORFs conserved in Arabidopsis</title> <p>It has been shown that the amino acid sequence of uORFs in monocot and dicot plants can be similar [##REF##12205086##46##]. Sequence similarity was observed at the amino acid level across the major taxonomic groups (e.g., Arabidopsis and rice) (Table ##TAB##2##3##). We identified 11 rice genes that contained uORFs conserved in Arabidopsis, of which nine were also conserved in additional cereals (at least two others). For example, a rice cDNA encoding Ankyrin-3 contains an uORF that is conserved in the cereals and Arabidopsis, but it contains a nested uORF that appears to be conserved only in rice and Arabidopsis. Therefore, it is likely that after the split between the two major groups of angiosperms (monocots and dicots), the rice gene has gained an additional in-frame and internal start codon, that is not present in the other cereals, making a nested uORF that is shorter by 33 nucleotides. It would be of interest to determine if the nested uORF is still functional.</p> <p>Conservation of uORF sequence within the cereals might simply reflect a relatively recent ancestor, rather than conservation of function, therefore it is difficult to predict whether these uORFs are likely to be sequence dependent or sequence independent uORFs [##REF##12890013##18##,##REF##12117416##19##]. However, uORFs that are conserved across both monocots and dicots suggest that these uORFs have a role in a sequence dependent manner. Indeed, six rice uORFs (out of 15, excluding nested uORFs, Table ##TAB##2##3##) that were conserved in Arabidopsis had a notable amino acid composition that was rich in serine or arginine (at least 20%). It has been suggested that uORF peptides that are rich in serine could either promote or inhibit ribosomal stalling through their phosphorylation [##REF##17663791##12##,##REF##18160716##54##], while arginine rich motifs can be involved in RNA binding [##REF##16314457##55##]. Interestingly, of these six rice uORFs two (AK101100 and AK067412) are found in genes involved in phosphorylation, a function that appears to be over-represented in this dataset (Table ##TAB##2##3##). We hypothesize that the main protein of these genes could have dual functions, the primary function is as a <italic>trans</italic>-acting factor in an unknown signalling cascade, and a secondary function as a regulator of mORF expression whereby the mORF protein phosphorylates the serine-rich uORF peptides, resulting in a conformational change that allows the uORF peptides to bind and stall ribosomes [##REF##11707416##16##].</p> <p>There are uORFs previously identified in Arabidopsis that were not identified in this study. For example, the Arabidopsis auxin response factor (ARF) genes [##REF##16227452##35##]<italic>ETTIN (ETT) </italic>and <italic>MONOPTEROS (MP) </italic>contain uORFs and while orthologues of these genes were found in the rice, sorghum and wheat assembled EST datasets, the uORFs showed no sequence similarity (by ClustalW) and were of different lengths (data not shown). Similarly, uORFs found in Arabidopsis genes <italic>AtMHX </italic>and <italic>AtNMT1 </italic>encoding encoding a tonoplast transporter [##REF##15710632##56##] and a phosphoethanolamine N-methyltransferase [##REF##16960350##57##] respectively were not identified because the uORFs were not conserved in rice and at least two other cereals. Finally, the gene containing the uORF in Arabidopsis <italic>sac51 </italic>encoding a bHLH-type transcription factor [##REF##16936072##58##] could not be identified in our rice dataset as we could not identify a clear orthologue. Therefore, it will be of interest to monitor new rice full-length cDNAs and high quality sequences for cereals as they become available to see if more conserved uORFs can be found.</p> <p>Recently, a pair-wise comparative approach was used to identify conserved uORFs within homology groups that also included paralogs and ohnologs (homologous genes arising by whole-genome duplication) using rice and Arabidopsis full-length cDNAs [##REF##17663791##12##]. Compared to the 11 genes we had identified Hayden and Jorgensen [##REF##17663791##12##] reported that 19 genes contained conserved uORFs between rice and Arabidopsis. Interestingly only four genes (S-Adenosylmethionine decarboxylase, Trehalose-6-phosphate phosphatase, Auxilin-like protein, and Ankyrin-3) were in common highlighting the benefits of complementary search methods. For example, we used the modified reciprocal best hit method to find putative orthologues. It is likely that some of the homologue groups identified by Hayden and Jorgensen [##REF##17663791##12##] may not be true orthologues. For example, homologue group 12 identified by Hayden and Jorgensen [##REF##17663791##12##] were not reciprocal best hit pairs according to our analysis, and therefore are likely to be paralogues. Our approach is deliberately conservative, eliminating paralogues, to maximise the finding of all conserved uORFs independent of their length.</p> <p>One possible criticism of our approach is that we have included uORFs as short as 9 nt. However, there are two independent reports that showed that the tiny uORF of SAMDC is functional [##REF##15821146##27##,##REF##16176926##59##], although there is controversy regarding the type of effect and conditions under which the tiny uORF of SAMDC exerts its affect on downstream translation. Therefore, there is insufficient data to conclude one way or the other, and as such we have elected to be conservative. This has allowed us to find several genes (e.g. protein phosphatase 2a, a protein containing a ww domain, and GAMYB-binding protein) that were not found by Hayden and Jorgensen's 'uORF-Finder' program [##REF##17663791##12##] because it only detected conserved uORFs greater than 63 codons.</p> <title>Better quality assembled EST data is needed</title> <p>One unavoidable limitation of using incomplete assembled EST data for orthology determination is that orthologues could be falsely assigned in situations where sequences have multiple protein domains. This will increase the number of putative orthologues identified prior to the prediction of uORFs, which is not necessarily harmful as these predictions are manually curated. However, to minimise this problem, we used a sequence coverage cutoff of at least 70% of any of the protein sequences in the alignment (see Methods). We also grouped the orthologues into several datasets representing the number of orthologues that could be found for each gene. For example, the datasets included rice genes that had orthologues in four other cereals (5 out of 5 dataset), in three other cereals (4 out of 5 dataset), and in two other cereals (3 out of 5 dataset). This grouping of orthologues will also help minimise the effects of missing, incomplete, or erroneous assembled EST data.</p> <p>There are reports of conserved uORFs in monocots and dicots that share high sequence similarity that were not found by our pipeline, due to either lack of sequence conservation or due to limitations in the assembled ESTs currently available. For example, the uORF found in the basic region leucine zipper (bZIP)-type transcription factor <italic>AtB2/AtbZIP11 </italic>was found to be conserved in rice and barley [##REF##15208401##32##], but not in the other cereals included in this study because the sequences are not represented in the other datasets. Current limitations include incomplete data (i.e., not all sequences are represented) and poor quality sequence data, leading to frame-shifts and incorrect prediction of uORFs. Therefore, it is possible to obtain higher numbers of conserved uORFs if the cluster size was relaxed to two out of five, but this approach would reduce the power of comparative R-nomics, and would require significant manual curation.</p> <title>Sequence dependent and independent uORF</title> <p>The cereal uORFs identified here are likely to encode bioactive peptides as selection has occurred at the peptide level. Those cereal uORFs that showed sequence conservation at the amino acid level with Arabidopsis are likely to be classified as sequence-dependent, as the encoded uORF peptide has remained conserved across the angiosperms, suggesting the peptide is directly involved in translational control [##REF##11139406##9##] or has some other biological activity [##REF##17663791##12##, ####REF##16438715##13##, ##REF##18237443##14##, ##REF##15777708##15####15777708##15##]. Some identified uORFs were conserved only within the cereals, indicating a relative recent origin or selective loss of the uORFs in Arabidopsis. We cannot rule out the possibility that some conserved cereal uORFs could also act in a sequence-independent manner, as a recent paper reported a conserved uORF in human and mouse ribosomal protein S6 kinase genes (the same finding by our analysis in cereals, Table ##TAB##0##1##), and suggested that the uORF translational control of the main ORF was through reinitiation [##REF##18358843##26##]. Experiments are needed to confirm the biological activity of the uORF in ribosomal protein S6 kinase gene.</p> <p>The sequence context surrounding an uORF (ignoring secondary structure) does not appear to play a major role in its recognition and initiation of translation at an uORF. We hypothesize that this sub-optimal uORF sequence context (compared to optimal Kozak consensus [##REF##9426620##47##] sequence for the main coding region) would allow for leaky scanning [##REF##14990743##48##,##REF##15741183##60##] of the uORF, and preferential initiation at the downstream main coding region. An optimal uORF sequence context would provide rigid control in the translational regulation of the main coding region, as initiation would predominantly start at the uORF resulting in reduced availability of initiation factors, such as eIF2, for re-initiation at the downstream main open reading frame.</p> <p>Sequence-independent uORFs allow for low-level translation of the downstream main coding region [##REF##10446211##61##]. Low-level translation is possible, as sequence-independent uORFs do not cause ribosomal stalling as seen in sequence-dependent uORFs. The regulatory mechanism of the sequence-independent uORF involves other factors (uORF recognition, length, stop codon environment, and the downstream intercistronic sequence) that influence reinitiation efficiency [##REF##12890013##18##,##REF##12117416##19##], and more recently leaky scanning [##REF##14990743##48##], to regulate downstream translation. We analysed the codon usage of conserved uORFs and found no preferential usage of rare codons in the uORFs. Therefore, it is unlikely that the uORF codon usage could contribute to low-level translation as seen for certain rare codons in <italic>Xenopus laevis </italic>[##REF##12799445##50##] and <italic>Eschericia coli </italic>[##REF##16703415##51##] that can reduce translation.</p> <p>In conclusion, this study showed that the uORFSCAN pipeline is a useful tool for identifying conserved uORFs in closely related species. This pipeline has allowed us to identify 29 conserved uORFs in cereals. Possibly more conserved uORFs will be identified once the cDNA and assembled EST datasets become more comprehensive. These conserved rice uORFs will be useful for future functional analyses that should provide some perspective into downstream translational regulation by uORFs.</p>
<title>Discussion and conclusion</title> <title>Conserved uORFs are rare</title> <p>This study provides a method to identify conserved uORFs from large assembled EST datasets. We developed a pipeline that used a modified reciprocal best hit method to identify putative orthologous sequences that were then analysed by a comparative R-nomics program called uORFSCAN to find conserved uORFs. We showed that this pipeline was successful in identifying 29 rice uORFs that are conserved at the amino acid level (median value of 36% sequence similarity) in wheat, barley, maize, sorghum, and in some cases (33%) Arabidopsis.</p> <p>The number of conserved uORFs that share sequence similarity in the transcriptome of cereals appears to be low. This is consistent with reports of conserved uORFs in distantly related plants (i.e., rice and Arabidopsis) [##REF##17663791##12##] and in <italic>Drosophila melanogaster </italic>[##REF##18237443##14##]. One explanation is that genes controlled at the level of translation by uORFs have low levels of transcription [##REF##15821146##27##] and therefore are under- represented in cDNA and assembled EST databases. Another explanation for the low numbers of conserved cereal uORFs is that the uORFs have evolved in both length and sequence such that they no longer share sequence similarity across minor taxonomic groups (i.e., within the cereals) (Table ##TAB##1##2## and ##TAB##3##4##). Furthermore, if the codon usage of cereal uORFs rather than the uORF-encoded peptide were a major controlling mechanism then amino acid sequence may not be conserved.</p> <title>Cereal uORFs conserved in Arabidopsis</title> <p>It has been shown that the amino acid sequence of uORFs in monocot and dicot plants can be similar [##REF##12205086##46##]. Sequence similarity was observed at the amino acid level across the major taxonomic groups (e.g., Arabidopsis and rice) (Table ##TAB##2##3##). We identified 11 rice genes that contained uORFs conserved in Arabidopsis, of which nine were also conserved in additional cereals (at least two others). For example, a rice cDNA encoding Ankyrin-3 contains an uORF that is conserved in the cereals and Arabidopsis, but it contains a nested uORF that appears to be conserved only in rice and Arabidopsis. Therefore, it is likely that after the split between the two major groups of angiosperms (monocots and dicots), the rice gene has gained an additional in-frame and internal start codon, that is not present in the other cereals, making a nested uORF that is shorter by 33 nucleotides. It would be of interest to determine if the nested uORF is still functional.</p> <p>Conservation of uORF sequence within the cereals might simply reflect a relatively recent ancestor, rather than conservation of function, therefore it is difficult to predict whether these uORFs are likely to be sequence dependent or sequence independent uORFs [##REF##12890013##18##,##REF##12117416##19##]. However, uORFs that are conserved across both monocots and dicots suggest that these uORFs have a role in a sequence dependent manner. Indeed, six rice uORFs (out of 15, excluding nested uORFs, Table ##TAB##2##3##) that were conserved in Arabidopsis had a notable amino acid composition that was rich in serine or arginine (at least 20%). It has been suggested that uORF peptides that are rich in serine could either promote or inhibit ribosomal stalling through their phosphorylation [##REF##17663791##12##,##REF##18160716##54##], while arginine rich motifs can be involved in RNA binding [##REF##16314457##55##]. Interestingly, of these six rice uORFs two (AK101100 and AK067412) are found in genes involved in phosphorylation, a function that appears to be over-represented in this dataset (Table ##TAB##2##3##). We hypothesize that the main protein of these genes could have dual functions, the primary function is as a <italic>trans</italic>-acting factor in an unknown signalling cascade, and a secondary function as a regulator of mORF expression whereby the mORF protein phosphorylates the serine-rich uORF peptides, resulting in a conformational change that allows the uORF peptides to bind and stall ribosomes [##REF##11707416##16##].</p> <p>There are uORFs previously identified in Arabidopsis that were not identified in this study. For example, the Arabidopsis auxin response factor (ARF) genes [##REF##16227452##35##]<italic>ETTIN (ETT) </italic>and <italic>MONOPTEROS (MP) </italic>contain uORFs and while orthologues of these genes were found in the rice, sorghum and wheat assembled EST datasets, the uORFs showed no sequence similarity (by ClustalW) and were of different lengths (data not shown). Similarly, uORFs found in Arabidopsis genes <italic>AtMHX </italic>and <italic>AtNMT1 </italic>encoding encoding a tonoplast transporter [##REF##15710632##56##] and a phosphoethanolamine N-methyltransferase [##REF##16960350##57##] respectively were not identified because the uORFs were not conserved in rice and at least two other cereals. Finally, the gene containing the uORF in Arabidopsis <italic>sac51 </italic>encoding a bHLH-type transcription factor [##REF##16936072##58##] could not be identified in our rice dataset as we could not identify a clear orthologue. Therefore, it will be of interest to monitor new rice full-length cDNAs and high quality sequences for cereals as they become available to see if more conserved uORFs can be found.</p> <p>Recently, a pair-wise comparative approach was used to identify conserved uORFs within homology groups that also included paralogs and ohnologs (homologous genes arising by whole-genome duplication) using rice and Arabidopsis full-length cDNAs [##REF##17663791##12##]. Compared to the 11 genes we had identified Hayden and Jorgensen [##REF##17663791##12##] reported that 19 genes contained conserved uORFs between rice and Arabidopsis. Interestingly only four genes (S-Adenosylmethionine decarboxylase, Trehalose-6-phosphate phosphatase, Auxilin-like protein, and Ankyrin-3) were in common highlighting the benefits of complementary search methods. For example, we used the modified reciprocal best hit method to find putative orthologues. It is likely that some of the homologue groups identified by Hayden and Jorgensen [##REF##17663791##12##] may not be true orthologues. For example, homologue group 12 identified by Hayden and Jorgensen [##REF##17663791##12##] were not reciprocal best hit pairs according to our analysis, and therefore are likely to be paralogues. Our approach is deliberately conservative, eliminating paralogues, to maximise the finding of all conserved uORFs independent of their length.</p> <p>One possible criticism of our approach is that we have included uORFs as short as 9 nt. However, there are two independent reports that showed that the tiny uORF of SAMDC is functional [##REF##15821146##27##,##REF##16176926##59##], although there is controversy regarding the type of effect and conditions under which the tiny uORF of SAMDC exerts its affect on downstream translation. Therefore, there is insufficient data to conclude one way or the other, and as such we have elected to be conservative. This has allowed us to find several genes (e.g. protein phosphatase 2a, a protein containing a ww domain, and GAMYB-binding protein) that were not found by Hayden and Jorgensen's 'uORF-Finder' program [##REF##17663791##12##] because it only detected conserved uORFs greater than 63 codons.</p> <title>Better quality assembled EST data is needed</title> <p>One unavoidable limitation of using incomplete assembled EST data for orthology determination is that orthologues could be falsely assigned in situations where sequences have multiple protein domains. This will increase the number of putative orthologues identified prior to the prediction of uORFs, which is not necessarily harmful as these predictions are manually curated. However, to minimise this problem, we used a sequence coverage cutoff of at least 70% of any of the protein sequences in the alignment (see Methods). We also grouped the orthologues into several datasets representing the number of orthologues that could be found for each gene. For example, the datasets included rice genes that had orthologues in four other cereals (5 out of 5 dataset), in three other cereals (4 out of 5 dataset), and in two other cereals (3 out of 5 dataset). This grouping of orthologues will also help minimise the effects of missing, incomplete, or erroneous assembled EST data.</p> <p>There are reports of conserved uORFs in monocots and dicots that share high sequence similarity that were not found by our pipeline, due to either lack of sequence conservation or due to limitations in the assembled ESTs currently available. For example, the uORF found in the basic region leucine zipper (bZIP)-type transcription factor <italic>AtB2/AtbZIP11 </italic>was found to be conserved in rice and barley [##REF##15208401##32##], but not in the other cereals included in this study because the sequences are not represented in the other datasets. Current limitations include incomplete data (i.e., not all sequences are represented) and poor quality sequence data, leading to frame-shifts and incorrect prediction of uORFs. Therefore, it is possible to obtain higher numbers of conserved uORFs if the cluster size was relaxed to two out of five, but this approach would reduce the power of comparative R-nomics, and would require significant manual curation.</p> <title>Sequence dependent and independent uORF</title> <p>The cereal uORFs identified here are likely to encode bioactive peptides as selection has occurred at the peptide level. Those cereal uORFs that showed sequence conservation at the amino acid level with Arabidopsis are likely to be classified as sequence-dependent, as the encoded uORF peptide has remained conserved across the angiosperms, suggesting the peptide is directly involved in translational control [##REF##11139406##9##] or has some other biological activity [##REF##17663791##12##, ####REF##16438715##13##, ##REF##18237443##14##, ##REF##15777708##15####15777708##15##]. Some identified uORFs were conserved only within the cereals, indicating a relative recent origin or selective loss of the uORFs in Arabidopsis. We cannot rule out the possibility that some conserved cereal uORFs could also act in a sequence-independent manner, as a recent paper reported a conserved uORF in human and mouse ribosomal protein S6 kinase genes (the same finding by our analysis in cereals, Table ##TAB##0##1##), and suggested that the uORF translational control of the main ORF was through reinitiation [##REF##18358843##26##]. Experiments are needed to confirm the biological activity of the uORF in ribosomal protein S6 kinase gene.</p> <p>The sequence context surrounding an uORF (ignoring secondary structure) does not appear to play a major role in its recognition and initiation of translation at an uORF. We hypothesize that this sub-optimal uORF sequence context (compared to optimal Kozak consensus [##REF##9426620##47##] sequence for the main coding region) would allow for leaky scanning [##REF##14990743##48##,##REF##15741183##60##] of the uORF, and preferential initiation at the downstream main coding region. An optimal uORF sequence context would provide rigid control in the translational regulation of the main coding region, as initiation would predominantly start at the uORF resulting in reduced availability of initiation factors, such as eIF2, for re-initiation at the downstream main open reading frame.</p> <p>Sequence-independent uORFs allow for low-level translation of the downstream main coding region [##REF##10446211##61##]. Low-level translation is possible, as sequence-independent uORFs do not cause ribosomal stalling as seen in sequence-dependent uORFs. The regulatory mechanism of the sequence-independent uORF involves other factors (uORF recognition, length, stop codon environment, and the downstream intercistronic sequence) that influence reinitiation efficiency [##REF##12890013##18##,##REF##12117416##19##], and more recently leaky scanning [##REF##14990743##48##], to regulate downstream translation. We analysed the codon usage of conserved uORFs and found no preferential usage of rare codons in the uORFs. Therefore, it is unlikely that the uORF codon usage could contribute to low-level translation as seen for certain rare codons in <italic>Xenopus laevis </italic>[##REF##12799445##50##] and <italic>Eschericia coli </italic>[##REF##16703415##51##] that can reduce translation.</p> <p>In conclusion, this study showed that the uORFSCAN pipeline is a useful tool for identifying conserved uORFs in closely related species. This pipeline has allowed us to identify 29 conserved uORFs in cereals. Possibly more conserved uORFs will be identified once the cDNA and assembled EST datasets become more comprehensive. These conserved rice uORFs will be useful for future functional analyses that should provide some perspective into downstream translational regulation by uORFs.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Upstream open reading frames (uORFs) can down-regulate the translation of the main open reading frame (mORF) through two broad mechanisms: ribosomal stalling and reducing reinitiation efficiency. In distantly related plants, such as rice and Arabidopsis, it has been found that conserved uORFs are rare in these transcriptomes with approximately 100 loci. It is unclear how prevalent conserved uORFs are in closely related plants.</p> <title>Results</title> <p>We used a homology-based approach to identify conserved uORFs in five cereals (monocots) that could potentially regulate translation. Our approach used a modified reciprocal best hit method to identify putative orthologous sequences that were then analysed by a comparative R-nomics program called uORFSCAN to find conserved uORFs.</p> <title>Conclusion</title> <p>This research identified new genes that may be controlled at the level of translation by conserved uORFs. We report that conserved uORFs are rare (&lt;150 loci contain them) in cereal transcriptomes, are generally short (less than 100 nt), highly conserved (50% median amino acid sequence similarity), position independent in their 5'-UTRs, and their start codon context and the usage of rare codons for translation does not appear to be important.</p>
<title>List of abbreviations</title> <p>5'-UTR: five prime untranslated region; uORF: upstream open reading frame; EST: expressed sequence tag; RBH: reciprocal best hit; CDS: coding sequence.</p> <title>Authors' contributions</title> <p>MT conducted the research, analysed the data and drafted the manuscript. UB and CJS designed the research, participated in the study design, coordinated the study and contributed to the final manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>The authors thank Julian Schwerdt for high performance computing support at the South Australian Partnership for Advanced Computing (SAPAC), Andreas Schreiber for statistical support, and Rodney Davies for helpful discussions. This work was supported by the Australian Centre for Plant Functional Genomics (ACPFG) funded by Grains Research and Development Corporation (GRDC), Australian Research Council (ARC), the University of Adelaide, and the Government of South Australia. Michael Tran is the recipient of an ACPFG Postgraduate Scholarship.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Tran_Figure1.eps</bold>. 'Overview of the uORFSCAN pipeline'. The pipeline consists of four steps: 1) Identifying putative orthologues using a modified reciprocal best hit (rbh) method, 2) Clustering of orthologues according to how many cereal species they are found in, 3) Using uORFSCAN program to find conserved uORFs using a comparative approach, and 4) manual curation of predicted conserved cereal and Arabidopsis uORFs.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Tran_Figure2.eps</bold>. 'The position of conserved uORFs within their 5'-UTRs'. It contains rice uORFs conserved in four other cereals and in Arabidopsis.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Tran_Figure3.eps</bold>. 'A frequency distribution of the length (nt) of rice uORFs conserved in four other cereals and in Arabidopsis.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>The uORFs predicted by uORFSCAN in 5 out of 5</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left" colspan="2">Rice</td><td align="center" colspan="2">Wheat</td><td align="center" colspan="2">Barley</td><td align="center" colspan="2">Maize</td><td align="center" colspan="2">Sorghum</td><td align="left">Avg. A.A. similarity (%)</td><td align="left">Putative function<sup>b</sup></td></tr><tr><td colspan="10"><hr/></td><td/><td/></tr><tr><td align="left">Identifier</td><td align="left">5'-UTR</td><td align="left">Identifier</td><td align="left">5'-UTR<sup>a</sup></td><td align="left">Identifier</td><td align="left">5'-UTR<sup>a</sup></td><td align="left">Identifier</td><td align="left">5'-UTR<sup>a</sup></td><td align="left">Identifier</td><td align="left">5'-UTR<sup>a</sup></td><td/><td/></tr></thead><tbody><tr><td align="left">AK106095</td><td align="left">131_9_17</td><td align="left">TC265929</td><td align="left">113_9_16</td><td align="left">TC148181</td><td align="left">67_9_16</td><td align="left">TC288369</td><td align="left">131_9_17</td><td align="left">TC102998</td><td align="left">149_9_17</td><td align="left">100</td><td align="left">Oxidoreductase</td></tr><tr><td align="left">AK103391<sup>c</sup></td><td align="left">205_75_74</td><td align="left">TC269775</td><td align="left">251_75_62</td><td align="left">TC134190</td><td align="left">204_75_62</td><td align="left">TC294011</td><td align="left">215_75_75</td><td align="left">TC103599</td><td align="left">106_75_378</td><td align="left">88</td><td align="left">Trehalose-6-phosphate phosphatase</td></tr><tr><td align="left">AK100589 <sup>d,e</sup></td><td align="left">240_9_334</td><td align="left">TC264559</td><td align="left">201_9_317</td><td align="left">TC130707</td><td align="left">228_9_318</td><td align="left">TC292591</td><td align="left">286_9_320</td><td align="left">TC91317</td><td align="left">260_9_329</td><td align="left">50</td><td align="left">S-adenosylmethionine decarboxylase (AdoMetDC)</td></tr><tr><td/><td align="left">248_156_179</td><td/><td align="left">209_150_168</td><td/><td align="left">236_150_169</td><td/><td align="left">294_153_168</td><td/><td align="left">268_153_177</td><td align="left">90</td><td/></tr><tr><td/><td align="left">296_108_179</td><td/><td align="left">254_105_168</td><td/><td align="left">281_105_169</td><td/><td align="left">336_111_168</td><td/><td align="left">310_111_177</td><td align="left">92</td><td/></tr><tr><td align="left">AK073303</td><td align="left">67_9_142</td><td align="left">TC237149</td><td align="left">75_9_113</td><td align="left">TC132556</td><td align="left">81_9_139</td><td align="left">TC305609</td><td align="left">127_9_69</td><td align="left">TC102988</td><td align="left">222_9_69</td><td align="left">50</td><td align="left">Alkaline phytoceramidase</td></tr><tr><td/><td align="left">135_9_74</td><td/><td align="left">75_9_113</td><td/><td align="left">81_9_139</td><td/><td align="left">127_9_69</td><td/><td align="left">222_9_69</td><td align="left">50</td><td/></tr><tr><td align="left">AK072868 <sup>e</sup></td><td align="left">249_27_248</td><td align="left">TC247418</td><td align="left">258_27_266</td><td align="left">TC139536</td><td align="left">298_27_272</td><td align="left">TC306591</td><td align="left">444_27_564</td><td align="left">TC102544</td><td align="left">331_27_265</td><td align="left">11</td><td align="left">CBL-interacting protein kinase</td></tr><tr><td/><td align="left">259_195_70</td><td/><td align="left">268_198_85</td><td/><td align="left">308_198_91</td><td/><td align="left">260_192_583</td><td/><td align="left">341_195_87</td><td align="left">29</td><td/></tr><tr><td/><td align="left">269_39_216</td><td/><td align="left">278_39_234</td><td/><td align="left">318_39_240</td><td/><td align="left">768_39_228</td><td/><td align="left">351_39_233</td><td align="left">8</td><td/></tr><tr><td/><td align="left">338_90_96</td><td/><td align="left">347_93_111</td><td/><td align="left">387_93_117</td><td/><td align="left">576_93_366</td><td/><td align="left">420_90_113</td><td align="left">10</td><td/></tr><tr><td/><td align="left">392_36_96</td><td/><td align="left">404_36_111</td><td/><td align="left">444_36_117</td><td/><td align="left">633_36_366</td><td/><td align="left">474_36_113</td><td align="left">8</td><td/></tr><tr><td align="left">AK072649</td><td align="left">100_192_117</td><td align="left">TC236348</td><td align="left">79_192_117</td><td align="left">TC133316</td><td align="left">76_192_93</td><td align="left">TC305793</td><td align="left">180_192_116</td><td align="left">TC93140</td><td align="left">168_192_116</td><td align="left">78</td><td align="left">Ribosomal protein S6 kinase</td></tr><tr><td align="left">AK066145</td><td align="left">178_12_58</td><td align="left">TC266262</td><td align="left">149_12_73</td><td align="left">TC134484</td><td align="left">154_12_231</td><td align="left">TC286452</td><td align="left">224_12_70</td><td align="left">TC94546</td><td align="left">187_12_69</td><td align="left">33</td><td align="left">Ubiquitin-fold protein</td></tr><tr><td align="left">AK064792</td><td align="left">276_15_187</td><td align="left">TC267323</td><td align="left">254_15_188</td><td align="left">TC132983</td><td align="left">253_15_-9</td><td align="left">TC306152</td><td align="left">263_15_170</td><td align="left">TC107743</td><td align="left">230_15_150</td><td align="left">87</td><td align="left">F9L1.29 protein</td></tr><tr><td align="left">AK060523</td><td align="left">173_123_185</td><td align="left">TC235416</td><td align="left">201_126_157</td><td align="left">TC148319</td><td align="left">211_120_163</td><td align="left">TC305149</td><td align="left">255_129_195</td><td align="left">TC103609</td><td align="left">240_129_212</td><td align="left">58</td><td align="left">Ankyrin-3</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Criteria for verifying rice uORFs in 5 out of 5</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Accession</td><td align="left">FL- cDNA<sup>a</sup></td><td align="center">Upstream &amp; In-frame stop codon</td><td align="center">Agreement with genome annotation<sup>b</sup></td><td align="center" colspan="6">Alignment of uORFSCAN identified main proteins with UniProt proteins<sup>c</sup></td><td align="left">uORF valid</td></tr><tr><td/><td/><td/><td/><td colspan="6"><hr/></td><td/></tr><tr><td/><td/><td/><td/><td align="center">UniProt protein length (AA)</td><td align="center">Align length (AA)</td><td align="center">Identities (%)</td><td align="center">Expect</td><td align="center">Annotation</td><td align="center">GO classication <italic>(Arabidopsis thaliana)</italic></td><td/></tr></thead><tbody><tr><td align="left">AK106095</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">392</td><td align="left">392</td><td align="left">100</td><td align="left">2.2e-217</td><td align="center">Oxidoreductase</td><td align="center"><underline>[go:19538]</underline> protein metabolism<break/><underline>[go:16706]</underline> oxidoreductase activity</td><td align="center">Yes</td></tr><tr><td align="left">AK103391</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">371</td><td align="left">371</td><td align="left">100</td><td align="left">3.4e-194</td><td align="center">Trehalose-6-phosphate phosphatase</td><td align="center"><underline>[go:5992]</underline> trehalose biosynthesis<break/><underline>[go:9507]</underline> chloroplast<break/><underline>[go:4805]</underline> trehalose-phosphatase activity</td><td align="center">Yes</td></tr><tr><td align="left">AK100589</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">398</td><td align="left">398</td><td align="left">100</td><td align="left">1.1e-215</td><td align="center">AdoMetDC</td><td align="center"><underline>[go:6596]</underline> polyamine biosynthesis<break/><underline>[go:5694]</underline> chromosome</td><td align="center">Yes</td></tr><tr><td align="left">AK073303</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">257</td><td align="left">257</td><td align="left">100</td><td align="left">1.6e-141</td><td align="center">Acyl-CoA independent ceramide synthase</td><td align="center"><underline>[go:6672]</underline> ceramide metabolism<break/><underline>[go:16020]</underline> membrane<break/><underline>[go:3824]</underline> catalytic activity<break/><underline>[go:16811]</underline> hydrolase activity</td><td align="center">Yes</td></tr><tr><td align="left">AK072868</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">443</td><td align="left">443</td><td align="left">100</td><td align="left">3.6e-238</td><td align="center">uncharacterized protein (probable CBL-interacting serine/threonine-protein kinase 15)</td><td align="center"><underline>[go:6468]</underline> protein phosphorylation<break/><underline>[go:7165]</underline> signal transduction<break/><underline>[go:5524]</underline> ATP binding<break/><underline>[go:4672]</underline> protein kinase activity</td><td align="center">Yes</td></tr><tr><td align="left">AK072649</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">480</td><td align="left">488</td><td align="left">76</td><td align="left">9.6e-199</td><td align="center">Ribosomal protein S6 Kinase</td><td align="center"><underline>[go:45946]</underline> positive translation<break/><underline>[go:6468]</underline> protein phosphorylation<break/><underline>[go:9507]</underline> chloroplast<break/><underline>[go:16301]</underline> kinase activity</td><td align="center">Yes</td></tr><tr><td align="left">AK066145</td><td align="left">No</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">119</td><td align="left">119</td><td align="left">100</td><td align="left">1.3e-59</td><td align="center">Membrane-anchored ubiquitin-fold protein</td><td align="center"><underline>[go:6464]</underline> protein modification</td><td align="center">Yes</td></tr><tr><td align="left">AK064792</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">197<sup>d</sup></td><td align="left">109</td><td align="left">108</td><td align="left">57</td><td align="left">8.4e-26</td><td align="center">F9L1.29 protein</td><td align="center">Not available</td><td align="center">Yes</td></tr><tr><td align="left">AK060523</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">166</td><td align="left">166</td><td align="left">100</td><td align="left">1.9e-88</td><td align="center">uncharacterized protein (probable ankyrin-3)</td><td align="center"><underline>[go:9507]</underline> chloroplast<break/><underline>[go:5515]</underline> protein binding</td><td align="center">Yes</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Rice uORFs predicted by uORFSCAN that are conserved in Arabidopsis</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center" colspan="2">Rice</td><td align="center" colspan="2">Arabidopsis</td><td align="left">Avg. A.A. similarity (%)</td><td align="left">Putative function<sup>b</sup></td></tr><tr><td colspan="2"><hr/></td><td colspan="2"><hr/></td><td/><td/></tr><tr><td align="left">Identifier</td><td align="left">5'-UTR<sup>a</sup></td><td align="left">Identifier</td><td align="left">5'-UTR<sup>a</sup></td><td/><td/></tr></thead><tbody><tr><td align="left">AK101100</td><td align="left">142_12_21<sup>c,d</sup></td><td align="left">AT1G51690.1</td><td align="left">555_12_1160</td><td align="left">33</td><td align="left">Protein phosphatase 2a</td></tr><tr><td align="left">AK066952</td><td align="left">365_66_182</td><td align="left">AT3G13225.1</td><td align="left">364_63_431</td><td align="left">27</td><td align="left">WW domain containing protein</td></tr><tr><td/><td align="left">368_63_182<sup>e</sup></td><td/><td align="left">364_63_431</td><td align="left">29</td><td/></tr><tr><td/><td align="left">503_51_59</td><td/><td align="left">553_51_254</td><td align="left">1</td><td/></tr><tr><td align="left">AK119592</td><td align="left">304_90_148<sup>c,d</sup></td><td align="left">AT3G01470.1</td><td align="left">162_87_120</td><td align="left">36</td><td align="left">Homeodomain leucine zipper protein</td></tr><tr><td align="left">AK100589</td><td align="left">248_156_179<sup>c,d</sup></td><td align="left">AT3G02470.3</td><td align="left">222_156_154</td><td align="left">82</td><td align="left">S-Adenosylmethionine decarboxylase</td></tr><tr><td align="left">AK103391</td><td align="left">176_30_148<sup>c,d,f</sup></td><td align="left">AT4G22590.1</td><td align="left">254_30_137</td><td align="left">44</td><td align="left">Trehalose-6-phosphate phosphatase</td></tr><tr><td/><td align="left">205_75_74<sup>g</sup></td><td/><td align="left">283_75_63<sup>g</sup></td><td align="left">71</td><td/></tr><tr><td align="left">AK069534</td><td align="left">813_9_432</td><td align="left">AT4G12770.1</td><td align="left">41_9_108</td><td align="left">50</td><td align="left">Auxilin-like protein</td></tr><tr><td align="left">AK069526</td><td align="left">214_126_544<sup>c</sup></td><td align="left">AT4G19110.2</td><td align="left">255_126_527</td><td align="left">44</td><td align="left">GAMYB-binding protein</td></tr><tr><td/><td align="left">690_9_185<sup>c</sup></td><td/><td align="left">603_9_296</td><td align="left">50</td><td/></tr><tr><td/><td align="left">820_36_28</td><td/><td align="left">398_36_474</td><td align="left">17</td><td/></tr><tr><td align="left">AK072868</td><td align="left">338_90_96<sup>c,d</sup></td><td align="left">AT5G58380.1</td><td align="left">11_87_295</td><td align="left">17</td><td align="left">CBL-interacting protein kinase</td></tr><tr><td align="left">AK060523</td><td align="left">173_123_185<sup>c</sup></td><td align="left">AT5G07840.1</td><td align="left">289_117_250</td><td align="left">36</td><td align="left">Ankyrin-3</td></tr><tr><td/><td/><td/><td align="left">313_93_250<sup>e</sup></td><td align="left">44</td><td/></tr><tr><td/><td align="left">206_90_185<sup>e</sup></td><td/><td align="left">313_93_250<sup>e</sup></td><td align="left">33</td><td/></tr><tr><td align="left">AK067412</td><td align="left">222_84_49<sup>c,h</sup></td><td align="left">AT5G50180.1</td><td align="left">357_84_79</td><td align="left">4</td><td align="left">Protein kinase ATN1</td></tr><tr><td align="left">AK102277</td><td align="left">228_117_150<sup>c</sup></td><td align="left">AT1G68550.1</td><td align="left">309_96_95</td><td align="left">21</td><td align="left">Hypothetical protein</td></tr><tr><td align="left"><italic>AK100332</italic></td><td align="left"><italic>1174_21_883</italic></td><td align="left"><italic>AT5G44800.1</italic></td><td align="left"><italic>359_21_3</italic></td><td align="left"><italic>14</italic></td><td align="left"><italic>Helicase</italic></td></tr><tr><td/><td align="left"><italic>1618_21_439</italic></td><td/><td align="left"><italic>359_21_3</italic></td><td align="left"><italic>17</italic></td><td/></tr><tr><td/><td align="left"><italic>1810_21_247</italic></td><td/><td align="left"><italic>359_21_3</italic></td><td align="left"><italic>17</italic></td><td/></tr><tr><td align="left"><italic>AK059639</italic></td><td align="left"><italic>1_45_784</italic><sup>c</sup></td><td align="left"><italic>ATCG00920.1</italic></td><td align="left"><italic>55_45_844</italic></td><td align="left"><italic>86</italic></td><td align="left"><italic>40S ribosomal protein S15</italic></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Criteria for verifying rice uORFs that are conserved in Arabidopsis</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Accession</td><td align="left">FL- cDNA<sup>a</sup></td><td align="left">Upstream &amp; In-frame stop codon</td><td align="left">Agreement with genome annotation<sup>b</sup></td><td align="left" colspan="6">Alignment of uORFSCAN identified main proteins with UniProt proteins<sup>c</sup></td><td align="left">uORF valid</td></tr><tr><td/><td/><td/><td/><td colspan="6"><hr/></td><td/></tr><tr><td/><td/><td/><td/><td align="left">UniProt protein length (AA)</td><td align="left">Align length (AA)</td><td align="left">Identities (%)</td><td align="left">Expect</td><td align="left">Annotation</td><td align="left">GO classication <italic>(Arabidopsis thaliana)</italic></td><td/></tr></thead><tbody><tr><td align="left">AK101100</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">525</td><td align="left">525</td><td align="left">100</td><td align="left">5.0e-287</td><td align="left">Protein phosphatase 2A</td><td align="left"><underline>[go:6470]</underline> protein dephosphorylation<break/><underline>[go:166]</underline> nucleotide binding</td><td align="left">Yes</td></tr><tr><td align="left">AK066952</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">860</td><td align="left">694</td><td align="left">99</td><td align="left">0</td><td align="left">WW domain containing protein</td><td align="left">Not available</td><td align="left">Yes<sup>d</sup></td></tr><tr><td align="left">AK119592</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">343</td><td align="left">343</td><td align="left">100</td><td align="left">6.8e-187</td><td align="left">Homeodomain leucine zipper protein</td><td align="left"><underline>[go:6355]</underline> regulation of transcription<break/><underline>[go:3677]</underline> DNA binding</td><td align="left">Yes</td></tr><tr><td align="left">AK100589</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">398</td><td align="left">398</td><td align="left">100</td><td align="left">1.1e-215</td><td align="left">S-Adenosylmethionine decarboxylase</td><td align="left"><underline>[go:6596]</underline> polyamine biosynthesis<break/><underline>[go:5694]</underline> chromosome</td><td align="left">Yes</td></tr><tr><td align="left">AK103391</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">371</td><td align="left">371</td><td align="left">100</td><td align="left">3.3e-194</td><td align="left">Trehalose-6-phosphate phosphatase</td><td align="left"><underline>[go:5992]</underline> trehalose biosynthesis<break/><underline>[go:9507]</underline> chloroplast</td><td align="left">Yes</td></tr><tr><td align="left">AK069534</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">1066<sup>e</sup></td><td align="left">485</td><td align="left">413</td><td align="left">61</td><td align="left">7.6e-117</td><td align="left">Auxilin-like protein</td><td align="left">Not available</td><td align="left">Yes<sup>f</sup></td></tr><tr><td align="left">AK069526</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">483</td><td align="left">483</td><td align="left">83</td><td align="left">5.8e-256</td><td align="left">GAMYB-binding protein</td><td align="left"><underline>[go:6468]</underline> protein phosphorylation<break/><underline>[go:5524]</underline> ATP binding<break/><underline>[go:16301]</underline> kinase activity</td><td align="left">Yes</td></tr><tr><td align="left">AK072868</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">443</td><td align="left">443</td><td align="left">100</td><td align="left">3.5e-238</td><td align="left">CBL-interacting kinase 15</td><td align="left"><underline>[go:6468]</underline> protein phosphorylation<break/><underline>[go:5524]</underline> ATP binding<break/><underline>[go:16301]</underline> kinase activity</td><td align="left">Yes</td></tr><tr><td align="left">AK060523</td><td align="left">No</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">166</td><td align="left">166</td><td align="left">99</td><td align="left">8.2e-88</td><td align="left">Ankyrin-3</td><td align="left"><underline>[go:5515]</underline> protein binding</td><td align="left">Yes</td></tr><tr><td align="left">AK067412</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">353</td><td align="left">353</td><td align="left">72</td><td align="left">1.2e-136</td><td align="left">Protein kinase ATN1</td><td align="left"><underline>[go:6468]</underline> protein phosphorylation<break/><underline>[go:5524]</underline> ATP binding<break/><underline>[go:16301]</underline> kinase activity</td><td align="left">Yes</td></tr><tr><td align="left">AK102277</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">338</td><td align="left">338</td><td align="left">99</td><td align="left">4.9e-179</td><td align="left">Hypothetical protein</td><td align="left">Not available</td><td align="left">Yes</td></tr><tr><td align="left">AK100332</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">4092<sup>e</sup></td><td align="left">2192</td><td align="left">872</td><td align="left">30</td><td align="left">5.3e-28</td><td align="left">Helicase</td><td align="left"><underline>[go:3676]</underline> nucleic acid binding<break/><underline>[go:6355]</underline> regulation of transcription<break/><underline>[go:5515]</underline> protein binding</td><td align="left">No<sup>g</sup></td></tr><tr><td align="left">AK059639</td><td align="left">No</td><td align="left">Yes</td><td align="left">Yes</td><td align="left">154</td><td align="left">154</td><td align="left">100</td><td align="left">2.6e-77</td><td align="left">40S ribosomal s15 protein</td><td align="left"><underline>[go:3735]</underline> structural part of ribosome<break/><underline>[go:6412]</underline> protein biosynthesis</td><td align="left">No<sup>h</sup></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T5"><label>Table 5</label><caption><p>Comparison of conserved cereal uORFs and their main ORF start context'</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center" colspan="7">In five cereals</td></tr></thead><tbody><tr><td align="left">Identifier</td><td align="left">uORF1</td><td align="left">uORF2</td><td align="left">uORF3</td><td align="left">uORF4</td><td align="left">uORF5</td><td align="left">Main ORF</td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="left">AK106095</td><td align="left">131_9_17<sup>a </sup>CCGATGC<sup>b</sup></td><td/><td/><td/><td/><td align="left">157_1179 CCCATGG</td></tr><tr><td align="left">AK103391</td><td align="left">205_75_74 TTGATGA</td><td/><td/><td/><td/><td align="left">354_1116 CAAATGG</td></tr><tr><td align="left">AK100589</td><td align="left">240_9_334 TGGATGT</td><td align="left">248_156_179 CTAATGG</td><td align="left">296_108_179 TTGATGT</td><td/><td/><td align="left">583_1197 CCAATGG</td></tr><tr><td align="left">AK073303</td><td align="left">67_9_142 TCCATGC</td><td align="left">135_9_74 CTCATGA</td><td/><td/><td/><td align="left">218_774 AGCATGG</td></tr><tr><td align="left">AK072868</td><td align="left">249_27_248 GGAATGC</td><td align="left">259_195_70 AAGATGT</td><td align="left">269_39_216 TGCATGC</td><td align="left">338_90_96 TTCATGA</td><td align="left">392_36_96 ACTATGG</td><td align="left">524_1332 GTGATGG</td></tr><tr><td align="left">AK072649</td><td align="left">100_192_117 CTCATGA</td><td/><td/><td/><td/><td align="left">409_1443 AAGATGG</td></tr><tr><td align="left">AK066145</td><td align="left">178_12_58 GCTATGG</td><td/><td/><td/><td/><td align="left">248_360 GAGATGG</td></tr><tr><td align="left">AK064792</td><td align="left">276_15_187 CGGATGC</td><td/><td/><td/><td/><td align="left">478_330 GGAATGG</td></tr><tr><td align="left">AK060523</td><td align="left">173_123_185 ACTATGG</td><td/><td/><td/><td/><td align="left">481_501 CGGATGG</td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="center" colspan="7">In rice and arabidopsis</td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="left">Identifier</td><td align="left">uORF1</td><td align="left">uORF2</td><td align="left">uORF3</td><td align="left">uORF4</td><td align="left">uORF5</td><td align="left">Main ORF</td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="left">AK101100</td><td align="left">142_12_21 GCCATGG</td><td/><td/><td/><td/><td align="left">175_1578 AAGATGG</td></tr><tr><td align="left">AK066952</td><td align="left">365_66_182 CCAATGA</td><td align="left">368_63_182 ATGATGA</td><td align="left">503_51_59 CTGATGA</td><td/><td/><td align="left">613_2085 GGGATGC</td></tr><tr><td align="left">AK119592</td><td align="left">304_90_148 CCGATGA</td><td/><td/><td/><td/><td align="left">542_1032 GCGATGG</td></tr><tr><td align="left">AK100589</td><td align="left">248_156_179 CTAATGG</td><td/><td/><td/><td/><td align="left">583_1197 CCAATGG</td></tr><tr><td align="left">AK103391</td><td align="left">176_30_148 AACATGA</td><td align="left">205_75_74 TTGATGA</td><td/><td/><td/><td align="left">354_1116 CAAATGG</td></tr><tr><td align="left">AK069534</td><td align="left">813_9_432 TCGATGA</td><td/><td/><td/><td/><td align="left">1254_1602 GAGATGC</td></tr><tr><td align="left">AK069526</td><td align="left">214_126_544 GATATGG</td><td align="left">690_9_185 TTGATGG</td><td align="left">820_36_28 CATATGA</td><td/><td/><td align="left">884_1455 AAAATGG</td></tr><tr><td align="left">AK072868</td><td align="left">338_90_96 TTCATGA</td><td/><td/><td/><td/><td align="left">524_1332 GTGATGG</td></tr><tr><td align="left">AK060523</td><td align="left">173_123_185 ACTATGG</td><td align="left">206_90_185 CCGATGC</td><td/><td/><td/><td align="left">481_501 CGGATGG</td></tr><tr><td align="left">AK067412</td><td align="left">222_84_49 CTGATGC</td><td/><td/><td/><td/><td align="left">355_1059 GGGATGG</td></tr><tr><td align="left">AK102277</td><td align="left">228_117_150 TCTATGC</td><td/><td/><td/><td/><td align="left">495_1017 GAAATGG</td></tr></tbody></table></table-wrap>
<disp-formula id="bmcM1"><label>(1)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1" name="1471-2164-9-361-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mtext>n</mml:mtext><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mtext>n</mml:mtext><mml:mrow><mml:mtext>obs</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mtext>N</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>n</mml:mtext></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mtext>n</mml:mtext></mml:msup><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>P</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mtext>N-n</mml:mtext></mml:mrow></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM2"><label>(2)</label><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2" name="1471-2164-9-361-i2" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mtext>n</mml:mtext><mml:mo>=</mml:mo><mml:msub><mml:mtext>n</mml:mtext><mml:mrow><mml:mtext>obs</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mtext>N</mml:mtext></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mtext>N</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>n</mml:mtext></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>P</mml:mi><mml:mtext>n</mml:mtext></mml:msup><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>P</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mtext>N-n</mml:mtext></mml:mrow></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></disp-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>TRAN_TableS1</bold>. 'The uORFs predicted by uORFSCAN in 4 out of 5'.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p><bold>TRAN_TableS2</bold>. 'The uORFs predicted by uORFSCAN in 3 out of 5'</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p><bold>TRAN_TableS3</bold>. 'The uORFs predicted by uORFSCAN in 4 out of 4'.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p><bold>TRAN_TableS4</bold>. 'The uORFs predicted by uORFSCAN in 3 out of 4'.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p><bold>TRAN_TableS5</bold>. 'The uORFs predicted by uORFSCAN in 3 out of 3'.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p><bold>Tran_FigureS1</bold>. 'The pattern of nucleotide sequence conservation calculated for the decanucleotide surrounding the uORF AUG triplet using WebLogo [##REF##15173120##62##]'. The overall height of each stack indicates the nucleotide sequence conservation at that position (measured in bits), whereas the height of nucleotide symbols (A, T, G, C) within the stack reflects the relative frequency of the corresponding nucleotide at that position. (B) Positions showing detectable nucleotide sequence conservation were magnified.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p><bold>Tran_FigureS2</bold>. 'Relative frequencies of codons showing significant deviation in codon usage between rice uORFs and rice main coding regions'. Rice uORF codon usage calculated from <ext-link ext-link-type="uri" xlink:href="http://www.bioinformatics.org/sms/codon_usage.html"/>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p><bold>TRAN_TableS6</bold>. 'ClustalW alignment of uORFs identified by uORFSCAN in 5 out of 5 cereals and in Arabidopsis'.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p><bold>uORFSCAN</bold>. 'uORFSCAN program'.</p></caption></supplementary-material>
<table-wrap-foot><p><sup>a </sup>Pre-orf distance_uORF length_intercistronic distance.</p><p><sup>b </sup>Functional annotation based on "The UniProt Knowledgebase (UniProt)" database.</p><p><sup>c </sup>One of several genes (identifiers) that are in multiple tables because different conserved uORFs were identified in the different datasets.</p><p><sup>d </sup>Previously reported upstream open reading frames (see 'Introduction' section on AdoMetDC).</p><p><sup>e </sup>Contain one or more nested uORFs.</p><p>Ribosomal rRNA genes have been removed.</p><p>See Table ##TAB##1##2## for criteria for verifying rice uORFs in 5 out of 5.</p></table-wrap-foot> <table-wrap-foot><p><sup>a </sup>Used rice cDNA in blastn search against "NCBI EST_Others" database (rice) to search for longer 5' ESTs.</p><p><sup>b </sup>Used rice cDNA in blastn search against "TIGR Rice Genome Annotation DB: Coding Sequences" database to verify the cDNA ORF.</p><p><sup>c </sup>Translated the rice cDNA in the same frame as the main open reading frame identified by uORFSCAN (include translations upstream of predicted start Methionine). The resulting protein sequence was used in a blastp search against "The UniProt Knowledgebase (UniProt)" database.</p><p><sup>d </sup>The genome annotation for the CDS is longer by the indicated number of base pairs.</p></table-wrap-foot> <table-wrap-foot><p><sup>a </sup>Pre-uORF distance_uORF length_intercistronic distance.</p><p><sup>b </sup>Functional annotation based on "the UniProt Knowledgebase (uniProt)" database.</p><p><sup>c </sup>Rice uORF is conserved in at least two orthologous cereal and Arabidopsis genes.</p><p><sup>d </sup>Rich in serine (at least 20%).</p><p><sup>e </sup>Nested uORF.</p><p><sup>f </sup>One of several genes (identifiers) that are in multiple tables because different conserved uORFs were identified in the different datasets.</p><p><sup>g </sup>Overlapping uORF.</p><p><sup>h </sup>Rice in arginine (at least 20%).</p><p>Ribosomal rRNA genes have been removed.</p><p>Rows in italics are false positive predictions (see Table ##TAB##3##4##. Criteria for verifying rice uORFs that are conserved in Arabidopsis)</p></table-wrap-foot> <table-wrap-foot><p><sup>a </sup>Used rice cDNA in blastn search against "NCBI EST_Others" database (rice) to search for longer 5' ESTs.</p><p><sup>b </sup>Used rice cDNA in blastn search against "TIGR Rice Genome Annotation DB: Coding Sequences" database to verify the cDNA ORF.</p><p><sup>c </sup>Translated the rice cDNA in the same frame as the main open reading frame identified by uORFSCAN (include translations upstream of predicted start Methionine). The resulting protein sequence was used in a blastp search against "The UniProt Knowledgebase (UniProt)" database.</p><p><sup>d </sup>The protein data suggests that the main open reading frame predicted by uORFSCAN extends further upstream, but does not overlap the predicted uORFs and so the uORFs are still valid.</p><p><sup>e </sup>The genome annotation for the CDS is longer by the indicated number of base pairs.</p><p><sup>f </sup>A shorter protein was identified, but does not overlap the predicted uORFs and so the uORFs are still valid.</p><p><sup>g </sup>A longer protein was identified indicating the main open reading frame extends further upstream, and does overlap the predicted uORFs and so the uORFs are not valid.</p><p><sup>h </sup>Possibly not functional because pre-orf distance is less than 20 nucleotides that is thought to be required for translation initiation.</p></table-wrap-foot> <table-wrap-foot><p><sup>a </sup>Pre-orf distance_uORF length_intercistronic distance.</p><p><sup>b </sup>uORF or mainORF sequence context from -3 position to +4.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-361-1"/> <graphic xlink:href="1471-2164-9-361-2"/> <graphic xlink:href="1471-2164-9-361-3"/>
<media xlink:href="1471-2164-9-361-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S4.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S5.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S6.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S7.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S8.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-361-S9.zip" mimetype="application" mime-subtype="x-zip-compressed"><caption><p>Click here for file</p></caption></media>
[{"collab": ["Selpi"], "surname": ["Bryant", "Kemp", "Cvijovic"], "given-names": ["CH", "GJL", "M"], "article-title": ["A first step towards learning which uORFs regulate gene expression"], "source": ["J Int Bioinformatics"], "year": ["2006"], "volume": ["3"], "fpage": ["31"], "lpage": ["44"]}]
BMC Genomics. 2008 Jul 31; 9:361