Search

Motivation: RNA editing is a phenomenon, which is responsible for the alteration of particular nucleotides in RNA sequences relative to their genomic templates. Recently, a large number of RNA editing instances in humans have been identified using bioinformatic screens and high-throughput experimental investigations utilizing next-generation sequencing technologies. However, the available data on RNA editing are not uniform and difficult to access., Results: Here, we describe a new database DARNED (DAtabase of RNa EDiting) that provides centralized access to available published data related to RNA editing. RNA editing locations are mapped on the reference human genome. The current release of the database contains information on approximately 42 000 human genome coordinates corresponding to RNA locations that undergo RNA editing, mostly involving adenosine-to-inosine (A-to-I) substitutions. The data can be queried using a range of genomic coordinates, their corresponding functional localization in RNA molecules [Exons, Introns, CoDing Sequence (CDS) and UnTranslated Regions (UTRs)] and information regarding tissue/organ/cell sources where RNA editing has been observed. It is also possible to obtain RNA editing information for a specific gene or an RNA molecule using corresponding accession numbers. Search results provide information on the number of expressed sequence tags (ESTs) supporting edited and genomic bases, functional localization of RNA editing and existence of known single nucleotide polymorphisms (SNPs). Editing data can be explored in UCSC and Ensembl genome browsers, in conjunction with additional data provided by these popular genome browsers. DARNED has been designed for researchers seeking information on RNA editing and for the developers of novel algorithms for its prediction., Availability: DARNED is accessible at http://darned.ucc.ie

'Recoding' is a term used to describe non-standard read-out of the genetic code, and encompasses such phenomena as programmed ribosomal frameshifting, stop codon readthrough, selenocysteine insertion and translational bypassing. Although only a small proportion of genes utilize recoding in protein synthesis, accurate annotation of 'recoded' genes lags far behind annotation of 'standard' genes. In order to address this issue, provide a service to researchers in the field, and offer training data for developers of gene-annotation software, we have gathered together known cases of recoding within the Recode database. Recode-2 is an improved and updated version of the database. It provides access to detailed information on genes known to utilize translational recoding and allows complex search queries, browsing of recoding data and enhanced visualization of annotated sequence elements. At present, the Recode-2 database stores information on approximately 1500 genes that are known to utilize recoding in their expression--a factor of approximately three increase over the previous version of the database. Recode-2 is available at http://recode.ucc.ie.

The genetic code appears to be optimized in its robustness to missense errors and frameshift errors. In addition, the genetic code is near-optimal in terms of its ability to carry information in addition to the sequences of encoded proteins. As evolution has no foresight, optimality of the modern genetic code suggests that it evolved from less optimal code variants. The length of codons in the genetic code is also optimal, as three is the minimal nucleotide combination that can encode the twenty standard amino acids. The apparent impossibility of transitions between codon sizes in a discontinuous manner during evolution has resulted in an unbending view that the genetic code was always triplet. Yet, recent experimental evidence on quadruplet decoding, as well as the discovery of organisms with ambiguous and dual decoding, suggest that the possibility of the evolution of triplet decoding from living systems with non-triplet decoding merits reconsideration and further exploration. To explore this possibility we designed a mathematical model of the evolution of primitive digital coding systems which can decode nucleotide sequences into protein sequences. These coding systems can evolve their nucleotide sequences via genetic events of Darwinian evolution, such as point-mutations. The replication rates of such coding systems depend on the accuracy of the generated protein sequences. Computer simulations based on our model show that decoding systems with codons of length greater than three spontaneously evolve into predominantly triplet decoding systems. Our findings suggest a plausible scenario for the evolution of the triplet genetic code in a continuous manner. This scenario suggests an explanation of how protein synthesis could be accomplished by means of long RNA-RNA interactions prior to the emergence of the complex decoding machinery, such as the ribosome, that is required for stabilization and discrimination of otherwise weak triplet codon-anticodon interactions.

