Secondary structure model for the last two domains of single-stranded RNA phage Q beta

Genome-regulated Assembly of a ssRNA Virus May Also Prepare It for Infection

Many single-stranded, positive-sense RNA viruses regulate assembly of their infectious virions by forming multiple, cognate coat protein (CP)-genome contacts at sites termed Packaging Signals (PSs). We have determined the secondary structures of the bacteriophage MS2 ssRNA genome (gRNA) frozen in defined states using constraints from X-ray synchrotron footprinting (XRF). Comparison of the footprints from phage and transcript confirms the presence of multiple PSs in contact with CP dimers in the former. This is also true for a virus-like particle (VLP) assembled around the gRNA in vitro in the absence of the single-copy Maturation Protein (MP) found in phage. Since PS folds are present at many sites across gRNA transcripts, it appears that this genome has evolved to facilitate this mechanism of assembly regulation. There are striking differences between the gRNA-CP contacts seen in phage and the VLP, suggesting that the latter are inappropriate surrogates for aspects of phage structure/function. Roughly 50% of potential PS sites in the gRNA are not in contact with the protein shell of phage. However, many of these sit adjacent to, albeit not in contact with, PS-binding sites on CP dimers. We hypothesize that these act as PSs transiently during assembly but subsequently dissociate. Combining the XRF data with PS locations from an asymmetric cryo-EM reconstruction suggests that the genome positions of such dissociations are non-random and may facilitate infection. The loss of many PS-CP interactions towards the 3' end of the gRNA would allow this part of the genome to transit more easily through the narrow basal body of the pilus extruding machinery. This is the known first step in phage infection. In addition, each PS-CP dissociation event leaves the protein partner trapped in a non-lowest free-energy conformation. This destabilizes the protein shell which must disassemble during infection, further facilitating this stage of the life-cycle.

Structural Assembly of Qβ Virion and Its Diverse Forms of Virus-like Particles

The coat proteins (CPs) of single-stranded RNA bacteriophages (ssRNA phages) directly assemble around the genomic RNA (gRNA) to form a near-icosahedral capsid with a single maturation protein (Mat) that binds the gRNA and interacts with the retractile pilus during infection of the host. Understanding the assembly of ssRNA phages is essential for their use in biotechnology, such as RNA protection and delivery. Here, we present the complete gRNA model of the ssRNA phage Qβ, revealing that the 3' untranslated region binds to the Mat and the 4127 nucleotides fold domain-by-domain, and is connected through long-range RNA-RNA interactions, such as kissing loops. Thirty-three operator-like RNA stem-loops are located and primarily interact with the asymmetric A/B CP-dimers, suggesting a pathway for the assembly of the virions. Additionally, we have discovered various forms of the virus-like particles (VLPs), including the canonical T = 3 icosahedral, larger T = 4 icosahedral, prolate, oblate forms, and a small prolate form elongated along the 3-fold axis. These particles are all produced during a normal infection, as well as when overexpressing the CPs. When overexpressing the shorter RNA fragments encoding only the CPs, we observed an increased percentage of the smaller VLPs, which may be sufficient to encapsidate a shorter RNA.

Expansion of known ssRNA phage genomes: From tens to over a thousand

The first sequenced genome was that of the 3569-nucleotide single-stranded RNA (ssRNA) bacteriophage MS2. Despite the recent accumulation of vast amounts of DNA and RNA sequence data, only 12 representative ssRNA phage genome sequences are available from the NCBI Genome database (June 2019). The difficulty in detecting RNA phages in metagenomic datasets raises questions as to their abundance, taxonomic structure, and ecological importance. In this study, we iteratively applied profile hidden Markov models to detect conserved ssRNA phage proteins in 82 publicly available metatranscriptomic datasets generated from activated sludge and aquatic environments. We identified 15,611 nonredundant ssRNA phage sequences, including 1015 near-complete genomes. This expansion in the number of known sequences enabled us to complete a phylogenetic assessment of both sequences identified in this study and known ssRNA phage genomes. Our expansion of these viruses from two environments suggests that they have been overlooked within microbiome studies.

Asymmetric cryo-EM structure of the canonical Allolevivirus Qβ reveals a single maturation protein and the genomic ssRNA in situ

Single-stranded (ss) RNA viruses infect all domains of life. To date, for most ssRNA virions, only the structures of the capsids and their associated protein components have been resolved to high resolution. Qβ, an ssRNA phage specific for the conjugative F-pilus, has a T = 3 icosahedral lattice of coat proteins assembled around its 4,217 nucleotides of genomic RNA (gRNA). In the mature virion, the maturation protein, A2, binds to the gRNA and is required for adsorption to the F-pilus. Here, we report the cryo-electron microscopy (cryo-EM) structures of Qβ with and without symmetry applied. The icosahedral structure, at 3.7-Å resolution, resolves loops not previously seen in the published X-ray structure, whereas the asymmetric structure, at 7-Å resolution, reveals A2 and the gRNA. A2 contains a bundle of α-helices and replaces one dimer of coat proteins at a twofold axis. The helix bundle binds gRNA, causing denser packing of RNA in its proximity, which asymmetrically expands the surrounding coat protein shell to potentially facilitate RNA release during infection. We observe a fixed pattern of gRNA organization among all viral particles, with the major and minor grooves of RNA helices clearly visible. A single layer of RNA directly contacts every copy of the coat protein, with one-third of the interactions occurring at operator-like RNA hairpins. These RNA-coat interactions stabilize the tertiary structure of gRNA within the virion, which could further provide a roadmap for capsid assembly.

Exploring the architecture of viral RNA genomes

The genomes of RNA viruses contain local structural elements and long-range interactions that control various steps in virus replication. While many individual RNA elements have been characterized, it remains less clear how the structure and activity of such elements are integrated and regulated within the complex context of complete viral genomes. Recent technical advances, particularly the development of high-throughput solution structure mapping methods, have made secondary structural analysis of entire viral RNA genomes feasible. As a consequence, whole-genome structural models have been deduced for a number of plus-strand RNA viruses and retroviruses and these structures have provided intriguing functional and evolutionary insights into global genome architecture.

