Studies on the bacteriophages MS2. XXXIII. Comparison of the nucleotide sequences in related bacteriophage RNAs

Comparing human norovirus surrogates: murine norovirus and Tulane virus

Viral surrogates are widely used by researchers to predict human norovirus behavior. Murine norovirus (MNV) is currently accepted as the best surrogate and is assumed to mimic the survival and inactivation of human noroviruses. Recently, a new calicivirus, the Tulane virus (TV), was discovered, and its potential as a human norovirus surrogate is being explored. This study aimed to compare the behavior of the two potential surrogates under varying treatments of pH (2.0 to 10.0), chlorine (0.2 to 2,000 ppm), heat (50 to 75°C), and survival in tap water at room (20°C) and refrigeration (4°C) temperatures for up to 30 days. Viral infectivity was determined by the plaque assay for both MNV and TV. There was no significant difference between the inactivation of MNV and TV in all heat treatments, and for both MNV and TV survival in tap water at 20°C over 30 days. At 4°C, MNV remained infectious over 30 days at a titer of approximately 5 log PFU/ml, whereas TV titers decreased significantly by 5 days. MNV was more pH stable, as TV titers were reduced significantly at pH 2.0, 9.0, and 10.0, as compared with pH 7.0, whereas MNV titers were only significantly reduced at pH 10.0. After chlorine treatment, there was no significant difference in virus with the exception of at 2 ppm, where TV decreased significantly compared with MNV. Compared with TV, MNV is likely a better surrogate for human noroviruses, as MNV persisted over a wider range of pH values, at 2 ppm of chlorine, and without a loss of titer at 4°C.

Criteria for selection of surrogates used to study the fate and control of pathogens in the environment

This article defines the term surrogate as an organism, particle, or substance used to study the fate of a pathogen in a specific environment. Pathogenic organisms, nonpathogenic organisms, and innocuous particles have been used as surrogates for a variety of purposes, including studies on survival and transport as well as for method development and as "indicators" of certain conditions. This article develops a qualitative surrogate attribute prioritization process and allows investigators to select a surrogate by systematically detailing the experimental process and prioritizing attributes. The results are described through the use of case studies of various laboratories that have used this process. This article also discusses the history of surrogate and microbial indicator use and outlines the method by which surrogates can be used when conducting a quantitative microbial risk assessment. The ultimate goal of selecting a sufficiently representative surrogate is to improve public health through a health-based risk assessment framework. Under- or overestimating the resistance, inactivation, or movement may negatively impact risk assessments that, in turn, will impact health assessments and estimated safety levels. Reducing uncertainty in a risk assessment is one of the objectives of using surrogates and the ultimate motive for any experiment investigating potential exposure of a pathogen.

Reciprocal relationship between stem-loop potential and substitution density in retroviral quasispecies under positive Darwinian selection

Nucleic acids have the potential to form intrastrand stem-loops if complementary bases are suitably located. Computer analyses of poliovirus and retroviral RNAs have revealed a reciprocal relationship between "statistically significant" stem-loop potential and "sequence variability." The statistically significant stem-loop potential of a nucleic acid segment has been defined as a function of the difference between the folding energy of the natural segment (FONS) and the mean folding energy of a set of randomized (shuffled) versions of the natural segment (FORS-M). Since FONS is dependent on both base composition and base order, whereas FORS-M is solely dependent on base composition (a genomic characteristic), it follows that statistically significant stem-loop potential (FORS-D) is a function of base order (a local characteristic). In retroviral genomes, as in all DNA genomes studied, positive FORS-D values are widely distributed. Thus there have been pressures on base order both to encode specific functions and to encode stem-loops. As in the case of DNA genomes under positive Darwinian selection pressure, in HIV-1 specific function appears to dominate in rapidly evolving regions. Here high sequence variability, expressed as substitution density (not indel density), is associated with negative FORS-D values (impaired base-order-dependent stem-loop potential). This suggests that in these regions HIV-1 genomes are under positive selection pressure by host defenses. The general function of stem-loops is recombination. This is a vital process if, from among members of viral "quasispecies," functional genomes are to be salvaged. Thus, for rapidly evolving RNA genomes, it is as important to conserve base-order-dependent stem-loop potential as to conserve other functions.

