RECODE 2003

High-Throughput Sequencing of Phage Display Libraries Reveals Parasitic Enrichment of Indel Mutants Caused by Amplification Bias

The combination of phage display technology with high-throughput sequencing enables in-depth analysis of library diversity and selection-driven dynamics. We applied short-read sequencing of the mutagenized region on focused display libraries of two homologous nucleic acid modification eraser proteins-AlkB and FTO-biopanned against methylated DNA. This revealed enriched genotypes with small indels and concomitant doubtful amino acid motifs within the FTO library. Nanopore sequencing of the entire display vector showed additional enrichment of large deletions overlooked by region-specific sequencing, and further impacted the interpretation of the obtained amino acid motifs. We could attribute enrichment of these corrupted clones to amplification bias due to arduous FTO display slowing down host cell growth as well as phage production. This amplification bias appeared to be stronger than affinity-based target selection. Recommendations are provided for proper sequence analysis of phage display data, which can improve motive discovery in libraries of proteins that are difficult to display.

Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

Structural and Functional Characterization of Programmed Ribosomal Frameshift Signals in West Nile Virus Strains Reveals High Structural Plasticity Among cis-Acting RNA Elements

West Nile virus (WNV) is a prototypical emerging virus for which no effective therapeutics currently exist. WNV uses programmed -1 ribosomal frameshifting (-1 PRF) to synthesize the NS1' protein, a C terminally extended version of its non-structural protein 1, the expression of which enhances neuro-invasiveness and viral RNA abundance. Here, the NS1' frameshift signals derived from four WNV strains were investigated to better understand -1 PRF in this quasispecies. Sequences previously predicted to promote -1 PRF strongly promote this activity, but frameshifting was significantly more efficient upon inclusion of additional 3' sequence information. The observation of different rates of -1 PRF, and by inference differences in the expression of NS1', may account for the greater degrees of pathogenesis associated with specific WNV strains. Chemical modification and mutational analyses of the longer and shorter forms of the -1 PRF signals suggests dynamic structural rearrangements between tandem stem-loop and mRNA pseudoknot structures in two of the strains. A model is suggested in which this is employed as a molecular switch to fine tune the relative expression of structural to non-structural proteins during different phases of the viral replication cycle.

Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

We have identified conserved orthologs in completely sequenced genomes of double-strand DNA phages and arranged them into evolutionary families (phage orthologous groups [POGs]). Using this resource to analyze the collection of known phage genomes, we find that most orthologs are unique in their genomes (having no diverged duplicates [paralogs]), and while many proteins contain multiple domains, the evolutionary recombination of these domains does not appear to be a major factor in evolution of these orthologous families. The number of POGs has been rapidly increasing over the past decade, the percentage of genes in phage genomes that have orthologs in other phages has also been increasing, and the percentage of unknown "ORFans" is decreasing as more proteins find homologs and establish a family. Other properties of phage genomes have remained relatively stable over time, most notably the high fraction of genes that are never or only rarely observed in their cellular hosts. This suggests that despite the renowned ability of phages to transduce cellular genes, these cellular "hitchhiker" genes do not dominate the phage genomic landscape, and a large fraction of the genes in phage genomes maintain an evolutionary trajectory that is distinct from that of the host genes.

Mining Functional Elements in Messenger RNAs: Overview, Challenges, and Perspectives

Eukaryotic messenger RNA (mRNA) contains not only protein-coding regions but also a plethora of functional cis-elements that influence or coordinate a number of regulatory aspects of gene expression, such as mRNA stability, splicing forms, and translation rates. Understanding the rules that apply to each of these element types (e.g., whether the element is defined by primary or higher-order structure) allows for the discovery of novel mechanisms of gene expression as well as the design of transcripts with controlled expression. Bioinformatics plays a major role in creating databases and finding non-evident patterns governing each type of eukaryotic functional element. Much of what we currently know about mRNA regulatory elements in eukaryotes is derived from microorganism and animal systems, with the particularities of plant systems lagging behind. In this review, we provide a general introduction to the most well-known eukaryotic mRNA regulatory motifs (splicing regulatory elements, internal ribosome entry sites, iron-responsive elements, AU-rich elements, zipcodes, and polyadenylation signals) and describe available bioinformatics resources (databases and analysis tools) to analyze eukaryotic transcripts in search of functional elements, focusing on recent trends in bioinformatics methods and tool development. We also discuss future directions in the development of better computational tools based upon current knowledge of these functional elements. Improved computational tools would advance our understanding of the processes underlying gene regulations. We encourage plant bioinformaticians to turn their attention to this subject to help identify novel mechanisms of gene expression regulation using RNA motifs that have potentially evolved or diverged in plant species.

FSscan: a mechanism-based program to identify +1 ribosomal frameshift hotspots