Background: Ornithine decarboxylase antizymes are proteins which negatively regulate cellular polyamine levels via their affects on polyamine synthesis and cellular uptake. In virtually all organisms from yeast to mammals, antizymes are encoded by two partially overlapping open reading frames (ORFs). A +1 frameshift between frames is required for the synthesis of antizyme. Ribosomes change translation phase at the end of the first ORF in response to stimulatory signals embedded in mRNA. Since standard sequence analysis pipelines are currently unable to recognise sites of programmed ribosomal frameshifting, proper detection of full length antizyme coding sequences (CDS) requires conscientious manual evaluation by a human expert. The rapid growth of sequence information demands less laborious and more cost efficient solutions for this problem. This manuscript describes a rapid and accurate computer tool for antizyme CDS detection that requires minimal human involvement., Results: We have developed a computer tool, OAF (ODC antizyme finder) for identifying antizyme encoding sequences in spliced or intronless nucleic acid sequenes. OAF utilizes a combination of profile hidden Markov models (HMM) built separately for the products of each open reading frame constituting the entire antizyme coding sequence. Profile HMMs are based on a set of 218 manually assembled antizyme sequences. To distinguish between antizyme paralogs and orthologs from major phyla, antizyme sequences were clustered into twelve groups and specific combinations of profile HMMs were designed for each group. OAF has been tested on the current version of dbEST, where it identified over six thousand Expressed Sequence Tags (EST) sequences encoding antizyme proteins (over two thousand antizyme CDS in these ESTs are non redundant)., Conclusion: OAF performs well on raw EST sequences and mRNA sequences derived from genomic annotations. OAF will be used for the future updates of the RECODE database. OAF can also be useful for identifying novel antizyme sequences when run with relaxed parameters. It is anticipated that OAF will be used for EST and genome annotation purposes. OAF outputs sequence annotations in fasta, genbank flat file or XML format. The OAF web interface and the source code are freely available at http://recode.ucc.ie/oaf/ and at a mirror site http://recode.genetics.utah.edu/oaf/.

Unlabelled: Correct annotation of genes encoding release factors in bacterial genomes is often complicated by utilization of +1 programmed ribosomal frameshifting during synthesis of release factor 2, RF2. In the absence of robust computational approaches for predicting ribosomal frameshifting, the success of proper annotation depends on annotators' familiarity with this phenomenon. Here we describe a novel computer tool that allows automatic discrimination of genes encoding class-I bacterial release factors, RF1, RF2 and RFH. Most usefully, this program identifies and automatically annotates +1 frameshifting in RF2 encoding genes. Comparison of ARFA performance with existing annotations of bacterial genomes revealed that only 20% of RF2 genes utilizing ribosomal frameshifting during their expression are annotated correctly., Availability: The PHP based web interface of ARFA and the source code are located at http://recode.genetics.utah.edu/arfa

Background: While all codons that specify amino acids are universally recognized by tRNA molecules, codons signaling termination of translation are recognized by proteins known as class-I release factors (RF). In most eukaryotes and archaea a single RF accomplishes termination at all three stop codons. In most bacteria, there are two RFs with overlapping specificity, RF1 recognizes UA(A/G) and RF2 recognizes U(A/G)A., The Hypothesis: First, we hypothesize that orthologues of the E. coli K12 pseudogene prfH encode a third class-I RF that we designate RFH. Second, it is likely that RFH responds to signals other than conventional stop codons. Supporting evidence comes from the following facts: (i) A number of bacterial genomes contain prfH orthologues with no discernable interruptions in their ORFs. (ii) RFH shares strong sequence similarity with other class-I bacterial RFs. (iii) RFH contains a highly conserved GGQ motif associated with peptidyl hydrolysis activity (iv) residues located in the areas supposedly interacting with mRNA and the ribosomal decoding center are highly conserved in RFH, but different from other RFs. RFH lacks the functional, but non-essential domain 1. Yet, RFH-encoding genes are invariably accompanied by a highly conserved gene of unknown function, which is absent in genomes that lack a gene for RFH. The accompanying gene is always located upstream of the RFH gene and with the same orientation. The proximity of the 3' end of the former with the 5' end of the RFH gene makes it likely that their expression is co-regulated via translational coupling. In summary, RFH has the characteristics expected for a class-I RF, but likely with different specificity than RF1 and RF2., Testing the Hypothesis: The most puzzling question is which signals RFH recognizes to trigger its release function. Genetic swapping of RFH mRNA recognition components with its RF1 or RF2 counterparts may reveal the nature of RFH signals., Implications of the Hypothesis: The hypothesis implies a greater versatility of release-factor like activity in the ribosomal A-site than previously appreciated. A closer study of RFH may provide insight into the evolution of the genetic code and of the translational machinery responsible for termination of translation., Reviewers: This article was reviewed by Daniel Wilson (nominated by Eugene Koonin), Warren Tate (nominated by Eugene Koonin), Yoshikazu Nakamura (nominated by Eugene Koonin) and Eugene Koonin.

