Horizontal gene transfer in microbial genome evolution

Horizontal gene transfer is the collective name for processes that permit the exchange of DNA among organisms of different species. Only recently has it been recognized as a significant contribution to inter-organismal gene exchange. Traditionally, it was thought that microorganisms evolved clonally, passing genes from mother to daughter cells with little or no exchange of DNA among diverse species. Studies of microbial genomes, however, have shown that genomes contain genes that are closely related to a number of different prokaryotes, sometimes to phylogenetically very distantly related ones. (Doolittle et al., 1990, J. Mol. Evol. 31, 383-388; Karlin et al., 1997, J. Bacteriol. 179, 3899-3913; Karlin et al., 1998, Annu. Rev. Genet. 32, 185-225; Lawrence and Ochman, 1998, Proc. Natl. Acad. Sci. USA 95, 9413-9417; Rivera et al., 1998, Proc. Natl. Acad. Sci. USA 95, 6239-6244; Campbell, 2000, Theor. Popul. Biol. 57 71-77; Doolittle, 2000, Sci. Am. 282, 90-95; Ochman and Jones, 2000, Embo. J. 19, 6637-6643; Boucher et al. 2001, Curr. Opin., Microbiol. 4, 285-289; Wang et al., 2001, Mol. Biol. Evol. 18, 792-800). Whereas prokaryotic and eukaryotic evolution was once reconstructed from a single 16S ribosomal RNA (rRNA) gene, the analysis of complete genomes is beginning to yield a different picture of microbial evolution, one that is wrought with the lateral movement of genes across vast phylogenetic distances. (Lane et al., 1988, Methods Enzymol. 167, 138-144; Lake and Rivera, 1996, Proc. Natl. Acad. Sci. USA 91, 2880-2881; Lake et al., 1999, Science 283, 2027-2028).     

Comparison of analyses of DNA curvature

