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. ^

Characterization of osteopontin charge forms produced by osteoblasts

Osteopontin (OPN) is a highly-phosphorylated extracellular matrix protein localized in bone, kidney, placenta, T-lymphocytes, macrophages, smooth muscle of the vascular system, milk, urine, and plasma. In ROS 17/2.8 osteoblast-like osteosarcoma cells, 1,25-dihydroxyvitamin D3 [1,25(OH)2D 3] regulates OPN at the transcriptional level resulting in increased steady state mRNA levels and increased production of OPN protein, maximal at 48 hours. Using ROS 17/2.8 cells as an osteoblast model, OPN was purified from culture medium after three hour treatments of either vehicle (ethanol) or 1,25(OH)2D3 via barium citrate precipitation followed by immunoaffinity chromatography. ^ Here, further evidence of regulation of OPN by 1,25(OH)2D 3 at the posttranslational level is presented. Prior to the up-regulation of OPN at the transcriptional level, 1,25(OH)2D3 induces a shift in OPN isoelectric point (pI) detected on two-dimensional gels from pI 4.6 to pI 5.1. Loading equal amounts of [32P]-labeled OPN recovered from ROS 17/2.8 cells exposed to 1,25(OH)2D3 or vehicle alone for three hours reveals that the shift from pI 4.6 to 5.1 is the result of reduced phosphorylation. Using structural analogs to 1,25(OH) 2D3, analog AT [25-(OH)-16-ene-23-yne-D3], which triggers Ca2+ influx through voltage sensitive Ca2+ channels but does not bind to the vitamin D receptor, mimicked the OPN pI shift while analog BT [1,25(OH)2-22-ene-24-cyclopropyl-D 3], which binds to the vitamin D receptor but does not allow Ca 2+ influx, did not. Inclusion of the Ca2+ channel blocker nifedipine also blocks the charge shift conversion of OPN. Further analysis of the signaling pathway initiated by 1,25(OH)2D3 reveals that inhibition of the cyclic 3′,5′ -adenosine monophosphate-dependent kinase, protein kinase A, or inhibition of the cyclic 3′,5′-guanine monophosphate-dependent kinase, protein kinase G, also prevents the charge shift conversion. ^ Isolation of OPN from rat femurs and tibiae provides evidence for the existence of these two OPN charge forms in vivo, evidenced by differential migration on isoelectric focusing gels and sodium dodecyl sulfate-polyacrylamide gels. Peptide sequencing of rat long bone fractions revealed the presence of a presumed dentin specific protein, dentin matrix protein-1 (DMP-1). Western blot analysis confirmed the existence of DMP-1 in these fractions. ^ Using the OPN charge forms in functional assays, it was determined that the charge forms have differential roles in both cell surface and mineralization functions. In cell attachment assays and Ca2+ influx assays using PC-3 prostate cancer cells, the pI 5.1 charge form of OPN was found to permit binding and increase intracellular Ca2+ concentrations of PC-3 cells. The increase in intracellular Ca2+ concentration was found to be integrin αvβ3-dependent. In mineralization assays, the pI 4.6 charge form of OPN promoted hydroxyapatite formation, while the pI 5.1 charge form had improved Ca2+ binding ability. ^ In conclusion, these findings suggest that 1,25(OH) 2D3 regulates OPN not only at the transcriptional level, but also plays a role in determination of the OPN phosphorylation state. The latter involves a short term (less than three hours) treatment and is associated with membrane-initiated Ca2+ influx. Functional assays utilizing the two OPN charge forms reveal the dependence of OPN post-translational state on its function. ^

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. ^

IN VITRO METABOLISM OF 6-THIOGUANINE IN RAT SUBCELLULAR FRACTIONS

The metabolism of the antitumor agent 6-thioguanine (TG, NSC-752) by rat liver was studied in vitro. Livers from adult male Sprague-Dawley rats were homogenized and the "liver homogenate" was subjected to differential centrifugation to obtain the "10,000 x g pellet", the "post-mitochondrial fraction", the "cytosol fraction", and the "microsomes". The homogenity of each fraction was estimated by appropriate marker enzyme assays. To delineate the in vitro metabolism of TG by rat liver, 0.2 mM of {8-('14)C}TG was incubated with different subcellular fractions in KCl-Tris-MgCl(,2) buffer, pH 7.4 at 37(DEGREES). The metabolites formed were identified by chromatography, UV spectrometry, as well as mass spectrometry. After a 1 hr incubation, TG was metabolized by the liver homogenate, the 10,000 x g pellet and the post-mitochondrial fraction mainly to 6-thioguanosine (TGR), accompanied by varying lesser amounts of 6-thiouric acid (TUA), allantoin, guanine-6-sulfinic acid (G-SO(,2)H) and an unknown product. In comparison, the cytosal fraction converted TG almost entirely to TGR and TUA in equal amounts. The formation of TGR from TG was limited by the endogenous supply of ribose-1-phosphate. With the microsomal fraction, however, TG was metabolized significantly to G-SO(,2)H and the unknown, accompanied with some TGR. After a 5 hr incubation the metabolism of TG was changed to favor the catabolic route, yielding mostly TUA in the post-mitochondrial and cytosol fractions; but mainly allantoin in the liver homogenate fraction. The kinetic studies of TG metabolism by the subcellar fractions indicated that the formation of TGR served as a depot form of TG. The level of TGR decreased when the catabolism of TG became prominent. The oxidation of TG to GSO(,2)H mediated by the hepatic microsomes represented a new catabolic pathway of TG. This GSO(,2)H, under acidic conditions, readily decomposes to guanine and inorganic sulfate. In the presence of reduced glutathione in Tris buffer, pH 7.8 at 25(DEGREES), GSO(,2)H is adducted to glutathione chemically to form S-(2-amino-purin-6-yl) glutathione and conceivably, inorganic sulfate. Therefore, the formation of GSO(,2)H from TG might have implication in the desulfuration mechanism of TG. On the other hand, the unknown formed from TG by the action of the microsomal enzymes appeared to be a TG conjugate. However, it is neither a glutathione, a glucuronide, nor a ribose conjugate. Additionally, the deamination of TG by guanine deaminase (E.C.3.5.4.3) isolated from rat liver was also investigated. TG is a poorer substrate (Km = 4.8 x 10('-3)M) for guanine deaminase than that of guanine (Km = 4.7 x 10('-6)M) at pH 7.25, optimal pH for TG as a substrate. TG is also a competitive inhibitor of guanine for guanine deaminase, with a ki of 2.2 x 10('-4)M. ^

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.) ^

