Systematic analysis of coding and noncoding DNA sequences using methods of statistical linguistics

We compare the statistical properties of coding and noncoding regions in eukaryotic and viral DNA sequences by adapting two tests developed for the analysis of natural languages and symbolic sequences. The data set comprises all 30 sequences of length above 50 000 base pairs in GenBank Release No. 81.0, as well as the recently published sequences of C. elegans chromosome III (2.2 Mbp) and yeast chromosome XI (661 Kbp). We find that for the three chromosomes we studied the statistical properties of noncoding regions appear to be closer to those observed in natural languages than those of coding regions. In particular, (i) a n-tuple Zipf analysis of noncoding regions reveals a regime close to power-law behavior while the coding regions show logarithmic behavior over a wide interval, while (ii) an n-gram entropy measurement shows that the noncoding regions have a lower n-gram entropy (and hence a larger "n-gram redundancy") than the coding regions. In contrast to the three chromosomes, we find that for vertebrates such as primates and rodents and for viral DNA, the difference between the statistical properties of coding and noncoding regions is not pronounced and therefore the results of the analyses of the investigated sequences are less conclusive. After noting the intrinsic limitations of the n-gram redundancy analysis, we also briefly discuss the failure of the zeroth- and first-order Markovian models or simple nucleotide repeats to account fully for these "linguistic" features of DNA. Finally, we emphasize that our results by no means prove the existence of a "language" in noncoding DNA.     

Asymmetry of coding versus noncoding strand in coding sequences of different genomes

We have used the asymmetry between the coding and noncoding strands in different codon positions of coding sequences of DNA as a parameter to evaluate the coding probability for open reading frames (ORFs). The method enables an approximation of the total number of coding ORFs in the set of analyzed sequences as well as an estimation of the coding probability for the ORFs. The asymmetry observed in the nucleotide composition of codons in coding sequences has been used successfully for analysis of the genomes completed at the time of this analysis.     

Triplet correlation in DNA sequences and stability of heteroduplexes

Rhodococcus equi was isolated from the lungs of six foals with bronchopneumonia. All isolates expressed 15-17-kd antigens by immunoblot analysis and contained a virulence-associated plasmid of 85 or 90 kb. Immunohistochemically, R. equi from all pulmonary lesions showed the expression of 15-17-kd antigens mainly in the phagocytic cells. The specific monoclonal antibody to 15-17-kd antigens of R. equi (MAb 10G5) may be an aid in the diagnosis of R. equi-induced pneumonia.     

DNA sequences coding for the H2B histone of Psammechinus miliaris

Starting from restriction enzyme cleavage sites of known topologies the DNA sequences coding for two thirds of the H2B histone protein, together with some 3 inch extracistronic sequences, have been determined for the sea urchin Psammechinus miliaris. This unambiguously identifies this gene and further reveals its 5 inch-3inch polarity and accurate map position in a cloned histone DNA repeat unit.     

Ribosomal S27a coding sequences upstream of ubiquitin coding sequences in the genome of a pestivirus

Molecular characterization of cytopathogenic (cp) bovine viral diarrhea virus (BVDV) strain CP Rit, a temperature-sensitive strain widely used for vaccination, revealed that the viral genomic RNA is about 15.2 kb long, which is about 2.9 kb longer than the one of noncytopathogenic (noncp) BVDV strains. Molecular cloning and nucleotide sequencing of parts of the genome resulted in the identification of a duplication of the genomic region encoding nonstructural proteins NS3, NS4A, and part of NS4B. In addition, a nonviral sequence was found directly upstream of the second copy of the NS3 gene. The 3' part of this inserted sequence encodes an N-terminally truncated ubiquitin monomer. This is remarkable since all described cp BVDV strains with ubiquitin coding sequences contain at least one complete ubiquitin monomer. The 5' region of the nonviral sequence did not show any homology to cellular sequences identified thus far in cp BVDV strains. Databank searches revealed that this second cellular insertion encodes part of ribosomal protein S27a. Further analyses included molecular cloning and nucleotide sequencing of the cellular recombination partner. Sequence comparisons strongly suggest that the S27a and the ubiquitin coding sequences found in the genome of CP Rit were both derived from a bovine mRNA encoding a hybrid protein with the structure NH2-ubiquitin-S27a-COOH. Polyprotein processing in the genomic region encoding the N-terminal part of NS4B, the two cellular insertions, and NS3 was studied by a transient-expression assay. The respective analyses showed that the S27a-derived polypeptide, together with the truncated ubiquitin, served as processing signal to yield NS3, whereas the truncated ubiquitin alone was not capable of mediating the cleavage. Since the expression of NS3 is strictly correlated with the cp phenotype of BVDV, the altered genome organization leading to expression of NS3 most probably represents the genetic basis of cytopathogenicity of CP Rit.     