In vitro evolution and affinity-maturation with Coliphage qβ display

The Escherichia coli bacteriophage, Qβ (Coliphage Qβ), offers a favorable alternative to M13 for in vitro evolution of displayed peptides and proteins due to high mutagenesis rates in Qβ RNA replication that better simulate the affinity maturation processes of the immune response. We describe a benchtop in vitro evolution system using Qβ display of the VP1 G-H loop peptide of foot-and-mouth disease virus (FMDV). DNA encoding the G-H loop was fused to the A1 minor coat protein of Qβ resulting in a replication-competent hybrid phage that efficiently displayed the FMDV peptide. The surface-localized FMDV VP1 G-H loop cross-reacted with the anti-FMDV monoclonal antibody (mAb) SD6 and was found to decorate the corners of the Qβ icosahedral shell by electron microscopy. Evolution of Qβ-displayed peptides, starting from fully degenerate coding sequences corresponding to the immunodominant region of VP1, allowed rapid in vitro affinity maturation to SD6 mAb. Qβ selected under evolutionary pressure revealed a non-canonical, but essential epitope for mAb SD6 recognition consisting of an Arg-Gly tandem pair. Finally, the selected hybrid phages induced polyclonal antibodies in guinea pigs with good affinity to both FMDV and hybrid Qβ-G-H loop, validating the requirement of the tandem pair epitope. Qβ-display emerges as a novel framework for rapid in vitro evolution with affinity-maturation to molecular targets.

Functional long-range RNA-RNA interactions in positive-strand RNA viruses

Long-range RNA–RNA interactions, many of which span several thousands of nucleotides, have been discovered within the genomes of positive-strand RNA viruses. These interactions mediate fundamental viral processes, including translation, replication and transcription.

In certain plant viruses that have uncapped, non-polyadenylated RNA genomes, translation initiation is facilitated by 3′ cap-independent translational enhancers (3′ CITEs) that are located in or near to their 3′ UTRs. These RNA elements function by binding to either the ribosome-recruiting eukaryotic translation initiation factor 4F (eIF4F) complex or ribosomal subunits, and they enhance translation initiation by engaging the 5′ end of the genome via a 5′-to-3′ RNA-based bridge.

The activities of the internal ribosome entry sites (IRESs) in the 5′ UTRs of various viruses are modulated by RNA-based interactions between the IRESs and elements near to the 3′ ends of their genomes.

In several plant viruses, translational recoding events, including ribosomal frameshifting and stop codon readthrough, have been found to rely on long-range RNA–RNA interactions.

Multiple 5′-to-3′ base-pairing interactions facilitate genome circularization in flaviviruses, which has been proposed to reposition the 5′-bound RNA-dependent RNA polymerase (RdRp) to the initiation site of negative-strand synthesis at the 3′ terminus.

The long-distance interaction between two cis-acting replication elements in tombusviruses generates a bipartite RNA platform for the assembly of the replicase complex and repositions the internally bound RdRp to the 3′ terminus.

Tombusviruses also rely on several long-range interactions that mediate the premature termination of the RdRp during negative-strand synthesis that leads to transcription of subgenomic mRNAs (sgmRNAs).

In a coronavirus, an exceptionally long-range interaction, which spans ∼26,000 nucleotides, promotes polymerase repriming during the discontinuous template synthesis step of sgmRNA-N transcription.

A challenge for the future will be to determine how these long-range interactions are integrated and regulated in the complex context of viral RNA genomes.

Long-range intragenomic RNA–RNA interactions in the genomes of positive-strand RNA viruses involve direct nucleotide base pairing and can span distances of thousands of nucleotides. In this Review, Nicholson and White discuss recent insights into the structure and function of these genomic features and highlight their diverse roles in the gene expression and genome replication of positive-strand RNA viruses.

Positive-strand RNA viruses are important human, animal and plant pathogens that are defined by their single-stranded positive-sense RNA genomes. In recent years, it has become increasingly evident that interactions that occur between distantly positioned RNA sequences within these genomes can mediate important viral activities. These long-range intragenomic RNA–RNA interactions involve direct nucleotide base pairing and can span distances of thousands of nucleotides. In this Review, we discuss recent insights into the structure and function of these intriguing genomic features and highlight their diverse roles in the gene expression and genome replication of positive-strand RNA viruses.

On the role of four small hairpins in the HIV-1 RNA genome

An RNA secondary structure model for the complete HIV-1 genome has recently been published based on SHAPE technology. Several well-known RNA motifs such as TAR and RRE were confirmed and numerous new structured motifs were described that may play important roles in virus replication. The 9 kb viral RNA genome is densely packed with many RNA hairpin motifs and the collective fold may play an important role in HIV-1 biology. We initially focused on 16 RNA hairpin motifs scattered along the viral genome. We considered conservation of these structures, despite sequence variation among virus isolates, as a first indication for a significant function. Four relatively small hairpins exhibited considerable structural conservation and were selected for experimental validation in virus replication assays. Mutations were introduced into the HIV-1 RNA genome to destabilize individual RNA structures without affecting the protein-coding properties (silent codon changes). No major virus replication defects were scored, suggesting that these four hairpin structures do not play essential roles in HIV-1 replication.

Diversity of pili-specific bacteriophages: genome sequence of IncM plasmid-dependent RNA phage M

Background

Bacteriophages of the Leviviridae family are small RNA viruses with linear, positive-sense, single-stranded RNA genomes that encode only four proteins. All phages of this family require bacterial pili to attach to and infect cells. Leviviridae phages utilizing F-pili for this purpose have been extensively studied. RNA phages specific for conjugative plasmid-encoded pili other than that of plasmid F have been isolated, but are much less understood and their relation to the F-pili-specific phages in many cases is not known.

Results

