Analysis of five presumptive protein-coding sequences clustered between the primosome genes, 41 and 61, of bacteriophages T4, T2, and T6.

In bacteriophage T4, there is a strong tendency for genes that encode interacting proteins to be clustered on the chromosome. There is 1.6 kb of DNA between the DNA helicase (gene 41) and the DNA primase (gene 61) genes of this virus. The DNA sequence of this region suggests that it contains five genes, designated as open reading frames (ORFs) 61.1 to 61.5, predicted to encode proteins ranging in size from 5.94 to 22.88 kDa. Are these ORFs actually genes? As one test, we compared the DNA sequence of this region in bacteriophages T2, T4, and T6 and found that ORFs 61.1, 61.3, 61.4, and 61.5 are highly conserved among the three closely related viruses. In contrast, ORF 61.2 is conserved between phages T4 and T6 yet is absent from phage T2, where it is replaced by another ORF, T2 ORF 61.2, which is not found in the T4 and T6 genomes. As a second, independent test for coding sequences, we calculated the codon base position preferences for all ORFs in this region that could encode proteins that contain at least 30 amino acids. Both the T4/T6 and T2 versions of ORF 61.2, as well as the other ORFs, have codon base position preferences that are indistinguishable from those of known T4 genes (coefficients of 0.81 to 0.94); the six other possible ORFs of at least 90 bp in this region are ruled out as genes by this test (coefficients less than zero). Thus, both evolutionary conservation and codon usage patterns lead us to conclude that ORFs 61.1 to 61.5 represent important protein-coding sequences for this family of bacteriophages. Because they are located between the genes that encode the two interacting proteins of the T4 primosome (DNA helicase plus DNA primase), one or more may function in DNA replication by modulating primosome function.     

Identifying constraints on the higher-order structure of RNA: continued development and application of comparative sequence analysis methods.

Comparative sequence analysis addresses the problem of RNA folding and RNA structural diversity, and is responsible for determining the folding of many RNA molecules, including 5S, 16S, and 23S rRNAs, tRNA, RNAse P RNA, and Group I and II introns. Initially this method was utilized to fold these sequences into their secondary structures. More recently, this method has revealed numerous tertiary correlations, elucidating novel RNA structural motifs, several of which have been experimentally tested and verified, substantiating the general application of this approach. As successful as the comparative methods have been in elucidating higher-order structure, it is clear that additional structure constraints remain to be found. Deciphering such constraints requires more sensitive and rigorous protocols, in addition to RNA sequence datasets that contain additional phylogenetic diversity and an overall increase in the number of sequences. Various RNA databases, including the tRNA and rRNA sequence datasets, continue to grow in number as well as diversity. Described herein is the development of more rigorous comparative analysis protocols. Our initial development and applications on different RNA datasets have been very encouraging. Such analyses on tRNA, 16S and 23S rRNA are substantiating previously proposed associations and are now beginning to reveal additional constraints on these molecules. A subset of these involve several positions that correlate simultaneously with one another, implying units larger than a basepair can be under a phylogenetic constraint.     

Intraspecific DNA divergence in Drosophila: a study on parthenogenetic D. mercatorum

Drosophila mercatorum is a species that can give rise to totally homozygous parthenogenetic strains. Using the technique of DNA-DNA hybridization, we have assessed the overall single-copy DNA differences among three independently derived strains that represent three independent genomes. Among strains, the average difference between homoduplex and heteroduplex median melting temperatures is 1.3 degrees C. This represents greater than or equal to 1.3% base-pair mismatch. Normalized percent of reassociation indicates further genetic differences, probably reflecting insertion/deletion differences and/or regions of the genome that are highly variable. This overall intraspecific genetic variation is higher than generally is thought to exist but is consistent with growing evidence of extensive DNA diversity within species of invertebrates. High intraspecific DNA variation may be correlated with rapid phyletic rates of evolution. Because of this high level of variation, the technique of DNA-DNA hybridization may be used to study intraspecific variation in invertebrates but is limited in its usefulness for higher systematic studies.      

Quantitative analysis of ribosome binding sites in E.coli.

185 clones with randomized ribosome binding sites, from position -11 to 0 preceding the coding region of beta-galactosidase, were selected and sequenced. The translational yield of each clone was determined; they varied by more than 3000-fold. Multiple linear regression analysis was used to determine the contribution to translation initiation activity of each base at each position. Features known to be important for translation initiation, such as the initiation codon, the Shine/Dalgarno sequence, the identity of the base at position -3 and the occurrence of alternative ATGs, are all found to be important quantitatively for activity. No other features are found to be of general significance, although the effects of secondary structure can be seen as outliers. A comparison to a large number of natural E.coli translation initiation sites shows the information profile to be qualitatively similar although differing quantitatively. This is probably due to the selection for good translation initiation sites in the natural set compared to the low average activity of the randomized set.     

Autogenous regulatory site on the bacteriophage T4 gene 32 messenger RNA

We have identified the binding site on the bacteriophage T4 gene 32 mRNA responsible for autogenous translational regulation. We demonstrate that this site is largely unstructured and overlaps the initiation codon of gene 32 as previously predicted. Co-operative binding of gene 32 protein to this site specifically blocks the formation of 30 S-tRNA(fMet)-gene 32 mRNA ternary complexes and initiation of translation. The translational operator is bound co-operatively by gene 32 protein and this binding is facilitated by a nucleation site far upstream from the initiation codon. A similar unstructured mRNA lacking this nucleation site is also bound co-operatively, but only at concentrations of gene 32 protein higher than those needed to repress binding of ribosomes to the gene 32 mRNA. Some sequence-specific interactions may also influence this binding. Comparison of the bacteriophage T2, T4 and T6 gene 32 operator sequences leads us to propose that the nucleation site is a pseudoknot.     

Translation initiation in Escherichia coli: sequences within the ribosome-binding site

The translational roles of the Shine-Dalgarno sequence, the initiation codon, the space between them, and the second codon have been studied. The Shine-Dalgarno sequence UAAGGAGG initiated translation roughly four times more efficiently than did the shorter AAGGA sequence. Each Shine-Dalgarno sequence required a minimum distance to the initiation codon in order to drive translation; spacing, however, could be rather long. Initiation at AUG was more efficient than at GUG or UUG at each spacing examined; initiation at GUG was only slightly better than UUG. Translation was also affected by residues 3' to the initiation codon. The second codon can influence the rate of initiation, with the magnitude depending on the initiation codon. The data are consistent with a simple kinetic model in which a variety of rate constants contribute to the process of translation initiation.