Linguistics of Nucleotide Sequences II: Stationary Words in Genetic Texts and the Zonal Structure of DNA

Abstract

Words are irregularly distributed in genetic texts. The analysis of this irregularity leads to the notion of stationary and non-stationary words. The polyW and polyS tracts are shown to be the most non-stationary words in genetic texts (here W-(AT), S-{G,C}, a polyW tract is a sequence of AT nucleotides and a polyS tract is a sequence of G,C nucleotides. The distribution of stationary words suggests a method for partitioning DNA into zones. The zones obtained in the case of the phage are interpreted in the light of the Dowe hypothesis of the modular structure of bacteriophage genomes.     

Genome inhomogeneity is determined mainly by WW and SS dinucleotides

According to the hypothesis of the modular structure of DNA,genomes consist of modules of various nature which may differin statistical characteristics. Statistical analysis helps inrevealing the differences in statistical characteristics andpredicting the modular structure. In this connection the questionabout the contribution of each word of length l (l-tuple) tothe inhomogeneity of genetic text arises. The notion of stationary(i.e. relatively evenly distributed over a genome) versus non-stationaryl-tuples has been introduced previously. In this paper, thedinucleotide distributions for all long sequences from GenBankwere analyzed and it was shown that non-stationary dinucleotidesare closely associated with polyW and polyS tracts (W denotes‘weak’ nucleotides A or T, while S stands for the‘strong’ nucleotides G or C). Thus, genome inhomogeneityis shown to be determined mainly by AA, TT, GG, CC, AT, TA,GC and CG dinucleotides. It has been demonstrated that neither‘codon usage’ nor the ‘isochore model’can account for this phenomenon.     

Co-polymer tracts in eukaryotic, prokaryotic, and organellar DNA.

Large variations in DNA base composition and noticeable strand asymmetries are known to occur between different organisms and within different regions of the genomes of single organisms. Apparently such composition and sequence biases occur to fulfill structural rather than informational requirements. Here we report the wide occurrence of a more subtle biasing of DNA sequence that can have structural consequences: an increase or a suppression of the number of long tracts of two-base co-polymers. Strong biases were observed when the DNA sequences of the longest eukaryotic, prokaryotic, and organellar entries in the GenBank data base (totaling 773 kilobases) were analyzed for the number of occurrences of tracts of the two-base co-polymers (A,T)n, (G,C)n, and (A,C)n as a function of tract length. (The expression (A,T)n is used here to denote an uninterrupted tract, n nucleotides in length, of A and T bases in any proportion or order, terminated at each end by a G or C residue.) Characteristic differences are also observed in tract biases of eukaryotic vs. prokaryotic organisms.     

On the question of DNA bending: two-dimensional NMR studies on d(GTTTTAAAAC)2 in solution

It is very well documented that the presence of an An.Tn tract causes intrinsic DNA bending. Hagerman demonstrated that the sequence in which the An.Tn tracts are joined plays a very crucial role in determining DNA bending. For example, Hagerman showed that the polymer with a repeat of d(GA4T4C)n greater than or equal to 10 is bent but the polymer with a repeat of d(GT4A4C)n greater than or equal to 10 is not bent [Hagerman, P. J. (1986) Nature (London) 326, 720-722]. Earlier we have shown that the decamer repeat d(GA4T4C)2 is itself bent with a finite structural discontinuity at the A----T sequence [Sarma, M. H., Gupta, G., & Sarma, R. H. (1988) Biochemistry 27, 3423-3432]. In the present article, we summarize our studies on the decamer repeat d(GT4A4C)2 structure in solution. By employment of 1D and 2D 1H NMR studies at 500 MHz a complete sequential assignment has been made for the exchangeable and nonexchangeable protons belonging to the ten nucleotides. NOESY data were collected for d(GT4A4C)2 at 17 degrees C in D2O for three mixing times, 150, 100, and 50 ms. A quantitative NOESY simulation technique was employed to arrive at a structural model of d(GT4A4C)2 in solution. Our detailed analyses revealed the following structural features: (i) The duplex adopts the gross morphology of a B-DNA. (ii) All the A.T pairs are propeller twisted (less than or equal to -15 degrees). (iii) Although both A and T nucleotides belong to the C2'-endo,anticonformational domain, there is a mild variation in the actual conformation of the A and T residues. (iv) Even though there is a subtle conformational difference in the A and T nucleotides, two structural frames of T4.A4 segments are joined at the T----A sequence in such a way that there is no finite discontinuity at the junction; i.e., two neighboring frames exactly coincide at the T----A junction. Thus, our studies on d(GA4T4C)2 (Sarma et al., 1988) and on d(GT4A4C)2 (this article) reveal the structural peculiarity of the An.Tn tract and the effect of A----T/T----A sequence in causing DNA bending.     

Xrep, a plasmid-stimulating X chromosomal sequence bearing similarities to the BK virus replication origin and viral enhancers.

The human X chromosome-linked fragment, "Xrep," was sequenced because it exerts a positive effect on plasmid growth in both E. coli and Saccharomyces cerevisiae. The sequence revealed three features similar to the human BK virus replication origin: Xrep has a true palindrome, CCTCC(T)3CCTCC, which is similar to "true" palindrome-like sequences found at the replication origins of polyoma [CCTC(T/C)10CTCC], BK [CCTC(A/G)8CCTCC] and SV40 [CCTCC(A)6GCCTCC] viruses. Twenty nucleotides away from the true palindrome, Xrep has the sequence GAATCCTATTCACTTTT while BK virus, the human analogue of SV40, has GAAATCCCTATTCTTTT in exactly the same position relative to the true palindrome. These two 17-mers differ only in the positions of two nucleotides comparing Xrep and BK virus. Also similar to the replication origins of DNA viruses, Xrep appears to have a cluster of enhancers adjacent to the origin-like sequences. Potent enhancer-like activity was detected in pSV1 X CAT/Xrep constructs. Xrep may originate from an endogenous virus, or from an X chromosomal replication origin.     

Effect of the methyl group on DNA bending and curvature: structure of d(GA4U4C)2 in solution

NMR studies on d(GA4T4C)2 and d(GT4A4C)2 indicated two important factors that contribute to intrinsic DNA bending in polymers containing A/T tracts [Sarma, M.H., et al. (1988) Biochemistry 27, 3423; Gupta, G., et al. (1988) Biochemistry 27, 7909]. They are (i) propeller-twisted A.T pairs with associated bifurcated H bonds inside the A/T tract and (ii) the base sequence that joins the two neighboring A/T tracts. As an extension of our bending project, we carried out quantitative NMR studies on the decamer d(GA4U4C)2, a structural analogue of d(GA4T4C)2, to examine the effect of the methyl group on DNA bending. On the basis of quantitative NMR analysis, we arrive at the following results. (i) The decamer d(GA4U4C)2 adopts the gross morphology of a right-handed B-DNA duplex with A and U nucleotides belonging to C2'-endo,anti domain. (ii) A.U pairs are propeller twisted and hence can result in an array of interstrand bifurcated H bonds involving N6 of A and O4 of U (one base pair apart) inside the A/U tract. (iii) The orientations of A and U with respect to the long axis of the molecule are different; as a result, at the A5-U6 sequence that joins the two A/U tracts, two neighboring frames of reference do not exactly coincide in space and a junction is created at A5-U6. (iv) Inside the A/U tract, intrastrand stacking is more compact (average separation between secessive base planes being 3.2 A) than at the A5-U6 junction, where average separation between the base planes of A5 and U6 is 3.6 A.(ABSTRACT TRUNCATED AT 250 WORDS)     

Stability and mismatch discrimination of locked nucleic acid-DNA duplexes

Locked nucleic acids (LNA; symbols of bases, +A, +C, +G, and +T) are introduced into chemically synthesized oligonucleotides to increase duplex stability and specificity. To understand these effects, we have determined thermodynamic parameters of consecutive LNA nucleotides. We present guidelines for the design of LNA oligonucleotides and introduce free online software that predicts the stability of any LNA duplex oligomer. Thermodynamic analysis shows that the single strand-duplex transition is characterized by a favorable enthalpic change and by an unfavorable loss of entropy. A single LNA modification confines the local conformation of nucleotides, causing a smaller, less unfavorable entropic loss when the single strand is restricted to the rigid duplex structure. Additional LNAs adjacent to the initial modification appear to enhance stacking and H-bonding interactions because they increase the enthalpic contributions to duplex stabilization. New nearest-neighbor parameters correctly forecast the positive and negative effects of LNAs on mismatch discrimination. Specificity is enhanced in a majority of sequences and is dependent on mismatch type and adjacent base pairs; the largest discriminatory boost occurs for the central +C·C mismatch within the +T+C+C sequence and the +A·G mismatch within the +T+A+G sequence. LNAs do not affect specificity in some sequences and even impair it for many +G·T and +C·A mismatches. The level of mismatch discrimination decreases the most for the central +G·T mismatch within the +G+G+C sequence and the +C·A mismatch within the +G+C+G sequence. We hypothesize that these discrimination changes are not unique features of LNAs but originate from the shift of the duplex conformation from B-form to A-form.     

Total synthesis of the structural gene for the precursor of a tyrosine suppressor transfer RNA from Escherichia coli. 11. Enzymatic joining to form the total DNA duplex.

The DNA duplex corresponding to the entire length (126 nucleotides) of the precursor for an Escherichia coli tyrosine tRNA has been synthesized. Duplex [I] (Sekiya, T., Besmer, P., Takeya, T., and Khorana, H. G.(1976) J. Biol. Chem. 251, 634-641), corresponding to the nucleotide sequence 1-26, containing single-stranded ends and carrying one appropriately labeled 5'-phosphate group, was joined to duplex [II] (Loewen, P. C., Miller, R. C., Panet, A., Sekiya, T., and Khorana, H. G. (1976) J. Biol. Chem. 251, 642-650) (nucleotide sequence 23-66 or 23-60) was phosphorylated with [gamma-33P]ATP at the 5'-OH ends. Duplex [III] (Panet, A., Kleppe, R., Kleppe, K., and Khorana, H. G. (1976) J. Biol. Chem. 251, 651-657) (nucleotide sequence 57-94 (Fig. 2)) was also phosphorylated at 5'-ends with [gamma-33P]ATP and was joined to duplex [IV] (Caruthers, M. H., Kleppe, R., Kleppe, K., and Khorana, H. G. (1976) J. Biol. Chem. 251, 658-666) (nucleotide sequence 90-126) which carried a 33P-labeled phosphate group on nucleotide 90. The joined product, duplex [III + IV] (nucleotide sequence 57-126) was characterized. The latter duplex was joined to the duplex [I + II] to give the total duplex. The latter contains singlestranded ends (nucleotides 1 to 6 and 121 to 126) which can either be "filled in" to produce the completely base-paired duplex or may be used to add the promoter and terminator regions at the appropriate ends.     

Effect of strand orientation on the interaction of berenil with DNA triple helices

Using circular dichroism spectroscopy the ability of berenil, a minor groove binding drug, to induce triple helix formation was investigated with two oligonucleotides designed to form two intramolecular triplexes containing T*A:T and G*G:C triplets, which differ only by the orientation of their third strand: 5'-d(G4A4G4-[T4]-C4T4C4-[T4]-G4T4G4), and 5'-d(G4T4G4-[T4]-G4A4G4-[T4]-C4T4C4), where [T4] represents a stretch of four thymine residues. We demonstrate that when added to the duplex form of these oligonucleotides, berenil induces triplex structure formation only if the orientation of third strand is anti-parallel to the purine strand.     

Identification of sequence polymorphisms in the D-loop region of mitochondrial DNA as a risk factor for lung cancer

Accumulation of single nucleotide polymorphisms (SNPs) in the displacement loop (D-loop) of mitochondrial DNA (mtDNA) may be associated with an increased cancer risk. We investigated the lung cancer risk profile of D-loop SNPs in a case-controlled study. The minor alleles of nucleotides 235A/G and 324A/G were associated with an increased risk for lung cancer patients. The minor alleles of the nucleotides 151C/T, 200A/G, 524C/CA, and 16274G/A were specifically associated with the cancer risk of squamous cell carcinoma, whereas the minor allele of nucleotide 16298T/C was specifically associated with the risk of small cell lung cancer. In conclusion, SNPs in mtDNA are potential modifiers of lung cancer risk. The analysis of genetic polymorphisms in the mitochondrial D-loop can help identify subgroups of patients who are at a high risk of developing lung cancer.     

'Smoking genes': a genetic association study

Some controversy exists on the specific genetic variants that are associated with nicotine dependence and smoking-related phenotypes. The purpose of this study was to analyse the association of smoking status and smoking-related phenotypes (included nicotine dependence) with 17 candidate genetic variants: CYP2A6*1×2, CYP2A6*2 (1799T>A) [rs1801272], CYP2A6*9 (-48T>G) [rs28399433], CYP2A6*12, CYP2A13*2 (3375C>T) [rs8192789], CYP2A13*3 (7520C>G), CYP2A13*4 (579G>A), CYP2A13*7 (578C>T) [rs72552266], CYP2B6*4 (785A>G), CYP2B6*9 (516G>T), CHRNA3 546C>T [rs578776], CHRNA5 1192G>A [rs16969968], CNR1 3764C>G [rs6928499], DRD2-ANKK1 2137G>A (Taq1A) [rs1800497], 5HTT LPR, HTR2A -1438A>G [rs6311] and OPRM1 118A>G [rs1799971]. We studied the genotypes of the aforementioned polymorphisms in a cohort of Spanish smokers (cases, N = 126) and ethnically matched never smokers (controls, N = 80). The results showed significant between-group differences for CYP2A6*2 and CYP2A6*12 (both P<0.001). Compared with carriers of variant alleles, the odds ratio (OR) for being a non-smoker in individuals with the wild-type genotype of CYP2A6*12 and DRD2-ANKK1 2137G>A (Taq1A) polymorphisms was 3.60 (95%CI: 1.75, 7.44) and 2.63 (95%CI: 1.41, 4.89) respectively. Compared with the wild-type genotype, the OR for being a non-smoker in carriers of the minor CYP2A6*2 allele was 1.80 (95%CI: 1.24, 2.65). We found a significant genotype effect (all P≤0.017) for the following smoking-related phenotypes: (i) cigarettes smoked per day and CYP2A13*3; (ii) pack years smoked and CYP2A6*2, CYP2A6*1×2, CYP2A13*7, CYP2B6*4 and DRD2-ANKK1 2137G>A (Taq1A); (iii) nicotine dependence (assessed with the Fagestrom test) and CYP2A6*9. Overall, our results suggest that genetic variants potentially involved in nicotine metabolization (mainly, CYP2A6 polymorphisms) are those showing the strongest association with smoking-related phenotypes, as opposed to genetic variants influencing the brain effects of nicotine, e.g., through nicotinic acetylcholine (CHRNA5), serotoninergic (HTR2A), opioid (OPRM1) or cannabinoid receptors (CNR1).     

Solution structure of an LNA hybridized to DNA: NMR study of the d(CT(L)GCT(L)T(L)CT(L)GC):d(GCAGAAGCAG) duplex containing four locked nucleotides

We have used two-dimensional (1)H NMR spectroscopy at 750 MHz to determine a high-resolution solution structure of an oligonucleotide containing restricted nucleotides with a 2'-O, 4'-C-methylene bridge (LNA) hybridized to the complementary DNA strand. The LNA:DNA duplex examined contained four thymidine LNA modifications (T(L), d(C1T(L)2G3C4T(L)5T(L)6C7T(L)8G9C10):d( G11C12A13G14A15A16G17C 18A19G20). A total relaxation matrix approach was used to obtain interproton distance bounds from NOESY cross-peak intensities. These distance bounds were used as restraints in molecular dynamics (rMD) calculations. Forty final structures were generated for the duplex from A-form and B-form DNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the 40 structures of the complex was 0.6 A. The sugar puckerings are averaged values of a dynamic interchange between N- and S-type conformation except in case of the locked nucleotides that were found to be fixed in the C3'-endo conformation. Among the other nucleotides in the modified strand, the furanose ring of C7 and G9 is predominantly in the N-type conformation whereas that of G3 is in a mixed conformation. The furanose rings of the nucleotides in the unmodified complementary strand are almost exclusively in the S-type conformation. Due to these different conformations of the sugars in the two strands, there is a structural strain between the A-type modified strand and the B-type unmodified complementary strand. This strain is relaxed by decreasing the value of rise and compensating with tip, buckle, and propeller twist. The values of twist vary along the strand but for a majority of the base pairs a value even lower than that of A-DNA is observed. The average twist over the sequence is 32+/-1 degrees. On the basis of the structure, we conclude that the high stability of LNA:DNA duplexes is caused by a local change of the phosphate backbone geometry that favors a higher degree of stacking.     

Genetic variations of the NR3C1 gene in children with sporadic nephrotic syndrome

Previous studies have demonstrated that the genetic variations of glucocorticoid receptor gene (NR3C1) are associated with both familial steroid resistance and acquired steroid resistance in some diseases, such as Cushing's disease, leukemia, lupus nephritis, and female pseudohermaphroditism. In this study, we examined the genetic variations of NR3C1 in 35 children with sporadic steroid-resistant nephrotic syndrome (SRNS), and in 83 cases with sporadic steroid-sensitive NS (SSNS) using polymerase chain reaction, denaturing high-performance liquid chromatography and DNA sequencing, and analyzed possible associations between NR3C1 variants and steroid resistance in sporadic NS. No causative mutations were found; however, six previously identified and six novel polymorphisms, 1206C > T, 1374A > G, 2382C > T, 2193T > G, IVS7-68_-63delAAAAAA, and IVS8-9C > G, were detected. Two novel haplotypes, [1374A > G; IVS7-68_-63delAAAAAA; IVS8-9C > G; 2382C > T] and [1896C > T; 2166C > T; 2430T > C], of NR3C1 were also identified in sporadic NS and controls. The odds ratios (95% Confidence Interval) for the two novel NR3C1 haplotypes in the sporadic nephrotic children at risk of steroid resistance were 4.970 (0.889-27.788) and 2.194 (0.764-6.306), respectively, but the association between NR3C1 haplotypes and steroid resistance was not significant. Further studies on the possible association between the two novel NR3C1 haplotypes and steroid resistance in sporadic NS in larger cohorts are required.     

Number of genes and base composition of mitochondrial tRNA from Saccharomyces cerevisiae.

Increasing amounts of mitochondrial [32P] tRNA (4S fraction), were hybridized with mitochondrial DNA OF Saccharomyces cerevisiae. At saturation, the calculated number of genes for 4S mitochondrial RNA was 20. Mitochondrial [32P] tRNA eluted from the hydrids obtained either with an excess of tRNA or an excess of DNA showed, after alkaline hydrolysis and chromatography, a G+C content of 28 and 35 p. cent respectively. This last value is similar to that found with the total 4S fraction. The odd nucleotides T (about 1T per sequence), U, hU are present in mitochondrial tRNA. Some sequence may begin with pG.     

The Gly-952 residue of Saccharomyces cerevisiae DNA polymerase alpha is important in discriminating correct deoxyribonucleotides from incorrect ones

Gly-952 is a conserved residue in Saccharomyces cerevisiae DNA polymerase alpha (pol alpha) that is strictly required for catalytic activity and for genetic complementation of a pol alpha-deficient yeast strain. This study analyzes the role of Gly-952 by characterizing the biochemical properties of Gly-952 mutants. Analysis of the nucleotide incorporation specificity of pol alpha G952A showed that this mutant incorporates nucleotides with extraordinarily low fidelity. In a steady-state kinetic assay to measure nucleotide misincorporation, pol alpha G952A incorporated incorrect nucleotides more efficiently than correct nucleotides opposite template C, G, and T. The fidelity of the G952A mutant polymerase was highest at template A, where the ratio of incorporation of dCMP to dTMP was as high as 0.37. Correct nucleotide insertion was 500- to 3500-fold lower for G952A than for wild type pol alpha, with up to 22-fold increase in pyrimidine misincorporation. The Km for G952A pol alpha bound to mismatched termini T:T, T:C, C:A, and A:C was 71- to 460-fold lower than to a matched terminus. Furthermore, pol alpha G952A preferentially incorporated pyrimidine instead of dAMP opposite an abasic site, cis-syn cyclobutane di-thymine, or (6-4) di-thymine photoproduct. These data demonstrate that Gly-952 is a critical residue for catalytic efficiency and error prevention in S. cerevisiae pol alpha.     

Association Between Single-Nucleotide Polymorphism in CISH Gene and Susceptibility to Tuberculosis in Chinese Han Population

The cytokine-inducible SRC homology 2 domain (CISH) gene is up-regulated by IL-2 in response to infection, and inhibits microbial infection. The objective of the present study was to examine whether genetic variants of CISH (SNPs) are associated with increased susceptibility to tuberculosis (TB) in individuals of Chinese Han ethnicity. We sequenced five previously identified SNPs of CISH in patients with TB or healthy controls. Three of the SNPs, rs148685070 [position ?639; C/C], rs414171 [position ?292; A/T], and rs6768300 [position ?163; C/G]) are located in the promoter region, while the fourth (rs2239751 [position +1320; A/C]) near the translation start site, and the fifth (rs622502 [position +3415; C/G]) in the third intron. The AA genotypes of the SNPs rs2239751 and rs414171 were significantly associated with TB. Multivariate logistic regression analysis demonstrated that subjects with the rs414171 AA genotype were more likely to have TB than those with the AT genotype. By contrast, we did not observe genetic variants of the rs148685070 SNP. In conclusion, two genetic variants in CISH gene appear to increase susceptibility to TB in Chinese Han population.     

Linguistics of nucleotide sequences. I: The significance of deviations from mean statistical characteristics and prediction of the frequencies of occurrence of words

Mathematical models of the generation of genetic texts appeared simultaneously with the first sequencing DNA. They are used to establish functional and evolutionary relations between genetic texts, to predict the number and distribution of specific sites in a sequence and to identify "meaningful" words. The present paper deals with two problems: 1) The significance of deviations from the mean statistical characteristics in a genetic text. Anyone who has addressed himself to the statistical analysis of sequenced DNA is familiar with the question: what deviations from the expected frequencies of occurrence of particular words testify to the "biological" significance of those words? We propose a formula for the variance of the number of word's occurrences in the text, with allowance for word overlaps, making it possible to assess the significance of the deviations from the expected statistical characteristics. 2) A new method for predicting the frequencies of occurrence of particular words in a genetic text using the statistical characteristics of "spaced" L-grams. The method can be used for predicting the number of restriction sites in human DNA and in planning experiments on the physical mapping and sequencing of the human genome.     