Nucleotide sequence and mutation rate of the H strain of hepatitis C virus

Patient H is an American patient who was infected with hepatitis C virus (HCV) in 1977. The patient became chronically infected and has remained so for the past 13 years. In this study, we compared the nucleotide and predicted amino acid sequences of the HCV genome obtained from plasma collected in 1977 with that collected in 1990. We find that the two HCV isolates differ at 123 of the 4923 (2.50%) nucleotides sequenced. We estimate that the mutation rate of the H strain of HCV is approximately 1.92 x 10(-3) base substitutions per genome site per year. The nucleotide changes were exclusively base substitutions and were unevenly distributed throughout the genome. A relatively high rate of change was observed in the region of the HCV genome that corresponds to the non-structural protein 1 gene region of flaviviruses, where 44 of 960 (4.6%) nucleotides were different. Within this region there was a 39-nucleotide domain in which 28.2% of the nucleotides differed between the two isolates. In contrast, relatively few nucleotide substitutions were observed in the 5' noncoding region, where only 2 of 276 (0.7%) nucleotides were different. Our results suggest that the mutation rate of the HCV genome is similar to that of other RNA viruses and that genes appear to be evolving at different rates within the virus genome.

Secondary structure of the central region of bacteriophage MS2 RNA. Conservation and biological significance

The RNA of the Escherichia coli RNA phages is highly structured with 75% of the nucleotides estimated to take part in base-pairing. We have used enzymatic and chemical sensitivity of nucleotides, phylogenetic sequence comparison and the phenotypes of constructed mutants to develop a secondary structure model for the central region (900 nucleotides) of the group I phage MS2. The RNA folds into a number of, mostly irregular, helices and is further condensed by several long-distance interactions. There is substantial conservation of helices between the related groups I and II, attesting to the relevance of discrete RNA folding. In general, the secondary structure is thought to be needed to prevent annealing of plus and minus strand and to confer protection against RNase. Superimposed, however, are features required to regulate translation and replication. The MS2 RNA section studied here contains three translational start sites, as well as the binding sites for the coat protein and the replicase enzyme. Considering the density of helices along the RNA, it is not unexpected to find that all these sites lie in helical regions. This fact, however, does not mean that these sites are recognized as secondary structure elements by their interaction partners. This holds true only for the coat protein binding site. The other four sites function in the unfolded state and the stability of the helix in which they are contained serves to negatively control their accessibility. Mutations that stabilize helices containing ribosomal binding sites reduce their efficiency and vice versa. Comparison of homologous helices in different phage RNAs indicates that base substitutions have occurred in such a way that the thermodynamic stability of the helix is maintained. The evolution of individual helices shows several distinct size-reduction patterns. We have observed codon deletions from loop areas and shortening of hairpins by base-pair deletions from either the bottom, the middle or the top of stem structures. Evidence for the coaxial stacking of some helical segments is discussed.

Determination of the RNA secondary structure that regulates lysis gene expression in bacteriophage MS2

The lysis gene of the RNA bacteriophage MS2 is not expressed unless translation of the overlapping coat gene takes place. To understand the molecular basis for this translational coupling the RNA secondary structure around the lysis gene start was analyzed with structure-specific enzymes and chemicals. The existence of a hairpin between nucleotides 1636 and 1707 is in agreement with the structural mapping data and also with the conservation of base-pairing in the related M12 phage. In this hairpin, the G residues in the Shine and Dalgarno region and start codon are inaccessible to RNase T1, which is consistent with the fact that ribosomal access to the lysis gene is blocked when there is no coat gene translation. Deletions or point mutations that are predicted to destabilize the hairpin give rise to lysis protein synthesis that is independent of coat gene translation. Base substitutions that are not expected to weaken the helix do not lead to independent lysis gene expression. Finally, nucleotide changes that strengthen the hairpin lead neither to uncoupled nor to coupled synthesis of the lysis protein. Structural analysis of mutant MS2 RNA shows that small changes in the stability of the secondary structure lead to substantial differences in translation initiation. The function of the hairpin structure in coupling lysis gene to coat gene translation requires that its stability is kept within narrow limits.

Is there a higher level genetic code that directs evolution?