Cloning, characterization, and properties of seven triplet repeat DNA sequences

Several neuromuscular and neurodegenerative diseases are caused by genetically unstable triplet repeat sequences (CTG.CAG, CGG.CCG, or AAG.CTT) in or near the responsible genes. We implemented novel cloning strategies with chemically synthesized oligonucleotides to clone seven of the triplet repeat sequences (GTA.TAC, GAT.ATC, GTT.AAC, CAC.GTG, AGG.CCT, TCG.CGA, and AAG.CTT), and the adjoining paper (Ohshima, K., Kang, S., Larson, J. E., and Wells, R. D.(1996) J. Biol. Chem. 271, 16784-16791) describes studies on TTA.TAA. This approach in conjunction with in vivo expansion studies in Escherichia coli enabled the preparation of at least 81 plasmids containing the repeat sequences with lengths of approximately 16 up to 158 triplets in both orientations with varying extents of polymorphisms. The inserts were characterized by DNA sequencing as well as DNA polymerase pausings, two-dimensional agarose gel electrophoresis, and chemical probe analyses to evaluate the capacity to adopt negative supercoil induced non-B DNA conformations. AAG.CTT and AGG.CCT form intramolecular triplexes, and the other five repeat sequences do not form any previously characterized non-B structures. However, long tracts of TCG.CGA showed strong inhibition of DNA synthesis at specific loci in the repeats as seen in the cases of CTG.CAG and CGG.CCG (Kang, S., Ohshima, K., Shimizu, M., Amirhaeri, S., and Wells, R. D.(1995) J. Biol. Chem. 270, 27014-27021). This work along with other studies (Wells, R. D.(1996) J. Biol. Chem. 271, 2875-2878) on CTG.CAG, CGG.CCG, and TTA.TAA makes available long inserts of all 10 triplet repeat sequences for a variety of physical, molecular biological, genetic, and medical investigations. A model to explain the reduction in mRNA abundance in Friedreich's ataxia based on intermolecular triplex formation is proposed.     

The coding and effector transfer of movement sequences.

Three experiments utilizing a 14-element arm movement sequence were designed to determine if reinstating the visual–spatial coordinates, which require movements to the same spatial locations utilized during acquisition, results in better effector transfer than reinstating the motor coordinates, which require the same pattern of homologous muscle activation. Results demonstrated better transfer when visual–spatial coordinates were reinstated than when motor coordinates where reinstated regardless of the amount of practice (1, 4, or 12 days; Experiments 1–3, respectively). Transfer (left to right and right to left) was symmetric when visual–spatial coordinates were reinstated but not when motor coordinates were reinstated. When motor coordinates were reinstated after 12 days of practice and vision occluded, transfer was better from right limb to left than vice versa. The data are also consistent with the notion that multiple codes (visual, spatial, and motor) are developed over practice, with each contributing to transfer performance when the respective coordinates are reinstated. Further, the results indicate a disruption of the linkage (concatenation) between subsequences when one or more coordinates are changed on the transfer test. (PsycINFO Database Record (c) 2010 APA, all rights reserved)     

Kinetic analysis of the coding properties of O6-methylguanine in DNA: the crucial role of the conformation of the phosphodiester bond