Dynamic shifts between open reading frames and the redefinition of codon meaning at specific sites, programmed by signals in mRNA, permits versatility of gene expression. Such alterations are characteristic of organisms in all domains of life and serve a variety of functional purposes. In this article, we concentrate on programmed ribosomal frameshifting, stop codon read-through and transcriptional slippage in the decoding of phage genes and bacterial mobile elements. Together with their eukaryotic counterparts, the genes encoding these elements are the richest known source of nonstandard decoding. Recent analyses revealed several novel sequences encoding programmed alterations in gene decoding and provide a glimpse of the emerging picture.

The rare codons AGG and AGA comprise 2% and 4%, respectively, of the arginine codons of Escherichia coli K-12, and their cognate tRNAs are sparse. At tandem occurrences of either rare codon, the paucity of cognate aminoacyl tRNAs for the second codon of the pair facilitates peptidyl-tRNA shifting to the +1 frame. However, AGG_AGG and AGA_AGA are not underrepresented and occur 4 and 42 times, respectively, in E. coli genes. Searches for corresponding occurrences in other bacteria provide no strong support for the functional utilization of frameshifting at these sequences. All sequences tested in their native context showed 1.5 to 11% frameshifting when expressed from multicopy plasmids. A cassette with one of these sequences singly integrated into the chromosome in stringent cells gave 0.9% frameshifting in contrast to two- to four-times-higher values obtained from multicopy plasmids in stringent cells and eight-times-higher values in relaxed cells. Thus, +1 frameshifting efficiency at AGG_AGG and AGA_AGA is influenced by the mRNA expression level. These tandem rare codons do not occur in highly expressed mRNAs.

Selenocysteine (Sec) and pyrrolysine (Pyl) are known as the 21st and 22nd amino acids in protein. Both are encoded by codons that normally function as stop signals. Sec specification by UGA codons requires the presence of a cis-acting selenocysteine insertion sequence (SECIS) element. Similarly, it is thought that Pyl is inserted by UAG codons with the help of a putative pyrrolysine insertion sequence (PYLIS) element. Herein, we analyzed the occurrence of Pyl-utilizing organisms, Pyl-associated genes, and Pyl-containing proteins. The Pyl trait is restricted to several microbes, and only one organism has both Pyl and Sec. We found that methanogenic archaea that utilize Pyl have few genes that contain in-frame UAG codons, and many of these are followed with nearby UAA or UGA codons. In addition, unambiguous UAG stop signals could not be identified. This bias was not observed in Sec-utilizing organisms and non-Pyl-utilizing archaea, as well as with other stop codons. These observations as well as analyses of the coding potential of UAG codons, overlapping genes, and release factor sequences suggest that UAG is not a typical stop signal in Pyl-utilizing archaea. On the other hand, searches for conserved Pyl-containing proteins revealed only four protein families, including methylamine methyltransferases and transposases. Only methylamine methyltransferases matched the Pyl trait and had conserved Pyl, suggesting that this amino acid is used primarily by these enzymes. These findings are best explained by a model wherein UAG codons may have ambiguous meaning and Pyl insertion can effectively compete with translation termination for UAG codons obviating the need for a specific PYLIS structure. Thus, Sec and Pyl follow dissimilar decoding and evolutionary strategies.