Because the genetic code is redundant for most amino acids, different codons can be used in a given position without altering the structure of the protein for which the gene codes. This flexibility permits information encoding structural, and therefore functional, properties of RNA and DNA to be transmitted simultaneously by a protein-coding sequence of DNA. Among the other messages that might be transmitted, it is proposed, is one modulating the evolution of the DNA itself.

Resolution of the discrepancy between a gene translation--termination codon and the deduced sequence for release of the encoded polypeptide

The translation-termination codon of the synthetase gene of the RNA phage MS2 has been determined, by nucleotide sequencing and suppression studies in vitro, to be UAG. However in one of the only two studies on the signals for polypeptide chain release at the end of genes, Capecchi and Klein [(1970) Nature (Lond.) 226, 1029-1033] deduced that the synthetase of an almost identical phage, R17, is released at UAA. Here we show that under certain conditions the synthetase is released at the UAG terminator but that this UAG is especially prone to read-through with resulting release at the downstream UAA codon. The possible significance of the UAG being in a context prone to leakiness is discussed but is unresolved.

Antigenic drift in influenza virus H3 hemagglutinin from 1968 to 1980: multiple evolutionary pathways and sequential amino acid changes at key antigenic sites

Surveys of the antigenic properties of a wide range of variants of the H3N2 (Hong Kong) influenza virus subtype have revealed complex patterns of variants cocirculating during each of the main epidemic eras of the subtype. We determined hemagglutinin (HA) gene sequences for 14 isolates chosen to give the wildest possible spread of variant types. The addition of these data to existing HA gene sequence information for other variants provides a comprehensive picture of HA gene evolution during antigenic drift among H3N2 subtype viruses. The data reveal the existence of multiple evolutionary pathways during at least one period of development of the subtype and strikingly demonstrate that amino acid changes are limited to a small number of locations on the HA molecule during antigenic drift. The occurrence of sequential amino acid changes at key positions within these variable regions suggests that the HA structure has remained constant during subtype evolution so that only limited possibilities remain for further antigenic drift among H3N2 viruses.

Bacteriophage MS2 lysis protein does not require coat protein to mediate cell lysis

We inserted all but the extreme 5' end of a DNA copy of the bacteriophage MS2 lysis gene downstream of a lac-induced promoter on a multicopy plasmid. Upon induction, cells harboring this plasmid began to lyse, showing that phage coat protein is not required for the lytic process itself.

Silent nucleotide substitutions and the molecular evolutionary clock

Half of the nucleotide substitutions during the evolutionary divergence of genes in animals, bacteria, and viruses are silent changes. These result from an inherent biochemical property of DNA and are fixed by genetic drift. Evolution may be viewed as a device for protecting DNA molecules from extinction.

Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application

A method for estimating the evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences is presented. This method is applied to genes of phi X174 and G4 genomes, histone genes and beta-globin genes, for which homologous nucleotide sequences are available for comparison to be made. It is shown that the rates of synonymous substitutions are quite uniform among the non-overlapping genes of phi X174 and G4 and among histone genes H4, H2B, H3 and H2A. A comparison between phi X174 and G4 reveals that, in the overlapping segments of the A-gene, the rate of synonymous substitution is reduced more significantly than the rate of amino acid substitution relative to the corresponding rate in the non-overlapping segment. It is also suggested that, in the coding region surrounding the splicing points of intervening sequences of beta-globin genes, there exist rigid secondary structures. It is in only these regions that the beta-globin genes show the slowing down of evolutionary rates of both synonymous and amino acid substitutions in the primate line.

Antigenic drift between the haemagglutinin of the Hong Kong influenza strains A/Aichi/2/68 and A/Victoria/3/75

A DNA copy of the gene coding for the influenza A/Aichi/2/68 haemagglutinin protein was cloned in the plasmid pBR322 and the complete nucleotide sequence determined. Comparison of this primary structure and the deduced amino acid sequence with the haemagglutinin gene and protein of strains belonging to the same (H3) subtype and to different subtypes, of both human (H2) and avian (Hav1) origin, documents further at the molecular level the two independent modes of antigenic variation of the virus--drift and shift.

Secondary structure of the 5' end of bacteriophage MS2 RNA Methoxyamine and kethoxal modification