Phage M has the smallest known Leviviridae genome to date and has the typical genome organization with maturation, coat and replicase genes in the 5′ to 3′ direction. The lysis gene is located in a different position than in other known Leviviridae phages and completely overlaps with the replicase gene in a different reading frame. It encodes a 37 residue long polypeptide that contains a transmembrane helix like the other known lysis proteins of leviviruses. Sequence identities of M proteins to those of other phages do not exceed 25% for maturation protein, 51% for coat protein and 41% for replicase. Similarities in protein sequences and RNA secondary structures at the 3′ untranslated region place phage M together with phages specific for IncP, IncC and IncH, but not IncF plasmid-encoded pili. Phylogenetic analysis using the complete genome sequences and replicase proteins suggests that phage M represents a lineage that branched off early in the course of RNA phage specialization on different conjugative plasmids.

Conclusions

The genome sequence of phage M shows that it is clearly related to other conjugative pili-specific leviviruses but has an atypical location of the lysis gene. It provides a better view on the remarkable diversification of the plasmid-specific RNA phages.

Rethinking the evolution of single-stranded RNA (ssRNA) bacteriophages based on genomic sequences and characterizations of two R-plasmid-dependent ssRNA phages, C-1 and Hgal1

We have sequenced and characterized two R-plasmid-dependent single-stranded RNA bacteriophages (RPD ssRNA phages), C-1 and Hagl1. Phage C-1 requires a conjugative plasmid of the IncC group, while Hgal1 requires the IncH group. Both the adsorption rate constants and one-step growth curves are determined for both phages. We also empirically confirmed the lysis function of the predicted lysis genes. Genomic sequencing and phylogenetic analyses showed that both phages belong to the Levivirus group and are most closely related to another IncP-plasmid-dependent ssRNA phage, PRR1. Furthermore, our result strongly suggests that the stereotypical bauplans of genome organization found in Levivirus and Allolevivirus predate phage specialization for conjugative plasmids, suggesting that the utilization of conjugative plasmids for cell attachment and entry comprises independent evolutionary events for these two main clades of ssRNA phages. Our result is also consistent with findings of a previous study, making the Levivirus-like genome organization ancestral and the Allolevivirus-like genome derived. To obtain a deeper insight into the evolution of ssRNA phages, more phages specializing for various conjugative plasmids and infecting different bacterial species would be needed.

The three faces of riboviral spontaneous mutation: spectrum, mode of genome replication, and mutation rate

Riboviruses (RNA viruses without DNA replication intermediates) are the most abundant pathogens infecting animals and plants. Only a few riboviral infections can be controlled with antiviral drugs, mainly because of the rapid appearance of resistance mutations. Little reliable information is available concerning i) kinds and relative frequencies of mutations (the mutational spectrum), ii) mode of genome replication and mutation accumulation, and iii) rates of spontaneous mutation. To illuminate these issues, we developed a model in vivo system based on phage Qß infecting its natural host, Escherichia coli. The Qß RT gene encoding the Read-Through protein was used as a mutation reporter. To reduce uncertainties in mutation frequencies due to selection, the experimental Qß populations were established after a single cycle of infection and selection against RT⁻ mutants during phage growth was ameliorated by plasmid-based RT complementation in trans. The dynamics of Qß genome replication were confirmed to reflect the linear process of iterative copying (the stamping-machine mode). A total of 32 RT mutants were detected among 7,517 Qß isolates. Sequencing analysis of 45 RT mutations revealed a spectrum dominated by 39 transitions, plus 4 transversions and 2 indels. A clear template•primer mismatch bias was observed: A•C>C•A>U•G>G•U> transversion mismatches. The average mutation rate per base replication was ≈9.1×10⁻⁶ for base substitutions and ≈2.3×10⁻⁷ for indels. The estimated mutation rate per genome replication, μ_g, was ≈0.04 (or, per phage generation, ≈0.08), although secondary RT mutations arose during the growth of some RT mutants at a rate about 7-fold higher, signaling the possible impact of transitory bouts of hypermutation. These results are contrasted with those previously reported for other riboviruses to depict the current state of the art in riboviral mutagenesis.

Viral disease is a subject of major concern in public health. Diseases produced by riboviruses (RNA viruses sensu stricto) represent a special urgency, because these viruses display an exceptional capability to generate resistance mutations against antiviral drugs. Unfortunately, little is known about the rate and nature of spontaneous mutation in riboviruses. Thus, characterization of their mutation process may be helpful in the development of improved ways to counteract riboviral diseases. In this study, we investigated the mutation process in vivo of a model ribovirus, the bacteriophage Qß, focusing on three key aspects: i) the kinds and relative frequencies of mutations, ii) the mode of genome replication, and iii) the rate of spontaneous mutation. Our results, combined with other information about riboviral mutagenesis, depict a ribovirus mutation spectrum largely dominated by transitions, a predominantly linear mode of genome replication, and a mutation rate per genome replication on the order of 0.04 for bacteriophages and plant viruses but perhaps an order of magnitude higher for mammalian riboviruses.

Genome structure of caulobacter phage phiCb5

The complete genome sequence of caulobacter phage phiCb5 has been determined, and four open reading frames (ORFs) have been identified and characterized. As for related phages, the ORFs code for maturation, coat, replicase, and lysis proteins, but unlike other Leviviridae members, the lysis protein gene of phiCb5 entirely overlaps with the replicase in a different reading frame. The lysis protein of phiCb5 is about two times longer than that of the distantly related MS2 phage and presumably contains two transmembrane helices. Analysis of the proposed genome secondary structure revealed a stable 5' stem-loop, similar to other phages, and a substantially shorter 3' untranslated region (UTR) structure with only three stem-loops.

Influenza virus RNA structure: unique and common features

The influenza A virus genome consists of eight negative-sense RNA segments. Here we review the currently available data on structure-function relationships in influenza virus RNAs. Various ideas and hypotheses about the roles of influenza virus RNA folding in the virus replication are also discussed in relation to other viruses.

Predicting the sizes of large RNA molecules

