A conserved motif at the 3$\sp\prime$ end of genomic RNA of mouse hepatitis virus required for host protein binding and viral RNA replication

The initial step in coronavirus-mouse hepatitis virus (MHV) replication is the synthesis of negative strand RNA from a positive strand genomic RNA template. Our approach to studying MHV RNA replication is to identify the cis-acting signals for RNA synthesis and the protein(s) which recognizes these signals at the 3$\sp\prime$ end of genomic RNA of MHV. To determine whether host cellular and/or virus-specific proteins interact with the 3$\sp\prime$ end of the coronavirus genome, an RNase T$\sb1$ protection/gel mobility shift electrophoresis assay was used to examine cytoplasmic extracts from either mock- or MHV-JHM-infected 17Cl-1 murine cells for the ability to form complexes with defined regions of the genomic RNA. A conserved 11 nucleotide sequence UGAAUGAAGUU at nucleotide positions 36 to 26 from the 3$\sp\prime$ end of genomic RNA was identified to be responsible for the specific binding of host proteins, by using a series of RNA probes with deletions and mutations in this region. The RNA probe containing the 11 nucleotide sequence bound approximately four host cellular proteins with a highly labeled 120 kDa and three minor species with sizes of 103, 81 and 55 kDa, assayed by UV-induced covalent cross-linking. Mutation of the 11 nucleotide motif strongly inhibited cellular protein binding, and decreased the amount of the 103 and 81 kDa proteins in the complex to undetectable levels and strongly reduced the binding of the 120 kDa protein. Less extensive mutations within this 11 nucleotide motif resulted in variable decreases in RNA-protein complex formation depending on each probe tested. The RNA-protein complexes observed with cytoplasmic extracts from MHV-JHM-infected cells in both RNase protection/gel mobility shift and UV cross-linking assays were indistinguishable to those observed with extracts from uninfected cells.^ To investigate the possible role of this 3$\sp\prime$ protein binding element in viral RNA replication in vivo, defective interfering RNA molecules with complete or partial mutations of the 11 nucleotide conserved sequence were transcribed in vitro, transfected to host 17Cl-1 cells in the presence of helper virus MHV-JHM and analyzed by agarose gel electrophoresis, competitive RT-PCR and direct sequencing of the RT-PCR products. Both negative strand synthesis and positive strand replication of DI RNA were affected by mutation that disrupts RNA-protein complex formation, even though the 11 mutated nucleotides were converted to wild type sequence, presumably by recombination with helper virus. Kinetic analysis indicated that recombination between DI RNA and helper virus occurred 5.5 to 7.5 hours post infection when replication of positive strand DI RNA was barely observed. Replication of positive strand DI RNAs carrying partial mutations within the 11 nucleotide motif was dependent upon recombination events after transfection. Replication was strongly inhibited when reversion to wild type sequence did not occur, and after recombination, reached similar levels as wild type DI RNA. A DI RNA with mutation upstream of the protein binding motif replicated as efficiently as wild type without undergoing recombination. Thus the conserved 11 nucleotide host protein binding motif appears to play an important role in viral RNA replication. ^

Molecular evolution of primate immunodeficiency viruses and hepatitis delta virus

Primate immunodeficiency viruses, or lentiviruses (HIV-1, HIV-2, and SIV), and hepatitis delta virus (HDV) are RNA viruses characterized by rapid evolution. Infection by primate immunodeficiency viruses usually results in the development of acquired immunodeficiency syndrome (AIDS) in humans and AIDS-like illnesses in Asian macaques. Similarly, hepatitis delta virus infection causes hepatitis and liver cancer in humans. These viruses are heterogeneous within an infected patient and among individuals. Substitution rates in the virus genomes are high and vary in different lineages and among sites. Methods of phylogenetic analysis were applied to study the evolution of primate lentiviruses and the hepatitis delta virus. The following results have been obtained: (1) The substitution rate varies among sites of primate lentivirus genes according to the two parameter gamma distribution, with the shape parameter $\alpha$ being close to 1. (2) Primate immunodeficiency viruses fall into species-specific lineages. Therefore, viral transmissions across primate species are not as frequent as suggested by previous authors. (3) Primate lentiviruses have acquired or lost their pathogenicity several times in the course of evolution. (4) Evidence was provided for multiple infections of a North American patient by distinct HIV-1 strains of the B subtype. (5) Computer simulations indicate that the probability of committing an error in testing HIV transmission depends on the number of virus sequences and their length, the divergence times among sequences, and the model of nucleotide substitution. (6) For future investigations of HIV-1 transmissions, using longer virus sequences and avoiding the use of distant outgroups is recommended. (7) Hepatitis delta virus strains are usually related according to the geographic region of isolation. (8) Evolution of HDV is characterized by the rate of synonymous substitution being lower than the nonsynonymous substitution rate and the rate of evolution of the noncoding region. (9) There is a strong preference for G and C nucleotides at the third codon positions of the HDV coding region. ^