Programmed ribosomal frameshifting is an essential mechanism used for the expression of orf1b in coronaviruses. Comparative analysis of the frameshift region reveals a universal shift site U_UUA_AAC, followed by a predicted downstream RNA structure in the form of either a pseudoknot or kissing stem loops. Frameshifting in SARS-CoV has been characterized in cultured mammalian cells using a dual luciferase reporter system and mass spectrometry. Mutagenic analysis of the SARS-CoV shift site and mass spectrometry of an affinity tagged frameshift product confirmed tandem tRNA slippage on the sequence U_UUA_AAC. Analysis of the downstream pseudoknot stimulator of frameshifting in SARS-CoV shows that a proposed RNA secondary structure in loop II and two unpaired nucleotides at the stem I-stem II junction in SARS-CoV are important for frameshift stimulation. These results demonstrate key sequences required for efficient frameshifting, and the utility of mass spectrometry to study ribosomal frameshifting.

Background: Transcription slippage occurs on certain patterns of repeat mononucleotides, resulting in synthesis of a heterogeneous population of mRNAs. Individual mRNA molecules within this population differ in the number of nucleotides they contain that are not specified by the template. When transcriptional slippage occurs in a coding sequence, translation of the resulting mRNAs yields more than one protein product. Except where the products of the resulting mRNAs have distinct functions, transcription slippage occurring in a coding region is expected to be disadvantageous. This probably leads to selection against most slippage-prone sequences in coding regions., Results: To find a length at which such selection is evident, we analyzed the distribution of repetitive runs of A and T of different lengths in 108 bacterial genomes. This length varies significantly among different bacteria, but in a large proportion of available genomes corresponds to nine nucleotides. Comparative sequence analysis of these genomes was used to identify occurrences of 9A and 9T transcriptional slippage-prone sequences used for gene expression., Conclusions: IS element genes are the largest group found to exploit this phenomenon. A number of genes with disrupted open reading frames (ORFs) have slippage-prone sequences at which transcriptional slippage would result in uninterrupted ORF restoration at the mRNA level. The ability of such genes to encode functional full-length protein products brings into question their annotation as pseudogenes and in these cases is pertinent to the significance of the term 'authentic frameshift' frequently assigned to such genes.

The expression of some genes requires a high proportion of ribosomes to shift at a specific site into one of the two alternative frames. This utilized frameshifting provides a unique tool for studying reading frame control. Peptidyl-tRNA slippage has been invoked to explain many cases of programmed frameshifting. The present work extends this to other cases. When the A-site is unoccupied, the P-site tRNA can be repositioned forward with respect to mRNA (although repositioning in the minus direction is also possible). A kinetic model is presented for the influence of both, the cognate tRNAs competing for overlapping codons in A-site, and the stabilities of P-site tRNA:mRNA complexes in the initial and new frames. When the A-site is occupied, the P-site tRNA can be repositioned backward. Whether frameshifting will happen depends on the ability of the A-site tRNA to subsequently be repositioned to maintain physical proximity of the tRNAs. This model offers an alternative explanation to previously published mechanisms of programmed frameshifting, such as out-of-frame tRNA binding, and a different perspective on simultaneous tandem tRNA slippage.

Because of their compact genomes, retroelements (including retrotransposons and retroviruses) employ a variety of translational recoding mechanisms to express Gag and Pol. To assess the diversity of recoding strategies, we surveyed gag/pol gene organization among retroelements from diverse host species, including elements exhaustively recovered from the genome sequences of Caenorhabditis elegans, Drosophila melanogaster, Schizosaccharomyces pombe, Candida albicans, and Arabidopsis thaliana. In contrast to the retroviruses, which typically encode pol in the -1 frame relative to gag, nearly half of the retroelements surveyed encode a single gag-pol open reading frame. This was particularly true for the Ty1/copia group retroelements. Most animal Ty3/gypsy retroelements, on the other hand, encode gag and pol in separate reading frames, and likely express Pol through +1 or -1 frameshifting. Conserved sequences conforming to slippery sites that specify viral ribosomal frameshifting were identified among retroelements with pol in the -1 frame. None of the plant retroelements encoded pol in the -1 frame relative to gag; however, two closely related plant Ty3/gypsy elements encode pol in the +1 frame. Interestingly, a group of plant Ty1/copia retroelements encode pol either in a +1 frame relative to gag or in two nonoverlapping reading frames. These retroelements have a conserved stem-loop at the end of gag, and likely express pol either by a novel means of internal ribosomal entry or by a bypass mechanism.