We present a theory of the dependence on sequence of the three-dimensional size of large single-stranded (ss) RNA molecules. The work is motivated by the fact that the genomes of many viruses are large ssRNA molecules-often several thousand nucleotides long-and that these RNAs are spontaneously packaged into small rigid protein shells. We argue that there has been evolutionary pressure for the genome to have overall spatial properties-including an appropriate radius of gyration, R(g)-that facilitate this assembly process. For an arbitrary RNA sequence, we introduce the (thermal) average maximum ladder distance (MLD) and use it as a measure of the "extendedness" of the RNA secondary structure. The MLD values of viral ssRNAs that package into capsids of fixed size are shown to be consistently smaller than those for randomly permuted sequences of the same length and base composition, and also smaller than those of natural ssRNAs that are not under evolutionary pressure to have a compact native form. By mapping these secondary structures onto a linear polymer model and by using MLD as a measure of effective contour length, we predict the R(g) values of viral ssRNAs are smaller than those of nonviral sequences. More generally, we predict the average MLD values of large nonviral ssRNAs scale as N(0.67+/-0.01), where N is the number of nucleotides, and that their R(g) values vary as MLD(0.5) in an ideal solvent, and hence as N(0.34). An alternative analysis, which explicitly includes all branches, is introduced and shown to yield consistent results.

Detection of a novel intergenogroup recombinant Norovirus from Kolkata, India

Mutation and recombination are recognized as important driving forces of evolution among RNA viruses. An intergenogroup recombinant norovirus strain [Hu/Kol/NLV/L8775/AB290150/2006/India] was detected in the faecal specimen of a 17 year old male, who had suffered from acute watery diarrhea and severe dehydration. Sequence analysis confirmed that this novel recombinant strain had a polymerase gene fragment that closely resembled a Norovirus (NoV) genogroup-I genotype-3 virus (HuCV/NLV/GI.3/VA98115/AY038598/1998/USA) and a capsid gene resembling NoV genogroup-II genotype-4 virus (NoV/Hu/GII.4/Terneuzen70/EF126964/2006/NL). The crossing over and recombination was observed at nucleotide (nt) 790 of NoV GI VA98115 strain and nt808 of NoV GII Terneuzen70 strain. In both parent strains conserved nucleotide sequence and hairpin structure (DNA secondary structure) were reported at the junction point of ORF1 and ORF2, exhibiting the mechanism of recombination in these viruses. Thus this novel recombinant NoV is another step in evolution among NoVs, indicating that constant surveillance is important to successfully monitor emergence of these strains.

Structure and Function of RNA Replication

Contrary to their host cells, many viruses contain RNA as genetic material and hence encode an RNA-dependent RNA polymerase to replicate their genomes. This review discusses the present status of our knowledge on the structure of these enzymes and the mechanisms of RNA replication. The simplest viruses encode only the catalytic subunit of the replication complex, but other viruses also contribute a variable number of ancillary factors. These and other factors provided by the host cell play roles in the specificity and affinity of template recognition and the assembly of the replication complex. Usually, these host factors are involved in protein synthesis or RNA modification in the host cell, but they play roles in remodeling RNA-RNA, RNA-protein, and protein-protein interactions during virus RNA replication. Furthermore, viruses take advantage of and modify previous cell structural elements, frequently membrane vesicles, for the formation of RNA replication complexes.

A cryptic lysis gene near the start of the Qbeta replicase gene in the +1 frame

The maturation/lysis (A2) protein encoded by the group B single-stranded RNA bacteriophage Qbeta mediates lysis of host Escherichia coli cells. We found a frameshift mutation in the replicase (beta-subunit) gene of Qbeta cDNA causes cell lysis. The mutant has a single base deletion 73 nucleotides (nt) 3' from the start of the replicase gene with consequent translation termination at a stop codon 129-131 nt further 3'. The 43-amino acid C-terminal part of the 67-amino acid product encoded by what in WT (wild-type) is the +1 frame, is rich in basic amino acids This 67-aa protein can mediate cell lysis whose characteristics indicate that the protein may cause lysis by a different mechanism and via a different target, than that caused by the A2 maturation/lysis protein. Synthesis of a counterpart of the newly discovered lysis product in wild-type phage infection would require a hypothetical ribosomal frameshifting event. The lysis gene of group A RNA phages is also short, 75 codons in MS2, and partially overlaps the first part of their equivalently located replicase gene, raising significant evolutionary implications for the present finding.

Nucleotide sequence of a mitochondrial RNA virus from the plant pathogenic fungus, Helicobasidium mompa Tanaka

A double-stranded (ds) RNA (2411 bp) from a strain V18 of the violet root rot basidiomycetous fungus, Helicobasidium mompa was sequenced. Using the fungal mitochondrial genetic code in which UGA codes for tryptophan, the positive strand of V18 dsRNA was found to contain a long open-reading frame with the potential to encode a protein of 700 amino acids (molecular mass 79,805 Da), including conserved motifs characteristic of RNA-dependent RNA polymerase (RDRP). This putative RDRP was shown to be related to putative RDRPs of several fungal mitochondrial viruses. It is proposed that V18 dsRNA is assigned to the genus Mitovirus in the family Narnaviridae and designated as H. mompa mitovirus 1-18 (HmMV1-18). Like other mitoviruses, HmMV1-18 RNA can be folded into potentially stable stem-loop structures at both the 5'- and 3'-termini, and both terminal sequences have inverted complementarity with the potential to form panhandle structure. BLAST analysis indicates that the RDRP encoded by HmMV1-18 is more closely related to those encoded by mitochondrial viruses of some ascomycetes than to that of the unassigned RsM2-1A1 dsRNA in the basidiomycetous Rhizoctonia solani. HmMV1-18 is the first member of the genus Mitovirus from basidiomycete fungi.

Molecular detection and genotyping of male-specific coliphages by reverse transcription-PCR and reverse line blot hybridization