In +1 programmed ribosomal frameshifting (PRF), ribosomes skip one nucleotide toward the 3'-end during translation. Most of the genes known to demonstrate +1 PRF have been discovered by chance or by searching homologous genes. Here, a bioinformatic framework called FSscan is developed to perform a systematic search for potential +1 frameshift sites in the Escherichia coli genome. Based on a current state of the art understanding of the mechanism of +1 PRF, FSscan calculates scores for a 16-nt window along a gene sequence according to different effects of the stimulatory signals, and ribosome E-, P- and A-site interactions. FSscan successfully identified the +1 PRF site in prfB and predicted yehP, pepP, nuoE and cheA as +1 frameshift candidates in the E. coli genome. Empirical results demonstrated that potential +1 frameshift sequences identified promoted significant levels of +1 frameshifting in vivo. Mass spectrometry analysis confirmed the presence of the frameshifted proteins expressed from a yehP-egfp fusion construct. FSscan allows a genome-wide and systematic search for +1 frameshift sites in E. coli. The results have implications for bioinformatic identification of novel frameshift proteins, ribosomal frameshifting, coding sequence detection and the application of mass spectrometry on studying frameshift proteins.

RNA modifications: a mechanism that modulates gene expression

Accuracy in the flow of genetic information from DNA to protein, or gene expression, is essential to the viability of an organisms. Pre-mRNA splicing and protein translation are two major steps in eukaryotic gene expression that necessitate the production of accurate gene products. Both processes occur in large complexes, consisting of both proteins and noncoding RNAs. Interestingly, the RNA components contain a large number of posttranscriptional modifications, including 2'-O-methylation and pseudouridylation, which are functionally important. In this chapter, we highlight the functional aspects of the modifications of spliceosomal snRNA and rRNA and provide a framework for understanding how posttranscriptional modifications are capable of influencing gene expression.

Recode-2: new design, new search tools, and many more genes

‘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

Selection for minimization of translational frameshifting errors as a factor in the evolution of codon usage

In a wide range of genomes, it was observed that the usage of synonymous codons is biased toward specific codons and codon patterns. Factors that are implicated in the selection for codon usage include facilitation of fast and accurate translation. There are two types of translational errors: missense errors and processivity errors. There is considerable evidence in support of the hypothesis that codon usage is optimized to minimize missense errors. In contrast, little is known about the relationship between codon usage and frameshifting errors, an important form of processivity errors, which appear to occur at frequencies comparable to the frequencies of missense errors. Based on the recently proposed pause-and-slip model of frameshifting, we developed Frameshifting Robustness Score (FRS). We used this measure to test if the pattern of codon usage indicates optimization against frameshifting errors. We found that the FRS values of protein-coding sequences from four analyzed genomes (the bacteria Bacillus subtilis and Escherichia coli, and the yeasts Saccharomyces cerevisiae and Schizosaccharomyce pombe) were typically higher than expected by chance. Other properties of FRS patterns observed in B. subtilis, S. cerevisiae and S. pombe, such as the tendency of FRS to increase from the 5′- to 3′-end of protein-coding sequences, were also consistent with the hypothesis of optimization against frameshifting errors in translation. For E. coli, the results of different tests were less consistent, suggestive of a much weaker optimization, if any. Collectively, the results fit the concept of selection against mistranslation-induced protein misfolding being one of the factors shaping the evolution of both coding and non-coding sequences.

Transterm: a database to aid the analysis of regulatory sequences in mRNAs