Phenotypic and genotypic analysis of mouse hepatitis virus (JHM) complementation group A temperature-sensitive mutants

Temperature sensitive (ts) mutant viruses have helped elucidate replication processes in many viral systems. Several panels of replication-defective ts mutants in which viral RNA synthesis is abolished at the nonpermissive temperature (RNA$\sp{-})$ have been isolated for Mouse Hepatitis Virus, MHV (Robb et al., 1979; Koolen et al., 1983; Martin et al., 1988; Schaad et al., 1990). However, no one had investigated genetic or phenotypic relationships between these different mutant panels. In order to determine how the panel of MHV-JHM RNA$\sp{-}$ ts mutants (Robb et al., 1979) were genetically related to other described MHV RNA$\sp{-}$ ts mutants, the MHV-JHM mutants were tested for complementation with representatives from two different sets of MHV-A59 ts mutants (Koolen et al., 1983; Schaad et al., 1990). The three ts mutant panels together were found to comprise eight genetically distinct complementation groups. Of these eight complementation groups, three complementation classes are unique to their particular mutant panel; genetically equivalent mutants were not observed within the other two mutant panels. Two complementation groups were common to all three mutant panels. The three remaining complementation groups overlapped two of the three mutant sets. Mutants MHV-JHM tsA204 and MHV-A59 ts261 were shown to be within one of these overlapping complementation groups. The phenotype of the MHV-JHM mutants within this complementation class has been previously characterized (Leibowitz et al., 1982; Leibowitz et al, 1990). When these mutants were grown at the permissive temperature, then shifted up to the nonpermissive temperature at the start of RNA synthesis, genome-length RNA and leader RNA fragments accumulated, but no subgenomic mRNA was synthesized. MHV-A59 ts261 produced leader RNA fragments identical to those observed with MHV-JHM tsA204. Thus, these two MHV RNA$\sp{-}$ ts mutants that were genetically equivalent by complementation testing were phenotypically similar as well. Recombination frequencies obtained from crosses of MHV-A59 ts261 with several of the gene 1 MHV-A59 mutants indicated that the causal mutation(s) of MHV-A59 ts261 was located near the overlapping junction of ORF1a and ORF1b, in the 3$\sp\prime$ end of ORF1a, or the 5$\sp\prime$ end of ORF1b. Sequence analysis of this junction and 1400 nucleotides into the 5$\sp\prime$ end of ORF1b of MHV-A59 ts261 revealed one nucleotide change from the wildtype MHV-A59. This substitution at nucleotide 13,598 (A to G) was a silent mutation in the ORF1a reading frame, but resulted in an amino acid change in ORF1b gene product (I to V). This amino acid change would be expressed only in the readthrough translation product produced upon successful ribosome frameshifting. A revertant of MHV-A59 ts261 (R2) also retained this guanidine residue, but had a second substitution at nucleotide 14,475 in ORF1b. This mutation results in the substitution of valine for an isoleucine.^ The data presented here suggest that the mutation in MHV-A59 ts261 (nucleotide 13,598) would be responsible for the MHV-JHM complementation group A phenotype. A second-site reversion at nucleotide 14,475 may correct this defect in the revertant. Sequencing of gene 1 immediately upstream of nucleotide 13,296 and downstream of nucleotide 15,010 must be conducted to test this hypothesis. ^

Involvement of a region of 23S ribosomal RNA of {\it E. coli\/} in decoding of termination codons

