Recent developments in methods for RNA sequencing using in vitro 32P-labeling

A variety of approaches that utilize in vitro 32P-labeling of RNA and of oligonucleotides in the sequence analysis of RNAs are described. These include 1) methods for 5'- and 3'- end labeling of RNAs; 2) end labeling and sequencing of oligonucleotides present in complete T1 RNase or pancreatic RNase digests of RNA; 3) use of random endonucleases, such as nuclease P1, for terminal sequence analysis of end labeled RNAs; and 4) use of base specific enzymes or chemical reagents in the sequence analysis of end-labeled RNAs. Also described is an approach to RNA sequencing, applied so far to tRNAs, which is based on partial and random alkaline cleavage of an RNA to generate a series of overlapping oligonucleotide fragments, all containing the original 3'-end of the RNA. Analysis of the 5'- end group of each of these oligonucleotides (following 5'-end labeling with 32P) provides the sequence of most of the tRNA. The above methods have been used to derive the sequences of several tRNAs, the ribosomal 5S and 5 x 8S RNAs, a viroid RNA, and large segments of both prokaryotic and eukaryotic ribosomal and messenger RNAs.     

A computational proposal for designing structured RNA pools for in vitro selection of RNAs

Although in vitro selection technology is a versatile experimental tool for discovering novel synthetic RNA molecules, finding complex RNA molecules is difficult because most RNAs identified from random sequence pools are simple motifs, consistent with recent computational analysis of such sequence pools. Thus, enriching in vitro selection pools with complex structures could increase the probability of discovering novel RNAs. Here we develop an approach for engineering sequence pools that links RNA sequence space regions with corresponding structural distributions via a "mixing matrix" approach combined with a graph theory analysis. We define five classes of mixing matrices motivated by covariance mutations in RNA; these constructs define nucleotide transition rates and are applied to chosen starting sequences to yield specific nonrandom pools. We examine the coverage of sequence space as a function of the mixing matrix and starting sequence via clustering analysis. We show that, in contrast to random sequences, which are associated only with a local region of sequence space, our designed pools, including a structured pool for GTP aptamers, can target specific motifs. It follows that experimental synthesis of designed pools can benefit from using optimized starting sequences, mixing matrices, and pool fractions associated with each of our constructed pools as a guide. Automation of our approach could provide practical tools for pool design applications for in vitro selection of RNAs and related problems.     

Microarrays for identifying binding sites and probing structure of RNAs

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure–function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.     

Physical map of the Kirsten sarcoma virus genome as determined by fingerprinting RNase T1-resistant oligonucleotides.

From analysis of the large RNase T1-resistant oligonucleotides of Kirsten sarcoma virus (Ki-SV), a physical map of the virus genome was deduced. Kirsten murine leukemia virus (Ki-MuLV) sequences were detected in T1 oligonucleotides closest to the 3' end of the viral RNA and extended approximately 1,000 nucleotides into the genome. The rat genetic sequences started at this point and extended all the way to the very 5' end of the RNA molecules, where a small stretch of Ki-MuLV sequence was detected. By comparison of the fingerprints of Ki-SV RNA and the RNA of the endogenous rat src genetic sequences, it was found that more than 50% of the T1 oligonucleotides were similar between Ki-SV and the endogenous rat src RNA, suggesting an identical primary nucleotide sequence in over 50% of the viral genomes. The results indicate that Ki-SV arose by recombination between the 5' and 3' ends of Ki-MuLV and a large portion of the homologous sequences of the endogenous rat src RNA.     

In vitro RNA random pools are not structurally diverse: a computational analysis

In vitro selection of functional RNAs from large random sequence pools has led to the identification of many ligand-binding and catalytic RNAs. However, the structural diversity in random pools is not well understood. Such an understanding is a prerequisite for designing sequence pools to increase the probability of finding complex functional RNA by in vitro selection techniques. Toward this goal, we have generated by computer five random pools of RNA sequences of length up to 100 nt to mimic experiments and characterized the distribution of associated secondary structural motifs using sets of possible RNA tree structures derived from graph theory techniques. Our results show that such random pools heavily favor simple topological structures: For example, linear stem-loop and low-branching motifs are favored rather than complex structures with high-order junctions, as confirmed by known aptamers. Moreover, we quantify the rise of structural complexity with sequence length and report the dominant class of tree motifs (characterized by vertex number) for each pool. These analyses show not only that random pools do not lead to a uniform distribution of possible RNA secondary topologies; they point to avenues for designing pools with specific simple and complex structures in equal abundance in the goal of broadening the range of functional RNAs discovered by in vitro selection. Specifically, the optimal RNA sequence pool length to identify a structure with x stems is 20x.     

Prediction of common folding structures of homologous RNAs.

We have developed an algorithm and a computer program for simultaneously folding homologous RNA sequences. Given an alignment of M homologous sequences of length N, the program performs phylogenetic comparative analysis and predicts a common secondary structure conserved in the sequences. When the structure is not uniquely determined, it infers multiple structures which appear most plausible. This method is superior to energy minimization methods in the sense that it is not sensitive to point mutation of a sequence. It is also superior to usual phylogenetic comparative methods in that it does not require manual scrutiny for covariation or secondary structures. The most plausible 1-5 structures are produced in O(MN2 + N3) time and O(N2) space, which are the same requirements as those of widely used dynamic programs based on energy minimization for folding a single sequence. This is the first algorithm probably practical both in terms of time and space for finding secondary structures of homologous RNA sequences. The algorithm has been implemented in C on a Sun SparcStation, and has been verified by testing on tRNAs, 5S rRNAs, 16S rRNAs, TAR RNAs of human immunodeficiency virus type 1 (HIV-1), and RRE RNAs of HIV-1. We have also applied the program to cis-acting packaging sequences of HIV-1, for which no generally accepted structures yet exist, and propose potentially stable structures. Simulation of the program with random sequences with the same base composition and the same degree of similarity as the above sequences shows that structures common to homologous sequences are very unlikely to occur by chance in random sequences.     

A computer method for finding common base paired helices in aligned sequences: application to the analysis of random sequences.

We describe a new computer program that identifies conserved secondary structures in aligned nucleotide sequences of related single-stranded RNAs. The program employs a series of hash tables to identify and sort common base paired helices that are located in identical positions in more than one sequence. The program gives information on the total number of base paired helices that are conserved between related sequences and provides detailed information about common helices that have a minimum of one or more compensating base changes. The program is useful in the analysis of large biological sequences. We have used it to examine the number and type of complementary segments (potential base paired helices) that can be found in common among related random sequences similar in base composition to 16S rRNA from Escherichia coli. Two types of random sequences were analyzed. One set consisted of sequences that were independent but they had the same mononucleotide composition as the 16S rRNA. The second set contained sequences that were 80% similar to one another. Different results were obtained in the analysis of these two types of random sequences. When 5 sequences that were 80% similar to one another were analyzed, significant numbers of potential helices with two or more independent base changes were observed. When 5 independent sequences were analyzed, no potential helices were found in common. The results of the analyses with random sequences were compared with the number and type of helices found in the phylogenetic model of the secondary structure of 16S ribosomal RNA. Many more helices are conserved among the ribosomal sequences than are found in common among similar random sequences. In addition, conserved helices in the 16S rRNAs are, on the average, longer than the complementary segments that are found in comparable random sequences. The significance of these results and their application in the analysis of long non-ribosomal nucleotide sequences is discussed.     

Deletion mutant of the Bratislava-77 strain of Rous sarcoma virus containing a fusion of the group-specific antigen and envelope genes.

The genetic compositions of two independently derived preparations of the Bratislava-77 strain (B77) of Rous sarcoma virus were analyzed after each was passaged seven or more times in duck embryo fibroblasts. RNase, T1-resistant oligonucleotide fingerprint analysis of virion RNA from both preparations of duck-passaged B77 revealed the presence of two large noncontiguous deletions. Approximately 75% of the RNAs contained a deletion which spans oligonucleotides 304 to 4 on the viral genome (about 3,500 nucleotides) and encompasses all of the B77 polymerase gene. More than 90% of the RNAs also contained a deletion which spans src-specific oligonucleotides 6 and 5(about 2,200 nucleotides) and is identical to the deletion observed in transformation-defective B77. Virion RNA from duck-passaged B77 also contained two oligonucleotides (D1 and D2) not observed in the RNA of B77 virus grown on chicken embryo fibroblasts. Analysis of the virion RNA of duck-passaged B77 by denaturing agarose gel electrophoresis revealed four major subunits with molecular weights of 3.40 x 10(6), 2.65 x 10(6), 2.25 x 10(6), and 1.55 x 10(6). Whereas the 3.40- and 2.65-megadalton (Mdal) RNA species comigrated with the nondefective and transformation-defective RNAs of B77 propagated on chicken embryo fibroblasts, no counterparts to the 2.25- and 1.55-Mdal RNAs were observed in the RNA of B77 grown on chicken embryo fibroblasts. Oligonucleotide fingerprint analysis of these RNA species revealed that the 2.65-Mdal RNA contains the src-specific deletion and that 2.25-Mdal RNA contains the polymerase region deletion; both of these deletions were observed in the 1.55-Mdal RNA, which was the major RNA subunit species detected in duck-passaged B77. The new oligonucleotides (D1 and D2) observed in the duck-passaged virus were present in the 2.25- and 1.55-Mdal RNA species in vitro and in vivo and directs the synthesis of a 130,000-dalton protein (p130). p130 contains antigenic determinants specific for p27 (gag gene) and gp85 (env gene) but does not contain sequences which cross-react with antisera directed against the alpha beta form of RNA-dependent DNA polymerase (pol gene). This RNA, therefore, is generated by a fusion of the gag and env genes of Rous sarcoma virus B77.     

Frequency of RNA-RNA interaction in a model of the RNA World

The RNA World model for prebiotic evolution posits the selection of catalytic/template RNAs from random populations. The mechanisms by which these random populations could be generated de novo are unclear. Non-enzymatic and RNA-catalyzed nucleic acid polymerizations are poorly processive, which means that the resulting short-chain RNA population could contain only limited diversity. Nonreciprocal recombination of smaller RNAs provides an alternative mechanism for the assembly of larger species with concomitantly greater structural diversity; however, the frequency of any specific recombination event in a random RNA population is limited by the low probability of an encounter between any two given molecules. This low probability could be overcome if the molecules capable of productive recombination were redundant, with many nonhomologous but functionally equivalent RNAs being present in a random population. Here we report fluctuation experiments to estimate the redundancy of the set of RNAs in a population of random sequences that are capable of non-Watson-Crick interaction with another RNA. Parallel SELEX experiments showed that at least one in 10(6) random 20-mers binds to the P5.1 stem-loop of Bacillus subtilis RNase P RNA with affinities equal to that of its naturally occurring partner. This high frequency predicts that a single RNA in an RNA World would encounter multiple interacting RNAs within its lifetime, supporting recombination as a plausible mechanism for prebiotic RNA evolution. The large number of equivalent species implies that the selection of any single interacting species in the RNA World would be a contingent event, i.e., one resulting from historical accident.     

Isolation and sequence analysis of sea urchin (Lytechinus pictus) histone H4 messenger RNA

Histone messenger RNAs isolated from early blastula stage Lytechinus pictus sea urchin embryos have been separated into discrete RNA bands on polyacrylamide gels. The most rapidly migrating of these molecules, the putative histone H4 mRNA, has been digested with T₁ ribonuclease to generate oligonucleotides for nucleotide sequence analysis. Many of these sequences are colinear with the highly conserved amino acid sequence of histone H4 protein as determined for both cows and peas.Histone H4 messenger RNA hybridizes in conditions of DNA excess to sea urchin DNA which is repeated approximately 470-fold. Despite this level of repetition the nucleotide sequence of the H4 messenger RNA reflects little evolutionary divergence within the H4 genes of L. pictus as judged by the stoichiometric yield of T₁ oligonucleotides and the hybridization and thermal stability of histone H4 mRNA-DNA hybrids.     

Genetic structure, transforming sequence, and gene product of avian sarcoma virus UR1

We analyzed the genetic structure and gene products of the newly isolated avian sarcoma virus UR1, which recently has been shown to be replication defective and to contain no sequences homologous to the src gene of Rous sarcoma virus. The sizes of the genomic RNAs of UR1 and its associated helper virus, UR1AV, were determined to be 29S and 35S (5.9 and 8.5 kilobases), respectively, by gel electrophoresis and sucrose gradient sedimentation. RNase T1 oligonucleotide mapping of purified viral RNAs indicated that UR1 RNA contains eight unique oligonucleotides in the middle of the genome and shares four 5'-terminal and three 3'-terminal oligonucleotides with UR1AV RNA. The unique sequences of UR1 and Fujinami sarcoma virus were found to be closely related to each other by molecular hybridization of UR1 RNA with DNA complementary to the unique sequence of Fujinami sarcoma virus RNA, but minor differences were found by oligonucleotides fingerprinting. In the regions flanking the unique sequences, UR1 and Fujinami sarcoma viral RNAs contain distinct oligonucleotides, which are shared with oligonucleotides of the respective helper viral RNAs. Cell transformed with UR1 produce a single 29S RNA species which contains a UR1 unique sequence; this species is most likely the mRNA coding for the transforming protein. In UR1-transformed cells, a phosphoprotein fo 150,000 daltons (p150) was detected by immunoprecipitation with antiserum against gag proteins. p150 was associated with a protein kinase activity that was capable of phosphorylating p150 itself, immunoglobulin G of antiserum, and a soluble substrate, alpha-casein. This enzyme transferred phosphate exclusively to tyrosine residues of substrates in vitro, but p 150 labeled in vivo with 32P contained both phosphoserine and phosphotyrosine. The in vitro kinase reaction was not affected by the presence of cyclic AMP or cyclic GMP and strongly preferred Mn2+ over Mg2+. Thus, the properties of UR1 protein are almost identical to those of Fujinami sarcoma virus protein.     

Small RNA sequences are readily stabilized by inclusion in a carrier rRNA

This laboratory previously showed that an RNA derived from 5S ribosomal RNA could be used as a carrier to harbor a nucleic acid "tag" for monitoring genetically engineered or naturally occurring bacteria. The prototype system expressed a specific tagged RNA that was stable and accumulated to high levels. For such a system to be useful there should, however, be little limitation on the sequence composition and length of the insert. To test these limitations, a collection of insertion sequences were created and introduced into the artificial 5S rRNA cassette. This library consisted of random 13- and 50-base oligonucleotides that were inserted into the carrier RNA. We report here that essentially all of the insert-containing RNAs are stable and accumulate to detectable levels. Tagged RNAs were produced by both plasmid-borne and chromosomally integrated expression systems in E. coli and several Pseudomonas strains without obvious effect on the host cell. It is anticipated that in addition to its intended use in environmental monitoring, this system can be used for in vivo selection of useful artificial RNAs. Because the carrier lends stability to the RNAs, the system may also be useful in RNA production.     

Computer analysis of the sequence relationships among 4.5S RNA molecular species from various sources

Nucleotide sequence homology among 4.5S RNAs from various organisms was examined by computer analysis to evaluate their sequence relationships. Chloroplast 4.5S rRNAs of wheat and tobacco were not significantly related to Escherichia coli 4.5S RNA, but were closely related to the 3'-terminus of bacterial 23S rRNA. Significant sequence homology was found between rat Novikoff hepatoma 4.5S RNAI and mouse and hamster 4.5S RNAs, suggesting that these RNAs are products of a family of genes with diverged sequences. E. coli 4.5S RNA had no significant sequence homology with any rodent 4.5S RNAs as a whole sequence. The E. coli, mouse and hamster 4.5S RNAs, however, were found to share a homologous 14-nucleotide sequence at the center of the molecules, which is known to exist as a conserved sequence in both Alu and Alu-equivalent sequences of mammalian DNAs.     

Comparison of ribonucleotide sequences from the genome of vesicular stomatitis virus and two of its defective interfering particles.

RNA genomes from standard vesicular stomatitis virus and two defective interfering (DI) particles dI 0.33 (DI-T) and DI 0.52, were purified and digested with RNase T1. The resulting oligonucleotides were labeled at the 5' end with [32P]ATP and separated by two-dimensional electrophoresis in polyacrylamide gels. All of the major oligonucleotides containing 20 or more nucleotides were sequenced. Those oligonucleotides that were thought to be in common by their migration on polyacrylamide gels actually did have identical sequences. Those oligonucleotides thought to be unique to the DI RNAs either differed by only one nucleotide from oligonucleotides of the standard RNA or contained new sequences which were complementary to known sequences at the 5' end. These data indicate that RNAs from DI particles are not simple deletions but contain point mutations and additional complementary sequences.     

Sequences of the 5S rRNAs of the thermo-acidophilic archaebacterium Sulfolobus solfataricus (Caldariella acidophila) and the thermophilic eubacteria Bacillus acidocaldarius and Thermus aquaticus.

We have determined the nucleotide sequences of the 5 S rRNAs of three thermophilic bacteria: the archaebacterium Sulfolobus solfataricus, also named Caldariella acidophila, and the eubacteria Bacillus acidocaldarius and Thermus aquaticus. A 5 S RNA sequence for the latter species had already been published, but it looked suspect on the basis of its alignment with other 5 S RNA sequences and its base-pairing pattern. The corrected sequence aligns much better and fits in the universal five helix secondary structure model, as do the sequences for the two other examined species. The sequence found for Sulfolobus solfataricus is identical to that determined by others for Sulfolobus acidocaldarius. The secondary structure of its 5 S RNA shows a number of exceptional features which distinguish it not only from eubacterial and eukaryotic 5 S RNAs, but also from the limited number of archaebacterial 5 S RNA structures hitherto published. The free energy change of secondary structure formation is large in the three examined 5 S RNAs.