Messenger RNAs, in addition to coding for proteins, may contain regulatory elements that affect how the protein is translated. These include protein and microRNA-binding sites. Transterm (http://mRNA.otago.ac.nz/Transterm.html) is a database of regions and elements that affect translation with two major unique components. The first is integrated results of analysis of general features that affect translation (initiation, elongation, termination) for species or strains in Genbank, processed through a standard pipeline. The second is curated descriptions of experimentally determined regulatory elements that function as translational control elements in mRNAs. Transterm focuses on protein binding sites, particularly those in 3′-untranslated regions (3′-UTR). For this release the interface has been extensively updated based on user feedback. The data is now accessible by strain rather than species, for example there are 10 Escherichia coli strains (genomes) analysed separately. In addition to providing a repository of data, the database also provides tools for users to query their own mRNA sequences. Users can search sequences for Transterm or user defined regulatory elements, including protein or miRNA targets. Transterm also provides a central core of links to related resources for complementary analyses.

KnotInFrame: prediction of -1 ribosomal frameshift events

Programmed −1 ribosomal frameshift (−1 PRF) allows for alternative reading frames within one mRNA. First found in several viruses, it is now believed to exist in all kingdoms of life. Strong stimulators for −1 PRF are a heptameric slippery site and an RNA pseudoknot. Here, we present a new algorithm KnotInFrame, for the automatic detection of −1 PRF signals from genomic sequences. It finds the frameshifting stimulators by means of a specialized RNA-pseudoknot folding program, fast enough for genome-wide analyses. Evaluations on known −1 PRF signals demonstrate a high sensitivity.

PRFdb: a database of computationally predicted eukaryotic programmed -1 ribosomal frameshift signals

BACKGROUND

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.

RESULTS

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 (DeltaG degrees ) of the native mRNA sequence against a distribution of DeltaG degrees values of a pool of randomized sequences derived from the original. The data have been compiled in a user-friendly, easily searchable relational database.

CONCLUSION

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.

Site-specific release of nascent chains from ribosomes at a sense codon

"2A" oligopeptides are autonomous elements containing a D(V/I)EXNPGP motif at the C terminus. Protein synthesis from an open reading frame containing an internal 2A coding sequence yields two separate polypeptides, corresponding to sequences up to and including 2A and those downstream. We show that the 2A reaction occurs in the ribosomal peptidyltransferase center. Ribosomes pause at the end of the 2A coding sequence, over the glycine and proline codons, and the nascent chain up to and including this glycine is released. Translation-terminating release factors eRF1 and eRF3 play key roles in the reaction. On the depletion of eRF1, a greater proportion of ribosomes extend through the 2A coding sequence, yielding the full-length protein. In contrast, impaired eRF3 GTPase activity leads to many ribosomes failing to translate beyond 2A. Further, high-level expression of a 2A peptide-containing protein inhibits the growth of cells compromised for release factor activity and leads to errors in stop codon recognition. We propose that the nascent 2A peptide interacts with ribosomes to drive a highly unusual and specific "termination" reaction, despite the presence of a proline codon in the A site. After this, the majority of ribosomes continue translation, generating the separate downstream product.

A self-referential model for the formation of the genetic code

A model for the formation of the genetic code is presented where protein synthesis is directed initially by tRNA dimers. Proteins that are resistant to degradation and efficient RNA-binders protect the RNAs. Replication becomes elongational producing poly-tRNAs from which the mRNAs and ribosomes are derived. Attributions are successively fixed to tRNAs paired through the perfect palindromic anticodons, with the same bases at the extremities (5'ANA: UNU 3'; GNG: CNC; principal dinucleotides, pDiN). The 5' degeneracy is then developed. The first pairs to be encoded correspond to the hydropathy correlation outliers (Gly-CC: Pro-GG and Ser-GA: Ser-CU) and to the sector of homogeneous pDiN, composed by two pyrimidines or two purines. These amino acids are preferred in the N-ends of proteins, stabilizers of proteins against catabolism and strong RNA-binders. The next pairs complete the sector of homogeneous pDiN (Asp, Glu-UC: Leu-AG and Asn, Lys-UU: Phe-AA). This set of nine amino acids forms the protein cores with the predominant aperiodic conformation. Next enter the pairs with mixed pDiN (one purine and one pyrimidine), the RY attributions composing the protein N-ends and the YR attributions the C-ends. The last pair contains the main punctuation signs (Ile, Met, iMet-AU: Tyr, Stop-UA). The model indicates that genetic information emerged during the process of formation of the coding/decoding system and that genes were defined by the proteins. Stable proteins constructed the nucleoprotein system by binding to the RNAs that produced them. In this circular rationale, genes are memories in a metabolic system for production of proteins that stabilize it. The simplicity and the highly deterministic character of the process suggest that the Last Universal Common Ancestor populations could be composed, in early stages, of lineages bearing similar genetic codes.

The major form of hepatitis C virus alternate reading frame protein is suppressed by core protein expression

Hepatitis C virus (HCV) is a human RNA virus encoding 10 proteins in a single open reading frame. In the +1 frame, an ‘alternate reading frame’ (ARF) overlaps with the core protein-encoding sequence and encodes the ARF protein (ARFP). Here, we investigated the molecular regulatory mechanisms of ARFP expression in HCV target cells. Chimeric HCV-luciferase reporter constructs derived from the infectious HCV prototype isolate H77 were transfected into hepatocyte-derived cell lines. Translation initiation was most efficient at the internal AUG codon at position 86/88, resulting in the synthesis of a truncated ARFP named _86/88ARFP. Interestingly, _86/88ARFP synthesis was markedly enhanced in constructs containing an inactivated core protein reading frame. This enhancement was reversed by co-expression of core protein in trans, demonstrating suppression of ARFP synthesis by HCV core protein. In conclusion, our results indicate that translation of ARFP occurs mainly by alternative internal initiation at position 86/88 and is regulated by HCV core protein expression. The suppression of ARFP translation by HCV core protein suggests that ARFP expression is inversely linked to the level of viral replication. These findings define key mechanisms regulating ARFP expression and set the stage for further studies addressing the function of ARFP within the viral life cycle.

Ornithine decarboxylase antizyme finder (OAF): fast and reliable detection of antizymes with frameshifts in mRNAs

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/.

SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteins usually misidentified and misannotated in sequence databases. The presence of an in-frame UGA (stop) codon in the coding sequence of selenoprotein genes precludes their identification and correct annotation. The in-frame UGA codons are recoded to cotranslationally incorporate selenocysteine, a rare selenium-containing amino acid. The development of ad hoc experimental and, more recently, computational approaches have allowed the efficient identification and characterization of the selenoproteomes of a growing number of species. Today, dozens of selenoprotein families have been described and more are being discovered in recently sequenced species, but the correct genomic annotation is not available for the majority of these genes. SelenoDB is a long-term project that aims to provide, through the collaborative effort of experimental and computational researchers, automatic and manually curated annotations of selenoprotein genes, proteins and SECIS elements. Version 1.0 of the database includes an initial set of eukaryotic genomic annotations, with special emphasis on the human selenoproteome, for immediate inspection by selenium researchers or incorporation into more general databases. SelenoDB is freely available at http://www.selenodb.org.

Whole proteome analysis of post-translational modifications: applications of mass-spectrometry for proteogenomic annotation

While bacterial genome annotations have significantly improved in recent years, techniques for bacterial proteome annotation (including post-translational chemical modifications, signal peptides, proteolytic events, etc.) are still in their infancy. At the same time, the number of sequenced bacterial genomes is rising sharply, far outpacing our ability to validate the predicted genes, let alone annotate bacterial proteomes. In this study, we use tandem mass spectrometry (MS/MS) to annotate the proteome of Shewanella oneidensis MR-1, an important microbe for bioremediation. In particular, we provide the first comprehensive map of post-translational modifications in a bacterial genome, including a large number of chemical modifications, signal peptide cleavages, and cleavages of N-terminal methionine residues. We also detect multiple genes that were missed or assigned incorrect start positions by gene prediction programs, and suggest corrections to improve the gene annotation. This study demonstrates that complementing every genome sequencing project by an MS/MS project would significantly improve both genome and proteome annotations for a reasonable cost.

Genetic analysis of the E site during RF2 programmed frameshifting

The roles of the ribosomal E site are not fully understood. Prior evidence suggests that deacyl-tRNA in the E site can prevent frameshifting. We hypothesized that if the E-site codon must dissociate from its tRNA to allow for frameshifting, then weak codon:anticodon duplexes should allow for greater frameshifting than stronger duplexes. Using the well-characterized Escherichia coli RF2 (prfB) programmed frameshift to study frameshifting, we mutagenized the E-site triplet to all Unn and Cnn codons. Those variants should represent a very wide range of duplex stability. Duplex stability was estimated using two different methods. Frameshifting is inversely correlated with stability, as estimated by either method. These findings indicate that pairing between the deacyl-tRNA and the E-site codon opposes frameshifting. We discuss the implications of these findings on frame maintenance and on the RF2 programmed frameshift mechanism.

The three transfer RNAs occupying the A, P and E sites on the ribosome are involved in viral programmed -1 ribosomal frameshift

The -1 programmed ribosomal frameshifts (PRF), which are used by many viruses, occur at a heptanucleotide slippery sequence and are currently thought to involve the tRNAs interacting with the ribosomal P- and A-site codons. We investigated here whether the tRNA occupying the ribosomal E site that precedes a slippery site influences -1 PRF. Using the human immunodeficiency virus type 1 (HIV-1) frameshift region, we found that mutating the E-site codon altered the -1 PRF efficiency. When the HIV-1 slippery sequence was replaced with other viral slippery sequences, mutating the E-site codon also altered the -1 PRF efficiency. Because HIV-1 -1 PRF can be recapitulated in bacteria, we used a bacterial ribosome system to select, by random mutagenesis, 16S ribosomal RNA (rRNA) mutations that modify the expression of a reporter requiring HIV-1 -1 PRF. Three mutants were isolated, which are located in helices 21 and 22 of 16S rRNA, a region involved in translocation and E-site tRNA binding. We propose a novel model where -1 PRF is triggered by an incomplete translocation and depends not only on the tRNAs interacting with the P- and A-site codons, but also on the tRNA occupying the E site.

FSDB: a frameshift signal database

Programmed frameshifting is a recoding event in which a ribosome shifts reading frame by one or more nucleotides at a specific mRNA signal between overlapping genes. Programmed frameshifting is involved in the expression of many genes in a wide range of organisms, especially in viruses and bacteria. The mechanism of programmed frameshifting is not fully understood despite many studies, and there are few databases available for detailed information on programmed frameshifting. We have developed a database called FSDB (Frameshift Signal Database), which is a comprehensive compilation of experimentally known or computationally predicted data about programmed ribosomal frameshifting. FSDB provides a graphical view of frameshift signals and the genes using programmed frameshifting for their expression. It also allows the user himself/herself to find programmed frameshift sites in genomic sequences using a program called FSFinder (http://wilab.inha.ac.kr/fsfinder2). We believe FSDB will be a valuable resource for scientists studying programmed ribosomal frameshifting. FSDB is freely accessible at http://wilab.inha.ac.kr/fsdb/.

Anticodons, frameshifts, and hidden periodicities in tRNA sequences

Fourier analysis of the short-range periodicities for the complete set of sequences coding for tRNA genes in genome of Bacillus subtilis proves that periodicities with periods p = 2, 3, 4, and 6 sites are the inherent properties of tRNAs. The related periodicities should be understood in a broad statistical sense and their identifying needs the elaborate statistical methods. To improve the statistics, the analysis of significant periodicities was performed for the binary R-Y, S-W, and K-M sequences. Generally, such short-range periodicities are produced via biased positioning of particular nucleotides rather than via the tandem multiplication and subsequent modifications of repeats, though the latter mechanism may also be realized. Quasi-coherently piercing long segments of tRNA, the short-range periodicities create the effective long-range structural coupling between the acceptor stem and the anticodon loop and may participate in the mechanisms of molecular recognition. The periodicities with p = 2 and 4 provide the natural ground for the translation with spontaneous or programmed frameshifting and are present in tRNAs decoding the most frameshift-prone codons. The observation of short-range periodicities suggests that the mechanisms of amino-acylation of tRNAs and codon-anticodon pairing are not independent. Their study may also provide the important information related to the origin and evolution of the genetic code.

Identification of functional, endogenous programmed -1 ribosomal frameshift signals in the genome of Saccharomyces cerevisiae

In viruses, programmed -1 ribosomal frameshifting (-1 PRF) signals direct the translation of alternative proteins from a single mRNA. Given that many basic regulatory mechanisms were first discovered in viral systems, the current study endeavored to: (i) identify -1 PRF signals in genomic databases, (ii) apply the protocol to the yeast genome and (iii) test selected candidates at the bench. Computational analyses revealed the presence of 10 340 consensus -1 PRF signals in the yeast genome. Of the 6353 yeast ORFs, 1275 contain at least one strong and statistically significant -1 PRF signal. Eight out of nine selected sequences promoted efficient levels of PRF in vivo. These findings provide a robust platform for high throughput computational and laboratory studies and demonstrate that functional -1 PRF signals are widespread in the genome of Saccharomyces cerevisiae. The data generated by this study have been deposited into a publicly available database called the PRFdb. The presence of stable mRNA pseudoknot structures in these -1 PRF signals, and the observation that the predicted outcomes of nearly all of these genomic frameshift signals would direct ribosomes to premature termination codons, suggest two possible mRNA destabilization pathways through which -1 PRF signals could post-transcriptionally regulate mRNA abundance.

Transcriptional slippage controls production of type III secretion apparatus components in Shigella flexneri

During transcription, series of approximately 9 As or Ts can direct RNA polymerase to incorporate into the mRNA nucleotides not encoded by the DNA, changing the reading frame downstream from the slippage site. We detected series of 9 or 10 As in spa13, spa33 and mxiA encoding type III secretion apparatus components. Analysis of cDNAs indicated that transcriptional slippage occurs in spa13, mxiA and spa33. Changes in the reading frame were confirmed by using plasmids carrying slippage sites in the 5' part of lacZ. Slippage is required for production of Spa13 from two overlapping reading frames and should lead to production of truncated MxiA and Spa33 proteins. Complementation of spa13 and mxiA mutants with plasmids carrying altered sites indicated that slippage in spa13 is required for assembly of the secretion apparatus and that slippage sites in spa13 and mxiA have not been selected to encode Lys residues or to produce two proteins endowed with different activities. The presence of slippage sites decreases production of Spa13 by 70%, of MxiA and Spa33 by 15% and of Spa32 (encoded downstream from spa13) by 50%. These results suggest that transcriptional slippage controls protein production by reducing the proportion of mRNA translated into functional proteins.

Fluorescent T7 display phages obtained by translational frameshift

Lytic phages form a powerful platform for the display of large cDNA libraries and offer the possibility to screen for interactions with almost any substrate. To visualize these interactions directly by fluorescence microscopy, we constructed fluorescent T7 phages by exploiting the flexibility of phages to incorporate modified versions of its capsid protein. By applying translational frameshift sequences, helper plasmids were constructed that expressed a fixed ratio of both wild-type capsid protein (gp10) and capsid protein fused to enhanced yellow fluorescent protein (EYFP). The frameshift sequences were inserted between the 3′ end of the capsid gene and the sequence encoding EYFP. Fluorescent fusion proteins are only formed when the ribosome makes a −1 shift in reading frame during translation. Using standard fluorescence microscopy, we could sensitively monitor the enrichment of specific binders in a cDNA library displayed on fluorescent T7 phages. The perspectives of fluorescent display phages in the fast emerging field of single molecule detection and sorting technologies are discussed.

ARFA: a program for annotating bacterial release factor genes, including prediction of programmed ribosomal frameshifting

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

Recoding in bacteriophages and bacterial IS elements

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.

Transterm--extended search facilities and improved integration with other databases

Transterm has now been publicly available for >10 years. Major changes have been made since its last description in this database issue in 2002. The current database provides data for key regions of mRNA sequences, a curated database of mRNA motifs and tools to allow users to investigate their own motifs or mRNA sequences. The key mRNA regions database is derived computationally from Genbank. It contains 3' and 5' flanking regions, the initiation and termination signal context and coding sequence for annotated CDS features from Genbank and RefSeq. The database is non-redundant, enabling summary files and statistics to be prepared for each species. Advances include providing extended search facilities, the database may now be searched by BLAST in addition to regular expressions (patterns) allowing users to search for motifs such as known miRNA sequences, and the inclusion of RefSeq data. The database contains >40 motifs or structural patterns important for translational control. In this release, patterns from UTRsite and Rfam are also incorporated with cross-referencing. Users may search their sequence data with Transterm or user-defined patterns. The system is accessible at http://uther.otago.ac.nz/Transterm.html.

Protein repertoire of double-stranded DNA bacteriophages

The complexity and diversity of phage gene sets, which are produced by rapid evolution of phage genomes and rampant gene exchanges among phages, hamper the efforts to decipher the evolutionary relationships between individual phage proteins and reconstruct the complete set of evolutionary events leading to the known phages. To start unraveling the natural history of phages, we built the phage orthologous groups (POGs), a natural system of phage protein families that includes 6378 genes from 164 complete genome sequences of double-stranded DNA bacteriophages. Phage proteomes have high POG coverage: on average, 39 genes per phage genome belong to POGs, which is close to half of all genes in most phages. In an agreement with the notion of phage role in horizontal gene transfer, we see many cases of likely gene exchange between phages and their microbial hosts. At the same time, about 80% of all POGs are highly specific to phage genomes and are not commonly found in microbial genomes, indicating coherence and large degree of evolutionary independence of phage gene sets. The information on orthologous genes is essential for evolutionary classification of known bacteriophages and for reconstruction of ancestral phage genomes.

Identification of programmed translational -1 frameshifting sites in the genome of Saccharomyces cerevisiae

Frameshifting is a recoding event that allows the expression of two polypeptides from the same mRNA molecule. Most recoding events described so far are used by viruses and transposons to express their replicase protein. The very few number of cellular proteins known to be expressed by a -1 ribosomal frameshifting has been identified by chance. The goal of the present work was to set up a systematic strategy, based on complementary bioinformatics, molecular biology, and functional approaches, without a priori knowledge of the mechanism involved. Two independent methods were devised. The first looks for genomic regions in which two ORFs, each carrying a protein pattern, are in a frameshifted arrangement. The second uses Hidden Markov Models and likelihood in a two-step approach. When this strategy was applied to the Saccharomyces cerevisiae genome, 189 candidate regions were found, of which 58 were further functionally investigated. Twenty-eight of them expressed a full-length mRNA covering the two ORFs, and 11 showed a -1 frameshift efficiency varying from 5% to 13% (50-fold higher than background), some of which corresponds to genes with known functions. From other ascomycetes, four frameshifted ORFs are found fully conserved. Strikingly, most of the candidates do not display a classical viral-like frameshift signal and would have escaped a search based on current models of frameshifting. These results strongly suggest that -1 frameshifting might be more widely distributed than previously thought.

Genes in the postgenomic era

We outline three very different concepts of the gene-instrumental, nominal, and postgenomic. The instrumental gene has a critical role in the construction and interpretation of experiments in which the relationship between genotype and phenotype is explored via hybridization between organisms or directly between nucleic acid molecules. It also plays an important theoretical role in the foundations of disciplines such as quantitative genetics and population genetics. The nominal gene is a critical practical tool, allowing stable communication between bioscientists in a wide range of fields grounded in well-defined sequences of nucleotides, but this concept does not embody major theoretical insights into genome structure or function. The post-genomic gene embodies the continuing project of understanding how genome structure supports genome function, but with a deflationary picture of the gene as a structural unit. This final concept of the gene poses a significant challenge to conventional assumptions about the relationship between genome structure and function, and between genotype and phenotype.

Genome sequence of Blochmannia pennsylvanicus indicates parallel evolutionary trends among bacterial mutualists of insects

The distinct lifestyle of obligately intracellular bacteria can alter fundamental forces that drive and constrain genome change. In this study, sequencing the 792-kb genome of Blochmannia pennsylvanicus, an obligate endosymbiont of Camponotus pennsylvanicus, enabled us to trace evolutionary changes that occurred in the context of a bacterial-ant association. Comparison to the genome of Blochmannia floridanus reveals differential loss of genes involved in cofactor biosynthesis, the composition and structure of the cell wall and membrane, gene regulation, and DNA replication. However, the two Blochmannia species show complete conservation in the order and strand orientation of shared genes. This finding of extreme stasis in genome architecture, also reported previously for the aphid endosymbiont Buchnera, suggests that genome stability characterizes long-term bacterial mutualists of insects and constrains their evolutionary potential. Genome-wide analyses of protein divergences reveal 10- to 50-fold faster amino acid substitution rates in Blochmannia compared to related bacteria. Despite these varying features of genome evolution, a striking correlation in the relative divergences of proteins indicates parallel functional constraints on gene functions across ecologically distinct bacterial groups. Furthermore, the increased rates of amino acid substitution and gene loss in Blochmannia have occurred in a lineage-specific fashion, which may reflect life history differences of their ant hosts.

UMD (Universal Mutation Database): 2005 update

With the completion of the Human Genome Project, our vision of human genetic diseases has changed. The cloning of new disease-causing genes can now be performed in silico, and thousands of mutations are being identified in diagnostic and research laboratories yearly. Knowledge about these mutations and their association with clinical and biological data is essential for clinicians, geneticists, and researchers. To collect and analyze these data, we developed a generic software called Universal Mutation Databases (UMD) to create locus-specific databases. Here we report the new release (September 2004) of this freely available tool (www.umd.be), which allows the creation of LSDBs for virtually any gene and includes a large set of new analysis tools. We have implemented new features to integrate noncoding sequences, clinical data, pictures, monoclonal antibodies, and polymorphic markers (SNPs). Today the UMD retains all specifically designed tools to analyze mutations at the molecular level, as well as new sets of routines to search for genotype-phenotype correlations. We also created specific tools for infrequent mutations such as gross deletions and duplications, and deep intronic mutations. A large set of dedicated tools are now available for intronic mutations, including methods to calculate the consensus values (CVs) of potential splice sites and to search for exonic splicing enhancer (ESE) motifs. In addition, we have created specific routines to help researchers design new therapeutic strategies, such as exon skipping, aminoglycoside read-through of stop codons, or monoclonal antibody selection and epitope scanning for gene therapy.

A computational method to predict genetically encoded rare amino acids in proteins

A new method for predicting recoding by rare amino acids such as selenocysteine and pyrrolysine was used to survey a set of microbial genomes.

In several natural settings, the standard genetic code is expanded to incorporate two additional amino acids with distinct functionality, selenocysteine and pyrrolysine. These rare amino acids can be overlooked inadvertently, however, as they arise by recoding at certain stop codons. We report a method for such recoding prediction from genomic data, using read-through similarity evaluation. A survey across a set of microbial genomes identifies almost all the known cases as well as a number of novel candidate proteins.

Expression levels influence ribosomal frameshifting at the tandem rare arginine codons AGG_AGG and AGA_AGA in Escherichia coli

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.

A three-stemmed mRNA pseudoknot in the SARS coronavirus frameshift signal

A wide range of RNA viruses use programmed -1 ribosomal frameshifting for the production of viral fusion proteins. Inspection of the overlap regions between ORF1a and ORF1b of the SARS-CoV genome revealed that, similar to all coronaviruses, a programmed -1 ribosomal frameshift could be used by the virus to produce a fusion protein. Computational analyses of the frameshift signal predicted the presence of an mRNA pseudoknot containing three double-stranded RNA stem structures rather than two. Phylogenetic analyses showed the conservation of potential three-stemmed pseudoknots in the frameshift signals of all other coronaviruses in the GenBank database. Though the presence of the three-stemmed structure is supported by nuclease mapping and two-dimensional nuclear magnetic resonance studies, our findings suggest that interactions between the stem structures may result in local distortions in the A-form RNA. These distortions are particularly evident in the vicinity of predicted A-bulges in stems 2 and 3. In vitro and in vivo frameshifting assays showed that the SARS-CoV frameshift signal is functionally similar to other viral frameshift signals: it promotes efficient frameshifting in all of the standard assay systems, and it is sensitive to a drug and a genetic mutation that are known to affect frameshifting efficiency of a yeast virus. Mutagenesis studies reveal that both the specific sequences and structures of stems 2 and 3 are important for efficient frameshifting. We have identified a new RNA structural motif that is capable of promoting efficient programmed ribosomal frameshifting. The high degree of conservation of three-stemmed mRNA pseudoknot structures among the coronaviruses suggests that this presents a novel target for antiviral therapeutics.

Programmed ribosomal frameshifting in decoding the SARS-CoV genome

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.

Validation of Shewanella oneidensis MR-1 small proteins by AMT tag-based proteome analysis

Using stringent criteria for protein identification by accurate mass and time (AMT) tag mass spectrometric methodology, we detected 36 proteins of <101 amino acids in length, including 10 that were annotated as hypothetical proteins, in 172 global tryptic digests of Shewanella oneidensis MR-1 proteins. Peptides that map to the conserved, but functionally uncharacterized proteins SO4134 and SO2787, were the most frequently detected peptides in these samples, while those that map to hypotheticals SO2669 and SO2063, conserved hypotheticals SO0335 and SO2176, and the SlyX protein (SO1063) were observed at frequencies similar to those from essential small proteins (ribosomal proteins and translation initiation factor IF-1), suggesting that they may function in similarly important cellular functions. In addition, peptides were detected that map to 30 genes predicted to encode frameshifts, point mutations, or recoding signals. Of these 30 genes, peptides that map to positions beyond internal stop codons were detected in 13 genes (SO0101, SO0419, SO0590, SO0738, SO1113, SO1211, SO3079, SO3130, SO3240, SO4231, SO4328, SO4422, and SO4657). While expression of the full-length formate dehydrogenase encoded by SO0101 can be explained by incorporation of selenocysteine at the internal stop codon, the mechanism of translating downstream sequences in the remaining genes remains unknown.

P-site pairing subtleties revealed by the effects of different tRNAs on programmed translational bypassing where anticodon re-pairing to mRNA is separated from dissociation

Programmed ribosomal bypassing occurs in decoding phage T4 gene 60 mRNA. Half the ribosomes bypass a 50 nucleotide gap between codons 46 and 47. Peptidyl-tRNA dissociates from the "take-off" GGA, codon 46, and re-pairs to mRNA at a matched GGA "landing site" codon directly 5' of codon 47 where translation resumes. The system described here allows the contribution of peptidyl-tRNA re-pairing to be measured independently of dissociation. The matched GGA codons have been replaced by 62 other matched codons, giving a wide range of bypassing efficiencies. Codons with G or C in either or both of the first two codon positions yielded high levels of bypassing. The results are compared with those from a complementary study of non-programmed bypassing, where the combined effects of peptidyl-tRNA dissociation and reassociation were measured. The wild-type, GGA, matched codons are the most efficient in their gene 60 context in contrast to the relatively low value in the non-programmed bypassing study.

Predicting genes expressed via -1 and +1 frameshifts

Computational identification of ribosomal frameshift sites in genomic sequences is difficult due to their diverse nature, yet it provides useful information for understanding the underlying mechanisms and discovering new genes. We have developed an algorithm that searches entire genomic or mRNA sequences for frameshifting sites, and implements the algorithm as a web-based program called FSFinder (Frameshift Signal Finder). The current version of FSFinder is capable of finding -1 frameshift sites on heptamer sequences X XXY YYZ, and +1 frameshift sites for two genes: protein chain release factor B (prfB) and ornithine decarboxylase antizyme (oaz). We tested FSFinder on approximately 190 genomic and partial DNA sequences from a number of organisms and found that it predicted frameshift sites efficiently and with greater sensitivity and specificity than existing approaches. It has improved sensitivity because it considers many known components of a frameshifting cassette and searches these components on both + and - strands, and its specificity is increased because it focuses on overlapping regions of open reading frames and prioritizes candidate frameshift sites. FSFinder is useful for discovering unknown genes that utilize alternative decoding, as well as for analyzing frameshift sites. It is freely accessible at http://wilab.inha.ac.kr/FSFinder/.

P-site tRNA is a crucial initiator of ribosomal frameshifting

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.

Activation of the Pyrrolysine Suppressor tRNA Requires Formation of a Ternary Complex with Class I and Class II Lysyl-tRNA Synthetases

Monomethylamine methyltransferase of the archaeon Methanosarcina barkeri contains a rare amino acid, pyrrolysine, encoded by the termination codon UAG. Translation of this UAG requires the aminoacylation of the corresponding amber suppressor tRNAPyl. Previous studies reported that tRNAPyl could be aminoacylated by the synthetase-like protein PylS. We now show that tRNAPyl is efficiently aminoacylated in the presence of both the class I LysRS and class II LysRS of M. barkeri, but not by either enzyme acting alone or by PylS. In vitro studies show that both the class I and II LysRS enzymes must bind tRNAPyl in order for the aminoacylation reaction to proceed. Structural modeling and selective inhibition experiments indicate that the class I and II LysRSs form a ternary complex with tRNAPyl, with the aminoacylation activity residing in the class II enzyme.

Maintenance of the correct open reading frame by the ribosome

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.