Nucleotide sequence of the genome of the bacteriophage alpha 3: interrelationship of the genome structure and the gene products with those of the phages, phi X174, G4 and phi K.

The complete nucleotide sequence of the genome of the circular single-stranded DNA (isometric) phage alpha 3 has been determined and compared with that of the related phages phi X174 and G4. The alpha 3 genome consists of 6087 nucleotides, which is 701 nucleotides longer than the nucleotide sequence of the phi X174 genome and 510 nucleotides more than that of the G4 genome. The results demonstrated that the three phage species have 11 homologous genes (A, A*, B, C, K, D, E, J, F, G and H), the order of which is fundamentally identical, suggesting that they have evolved from a common ancestor. The sequence of some genes and untranslated intergenic regions, however, differs significantly from phage to phage: for example, the degree of amino acid sequence homology of the gene product is averaged at 47.7% between alpha 3 and phi X174 and 46.9% between alpha 3 and G4, and alpha 3 has a remarkable longer intergenic region composed of 758 nucleotides between the genes H and A compared with the counterparts of phi X174 and G4. Meanwhile, in vivo experiments of genetic complementation showed that alpha 3 can use none of the gene products of phi X174 and G4, whereas the related phage phi K can rescue alpha 3 nonsense mutants of the genes B, C, D and J. These sequencing and in vivo rescue results indicated that alpha 3 is closely related to phi K, but distantly remote from phi X174 or G4, and supported an evolutional hypothesis which has been so far proposed that the isometric phages are classified into three main groups: the generic representatives are phi X174, G4 and alpha 3.