Production by N-nitroso compounds of O6-alkylguanine (O6-alkylG) in DNA directs the misincorporation of thymine during DNA replication, leading to G:C to A:T transition mutations, despite the fact that DNA containing O6-alkylG:T base pairs is less stable than that containing O6-alkylG:C pairs. We have examined the kinetics of incorporation by Klenow fragment (KF) of Escherichia coli DNA polymerase I of thymine (T) and of cytosine (C) opposite O6-MeG in the template DNA strand. Both T and C were incorporated opposite O6-MeG much slower than nucleotides forming regular A:T or G:C base pairs. Using various concentrations of dTTP, dCTP, or their phosphorothioate (Sp)-dNTP alpha S analogues, or a mixture of dTTP and dCTP, the progress of incorporation of a single nucleotide in a single catalytic cycle of a preformed KF-DNA complex was measured (pre-steady-state kinetics). The results were consistent with the kinetic scheme (Kuchta, R. D., Benkovic, P., & Benkovic, S. J. (1988) Biochemistry 27, 6716-6725): (1) binding of dNTP to polymerase-DNA; (2) conformational change in polymerase; (3) formation of phosphodiester between the dNTP and the 3'-OH of the primer; (4) conformational change of polymerase; (5) release of pyrophosphate. The results were analyzed mathematically to identify the steps at which the rate constants differ significantly between the incorporation of T and C. The only significant difference was the 5-fold difference in the rates of formation of the phosphodiester bond (for dTTP, kforward = 3.9 s-1 and kback = 1.9 s-1; for dCTP, kforward = 0.7 s-1 and kback = 0.9 s-1). These pre-steady-state progress curves were biphasic with a rapid initial burst followed by an apparently steady-state rise. Deconvolution of these curves gave direct evidence for the importance of the conformational change after polymerization by showing that the curves represented the sum of the rapid accumulation of the product of step 3 followed by the slow conversion of that to the product of step 5 (because of the rapidity of the release of pyrophosphate there was no significant accumulation of the product of step 4). The equilibrium constants for each step suggest that the greatest change in the Gibbs free energy occurs at the conformational change after polymerization and that while the formation of the phosphodiester bond to T is slightly exothermic, that to C is slightly endothermic.(ABSTRACT TRUNCATED AT 400 WORDS)     

Individual DNA bands obtained by RAPD analysis of canine genomic DNA often contain multiple DNA sequences

We present a case of Langerhans' cell histiocytosis (LCH) of the liver and spleen in an adult. The imaging features are different from those in the few previously reported cases of individual organ involvement by LCH.     

Endogenous oncornaviral DNA sequences: evidence for two classes of viral DNA sequences in guinea pig cells

The nature of the endogenous viral DNA sequences in guinea pig cells was studied by hybridization. A segment of the viral RNA (r-VRNA) hybridizing to abundant (or reiterated) DNA sequences (R-VDNA) was isolated by recycling to a Cot of 300. The hybridization of the recycled VRNA, as well as the total VRNA, was followed by determining their kinetics and by Wetmur-Davidson analysis. The kinetics of hybridization of total VRNA were complex, did not follow a second-order kinetics, and revealed two slopes by Wetmur-Davidson analysis. The recycled RNA, on the other hand, had a second-order reaction rate expected of the hybridization between a single species of RNA and DNA sequences and yielded a single straight line in a Wetmur-Davidson plot. The Cot1/2 and slope of the recycled r-VRNA was almost identical to that of the abundant VDNA sequences obtained from the hybridization data of the total VRNA. Guinea pig 28S rRNA with or without recycling was used in monitoring hybridization rate. The kinetics of hybridization of 28S RNA followed a second-order reaction and produced a single straight line by Wetmur-Davidson plot, with a second-order reassociation rate constant of 9.6 x 10(-3) liters/mol-s, a Cot1/2 of 104 mol-s/liter, and reiteration frequency of 146. There was no difference in the kinetics of hybridization of 28S RNA before and after recycling. These experiments showed that guinea pig cells contain two classes of VDNA sequences. (i) R-VDNA sequences with a second-order reassociation rate constant of 8.2 x 10(-4) liters/mol-s, a Cot1/2 of 1,219 mol-s/liter, and a reiteration frequency of 12 represent 37.5% of the viral genome. (ii) Unique VDNA sequences with a second-order reassociation rate constant of 1.2 x 10(-4) liters/mol-s, a Cot1/2 of 7,692 mol-s/liter, and a reiteration frequency of 2 represent 62.5% of the viral genome.     