Involvement of E. coli 23S ribosomal RNA (rRNA) in decoding of termination codons was first indicated by the characterization of a 23S rRNA mutant that causes UGA-specific nonsense suppression. The work described here was begun to test the hypothesis that more 23S rRNA suppressors of specific nonsense mutations can be isolated and that they would occur non-randomly in the rRNA genes and be clustered in specific, functionally significant regions of rRNA.^ Approximately 2 kilobases of the gene for 23S rRNA were subjected to PCR random mutagenesis and the amplified products screened for suppression of nonsense mutations in trpA. All of the suppressor mutations obtained were located in a thirty-nucleotide part of the GTPase center, a conserved rRNA sequence and structure, and they and others made in that region by site-directed mutagenesis were shown to be UGA-specific in their suppression of termination codon mutations. These results proved the initial hypothesis and demonstrated that a group of nucleotides in this region are involved in decoding of the UGA termination codon. Further, it was shown that limitation of cellular availability or synthesis of L11, a ribosomal protein that binds to the GTPase center rRNA, resulted in suppression of termination codon mutations, suggesting the direct involvement of L11 in termination in vivo.^ Finally, in vivo analysis of certain site-specific mutations made in the GTPase center RNA demonstrated that (a) the G$\cdot$A base pair closing the hexanucleotide hairpin loop was not essential for normal termination, (b) the "U-turn" structure in the 1093 to 1098 hexaloop is critical for normal termination, (c) nucleotides A1095 and A1067, necessary for the binding to ribosomes of thiostrepton, an antibiotic that inhibits polypeptide release factor binding to ribosomes in vitro, are also necessary for normal peptide chain termination in vivo, and (d) involvement of this region of rRNA in termination is determined by some unique subset structure that includes particular nucleotides rather than merely by a general structural feature of the GTPase center.^ This work advances the understanding of peptide chain termination by demonstrating that the GTPase region of 23S rRNA participates in recognition of termination codons, through an associated ribosomal protein and specific conserved nucleotides and structural motifs in its RNA. ^

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

GENE EXPRESSION IN CHRONIC MYELOGENOUS LEUKEMIA