To refine the secondary structure model of the 5' end of the bacteriophage MS2 genome, 32P-labeled MS2 RNA was partially digested with T1 RNase or with Cm-RNase and the 5'-end fragment was isolated, renatured and submitted to treatment with methoxyamine or kethoxal. The resulting modified RNA was digested with T1 RNase and the products were separated by minifingerprinting. Methoxyamine-induced modification of exposed cytidines was detected by differential mobility of modified oligonucleotides, while kethoxal-induced alteration of exposed guanosines was monitored by resistance to T1 ribonuclease digestion. The positions of the modified residues are discussed in terms of an improved secondary structure model proposed for the 5' end of the viral RNA. The structure itself is discussed in relation to sequence conservation and biological function.

Secondary structure of MS2 phage RNA and bias in code word usage

Based on the secondary structural model of MS2 RNA, it is shown that, in base-pairing regions of the RNA, there is a bias in the use of synonymous codons which favours C and/or G over U and/or A in the third codon positions, and that in non-pairing regions, there is an opposite bias which favours U and/or A over C and/or G. This nature is interpreted as a result of selective constraint which stabilises the secondary structure of the single-stranded RNA genome of the MS2 phage.

Normal tRNAs promote ribosomal frameshifting

The addition of Ser AGC AGU tRNA to an E. coli cell-free protein synthesizing system which contains the endogenous tRNA levels results in up to 100% of the ribosomes translating the MS2 coat gene shifting into the -1 reading frame. An analogous phenomenon is seen at a much lower level without the tRNA addition, where a shift into the +1 frame can also be detected. Thus translation with the endogenous tRNA levels yields proteins which have the amino terminus of the coat protein but which are substantially larger than the coat protein and comprise about 5% of the coat translation. Since the lysis gene overlaps the 3' end of the coat gene in the +1 frame, we conclude that the reading frame shift into the +1 frame yields a hybrid protein. Also, we present evidence that ribosomes translating the synthetase gene shift into the -1 frame near the distal end of the gene. This frameshifting is promoted by thrACU ACC tRNA. Specific competitor tRNAs for both Thr and Ser tRNA-promoted frameshifting have been characterized. The generality of this new mechanism for producing additional proteins is unclear, but it investigation should increase understanding of the coding mechanism and its origin.

Bacteriophage MS2 RNA: a correlation between the stability of the codon: anticodon interaction and the choice of code words

The non-random distribution of degenerate code words in Bacteriophage MS2 RNA can be explained partially by considerations of the stability of the codon-anticodon complex in prokaryotic systems. Supporting this hypothesis we note that wobble codons are positively selected in codons having G and/or C in the first two positions. In contrast, wobble codons are statistically less likely in codons composed of A and U in the first two positions. Analyses of nucleotides adjacent to 5' and 3' ends of codons indicate a nonrandom distribution as well. It is thus likely that some elements of RNA evolution are independent of the structural needs of the RNA itself and of the translated protein product.

Codons and nearest-neighbor nucleotide pairs in mammalian messenger RNA

The codons in four mammalian messenger RNAs (rabbit beta hemoglobin, rat pre-proinsulin, rat pre-growth hormone and human chorionic somatomammotropin) show a predominance of C and G in third nucleotide positions. The C:U ratio is about 2 to 1, and the G:A ratio is about 4 to 1. The possibility is discussed that this disproportionally resulted from DNA replicative errors that favor C.G pairs over A.T pairs, as found in the E. coli mut T strain. "Nearest neighbor" base pairs ("doublets") in the protein-coding regions of phiX174 and in four mammalian mRNAs have been compared. Mammalian mRNA has a low content of CpG in comparison with expectations from its C and G content.

Nearest-neighbor doublets in protein-coding regions of MS2 RNA

'Nearest neighbor' base pairs (doublets') in the protein-coding regions of MS2 RNA have been tabulated with respect to their positions in the first two bases of amino acid codons, in the second two bases, or paired by contact between adjoining codons. tconsiderable variation is evident between numbers of doublets in each of these three possible positions, but the totals of each of the 16 doublets in the coding regions of the MS2 RNA molecule show much less varation. Compila-acid composition of proteins coded by such strands.