In recent years, there has been increased interest in the use of male-specific or F+ coliphages as indicators of microbial inputs to source waters. Sero- or genotyping of these coliphages can also be used for microbial source tracking (MST). Among the male-specific coliphages, the F+ RNA (FRNA) viruses are well studied, while little is known about the F+ DNA (FDNA) viruses. We have developed a reverse line blot hybridization (RLB) assay which allows for the simultaneous detection and genotyping of both FRNA as well as FDNA coliphages. These assays included a novel generic duplex reverse transcription-PCR (RT-PCR) assay for FRNA viruses as well as a generic PCR for FDNA viruses. The RT-PCR assays were validated by using 190 field and prototype strains. Subsequent DNA sequencing and phylogenetic analyses of RT-PCR products revealed the classification of six different FRNA clusters, including the well-established subgroups I through IV, and three different FDNA clusters, including one (CH) not previously described. Within the leviviruses, a potentially new subgroup (called JS) including strains having more than 40% nucleotide sequence diversity with the known levivirus subgroups (MS2 and GA) was identified. We designed subgroup-specific oligonucleotides that were able to genotype all nine (six FRNA, three FDNA) different clusters. Application of the method to a panel of 351 enriched phage samples from animal feces and wastewater, including known prototype strains (MS2, GA, Q beta, M11, FI, and SP for FRNA and M13, f1, and fd for FDNA), resulted in successful genotyping of 348 (99%) of the samples. In summary, we developed a novel method for standardized genotyping of F+ coliphages as a useful tool for large-scale MST studies.

Evolution of phage with chemically ambiguous proteomes

BACKGROUND

The widespread introduction of amino acid substitutions into organismal proteomes has occurred during natural evolution, but has been difficult to achieve by directed evolution. The adaptation of the translation apparatus represents one barrier, but the multiple mutations that may be required throughout a proteome in order to accommodate an alternative amino acid or analogue is an even more daunting problem. The evolution of a small bacteriophage proteome to accommodate an unnatural amino acid analogue can provide insights into the number and type of substitutions that individual proteins will require to retain functionality.

RESULTS

The bacteriophage Qbeta initially grows poorly in the presence of the amino acid analogue 6-fluorotryptophan. After 25 serial passages, the fitness of the phage on the analogue was substantially increased; there was no loss of fitness when the evolved phage were passaged in the presence of tryptophan. Seven mutations were fixed throughout the phage in two independent lines of descent. None of the mutations changed a tryptophan residue.

CONCLUSIONS

A relatively small number of mutations allowed an unnatural amino acid to be functionally incorporated into a highly interdependent set of proteins. These results support the 'ambiguous intermediate' hypothesis for the emergence of divergent genetic codes, in which the adoption of a new genetic code is preceded by the evolution of proteins that can simultaneously accommodate more than one amino acid at a given codon. It may now be possible to direct the evolution of organisms with novel genetic codes using methods that promote ambiguous intermediates.

3'-Terminal RNA secondary structures are important for accumulation of tomato bushy stunt virus DI RNAs

The plus-strand RNA genome of tomato bushy stunt virus (TBSV) contains a 351-nucleotide (nt)-long 3'-untranslated region. We investigated the role of the 3'-proximal 130 nt of this sequence in viral RNA accumulation within the context of a TBSV defective interfering (DI) RNA. Sequence comparisons between different tombusviruses revealed that the 3' portion of the 130-nt sequence is highly conserved and deletion analysis confirmed that this segment is required for accumulation of DI RNAs in protoplasts. Computer-aided sequence analysis and in vitro solution structure probing indicated that the conserved sequence consists of three stem-loop (SL) structures (5'-SL3-SL2-SL1-3'). The existence of SLs 1 and 3 was also supported by comparative secondary structure analysis of sequenced tombusvirus genomes. Formation of the stem regions in all three SLs was found to be very important, and modification of the terminal loop sequences of SL1 and SL2, but not SL3, decreased DI RNA accumulation in vivo. For SL3, alterations to an internal loop resulted in significantly reduced DI RNA levels. Collectively, these data indicate that all three SLs are functionally relevant and contribute substantially to DI RNA accumulation. In addition, secondary structure analysis of other tombusvirus replicons and related virus genera revealed that a TBSV satellite RNA and members of the closely related genus Aureusvirus (family Tombusviridae) share fundamental elements of this general structural arrangement. Thus, this secondary structure model appears to extend beyond tombusvirus genomes. These conserved 3'-terminal RNA elements likely function in vivo by promoting and/or regulating minus-strand synthesis.

RNA recombination in brome mosaic virus: effects of strand-specific stem-loop inserts

A model system of a single-stranded trisegment Brome mosaic bromovirus (BMV) was used to analyze the mechanism of homologous RNA recombination. Elements capable of forming strand-specific stem-loop structures were inserted at the modified 3' noncoding regions of BMV RNA3 and RNA2 in either positive or negative orientations, and various combinations of parental RNAs were tested for patterns of the accumulating recombinant RNA3 components. The structured negative-strand stem-loops that were inserted in both RNA3 and RNA2 reduced the accumulation of RNA3-RNA2 recombinants to a much higher extent than those in positive strands or the unstructured stem-loop inserts in either positive or negative strands. The use of only one parental RNA carrying the stem-loop insert reduced the accumulation of RNA3-RNA2 recombinants even further, but only when the stem-loops were in negative strands of RNA2. We assume that the presence of a stable stem-loop downstream of the landing site on the acceptor strand (negative RNA2) hampers the reattachment and reinitiation processes. Besides RNA3-RNA2 recombinants, the accumulation of nontargeted RNA3-RNA1 and RNA3-RNA3 recombinants were observed. Our results provide experimental evidence that homologous recombination between BMV RNAs more likely occurs during positive- rather than negative-strand synthesis.

Mutational analysis of the replication signals in the 3'-nontranslated region of citrus tristeza virus

