PRFdb: A database of computationally predicted eukaryotic programmed -1 ribosomal frameshift signals
© Belew et al; licensee BioMed Central Ltd. 2008
Received: 24 January 2008
Accepted: 17 July 2008
Published: 17 July 2008
The Programmed Ribosomal Frameshift Database (PRFdb) provides an interface to help researchers identify potential programmed -1 ribosomal frameshift (-1 PRF) signals in eukaryotic genes or sequences of interest.
To identify putative -1 PRF signals, sequences are first imported from whole genomes or datasets, e.g. the yeast genome project and mammalian gene collection. They are then filtered through multiple algorithms to identify potential -1 PRF signals as defined by a heptameric slippery site followed by an mRNA pseudoknot. The significance of each candidate -1 PRF signal is evaluated by comparing the predicted thermodynamic stability (ΔG°) of the native mRNA sequence against a distribution of ΔG° values of a pool of randomized sequences derived from the original. The data have been compiled in a user-friendly, easily searchable relational database.
The PRFdB enables members of the research community to determine whether genes that they are investigating contain potential -1 PRF signals, and can be used as a metasource of information for cross referencing with other databases. It is available on the web at http://dinmanlab.umd.edu/prfdb.
Canonical decoding of the genetic code requires translating ribosomes to convert triplets of bases (codons) into amino acid sequences. Although this is algorithm is employed for translation of the vast majority of mRNA sequences, in some special cases cis-acting mRNA elements direct ribosomes into alternative reading frames, dynamically "recoding" their sequence information (reviewed in ). Programmed -1 Ribosomal Frameshifting (-1 PRF) was first discovered in RNA viruses where it enables viral genomes to encode multiple peptides from a single mRNA . An individual -1 PRF signal consists of a heptameric 'slippery site' usually followed by an mRNA pseudoknot secondary structure separated by a suitable spacer region (reviewed in [3–5]). Unlike their viral counterparts, eukaryotic genome-encoded -1 PRF signals are predicted to direct elongating ribosomes into premature termination codons (reviewed in ). Such events have been shown to initiate rapid mRNA degradation in yeast through the Nonsense Mediated Decay (NMD) pathway . As such, -1 PRF is hypothesized to add a novel modality for regulation of gene expression at the post-transcriptional level.
There are currently three databases serving the translational recoding community. RECODE (please see Availability and requirements for more information) is a browsable collection of all the published translational recoding signals [8, 9]. RECODE's strength is as central repository of all empirically proven translational recoding signals. FSDB (please see Availability and requirements for more information) contains a compilation of a handful of known and predicted viral, prokaryotic and eukaryotic -1 and +1 PRF signals, and also allows users to input their own sequences to search for frameshift signals using a program called FSFinder . Although incomplete, this site provides tools not available through RECODE, in particular it integration of PseudoViewer, a powerful graphics tool that simplifies visualization of H-type pseudoknots . MLOGD (please see Availability and requirements for more information) is a suite of software that allows detection of new protein-coding sequences by identifying overlapping open reading frames . While all three of these sites have their strengths, a common weakness it that they do not provide well catalogued, searchable databases of all potential recoding signals of any one kind. To fill this gap, we have created PRFdb as a database of predicted -1 PRF signals in eukaryotic genomes. The methods used to search for predicted -1 PRF signals have been previously described, and importantly, we have empirically demonstrated that a significant number of -1 PRF signals so identified actually promote significant levels of frameshifting . The strength of the PRFdb is that it provides a tool for researchers outside of the translational recoding field to use to quickly search for and identify potential -1 PRF signals in genes in which they are interested.
Construction and content
In the PRFdb, the predicted -1 PRF signals are represented by: 1) the genes in which they reside; 2) the identity and location of their slippery sites; 3) graphical representations of their predicted secondary structures; 4) computationally identified minimum free energies (MFE); and 5) the thermodynamic significance of these mRNA structures as compared to randomized variants. Currently completed genomes in the PRFdb are Saccharomyces cerevisiae (88,683 sequences windows from 4,238 of 6,352 genes); Homo sapiens (24,104 sequence windows from 6826 of 17,891 genes); Mus musculus (19,313 windows from 6,053 of 15,620 genes); Rattus norvegicus (5924 sequence windows comprising 1982 of 5341 genes); Bos Taurus (10,690 windows from comprising 3349 of 9187 genes); Danio rerio (7640 sequence windows comprising 2618 of 6197 genes); and Xenopus tropicalis (9776 sequence windows comprising 2712 of 5126 genes). As of this writing, sequences from Xenopus laevis are being evaluated and entered into the database. The Arabidopsis thaliana genome is queued next in the pipeline.
If a sequence of interest is not currently in the PRFdb, it can be imported via its NCBI accession number. Sequences added in this manner will be filtered within hours of import. Sequences imported into the PRFdb are also folded using sequential windows across the entire sequence in order to create a graphical minimum free energy 'landscape.' This enables users to submit longer or shorter sequence strings for computational folding, a particularly useful feature e.g. for eliminating extraneous sequence that may not be involved in actual RNA folding. For example, the computational analysis of the 100 nucleotide sequence downstream of the slippery site of the mouse Ma3 -1 PRF signal provided by the PRFdb predicts tandem stem loop structures. However, when only 55 nt of downstream sequence are provided, PRFdB predicts the empirically documented pseudoknot structure .
Sequences to be analyzed by the PRFdb are imported into the database, filtered using RNAMotif , folded with secondary structure prediction algorithms, randomized using one or more randomization methods, and refolded. To avoid complications of untranslated intronic sequences, the PRFdb contains only mature mRNA sequences (cDNA sequences primarily).
Sequences are imported into the PRFdb from a web interface using Genbank accession numbers, yeast genome accessions, or raw sequences. Each new sequence is first passed through a simple text filter that searches for slippery sites following the International Union of Pure and Applied Chemistry (IUPAC) pattern 'N NNW WWH' (or X XXY YYZ), where N (X) denotes any three identical bases, W (Y) denotes AAA or UUU, H (Z) ≠ G, and spaces indicate the incoming (zero) reading frame. Since the distance between the end of the slippery site and the downstream stimulatory sequence is important, a spacer of 1 to 8 nt was incorporated into the algorithm. Downstream sequences (100, 75 and 50 nt) are passed to RNAMotif with a descriptor looking for the potential to form an mRNA pseudoknot. Sequence windows passing this minimal test are passed to multiple pseudoknot predicting mRNA secondary structure prediction algorithms, including Pknots , Nupack , and HotKnots . Mfold  and the Iterative Loop Matching algorithm are not used because Mfold does not predict pseudoknots and ILM does not provide the minimum free energy values for predicted secondary structures. After folding, every sequence is randomized using one or more algorithms including: Fisher-Yates shuffling to maintain dinucleotide frequencies, codon frequencies, or nucleotide frequencies. The resulting random sequence windows are then refolded without searching for pseudoknots. This process is repeated a fixed number of times (100 by default) to create a distribution of sequence specific randomized MFEs. These resulting distribution of randomized MFEs is then compared the MFE of the original sequence window. These values are used to compute a z score, thus providing a measurement of the significance of the native sequence.
As with most bioinformatics projects, the generation of the PRFdB engenders more questions than answers. Most importantly, why is this database of useful to anyone outside of the frameshifting community? The response is that -1 PRF represents fundamental mechanism for the regulation of gene expression at post-transcriptional gene with far-reaching and broad impact. Thus, the PRFdB represents a resource that is useful and available to any researcher interested in how their gene or pathway of interest might be regulated at this level. Here, we have sought to simplify such inquiries by constructing a user-friendly web-based interface. The other major question concerns what the user might do following identification of a potential -1 PRF signal in a gene of interest. The classical "bench-based" approach that we are currently following is to clone the element into reporters and quantitatively determine whether it actually promotes -1 PRF (we typically use a cutoff level of ≥ 0.5% of a readthrough control), and whether it can function as an mRNA destabilizing element. Computational followup strategies include cross database referencing (e.g. assessing how -1 PRF in an mRNA or set of mRNAs might affect the transcriptome as monitored by microarray analysis), and phylogenetic analyses. In sum, the PRFdB provides yet another set of information for deep, one gene at a time mining, as well as for the broader approaches typified by the rapidly developing field of systematic biology.
Discussion and conclusion
With regard to future development, the PRFdb is currently processing the Xenopus laevis genome. The Arabidopsis thaliana and Caenorabditis elegans genomes are currently in the queue, and additional genomes will be processed in the near future. In addition, the BLAST interface to the PRFdb has been used to discover Genbank sequences similar to the most statistically significant sequences in the database, thus providing a means to expand the PRFdb in a depth first manner. As more similar sequences are completed, comparative genomics studies using sequence and/or mRNA structure alignments will be incorporated to enable identification of conserved -1 PRF signals across species and/or genes. As time progresses, additional computational and empirical information will allow for improved scoring, helping to increase the statistical relevance of the predicted secondary mRNA structures. These improvements will continue to make the PRFdb more useful and accessible to the research community, providing a resource allowing individual users to identify -1 PRF signals in genes of interest, and as a metasource of information for cross referencing with other databases, e.g. genomes and DNA microarray databases.
Availability and requirements
The PRFdb is freely available on the web at http://dinmanlab.umd.edu/prfdb and can be accessed though any standard web browser. Access and data downloading require no special requirements.
We would like to thank the National Cancer Institute for the generous use of their computing hardware, patience, and technical assistance. ATB was supported in part by a training grant from the NIH (T32 AI51967); NLH was supported in part by a grant to the University of Maryland from the Howard Hughes Medical Institute Undergraduate Science Education Program; JLJ was supported in part by fellowship to JLJ from the National Libraries of Medicine (F37 LM8333), and the project was supported by grants from the Public Health Service to JDD (R21 GM068123, R01 GM058859).
- Gesteland RF, Atkins JF: Recoding: Dynamic reprogramming of translation. Annu Rev Biochem. 1996, 65: 741-768. 10.1146/annurev.bi.65.070196.003521.PubMedView ArticleGoogle Scholar
- Jacks T, Varmus HE: Expression of the Rous Sarcoma Virus pol gene by ribosomal frameshifting. Science. 1985, 230: 1237-1242. 10.1126/science.2416054.PubMedView ArticleGoogle Scholar
- Brierley I: Ribosomal frameshifting on viral RNAs. J Gen Virol. 1995, 76: 1885-1892. 10.1099/0022-1317-76-8-1885.PubMedView ArticleGoogle Scholar
- Farabaugh PJ: Programmed translational frameshifting. Microbiol Rev. 1996, 60: 103-134.PubMedPubMed CentralGoogle Scholar
- Dinman JD, Berry MJ: Translational Control in Biology and Medicine. Edited by: Mathews MB, Sonenberg N and Hershey JWB. 2006, Cold Spring Harbor, NY, Cold Spring Harbor Press, 22: Regulation of Termination and Recoding. 3, Cold Spring Harbor Monograph Series 48Google Scholar
- Dinman JD: Programmed Ribosomal Frameshifting Goes Beyond Viruses: Organisms from all three kingdoms use frameshifting to regulate gene expression, perhaps signaling a paradigm shift. Microbe Wash DC. 2006, 1: 521-527.PubMedPubMed CentralGoogle Scholar
- Plant EP, Wang P, Jacobs JL, Dinman JD: A programmed -1 ribosomal frameshift signal can function as a cis-acting mRNA destabilizing element. Nucleic Acids Res. 2004, 32: 784-790. 10.1093/nar/gkh256.PubMedPubMed CentralView ArticleGoogle Scholar
- Baranov PV, Gurvich OL, Hammer AW, Gesteland RF, Atkins JF: Recode 2003. Nucleic Acids Research. 2003, 31: 87-89. 10.1093/nar/gkg024.PubMedPubMed CentralView ArticleGoogle Scholar
- Baranov PV, Gurvich OL, Fayet O, Prere MF, Miller WA, Gesteland RF, Atkins JF, Giddings MC: RECODE: a database of frameshifting, bypassing and codon redefinition utilized for gene expression. Nucleic Acids Research. 2001, 29: 264-267. 10.1093/nar/29.1.264.PubMedPubMed CentralView ArticleGoogle Scholar
- Moon S, Byun Y, Han K: FSDB: a frameshift signal database. Comput Biol Chem. 2007, 31: 298-302. 10.1016/j.compbiolchem.2007.05.004.PubMedView ArticleGoogle Scholar
- Han K, Byun Y: PSEUDOVIEWER2: Visualization of RNA pseudoknots of any type. Nucleic Acids Res. 2003, 31: 3432-3440. 10.1093/nar/gkg539.PubMedPubMed CentralView ArticleGoogle Scholar
- Firth AE, Brown CM: Detecting overlapping coding sequences in virus genomes. BMC Bioinformatics. 2006, 7: 75-10.1186/1471-2105-7-75.PubMedPubMed CentralView ArticleGoogle Scholar
- Jacobs JL, Belew AT, Rakauskaite R, Dinman JD: Identification of functional, endogenous programmed -1 ribosomal frameshift signals in the genome of Saccharomyces cerevisiae. Nucleic Acids Res. 2007, 35: 165-174. 10.1093/nar/gkl1033.PubMedPubMed CentralView ArticleGoogle Scholar
- Manktelow E, Shigemoto K, Brierley I: Characterization of the frameshift signal of Edr, a mammalian example of programmed -1 ribosomal frameshifting. Nucleic Acids Res. 2005, 33: 1553-1563. 10.1093/nar/gki299.PubMedPubMed CentralView ArticleGoogle Scholar
- Macke TJ, Ecker DJ, Gutell RR, Gautheret D, Case DA, Sampath R: RNAMotif, an RNA secondary structure definition and search algorithm. Nucleic Acids Research. 2001, 29: 4724-4735. 10.1093/nar/29.22.4724.PubMedPubMed CentralView ArticleGoogle Scholar
- Rivas E, Eddy SR: A dynamic programming algorithm for RNA structure prediction including pseudoknots. Journal of Molecular Biology. 1999, 285: 2053-2068. 10.1006/jmbi.1998.2436.PubMedView ArticleGoogle Scholar
- Dirks RM, Pierce NA: An algorithm for computing nucleic acid base-pairing probabilities including pseudoknots. J Comput Chem. 2004, 25: 1295-1304. 10.1002/jcc.20057.PubMedView ArticleGoogle Scholar
- Ren J, Rastegari B, Condon A, Hoos HH: HotKnots: heuristic prediction of RNA secondary structures including pseudoknots. RNA. 2005, 11: 1494-1504. 10.1261/rna.7284905.PubMedPubMed CentralView ArticleGoogle Scholar
- Mathews DH, Sabina J, Zuker M, Turner DH: Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J Mol Biol. 1999, 288: 911-940. 10.1006/jmbi.1999.2700.PubMedView ArticleGoogle Scholar