It is generally believed that significant ribosomal frameshifting during translation does not occur without a functional purpose. The distribution of two frameshift-prone sequences, A_AAA_AAG and CCC_TGA, in coding regions of Escherichia coli has been analyzed. Although a moderate level of selection against the first sequence is evident, 68 genes contain A_AAA_AAG and 19 contain CCC_TGA. The majority of those tested in their genomic context showed >1% frameshifting. Comparative sequence analysis was employed to assess a potential biological role for frameshifting in decoding these genes. Two new candidates, in pheL and ydaY, for utilized frameshifting have been identified in addition to those previously known in dnaX and nine insertion sequence elements. For the majority of the shift-prone sequences no functional role can be attributed to them, and the frameshifting is likely erroneous. However, none of frameshift sequences is in the 306 most highly expressed genes. The unexpected conclusion is that moderate frameshifting during expression of at least some other genes is not sufficiently harmful for cells to trigger strong negative evolutionary pressure.

During translation, a string of non-overlapping triplet codons in messenger RNA is decoded into protein. The ability of a ribosome to decode mRNA without shifting between reading frames is a strict requirement for accurate protein biosynthesis. Despite enormous progress in understanding the mechanism of transfer RNA selection, the mechanism by which the correct reading frame is maintained remains unclear. In this report, evidence is presented that supports the idea that the translational frame is controlled mainly by the stability of codon-anticodon interactions at the P site. The relative instability of such interactions may lead to dissociation of the P-site tRNA from its codon, and formation of a complex with an overlapping codon, the process known as P-site tRNA slippage. We propose that this process is central to all known cases of +1 ribosomal frameshifting, including that required for the decoding of the yeast transposable element Ty3. An earlier model for the decoding of this element proposed 'out-of-frame' binding of A-site tRNA without preceding P-site tRNA slippage.

The RECODE database is a compilation of translational recoding events (programmed ribosomal frameshifting, codon redefinition and translational bypass). The database provides information about the genes utilizing these events for their expression, recoding sites, stimulatory sequences and other relevant information. The Database is freely available at http://recode.genetics.utah.edu/.

The mRNA encoding Escherichia coli polypeptide chain release factor 2 (RF2) has two partially overlapping reading frames. Synthesis of RF2 involves ribosomes shifting to the +1 reading frame at the end of the first open reading frame (ORF). Frameshifting serves an autoregulatory function. The RF2 gene sequences from the 86 additional bacterial species now available have been analyzed. Thirty percent of them have a single ORF and their expression does not require frameshifting. In the approximately 70% that utilize frameshifting, the sequence cassette responsible for frameshifting is highly conserved. In the E. coli RF2 gene, an internal Shine-Dalgarno (SD) sequence just before the shift site was shown earlier to be important for frameshifting. Mutagenic data presented here show that the spacer region between the SD sequence and the shift site influences frameshifting, and possible mechanisms are discussed. Internal translation initiation occurs at the shift site, but any functional role is obscure.

During the expression of a certain genes standard decoding is over-ridden in a site or mRNA specific manner. This recoding occurs in response to special signals in mRNA and probably occurs in all organisms. This review deals with the function and distribution of recoding with a focus on the ribosomal frameshifting used for gene expression in bacteria.

The RECODE database is a compilation of 'programmed' translational recoding events taken from the scientific literature and personal communications. The database deals with programmed ribosomal frameshifting, codon redefinition and translational bypass occurring in a variety of organisms. The entries for each event include the sequences of the corresponding genes, their encoded proteins for both the normal and alternate decoding, the types of the recoding events involved, trans-factors and cis-elements that influence recoding. The database is freely available at http://recode.genetics. utah.edu/.