Popular programs for characterizing DNA structure include Curves 5.1 (Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 63-91, 1988; Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 655-67, 1989) and Freehelix98 (Dickerson, R. E., Nucleic Acids Res. 26, 1906-1926, 1998), along with the more recent 3DNA (X. J. Lu, Z. Shakked and W. K. Olson., J. Mol. Biol. 300, 819-840 (2000). Given input of structural coordinates, all of these programs return values of the local helical parameters, such as roll, tilt, twist, etc. The first two programs also provide characterization of global curvature. Madbend (Strahs, D. and Schlick, T., J. Mol. Biol. 301, 643-663, 2000), a program that computes global curvature from local roll, tilt, and twist parameters, can be applied to the output of all three structural programs. We have compared the curvature predicted by the three programs with and without the use of Madbend. Global bend magnitudes and directions as well as values of helical kinks were calculated for four high-resolution DNA structures and four model DNA helices. Global curvature determined by Curves 5.1 without Madbend was found to differ from values obtained using Freehelix98 with or without Madbend or 3DNA and Curves 5.1 with Madbend. Using model helices, this difference was attributed the fact that Curves 5.1 is the only program sensitive to changes in axial displacement, such as shift and slide. Madbend produced robust values of bend magnitude and direction, and displayed little sensitivity to axis displacement or the source of local helical parameters. Madbend also appears to be the method of choice for bending comparisons of high-resolution structures with results from cyclization kinetics, a method that measures DNA curvature as a vectorial sum of local roll and tilt angles.     

Xenopus egg lysates repair heat-generated DNA nicks with an average patch size of 36 nucleotides.

Base excision repair (BER) is an essential DNA repair pathway since it processes spontaneous (endogenous) DNA damage such as abasic sites, oxidized and alkylated bases, as well as mismatches arising from deamination of cytosine and 5-methylcytosine. Some of these lesions are repaired by the exchange of a single deoxynucleotide [Dianov, G. et al. (1992) Mol. Cell. Biol. 12, 1605-1612; Wiebauer, K. and Jiricny, J. (1990) Proc. Natl. Acad. Sci. USA, 87, 5842-5845] or a few deoxynucleotides [Matsumoto, Y. et al. (1994) Mol. Cell. Biol., 14 6187-6197]. Here we report that DNA single strand breaks induced by hyperthermic conditions are repaired with an average patch size of approximately 36 nt in Xenopus laevis egg lysates.     

A computer graphics study of sequence-directed bending in DNA

We present theoretical results to account for the unusual physical properties of a 423 bp DNA restriction fragment isolated from the kinetoplast of the trypanosomatid Leishmania tarentolae. This fragment has an anomalously low electrophoretic mobility in polyacrylamide gels and a rotational relaxation time smaller than that of normally-behaved control fragments of the same molecular weight. Our earlier work (Proc. Natl. Acad. Sci. USA 79, 7664, 1982) has attributed these anomalies to the highly periodic distribution of the dinucleotide ApA in the DNA sequence. As originally proposed by Trifonov and Sussman (Proc. Natl. Acad. Sci. USA 77, 3816, 1980) local features of the DNA structure such as a small bend at ApA, if repeated with the periodicity of the helix, will cause systematic bending of the molecule. Computer graphics representations of DNA chain trajectories are presented for different structural models. It is shown that the structural model of Calladine (J. Mol. Biol. 161, 343, 1982) which is based on crystallographic data, is unsuccessful in predicting the systematic bending of DNA in solution.     

The curvature vector in nucleosomal DNAs and theoretical prediction of nucleosome positioning.

Using our model for predicting DNA superstructures from the sequence, the average distribution of the phases of curvature along the sequences of the set of the 177 nucleosomal DNAs investigated by Satchwell et al. (J. Mol. Biol. 191 (1986) 659) was calculated. The diagram obtained shows very significant features which allow the visualization of the intrinsic nucleosomal superstructure characterized by two quasi-parallel tracts of a flat left-handed superhelical turn connected by a left-handed inflection in a perpendicular direction; such a superstructure appears to be closely related to the nucleosome model of Travers and Klug (Phil. Trans. R. Soc. Lond. 317 (1987) 537). The nucleosomal curvature phase diagram was then adopted as a sensitive determinant for the nucleosome virtual positioning in DNAs via correlation function, obtaining a good agreement with the experimental mapping of SV40 regulatory region as recently investigated by Ambrose et al. (J. Mol. Biol. 209 (1989) 255). This analysis shows also the presence of a constant phase relation between the virtual nucleosome positions which suggests its possible implication in the nucleosome condensation in chromatin.     

Signal sequence of alkaline phosphatase of Escherichia coli.

The amino acid sequence of the signal sequence of phoA was determined by DNA sequencing by using the dideoxy chain termination technique (Sanger et al., Proc. Natl. Acad. Sci. U.S.A. 74:5463-5467, 1977). The template used was single-stranded DNA obtained from M13 on f1 phage derivatives carrying phoA, constructed by in vitro recombination. The results confirm the sequence of the first five amino acids determined by Sarthy et al. (J. Bacteriol. 139:932-939, 1979) and extend the sequence in the same reading frame into the amino terminal region of the mature alkaline phosphatase (Bradshaw et al., Proc. Natl. Acad. Sci. U.S.A., 78:3473-3477, 1981). As was predicted (Inouye and Beckwith, Proc. Natl. Acad. Sci. U.S.A. 74:1440-1444, 1977), the signal sequence was highly hydrophobic. The alteration of DNA sequence was identified for a promoter mutation that results in the expression of phoA independent of the positive control gene phoB and in insensitivity to high phosphate.     

A Computer Graphics Study of Sequence-Directed Bending in DNA

Abstract

We present theoretical results to account for the unusual physical properties of a 423 bp DNA restriction fragment isolated from the kinetoplast of the trypanosomatid Leishmania tarentolae. This fragment has an anomalously low electrophoretic mobility in Polyacrylamide gels and a rotational relaxation time smaller than that of normally-behaved control fragments of the same molecular weight. Our earlier work (Proc. Natl. Acad. Sci. USA 79, 7664, 1982) has attributed these anomalies to the highly periodic distribution of the dinucleotide ApA in the DNA sequence. As originally proposed by Trifonov and Sussman (Proc. Natl. Acad. Sci. USA 77, 3816,1980) local features of the DNA structure such as a small bend at ApA, if repeated with the periodicity of the helix, will cause systematic bending of the molecule.

Computer graphics representations of DNA chain trajectories are presented for different structural models. It is shown that the structural model of Calladine (J. Mol. Biol. 161, 343, 1982) which is based on crystallographic data, is unsuccessful in predicting the systematic bending of DNA in solution.     

A-Tract bending: insights into experimental structures by computational models

While solution structures of adenine tract (A-tract) oligomers have indicated a unique bend direction equivalent to negative global roll (commonly termed "minor-groove bending"), crystallographic data have not unambiguously characterized the bend direction; nevertheless, many features are shared by all A-tract crystal and solution structures (e.g. propeller twisting, narrow minor grooves, and localized water spines). To examine the origin of bending and to relate findings to the crystallographic and solution data, we analyze molecular dynamics trajectories of two solvated A-tract dodecamers: 1D89, d(CGCGA(6)CG), and 1D98, d(CGCA(6)GCG), using a new general global bending framework for analyzing bent DNA and DNA/protein complexes. It is significant that the crystallographically-based initial structures are converted from dissimilar to similar bend directions equivalent to negative global roll, with the average helical-axis bend ranging from 10.5 degrees to 14.1 degrees. The largest bend occurs as positive roll of 12 degrees on the 5' side of the A-tracts (supporting a junction model) and is reinforced by gradual curvature at each A-tract base-pair (bp) step (supporting a wedge model). The precise magnitude of the bend is subtly sequence dependent (consistent with a curved general sequence model). The conversion to negative global roll only requires small local changes at each bp, accumulated over flexible moieties both outside and inside the A-tract. In contrast, the control sequence 1BNA, d(CGCGA(2)TTCGCG), bends marginally (only 6.9 degrees ) with no preferred direction. The molecular features that stabilize the bend direction in the A-tract dodecamers include propeller twisting of AT base-pairs, puckering differences between A and T deoxyriboses, a narrow minor groove, and a stable water spine (that extends slightly beyond the A-tract, with lifetimes approaching 0.2 ns). The sugar conformations, in particular, are proposed as important factors that support bent DNA. It is significant that all these curvature-stabilizing features are also observed in the crystallographic structures, but yield overall different bending paths, largely due to the effects of sequences outside the A-tract. These results merge structural details reported for A-tract structures by experiment and theory and lead to structural and dynamic insights into sequence-dependent DNA flexibility, as highlighted by the effect of an A-tract variant of a TATA-box element on bending and flexibility required for TBP binding.     

Probable cloning artefacts previously interpreted as unusual leader sequences of rodent ornithine decarboxylase mRNAs--a cautionary tale

Messenger RNAs that have structurally unusual 5' leaders attract interest and provoke conjecture. Cloning and sequencing of two rodent ornithine decarboxylase (ODC) cDNAs, those for mouse [Kahana and Nathans, Proc. Natl. Acad. Sci. USA 82 (1985) 1673-1677] and, recently as published in this journal, for rat [Van Kranen et al., Gene 60 (1987) 145-155], have indicated the presence of such features. In both cases, the leader is unusually long and contains multiple AUG start codons preceding that which encodes the N terminus of the protein. In addition, the leader of the rat clone contains a 54-nt perfect inverted repeat. Because ODC expression appears to be regulated translationally, functional implications immediately suggest themselves. Certain unusual features of the mouse cDNA have proven artefactual [Brabant et al., Proc. Natl. Acad. Sci. USA 85 (1988) 2200-2204; Katz and Kahana, J. Biol. Chem. 263 (1988) 7604-7609]. It is likely that the putative leader sequence of rat ODC cDNA also resulted from a cloning artefact.     

Solution structure of an A-tract DNA bend

The solution structure of a DNA dodecamer d(GGCAAAAAACGG)/d(CCGTTTTTTGCC) containing an A-tract has been determined by NMR spectroscopy with residual dipolar couplings. The structure shows an overall helix axis bend of 19 degrees in a geometry consistent with solution and gel electrophoresis experiments. Fourteen degrees of the bending occurs in the GC regions flanking the A-tract. The remaining 5 degrees is spread evenly over its six AT base-pairs. The A-tract is characterized by decreasing minor groove width from the 5' to the 3' direction along the A strand. This is a result of propeller twist in the AT pairs and the increasing negative inclination of the adenine bases at the 3' side of the run of adenine bases. The four central thymine bases all have negative inclination throughout the A-tract with an average value of -6.1 degrees. Although this negative inclination makes the geometry of the A-tract different from all X-ray structures, the proton on N6 of adenine and the O4 of thymine one step down the helix are within distance to form bifurcated hydrogen bonds. The 5' bend of 4 degrees occurs at the junction between the GC flank and the A-tract through a combination of tilt and roll. The larger 3' bend, 10 degrees, occurs in two base steps: the first composed of tilt, -4.1 degrees, and the second a combination of tilt, -4.2 degrees, and roll, 6.0 degrees. This second step is a direct consequence of the change in inclination between an adjacent cytosine base, which has an inclination of -12 degrees, and the next base, a guanine, which has 3 degrees inclination. This bend is a combination of tilt and roll. The large change in inclination allows the formation of a hydrogen bond between the protons of N4 of the 3' cytosine and the O6 of the next 3' base, a guanine, stabilizing the roll component in the bend. These structural features differ from existing models for A-tract bends.For comparison, we also determined the structure of the control sequence, d(GGCAAGAAACGG)/d(CCGTTTCTTGCC), with an AT to GC transition in the center of the A-tract. This structure has no negative inclination in most of the bases within the A-tract, resulting in a bend of only 9 degrees. When ligated in phase, the control sequence has nearly normal mobility in gel electrophoresis experiments.     

DNA folding: structural and mechanical properties of the two-angle model for chromatin

We present a theoretical analysis of the structural and mechanical properties of the 30-nm chromatin fiber. Our study is based on the two-angle model introduced by Woodcock et al. (Woodcock, C. L., S. A. Grigoryev, R. A. Horowitz, and N. Whitaker. 1993. Proc. Natl. Acad. Sci. USA. 90:9021-9025) that describes the chromatin fiber geometry in terms of the entry-exit angle of the nucleosomal DNA and the rotational setting of the neighboring nucleosomes with respect to each other. We analytically explore the different structures that arise from this building principle, and demonstrate that the geometry with the highest density is close to the one found in native chromatin fibers under physiological conditions. On the basis of this model we calculate mechanical properties of the fiber under stretching. We obtain expressions for the stress-strain characteristics that show good agreement with the results of recent stretching experiments (Cui, Y., and C. Bustamante. 2000. Proc. Natl. Acad. Sci. USA. 97:127-132) and computer simulations (Katritch, V., C. Bustamante, and W. K. Olson. 2000. J. Mol. Biol. 295:29-40), and which provide simple physical insights into correlations between the structural and elastic properties of chromatin.     

The Chinese hamster Alu-equivalent sequence: a conserved highly repetitious, interspersed deoxyribonucleic acid sequence in mammals has a structure suggestive of a transposable element.

A consensus sequence has been determined for a major interspersed deoxyribonucleic acid repeat in the genome of Chinese hamster ovary cells (CHO cells). This sequence is extensively homologous to (i) the human Alu sequence (P. L. Deininger et al., J. Mol. Biol., in press), (ii) the mouse B1 interspersed repetitious sequence (Krayev et al., Nucleic Acids Res. 8:1201-1215, 1980) (iii) an interspersed repetitious sequence from African green monkey deoxyribonucleic acid (Dhruva et al., Proc. Natl. Acad. Sci. U.S.A. 77:4514-4518, 1980) and (iv) the CHO and mouse 4.5S ribonucleic acid (this report; F. Harada and N. Kato, Nucleic Acids Res. 8:1273-1285, 1980). Because the CHO consensus sequence shows significant homology to the human Alu sequence it is termed the CHO Alu-equivalent sequence. A conserved structure surrounding CHO Alu-equivalent family members can be recognized. It is similar to that surrounding the human Alu and the mouse B1 sequences, and is represented as follows: direct repeat-CHO-Alu-A-rich sequence-direct repeat. A composite interspersed repetitious sequence has been identified. Its structure is represented as follows: direct repeat-residue 47 to 107 of CHO-Alu-non-Alu repetitious sequence-A-rich sequence-direct repeat. Because the Alu flanking sequences resemble those that flank known transposable elements, we think it likely that the Alu sequence dispersed throughout the mammalian genome by transposition.     

Reading the three-dimensional structure of lattice model-designed proteins from their amino acid sequence.

While all the information required for the folding of a protein is contained in its amino acid sequence, one has not yet learned how to extract this information to predict the detailed, biological active, three-dimensional structure of a protein whose sequence is known. Using insight obtained from lattice model simulations of the folding of small proteins (fewer than 100 residues), in particular of the fact that this phenomenon is essentially controlled by conserved contacts (Mirny et al., Proc Natl Acad Sci USA 1995;92:1282) among (few) strongly interacting ("hot") amino acids (Tiana et al., J Chem Phys 1998;108:757-761), which also stabilize local elementary structures formed early in the folding process and leading to the (postcritical) folding core when they assemble together (Broglia et al., Proc Natl Acad Sci USA 1998;95:12930, Broglia & Tiana, J Chem Phys 2001;114:7267), we have worked out a successful strategy for reading the three-dimensional structure of lattice model-designed proteins from the knowledge of only their amino acid sequence and of the contact energies among the amino acids.     

Thermodynamic properties of globular proteins and the principle of stabilization of their native structure

A semi-empirical method has been used to estimate the thermodynamic parameters of hydration of buried surface areas of ribonuclease S, lysozyme and myoglobin from the model of complete unfolding according to Ooi et al. ((1987) Proc. Natl. Acad. Sci. USA 84, 3086-3090). The buried surface area of proteins is considered as the difference between the accessible surface area of native protein and the completely extended polypeptide chain according to Lee and Richards ((1971) J. Mol. Biol. 55, 379-400). The contributions of nonpolar and polar protein groups to the general value of Gibbs energy, enthalpy, entropy and heat capacity of hydration have been determined. The obtained results on the thermodynamic behavior of proteins in the process of complete unfolding are in good agreement with the results of microcalorimetric studies of thermal denaturation.     

A simple and efficient method for the separation and detection of small DNA fragments by electrophoresis in formamide containing agarose gels and Southern blotting to DBM-paper.

We show in this report that DNA fragments smaller than 300 bp are separated with high resolution by electrophoresis in concentrated (up to 7%) agarose gels containing 50% formamide. The separated DNA fragments can subsequently be quantitatively transferred to DBM-paper [Alwine, J.C. et al., Proc. Natl. Acad. Sci. USA (1977) 74, 5350-54] using the Southern technique [Southern, E.M., J. Mol. Biol. (1975) 98, 503-517], while preserving the sharpness of the original gel pattern. Since thin (0.2-0.4 mm) and thick (up to 5 mm) agarose slab gels can be easily handled in vertical or horizontal apparatus, this method should prove to be a very useful extention of the Southern technique, applicable to a variety of analytical and preparative purposes.     

Sequence similarities among monkey ori-enriched (ors) fragments

Nucleotide sequences have been determined for eight ors (ori-enriched sequence) fragments isolated from monkey DNA by a method that was designed to enrich for origins of DNA replication [Kaufmann et al., Mol. Cell. Biol. 5 (1985) 721-727]. Evidence has been presented that some or possibly all of these sequences can serve, albeit inefficiently, as oris in vivo [Frappier and Zannis-Hadjopoulos, Proc. Natl. Acad. Sci. USA 84 (1987) 6668-6672]. Two of the fragments were found to contain the long terminal repeat-like elements of the 'O-family' of moderately repetitive sequences that are present in human DNA as a transposon-like element [Paulson et al., Nature 315 (1985) 359-361]. Extensive pair-wise comparisons of the sequences failed to detect any statistically significant common sequences, except for long asymmetrically distributed A + T-rich stretches. Nonetheless, when the ors fragments were examined for the presence of published consensus sequences, seven of eight were found to contain the control sequence described by Dierks et al. [Cell 32 (1983) 695-706], and the same seven of eight were found to contain both the scaffold attachment region T consensus [Gasser and Laemmli, Cell 46 (1986) 521-530] and the minimal Saccharomyces cerevisiae autonomously replicating sequence consensus [e.g., Palzkill and Newlon, Cell 53 (1988) 441-450].