Citrus tristeza virus (CTV), a member of the Closteroviridae, has a 19.3-kb messenger-sense RNA genome consisting of 12 open reading frames with nontranslated regions (NTR) at the 5' and 3' termini. The 273 nucleotide (nt) 3'-NTR is highly conserved ( approximately 95%) among the sequenced CTV isolates in contrast to the highly diverse 5'-NTR sequences. The 3' replication signals were mapped to the 3' 234 nts within the NTR. This region of CTV does not contain a poly-A tract nor does it appear to fold as a tRNA-mimic. Instead, a computer-predicted thermodynamically stable secondary structure comprised of 10 stem-and-loop (SL) structures, referred to as SL1 to SL10 (5' to 3'), was common to all CTV isolates. This putative structure was used as a guide to examine the 3' requirements for replication in vivo. The resulting data suggest that a complex 3' structure is required for those functions that provide for efficient replication of CTV in vivo such as minus-strand initiation, regulation of strand asymmetry, effective translation of the myriad of viral mRNAs, or stability of RNAs. Deletions into the 3'-NTR, up to 66 nts from the 5' direction and 11 nts from the 3' direction, deleting or disrupting putative SL1, SL2 and SL3, or SL10, resulted in continued replication, suggesting that these sequences are not essential for basal-level replication, but are required for efficient replication. Predicted stem loops 3 through 10 were examined by mutations designed to alter the primary structures while preserving the secondary structures. Mutations designed to disrupt the predicted stems of SL3, SL5, SL7, SL9, or SL10 resulted in substantially reduced levels of replication, while compensatory mutations resulted in partial restorations of replication, suggesting that these predicted secondary structures are involved in replication. Also, the putative loop sequences of SL5, SL6, SL7, and SL9 tolerated mutagenesis with continued but reduced levels of replication. In contrast, all mutations introduced into putative SL4, SL8, and the stem of SL6 prevented replication, suggesting that the primary structure of these regions make up the core of the 3' replication signal. The 3' triplet, CCA, was shown to be necessary for efficient replication, but deletion of eleven nts to expose an internal CCA resulted in continued replication.

Nucleotide sequence of a ssRNA phage from Acinetobacter: kinship to coliphages

The complete nucleotide sequence of ssRNA phage AP205 propagating in Acinetobacter species is reported. The RNA has three large ORFs, which code for the following homologues of the RNA coliphage proteins: the maturation, coat and replicase proteins. Their gene order is the same as that in coliphages. RNA coliphages or Leviviridae fall into two genera: the alloleviviruses, like Q(beta), which have a coat read-through protein, and the leviviruses, like MS2, which do not have this coat protein extension. AP205 has no read-through protein and may therefore be classified as a levivirus. A major digression from the known leviviruses is the apparent absence of a lysis gene in AP205 at the usual position, overlapping the coat and replicase proteins. Instead, two small ORFs are present at the 5' terminus, preceding the maturation gene. One of these might encode a lysis protein. The other is of unknown function. Other new features concern the 3'-terminal sequence. In all ssRNA coliphages, there are always three cytosine residues at the 3' end, but in AP205, there is only a single terminal cytosine. Distantly related viruses, like AP205 and the coliphages, do not have significant sequence identity; yet, important secondary structural features of the RNA are conserved. This is shown here for the 3' UTR and the replicase-operator hairpin. Interestingly, although AP205 has the genetic map of a levivirus, its 3' UTR has the length and RNA secondary structure of an allolevivirus. Sharing features with both MS2 and Q(beta) suggests that, in an evolutionary sense, AP205 should be placed between Q(beta) and MS2. A phylogenetic tree for the ssRNA phages is presented.

Autonomous role of 3'-terminal CCCA in directing transcription of RNAs by Qbeta replicase

We have studied transcription in vitro by Qbeta replicase to deduce the minimal features needed for efficient end-to-end copying of an RNA template. Our studies have used templates ca. 30 nucleotides long that are expected to be free of secondary structure, permitting unambiguous analysis of the role of template sequence in directing transcription. A 3'-terminal CCCA (3'-CCCA) directs transcriptional initiation to opposite the underlined C; the amount of transcription is comparable between RNAs possessing upstream (CCA)(n) tracts, A-rich sequences, or a highly folded domain and is also comparable in single-round transcription assays to transcription of two amplifiable RNAs. Predominant initiation occurs within the 3'-CCCA initiation box when a wide variety of sequences is present immediately upstream, but CCA or a closely similar sequence in that position results in significant internal initiation. Removal of the 3'-A from the 3'-CCCA results in 5- to 10-fold-lower transcription, emphasizing the importance of the nontemplated addition of 3'-A by Qbeta replicase during termination. In considering whether 3'-CCCA could provide sufficient specificity for viral transcription, and consequently amplification, in vivo, we note that tRNA(His) is the only stable Escherichia coli RNA with 3'-CCCA. In vitro-generated transcripts corresponding to tRNA(His) served as poor templates for Qbeta replicase; this was shown to be due to the inaccessibility of the partially base-paired CCCA. These studies demonstrate that 3'-CCCA plays a major role in the control of transcription by Qbeta replicase and that the abundant RNAs present in the host cell should not be efficient templates.

Alternative conformations of a nucleic acid four-way junction

A crystal structure of a 108 nucleotide RNA-DNA complex containing a four-way junction was solved at 3.1 A resolution. The structure of the junction differs substantially from the "stacked-X" conformation observed previously, due to a 135 degrees rotation of the branches. Comparison of the two conformers provides insight into the factors contributing to the flexibility of four-way junctions. The stacked-X conformation maximizes base-stacking but causes unfavorable repulsion between phosphate groups, whereas the 135 degrees -rotated "crossed" conformation minimizes electrostatic clashes at the expense of reduced base-stacking. Despite the large rotation of the branches, both junction structures exhibit an antiparallel arrangement of the continuous strands and opposite polarity of the crossover strands.

A long-range pseudoknot in Qbeta RNA is essential for replication

A puzzling aspect of replication of bacteriophage Qbeta RNA has always been that replicase binds at an internal segment, the M-site, some 1450 nt away from the 3' end. Here, we report on the existence of a long-range pseudoknot, base-pairing eight nt in the loop of the 3' terminal hairpin to a single-stranded interdomain sequence located about 1200 nt upstream, close to the internal replicase binding site. Introduction of a single mismatch into this pseudoknot is sufficient to abolish replication, but the inhibition is fully reversed by a second-site substitution that restores the pairing. The pseudoknot is part of an elaborate structure that seems to hold the 3' end in a fixed position vis a vis the replicase binding site. Our results imply that the shape of the RNA confers the functonality. We discuss the possible relevance of our findings for replication of other viral RNAs.