The Database of Ribosomal Cross-links (DRC) was created in 1997. Here we describe new data incorporated into this database and several new features of the DRC. The DRC is freely available via World Wide Web at http://visitweb.com/database/ or http://www. mpimg-berlin-dahlem.mpg.de/ approximately ag_ribo/ag_brimacombe/drc/

23S rRNA from Escherichia coli was cleaved at single internucleotide bonds using ribonuclease H in the presence of appropriate chimeric oligonucleotides; the individual cleavage sites were between residues 384 and 385, 867 and 868, 1045 and 1046, and 2510 and 2511, with an additional fortuitous cleavage at positions 1117 and 1118. In each case, the 3' terminus of the 5' fragment was ligated to radioactively labeled 4-thiouridine 5'-,3'-biphosphate ("psUp"), and the cleaved 23S rRNA carrying this label was reconstituted into 50S subunits. The 50S subunits were able to associate normally with 30S subunits to form 70S ribosomes. Intra-RNA crosslinks from the 4-thiouridine residues were induced by irradiation at 350 nm, and the crosslink sites within the 23S rRNA were analyzed. The rRNA molecules carrying psUp at positions 867 and 1117 showed crosslinks to nearby positions on the opposite strand of the same double helix where the cleavage was located, and no crosslinking was detected from position 2510. In contrast, the rRNA carrying psUp at position 384 showed crosslinking to nt 420 (and sometimes also to 416 and 425) in the neighboring helix in 23S rRNA, and the rRNA with psUp at position 1045 gave a crosslink to residue 993. The latter crosslink demonstrates that the long helix 41-42 of the 23S rRNA (which carries the region associated with GTPase activity) must double back on itself, forming a "U-turn" in the ribosome. This result is discussed in terms of the topography of the GTPase region in the 50S subunit, and its relation to the locations of the 5S rRNA and the peptidyl transferase center.

The Database of Ribosomal Cross-links (DRC) provides a complete collection of all the published data produced by cross-linking studies on the Escherichia coli ribosome, as well as on its components and functional ligands. The DRC currently includes data on 986 cross-links from >100 research papers, yielded by >40 different reagents. For each cross-link, information is given concerning its location in the ribosome, the chemical or photochemical reagent applied, a brief description of the method(s) used to locate the cross-link, and the literature reference. The DRC is freely available via the World Wide Web at: http://Ribosome.Genebee.MSU.SU/DRC/ or at http://WWW:MPIMG-Berlin-Dahlem.MPG.DE/[symbol: see text]baranov/DRC/

A new approach for inserting a photo-label at a selected position within the long ribosomal RNA molecules has been developed. The Escherichia coli 16S rRNA was cleaved at a single internucleotide bond, 1141-1142, with RNase H in the presence of a complementary chimeric oligonucleotide. 4-Thiouridine 5', 3'-diphosphate was ligated to the 3'-end of the 5'fragment at the cleavage site with T4 RNA ligase. The 16S rRNA fragments containing this added photo-reactive nucleotide were assembled together with total 30S ribosomal proteins into small ribosomal subunits. The ability of such 30S particles containing fragmented rRNA to form 70S ribosomes has been demonstrated previously. Crosslinks were induced within the 30S subunits by mild UV irradiation. The sites of crosslinking within the 16S rRNA were then analyzed using RNase H digestion and reverse transcription. Two crosslinks from the thio-nucleotide attached to nt C1141 of 16S rRNA were observed, namely to nt U1295 and G1272. These results are in agreement with the established proximity of helix 39 and 41 in the 3D structure of the 30S ribosomal subunit, as shown by other intra RNA crosslinking data. These data furthermore allow us to refine the structural arrangement of helices 41 and 39 relative to one another.

Cleavage of 16S ribosomal RNA (rRNA) from E. coli "hammerhead" type ribozymes as well as by RNAase iI in the presence of "hymeric" (2'-deoxy-F-thymidine containing) oligonucleotides has been studied. The conditions for the cleavage of a desired single internucleotide bond have been found for a large molecule with a very complicated secondary and three-dimensional structure.