The goal of the present work was to identify and characterize gene sequences that are preferentially expressed in CML in an effort to better understand the molecular basis of the disease. As high abundance mRNAs generally encode proteins that are phenotypically characteristic of cells, positive-negative screening of a CML cDNA library was used to identify cDNA clones containing sequences preferentially transcribed in CML. One cDNA sequence that fulfilled this criterion, C-A3, has been characterized in some detail. It represents a small mRNA ((TURN)496 nucleotides) that is highly abundant ((TURN)2% of the poly(A('+))RNA) in cells from the chronic phase of CML. In situ hybridization to whole cells indicates the principal leukocytes that express C-A3 sequences are eosinophils, basophils and immature myelocytes. Surprisingly, CML patients with high numbers of myeloblasts do not have an abundance of C-A3 transcripts, although transcript levels remain elevated in patients with lymphoblasts. In AML, high transcript levels are only found sporadically and occasionally different sized transcripts can be detected. Sequences from the 3' end of the C-A3 message are present in 2-5 copies per haploid genome. The 3' end of C-A3 localizes to bands 8q21.1 and 8q23 by in situ chromosomal hybridization. This is a region that is often involved in hematopoietic malignancies. Restriction digests of human genomic DNA show a correlation between the presence of a 2.3 kb Hind III fragment and certain types of leukemia. All of the leukemic DNAs tested had this fragment. In comparison, only one of five normal DNAs had a band this size. Analysis of the nucleotide sequence indicates that C-A3 probably encodes a small, hydrophobic peptide which may be part of a larger protein. ^

AN INVESTIGATION INTO THE MOLECULAR MECHANISMS OF ANDROGEN ACTION IN THE SERTOLI CELL

Steroid hormones regulate target cell function via quantitative and qualitative changes in RNA and protein synthesis. In the testis, androgens are known to play an important role in the regulation of spermatogenesis. The Sertoli cell (SC), whose function is thought to be supportive to the developing germ cell, has been implicated as an androgen target cell. Although cytoplasmic androgen receptors and chromatin acceptor sites for androgen-receptor complexes have been found in SC, effects on RNA synthesis have not previously been demonstrated. In this study, SC RNA synthetic activity was characterized and the effect of testosterone on SC nuclear transcriptional activity in vitro assessed. SC exhibited two fold increases in RNA and ribonucleotide pool concentrations during sexual maturation. These changes appeared to correlate with a previously observed increase in protein concentration per cell over an age span of 15-60 days. Following incubation with ('3)H-uridine, SC from older animals incorporated more label into RNA than SC from younger animals. Since the relative concentration of cytidine nucleotides was higher in SC from older rats, the age-related increase in tritium incorporation may reflect an associated increase in incorporation of ('3)H-CMP into RNA. Alternatively, the enhanced labeling may be the result of either a change in the base composition of the RNA resulting in a higher proportion of CMP and UMP in the RNA, or compartmentalization of the nucleotide pools. The physiologic consequences of these maturational alterations of nucleotide pools remains to be elucidated. RNA polymerase activities were characterized in intact nuclei obtained from cultured rat SC. (alpha)-Amanitin resistant RNA polymerase I+III activity was identical when measured in low or high ionic strength (0.05 M or 0.25 M ammonium sulfate (AS)) in the presence of MnCl(,2) or MgCl(,2), with a divalent cation optimum of 1.6 mM. RNA polymerase II was most active in 0.25 M AS and 1.6 mM MnCl(,2). The apparent Km of RNA polymerase II for UTP was 0.016 mM in 0.05 M AS and 0.037 mM in 0.25 M AS. The apparent Km values for total polymerase activity was 0.008 mM and 0.036 mM at low and high ionic strenghts, respectively. These data indicate that Sertoli cell RNA polymerase activities have catalytic properties characteristic of eukaryotic polymerase activities in general. In the presence of 21 (mu)M testosterone, RNA polymerase II activity increased two fold at 15 minutes, then declined but was still elevated over control values six hours after androgen addition. Polymerase I+III activity was not greatly affected by testosterone. The stimulation of polymerase II measured at 15 minutes was dose-dependent, with a maximum at 0.53 nM and no further stimulation up to 10('-5) M (ED(,50) = 0.25 nM testosterone), and was androgen specific. The results of preliminary RNA isolation and characterization experiments suggested that the synthesis of several species of RNA was enhanced by testosterone administration. These findings have great potential importance since they represent the first demonstration of a direct effect of androgens on the transcriptional process in the Sertoli cell. Furthermore, the results of these studies constitute further evidence that the Sertoli cell is a target for androgen action in the testis. ^

A PHARMACOLOGICAL EVALUATION OF THE MECHANISM OF CYTOTOXICITY OF 9-BETA-D- XYLOFURANOSYLADENINE IN CHINESE HAMSTER OVARY CELLS

The biochemical determinants of cytotoxicity of the purine nucleoside analog, 9-(beta)-D-xylofuranosyladenine (xyl-A) were studied in wild-type Chinese hamster ovary cells and in nucleoside kinase deficient mutants. It was found that {('3)H}xyl-A was readily phosphorylated to the triphosphate level in both the wild-type and deoxycytidine kinase deficient mutant, but not by the adenosine kinase deficient cells. Values for the apparent Km and Vmax of this uptake process were 43.9 (mu)M and 118.7 nmol/min/10('9) cells, respectively. Cloning procedures indicated that the viability of CHO cells was decreased 90 per cent by a 5-hr incubation with 10 (mu)M xyl-A. However, the toxicity of xyl-A was increased 100-fold by the addition of a nontoxic concentration (10 (mu)M) of the adenosine deaminase inhibitor erythro-9-(2-hydroxy-3-nonyl)adenine (EHNA) to the medium. High-pressure liquid chromatographic analysis indicated that after 5 hr, the concentration of 9-(beta)-D-xylofuranosyladenine 5'-triphosphate (xyl-ATP) in cells incubated with xyl-A plus EHNA was 2.0 mM, four times greater than in those cells incubated with xyl-A alone. Incubation with xyl-A plus EHNA had no significant effect on the cellular concentrations of 5-phosphoribosyl-1-pyrophosphate after 1 hr whereas, treatment with 3'-dexoyadenosine (cordycepin) decreased the concentration of this metabolite. Determinations of the cellular nucleoside triphosphates indicated that under conditions that resulted in an intracellular accumulation of 500 (mu)M xyl-ATP, the endogenous concentrations of neither the ribonucleoside triphosphates nor deoxyribonucleoside triphosphates were significantly different from those of control cells. The ID(,50) for {('3)H}thymidine incorporation into DNA, 105 (mu)M xyl-ATP, was four-fold less than the ID(,50) for {('3)H}uridine incorporation into RNA suggesting that the process of DNA synthesis is more sensitive to the presence of xyl-ATP. When removed from exogenous xyl-A, CHO cells failed to recover their ability to synthesize RNA and DNA, although the intracellular xyl-ATP concentration decreased to less than 35 (mu)M. The selective inhibition of RNA synthesis by 6-azauridine did not prevent the expression of toxicity by xyl-ATP. However, the selective inhibition of DNA synthesis by ara-C significantly spared toxicity in cells that had accumulated an otherwise lethal concentration of xyl-ATP. It is shown that in cells which had accumulated 1.27 mM {('3)H}xyl-ATP, {('3)H}xyl-A was found to terminate cellular RNA chains at a frequency of 1.42 (mu)mol of {('3)H}xyl-A 3' termini per mol of mononucleotide. These results indicate that a general mechanism for the toxicity of xyl-A to CHO cells includes the cellular accumulation of xyl-ATP, which serves as a substrate for RNA synthesizing enzymes and subsequently is incorporated into nascent RNA transcripts as a chain terminator. A specific mechanism involving the premature termination of RNA primers required for the initiation of DNA synthesis is proposed to account for the inhibitory action of xyl-ATP on DNA synthesis. ^

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^

GLUCAGON GENE EXPRESSION: ANALYSIS OF PREPROGLUCAGON

Glucagon is a 29 amino acid polypeptide hormone produced in the (alpha) cells of the pancreatic islets. The purpose of this research was to understand better the role of glucagon in the regulation of metabolic processes. As with other polypeptide hormones, the synthesis of glucagon is thought to involve a larger precursor, which is then enzymatically cleaved to the functional form. The specific research objectives were to obtain cloned copies of the messenger RNA (mRNA) for pancreatic glucagon, to determine their primary sequences, and from this coding information to deduce the amino acid sequence of the initial glucagon precursor. From this suggested preproglucagon sequence and prior information on possible proglucagon intermediate processing products, the overall objective of this research is to propose a possible pathway for the biosynthesis of pancreatic glucagon.^ Synthetic oligodeoxynucleotide probes of 14-nucleotides (14-mer) and 17-nucleotides (a 17-mer) complementary to codons specifying a unique sequence of mature glucagon were synthesized. The ('32)P-labeled-14-mer was hybridized with size-fractionated fetal bovine pancreatic poly(A('+))RNA bound to nitrocellulose. RNA fractions of (TURN)14S were found to hybridize specifically, resulting in an (TURN)10-fold enrichment for these sequences. These poly(A('+))RNAs were translated in a cell-free system and the products analyzed by gel electrophoresis. The translation products were found to be enriched for a protein of the putative size of mammalian preproglucagon ((TURN)21 kd). These enriched RNA fractions were used to construct a complementary DNA (cDNA) library is plasmid pBR322.^ Screening of duplicate colony filters with the ('32)P-labeled-17-mer and a ('32)P-labeled-17-mer-primed cDNA probe indicated 25 possible glucagon clones from 3100 colonies screened. Restriction mapping of 6 of these clones suggested that they represented a single mRNA species. Primary sequence analysis of one clone containing a 1200 base pair DNA insert revealed that it contained essentially a full-length copy of glucagon cDNA.^ Analaysis of the cDNA suggested that it encoded an initial translation product of 180 amino acids with an M(,r) = 21 kd. The first initiation codon (ATG, methionine) followed by the longest open reading frame of 540 nucleotides was preceded by a 5'-untranslated region of 90 nucleotides, and was followed by a longer 3'-untranslated region of 471 nucleotides, resulting in a total of 1101 nucleotides. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

GENETIC ORGANIZATION OF THE STRUCTURAL GENES FOR NITRATE REDUCTASE (TN5)

Nitrate reductase in Escherichia coli is a membrane-bound anaerobic enzyme that is repressed by oxygen and induced by nitrate. The genetic organization of the structural genes for the two larger subunits of nitrate reductase ((alpha) and (beta)) was determined by immunoprecipitation analysis of the formation of these proteins in nitrate reductase-deficient mutants resulting from transposon Tn5 mutagenesis. The results suggested that the genes encoding the (alpha) and (beta) subunits (narG and H) were arranged in an operon with transcription in the direction promoter(--->)(alpha)(--->)(beta). Segments of the chromosome containing the Tn5 inserts from several of the mutants were cloned into plasmid pBR322 and the positions of the transposons determined by restriction mapping. The Tn5 insertion sites were localized on two contiguous EcoRI fragments spanning about 6.6 kilobases of DNA. The narI gene (proposed to encode the (gamma) subunit) was positioned immediately downstream from the (beta)-gene (narH) by Southern analysis of Tn10 insertions into the narI locus. A Tn10 insertion into the narK locus, proposed to encode a nitrate-sensitive repressor of other anaerobic enzymes, was located about 1.5 kilobases upstream from the narGHI operon promoter. The narL locus, proposed to encode a nitrate-sensitive positive regulator of the narGHI operon and known to be genetically linked to the other nar genes, was demonstrated to lie outside a 19.3-kilobase region of the chromosome which encompasses the other nar genes. The physical limit of the narGHI promoter was defined by studying the effect of Tn5 insertions into a hybrid plasmid containing the functional operon. The points of origin of the coding regions for the (alpha) and (beta) genes were deduced by alignment of the chromosomal map of Tn5 insertion sites with the sizes of (alpha) and (beta) subunit fragments produced by plasmids carrying these Tn5 inserts in the nar operon. The coding region for the (alpha) subunit (143,000 daltons) begins about 250 nucleotides downstream from the deduced limit of the promoter region and includes about 4.0 kilobases of DNA; the region encoding (beta) (60,000 daltons) lies immediately downstream from the (alpha)-gene and is approximately 1.6 kilobases in length. The adjacent region encoding the (gamma) subunit (19,000 daltons) is approximately 0.5 kilobase in length. ^

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

Mapping and characterization of the Xlsirt P11 RNA cis-acting localization element

In many organisms, polarity of the oocyte is established post-transcriptionally via subcellular RNA localization. Many RNAs are localized during oogenesis in Xenopus laevis, including Xlsirts ( Xenopus laevis short interspersed repeat transcripts) [Kloc, 1993]. Xlsirts constitute a large family defined by highly homologous repeat units 79–81 nucleotides in length. Endogenous Xlsirt RNAs use the METRO (Message Transport Organizer) pathway of localization, where RNAs are transported from the nucleus to the mitochondrial cloud in stage I oocytes. Secondly, RNAs anchor at the vegetal pole in stage II oocytes. Exogenous Xlsirt RNAs can also utilize the Late pathway of localization, which involves localization to the vegetal cortex during stage III of oogenesis and results in RNAs anchored in the cortex of the entire vegetal hemisphere. ^ The Xlsirts localization signal is contained within the repeat region. This study was designed to test the hypothesis that there are cis -acting localization elements in Xlsirts, and that higher order structure plays a role. Results of experiments on Xlsirt P11, a 1700 basepair (bp) family member, led to the conclusion that a 137-bp fragment of the repetitive region is necessary and sufficient for METRO and Late pathway localization. This analysis definitively demonstrates that the Xlsirt localization signal for the METRO and Late pathways reside within the repetitive region and not within the flanking regions. Analysis of Xlsirt linker scanning mutations revealed two METRO-pathway specific subelements, and one Late-pathway specific subelement. Functional, computer, and biochemical evidence relates the higher order structure of this element to its ability to function as a localization element. ^ Xlsirt 137 is 99% identical to the Xlsirt consensus sequence identified in this study, suggesting that it is the localization element for all localized Xlsirt family members. The repeat unit was reframed based on function, rather than arbitrarily based on sequence. This work supports the hypothesis presented in 1981 by George Spohr, who originally isolated the Xlsirts, which stated that the highly conserved repetitive elements must be constrained from variability due to some unknown function of the repeats themselves. These studies shed light on the mechanism of RNA localization, linking structure and function. ^

Integrating sequence information in microarray data analysis by free energy modeling

Microarray technology is a high-throughput method for genotyping and gene expression profiling. Limited sensitivity and specificity are one of the essential problems for this technology. Most of existing methods of microarray data analysis have an apparent limitation for they merely deal with the numerical part of microarray data and have made little use of gene sequence information. Because it's the gene sequences that precisely define the physical objects being measured by a microarray, it is natural to make the gene sequences an essential part of the data analysis. This dissertation focused on the development of free energy models to integrate sequence information in microarray data analysis. The models were used to characterize the mechanism of hybridization on microarrays and enhance sensitivity and specificity of microarray measurements. ^ Cross-hybridization is a major obstacle factor for the sensitivity and specificity of microarray measurements. In this dissertation, we evaluated the scope of cross-hybridization problem on short-oligo microarrays. The results showed that cross hybridization on arrays is mostly caused by oligo fragments with a run of 10 to 16 nucleotides complementary to the probes. Furthermore, a free-energy based model was proposed to quantify the amount of cross-hybridization signal on each probe. This model treats cross-hybridization as an integral effect of the interactions between a probe and various off-target oligo fragments. Using public spike-in datasets, the model showed high accuracy in predicting the cross-hybridization signals on those probes whose intended targets are absent in the sample. ^ Several prospective models were proposed to improve Positional Dependent Nearest-Neighbor (PDNN) model for better quantification of gene expression and cross-hybridization. ^ The problem addressed in this dissertation is fundamental to the microarray technology. We expect that this study will help us to understand the detailed mechanism that determines sensitivity and specificity on the microarrays. Consequently, this research will have a wide impact on how microarrays are designed and how the data are interpreted. ^