ConStruct: a tool for thermodynamic controlled prediction of conserved secondary structure

A tool for prediction of conserved secondary structure of a set of homologous single-stranded RNAs is presented. For each RNA of the set the structure distribution is calculated and stored in a base pair probability matrix. Gaps, resulting from a multiple sequence alignment of the RNA set, are introduced into the individual probability matrices. These 'aligned' probability matrices are summed up to give a consensus probability matrix emphasizing the conserved structural elements of the RNA set. Because the multiple sequence alignment is independent of any structural constraints, such an alignment may result in introduction of gaps into the homologous probability matrices that disrupt a common consensus structure. By use of its graphical user interface the presented tool allows the removal of such misalignments, which are easily recognized, from the individual probability matrices by optimizing the sequence alignment with respect to a structural alignment. From the consensus probability matrix a consensus structure is extracted, which is viewable in three different graphical representations. The functionality of the tool is demonstrated using a small set of U7 RNAs, which are involved in 3'-end processing of histone mRNA precursors. Supplementary Material lists further results obtained. Advantages and drawbacks of the tool are discussed in comparison to several other algorithms.

FUNCTIONS OF THE 3'-UNTRANSLATED REGIONS OF POSITIVE STRAND RNA VIRAL GENOMES

Positive strand RNA viral genomes are unique in the viral world in serving a dual role as mRNA and replicon. Since the origin of the minus-strand RNA replication intermediate is at the 3'-end of the genome, the 3'-untranslated region (UTR) clearly plays a role in viral RNA replication. The messenger role of this same RNA likely places functional demands on the 3'-UTR to serve roles typical of cellular mRNAs, including the regulation of RNA stability and translation. Current understanding indicates varied roles for positive strand RNA viral 3'-UTRs, with the dominant roles differing between viruses. Three case studies are discussed: turnip yellow mosaic virus RNA, whose 3' tRNA mimicry is thought to negatively regulate minus strand synthesis; brome mosaic virus, whose 3'-UTR contains a unique promoter element directing minus strand synthesis; and tobacco mosaic virus, whose 3'-UTR contains an enhancer of translational expression.

Rotavirus RNA replication requires a single-stranded 3' end for efficient minus-strand synthesis

The segmented double-stranded (ds) RNA genome of the rotaviruses is replicated asymmetrically, with viral mRNA serving as the template for the synthesis of minus-strand RNA. Previous studies with cell-free replication systems have shown that the highly conserved termini of rotavirus gene 8 and 9 mRNAs contain cis-acting signals that promote the synthesis of dsRNA. Based on the location of the cis-acting signals and computer modeling of their secondary structure, the ends of the gene 8 or 9 mRNAs are proposed to interact in cis to form a modified panhandle structure that promotes the synthesis of dsRNA. In this structure, the last 11 to 12 nucleotides of the RNA, including the cis-acting signal that is essential for RNA replication, extend as a single-stranded tail from the panhandled region, and the 5' untranslated region folds to form a stem-loop motif. To understand the importance of the predicted secondary structure in minus-strand synthesis, mutations were introduced into viral RNAs which affected the 3' tail and the 5' stem-loop. Analysis of the RNAs with a cell-free replication system showed that, in contrast to mutations which altered the structure of the 5' stem-loop, mutations which caused complete or near-complete complementarity between the 5' end and the 3' tail significantly inhibited (>/=10-fold) minus-strand synthesis. Likewise, incubation of wild-type RNAs with oligonucleotides which were complementary to the 3' tail inhibited replication. Despite their replication-defective phenotype, mutant RNAs with complementary 5' and 3' termini were shown to competitively interfere with the replication of wild-type mRNA and to bind the viral RNA polymerase VP1 as efficiently as wild-type RNA. These results indicate that the single-strand nature of the 3' end of rotavirus mRNA is essential for efficient dsRNA synthesis and that the specific binding of the RNA polymerase to the mRNA template is required but not sufficient for the synthesis of minus-strand RNA.

The natural 6 S RNA found in Q beta-infected cells is derived from host and phage RNA

The RNA of Escherichia coli infected with RNA bacteriophage Q beta was isolated and screened for replicable short-chained RNA. In contrast to earlier assumptions we show that (i) short-chained replicable RNA is a very minor part of the RNA synthesized in the infection cycle, and (ii) that the replicable RNA isolated from infected cells is derived from cellular RNA, in particular 23 S rRNA and 10 Sa RNA, and from Q beta RNA itself. None of the many RNA species known from in vitro experiments was found. The RNA species isolated were all inefficient templates. No replicable RNA could be isolated from non-infected cells. Even in cells expressing high amounts of Q beta replicase very few RNA species could be isolated. RNA generated in vitro in template-free synthesis is therefore not derived from RNA species found in vivo, and replicable RNA found in vitro is generated by a mechanism fundamentally different from the one operating in vivo.

Rescue of the RNA phage genome from RNase III cleavage

The secondary structure of the RNA from the single-stranded RNA bacteriophages, like MS2 and Qb, has evolved to serve a variety of functions such as controlling gene expression, exposing binding sites for the replicase and capsid proteins, allowing strand separation and so forth. On the other hand, all of these foldings have to perform in bacterial cells in which various RNA splitting enzymes are present. We therefore examined whether phage RNA structure is under selective pressure by host RNases. Here we show this to be true for RNase III. A fully double-stranded hairpin of 17 bp, which is an RNase III target, was inserted into a non-coding region of the MS2 RNA genome. In an RNase III-host these phages survived but in wild-type bacteria they did not. Here the stem underwent Darwinian evolution to a structure that was no longer a substrate for RNase III. This was achieved in three different ways: (i) the perfect stem was maintained but shortened by removing all or most of the insert; (ii) the stem acquired suppressor mutations that replaced Watson-Crick base pairs by mismatches; (iii) the stem acquired small deletions or insertions that created bulges. These insertions consist of short stretches of non-templated A or U residues. Their origin is ascribed to polyadenylation at the site of the RNase III cut (in the + or - strand) either by Escherichia coli poly(A) polymerase or by idling MS2 replicase.