Mutation analysis of coding sequences of the entire transforming growth factor beta type II receptor gene in sporadic human colon cancer using genomic DNA and intron primers

Mutations in the transforming growth factor beta type II receptor (TGFbeta RII) gene have been detected in several human cancers exhibiting microsatellite instability. To extend analyses of this gene, we previously investigated the exon-intron organization of the TGFbeta RII gene and defined seven exons and flanking intron sequences. In this study, we further determined the nucleotide sequences surrounding these seven exons and designed eight sets of intron-based primers to examine the entire coding region of the TGFbeta RII gene. Using these primers, we screened genomic DNA sequences from 30 sporadic colorectal cancers for mutations of the TGFbeta RII gene. TGFbeta RII mutations were detected in two of 30 tumors and both displayed microsatellite instability. One had a deletion in a polyadenine tract in exon 3 and the other had a point mutation in the kinase domain located in exon 7. There were no mutations in exons 1, 2, 4, 5 and 6. These results further implicate the polyadenine tract and kinase domain as mutational hotspots in the TGFbeta RII gene in cells with genomic instability and suggest that TGFbeta RII gene mutations occur rarely in cells lacking genomic instability.     

Rapid genomic changes in newly synthesized amphiploids of Triticum and Aegilops. II. Changes in low-copy coding DNA sequences

Postischemic myocardium possesses considerable contractile and metabolic reserves, but their mobilization could result in increased cell death. We tested the hypothesis that beta-adrenergic stimulation of reperfused myocardium would increase segment work more than O2 consumption, thereby improving efficiency without increased cell death. In 16 open-chest anesthetized dogs, the left anterior descending coronary artery (LAD) was ligated for 2 h; during the reperfusion period, isoproterenol (ISO; 0.1 microg/kg/min, i.v.) was administered to nine of the animals. Regional myocardial segment length and force were measured in the anterior (LAD) and posterior circumflex coronary artery (CFX) regions of the left ventricular myocardium. Work was calculated as the integrated products of force and shortening for each region. Regional myocardial O2 consumption was obtained from LAD flow and arterial and local venous O2 saturations. Infarct size (tetrazolium) was measured in the treated and untreated hearts at the end of the experiment. In untreated hearts, the first derivative of left ventricular pressure, cardiac output, and external work were significantly depressed during reperfusion; ISO restored all values to preocclusion levels. Regional myocardial work in both LAD and CFX regions was significantly increased by ISO (from 564 +/- 207 to 1,635 +/- 543 g/mm/min in LAD, and from 753 +/- 90 to 1,426 +/- 245 g/mm/min in CFX). Efficiency (work/oxygen consumption) of the reperfused region was similarly increased. LAD flow was significantly increased by ISO, and O2 extraction was unchanged. Infarct size was 28.2 +/- 4.7% in untreated hearts and 29.0 +/- 3.5% in ISO hearts. Thus isoproterenol stimulation significantly improved both regional and global function without subsequent evidence of increased cell death.     

Long-range translational coupling in single-stranded RNA bacteriophages: an evolutionary analysis