REGULATION OF BRAIN NORADRENERGIC-RECEPTOR FUNCTION BY ADRENOCORTICOTROPHIN AND ANTIDEPRESSANT DRUGS (ADRENERGIC RECEPTOR, CYCLIC-AMP, ADENYLATE CYCLASE, IMIPRAMINE, STRESS)

Human behavior appears to be regulated in part by noradrenergic transmission since antidepressant drugs modify the number and function of (beta)-adrenergic receptors in the central nervous system. Affective illness is also known to be associated with the endocrine system, particularly the hypothalamic-pituitary-adrenal axis. The aim of the present study was to determine whether hormones, in particular adrencorticotrophin (ACTH) and corticosterone, may influence behavior by regulating brain noradrenergic receptor function.^ Chronic treatment with ACTH accelerated the increase or decrease in rat brain (beta)-adrenergic receptor number induced by a lesion of the dorsal noradrenergic bundle or treatment with the antidepressant imipramine. Chronic administration of ACTH alone had no effect on (beta)-receptor number although it reduced norepinephrine stimulated cyclic AMP accumulation in brain slices. Treatment with imipramine also reduced the cyclic AMP response to norepinephrine but was accompanied by a decrease in (beta)-adrenergic receptor number. Both the imipramine and ACTH treatments reduced the affinity of (beta)-adrenergic receptors for norepinephrine, but only the antidepressant modified the potency of the neurotransmitter to stimulate second messenger production. Neither ACTH nor imipramine treatment altered Gpp(NH)p- or fluoride-stimulated adenylate cyclase, cyclic AMP, cyclic GMP, or cyclic GMP-stimulated cyclic AMP phosphodiesterase, or the activity of the guanine nucleotide binding protein (Gs). These findings suggested that post-receptor components of the cyclic nucleotide generating system are not influenced by the hormone or antidepressant. This conclusion was verified by the finding that neither treatment altered adenosine-stimulated cyclic AMP accumulation in brain tissue.^ A detailed examination of the (alpha)- and (beta)-adrenergic receptor components of norepinephrine-stimulated cyclic AMP production revealed that ACTH, but not imipramine, administration reduced the contribution of the (alpha)-receptor mediated response. Like ACTH treatment, corticosterone diminished the (alpha)-adrenergic component indicating that adrenal steroids probably mediate the neurochemical responses to ACTH administration. The data indicate that adrenal steroids and antidepressants decrease noradrenergic receptor function by selectively modifying the (alpha)- and (beta)-receptor components. The functional similarity in the action of the steroid and antidepressants suggests that adrenal hormones normally contribute to the maintenance of receptor systems which regulate affective behavior in man. ^

Molecular alterations in nucleotide excision repair genes in malignant glioma and their relevance to malignant progression and drug resistance

Recently, it has become apparent that DNA repair mechanisms are involved in the malignant progression and resistance to therapy of gliomas. Many investigators have shown that increased levels of O6-methyl guanine DNA alkyltransferase, a DNA monoalkyl adduct repair enzyme, are correlated with resistance of malignant glioma cell lines to nitrosourea-based chemotherapy. Three important DNA excision repair genes ERCC1 (excision repair cross complementation group 1), ERCC2 (excision repair cross complementation group 2), and ERCC6 (excision repair cross complementation group 6) have been studied in human tumors. Gene copy number variation of ERCC1 and ERCC2 has been observed in primary glioma tissues. A number of reports describing a relationship between ERCC1 gene alterations and resistance to anti-cancer drugs have been also described. The levels of ERCC1 gene expression, however, have not been correlated with drug resistance in gliomas. The expression of ERCC6 gene transcribes has been shown to vary with tissue types and to be highest in the brain. There have been no comprehensive studies so far, however, of ERCC6 gene expression and molecular alterations in malignant glioma. This project examined the ERCC1 expression levels and correlated them with cisplatin resistance in malignant glioma cell lines. We also examined the molecular alterations of ERCC6 gene in primary glioma tissues and cells and analyzed whether these alterations are related to tumor progression and chemotherapy resistance. Our results indicate the presence of mutations and/or deletions in exons II and V of the ERCC6 gene, and these alterations are more frequent in exon II. Furthermore, the mutations and/or deletions in exon II were shown to be associated with increased malignant grade of gliomas. The results on the Levels of ERCC1 gene transcripts showed that expression levels correlate with cisplatin resistance. The increase in ERCC1 mRNA induced by cisplatin could be down-regulated by cyclosporin A and herbimycin A. The results of this study are likely to provide useful information for clinical treatment of human gliomas. ^