RNA folding kinetics regulates translation of phage MS2 maturation gene

The gene for the maturation protein of the single-stranded RNA coliphage MS2 is preceded by an untranslated leader of 130 nt, which folds into a cloverleaf, i.e., three stem-loop structures enclosed by a long distance interaction (LDI). This LDI prevents translation because its 3' moiety contains the Shine-Dalgarno sequence of the maturation gene. Previously, several observations suggested that folding of the cloverleaf is kinetically delayed, providing a time window for ribosomes to access the RNA. Here we present direct evidence for this model. In vitro experiments show that ribosome binding to the maturation gene is faster than refolding of the denatured cloverleaf. This folding delay appears related to special properties of the leader sequence. We have replaced the three stem-loop structures by a single five nt loop. This change does not affect the equilibrium structure of the LDI. Nevertheless, in this construct, the folding delay has virtually disappeared, suggesting that now the RNA folds faster than ribosomes can bind. Perturbation of the cloverleaf by an insertion makes the maturation start permanently accessible. A pseudorevertant that evolved from an infectious clone carrying the insertion had overcome this defect. It showed a wild-type folding delay before closing down the maturation gene. This experiment reveals the biological significance of retarded cloverleaf formation.

Altered 3'-terminal RNA structure in phage Qbeta adapted to host factor-less Escherichia coli

The RNA phage Qbeta requires for the replication of its genome an RNA binding protein called Qbeta host factor or Hfq protein. Our previous results suggested that this protein mediates the access of replicase to the 3'-end of the Qbeta plus strand RNA. Here we report the results of an evolutionary experiment in which phage Qbeta was adapted to an Escherichia coli Q13 host strain with an inactivated host factor (hfq) gene. This strain initially produced phage at a titer approximately 10,000-fold lower than the wild-type strain and with minute plaque morphology, but after 12 growth cycles, phage titer and plaque size had evolved to levels near those of the wild-type host. RNAs isolated from adapted Qbeta mutants were efficient templates for replicase without host factor in vitro. Electron microscopy showed that mutant RNAs, in contrast to wild-type RNA, efficiently interacted with replicase at the 3'-end in the absence of host factor. The same set of four mutations in the 3'-terminal third of the genome was found in several independently evolved phage clones. One mutation disrupts the base pairing of the 3'-terminal CCCOH sequence, suggesting that the host factor stimulates activity of the wild-type RNA template by melting out its 3'-end.

An oligonucleotide hybridization assay for the identification and enumeration of F-specific RNA phages in surface water

F-specific RNA phages can be used as model organisms for enteric viruses to monitor the effectiveness of sewage treatment, and to assess the potential contamination of surface water with these viruses. In this paper a method is described which identifies RNA phages quantitatively by a plaque hybridization assay. Oligonucleotide probes were developed that can assign phages to their phylogenetic subgroups. Such a distinction is important, since some subgroups preferentially occur in sewage of human origin, while others tend to be associated with animal wastewater. The method has been tested on a large number of isolates and represents an improvement in time and reliability over the previously used serological classification.

Random removal of inserts from an RNA genome: selection against single-stranded RNA

We have monitored the evolution of insertions in two MS2 RNA regions of known secondary structure where coding pressure is negligible or absent. Base changes and shortening of the inserts proceed until the excessive nucleotides can be accommodated in the original structure. The stems of hairpins can be dramatically extended but the loops cannot, revealing natural selection against single-stranded RNA. The 3' end of the MS2 A-protein gene forms a small hairpin with an XbaI sequence in the loop. This site was used to insert XbaI fragments of various sizes. Phages produced by these MS2 cDNA clones were not wild type, nor had they retained the full insert. Instead, every revertant phage had trimmed the insert in a different way to leave a four- to seven-membered loop to the now extended stem. Similar results were obtained with inserts in the 5' untranslated region. The great number of different revertants obtained from a single starting mutant as well as sequence inspection of the crossover points suggest that the removal of redundant RNA occurs randomly. The only common feature among all revertants appears the potential to form a hairpin with a short loop, suggesting that single-stranded RNA negatively affects the viability of the phage. To test this hypothesis, we introduced XbaI fragments of 34 nucleotides that could form either a long stem with a small loop or a short stem with a large loop (26 nucleotides). The base-paired inserts were perfectly maintained for many generations, whereas the unpaired versions were quickly trimmed back to reduce the size of the loop. These data confirm that single-stranded RNA adversely affects phage fitness and is strongly selected against. The repair of the RNA genome that we describe here appears as the result of random recombination. Of the plethora of recombinants, only those able to adopt a base-paired structure survive. The frequency with which our inserts are removed seems higher than measured by others for small inserts in a reading frame in Q beta RNA. To account for this higher frequency, we suggest models in which the single-stranded nature of our inserts induces random recombination at the site of the insertion.

The global folding of four-way helical junctions in RNA, including that in U1 snRNA

Helical junctions are important elements in the architecture of folded RNA molecules. The global geometry of fully base-paired four-way junctions between RNA helices has been analyzed by comparative gel electrophoresis. Junctions appear to fold by pairwise coaxial helical stacking in one of two possible stereochemically equivalent isomers based upon alternative selections of stacking partners. In the presence of 1 mM Mg2+, the two continuous helical axes are approximately at right angles to each other for all junctions studied, but the RNA junctions exhibit significant sequence-dependent differences in their structures as a function of ionic conditions. The four-way junction found in the U1 snRNA folded by coaxial helical stacking. It retained the 90 degrees crossed stacked structure under all ionic conditions tested, despite the presence of a G.A mismatch at the point of strand exchange.