In coliphage MS2 RNA a long-distance interaction (LDI) between an internal segment of the upstream coat gene and the start region of the replicase gene prevents initiation of replicase synthesis in the absence of coat gene translation. Elongating ribosomes break up the repressor LDI and thus activate the hidden initiation site. Expression studies on partial MS2 cDNA clones identified base pairing between 1427-1433 and 1738-1744, the so-called Min Jou (MJ) interaction, as the molecular basis for the long-range coupling mechanism. Here, we examine the biological significance of this interaction for the control of replicase gene translation. The LDI was disrupted by mutations in the 3'-side and the evolutionary adaptation was monitored upon phage passaging. Two categories of pseudorevertants emerged. The first type had restored the MJ interaction but not necessarily the native sequence. The pseudorevertants of the second type acquired a compensatory substitution some 80 nt downstream of the MJ interaction that stabilizes an adjacent LDI. In one examined case we confirmed that the second site mutations had restored coat-replicase translational coupling. Our results show the importance of translational control for fitness of the phage. They also reveal that the structure that buries the replicase start extends to structure elements bordering the MJ interaction.     

The "clustered structure" of the purines/pyrimidines distribution in DNA distinguishes systematically between coding and non-coding sequences

A method allowing to measure the inhomogeneous distribution of purines/pyrimidines in nucleotide sequences is developed. We show that this measure relates to the coding or non-coding character of the considered sequence. Coding sequences present a near to the random Pu or Py distribution. This property is shared by both protein-coding DNA and functional RNA-coding DNA. Non-coding sequences present a highly clustered inhomogeneity. We propose the hypothesis, corroborated with appropriate computer simulations, that this is due to the action of various transposition events accumulated for long time periods.     

Long-range and short-range oxidative damage to DNA: photoinduced damage to guanines in ethidium-DNA assemblies

Short-range and long-range photoreactions between ethidium and DNA have been characterized. While no DNA reaction is observed upon excitation into the visible absorption band of ethidium, higher-energy irradiation (313-340 nm) leads both to direct strand cleavage at the 5'-G of 5'-GG-3' doublets and to piperidine-sensitive lesions at guanine. This reactivity is not consistent with oxidation of guanine by either electron transfer or singlet oxygen as shown by comparison with reactions of a rhodium intercalator and methylene blue, respectively. By covalently tethering ethidium to one end of a DNA duplex, we demonstrate the presence of two distinct reactions, one short-range and the other long-range. The short-range reaction involves a covalent modification of guanine by ethidium, based upon HPLC analysis of the nucleoside products and studies with ethidium derivatives. The long-range reaction is entirely consistent with oxidation of guanine by DNA-mediated electron transfer. The yield of this electron-transfer reaction is not attenuated with distance; equal yields of guanine damage are observed at a proximal (17 A Et-GG separation) and distal (44 A Et-GG separation) site. These results are quite similar to those previously observed with a covalently tethered rhodium photooxidant and underscore the unique ability of the DNA base stack to facilitate long-range electron transfer so as to effect oxidative damage from a distance.     

Tn916 transposition in Haemophilus influenzae Rd: preferential insertion into noncoding DNA

OBJECTIVES: A low level of education is associated with an increased risk of developing a dementia disorder, as well as with a higher risk of cardiovascular disease. The aim of this study was to investigate the association between education and cardiovascular risk factors, and to study the relation between these factors and cognitive function in elderly men. DESIGN: Cross-sectional population-based study. SETTING: Uppsala, Sweden. SUBJECTS: 504 men aged 69-74 years, participants in a longitudinal health survey concerning cardiovascular risk factors. MAIN OUTCOME MEASURE: Cognitive function as measured by a composite score of 13 standard psychometric tests. RESULTS: A low level of education was associated with poorer cognitive performance, as well as with obesity, smoking, diabetes, high concentrations of serum triglycerides and plasma fibrinogen. In the entire cohort, subjects with obesity, smoking, diabetes or hypertriglyceridaemia showed impaired cognitive test results, independent of socio-economic factors. When stroke cases were excluded, obesity and smoking were still related to impaired cognitive function. CONCLUSIONS: Smoking and obesity with associated metabolic disturbances are inversely related both to educational level and to cognitive function. Cognitive decline of vascular origin is potentially preventable by treatment of risk factors. The question of whether the increased vascular risk contributes to the higher prevalence of cognitive disorders in individuals with low socio-economic status, needs to be further evaluated in longitudinal population-based studies.