The most frequent short sequences in non-coding DNA

The purpose of this work is to determine the most frequent short sequences in non-coding DNA. They may play a role in maintaining the structure and function of eukaryotic chromosomes. We present a simple method for the detection and analysis of such sequences in several genomes, including Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster and Homo sapiens. We also study two chromosomes of man and mouse with a length similar to the whole genomes of the other species. We provide a list of the most common sequences of 9–14 bases in each genome. As expected, they are present in human Alu sequences. Our programs may also give a graph and a list of their position in the genome. Detection of clusters is also possible. In most cases, these sequences contain few alternating regions. Their intrinsic structure and their influence on nucleosome formation are not known. In particular, we have found new features of short sequences in C. elegans, which are distributed in heterogeneous clusters. They appear as punctuation marks in the chromosomes. Such clusters are not found in either A. thaliana or D. melanogaster. We discuss the possibility that they play a role in centromere function and homolog recognition in meiosis.     

Extending assembly of short DNA sequences to handle error

Inexpensive de novo genome sequencing, particularly in organisms with small genomes, is now possible using several new sequencing technologies. Some of these technologies such as that from Illumina's Solexa Sequencing, produce high genomic coverage by generating a very large number of small reads ( approximately 30 bp). While prior work shows that partial assembly can be performed by k-mer extension in error-free reads, this algorithm is unsuccessful with the sequencing error rates found in practice. We present VCAKE (Verified Consensus Assembly by K-mer Extension), a modification of simple k-mer extension that overcomes error by using high depth coverage. Though it is a simple modification of a previous approach, we show significant improvements in assembly results on simulated and experimental datasets that include error. AVAILABILITY: http://152.2.15.114/~labweb/VCAKE     

Fluorescence-based detection of short DNA sequences under non-denaturing conditions

The ability of peptide nucleic acid (PNA) to open up duplex DNA in a highly sequence-specific manner makes it possible to detect short DNA sequences on the background of or within genomic DNA under non-denaturing conditions. To do so, chosen marker sites in double-stranded DNA are locally opened by a pair of PNA openers, thus transforming one strand within the target region (20-30 bp) into the single-stranded form. Onto this accessible DNA sequence a circular oligonucleotide probe is assembled, which serves as a template for rolling circle amplification (RCA). Both homogeneous and heterogeneous assay formats are investigated, as are different formats for fluorescence-based amplicon detection. Our recent data with immobilized analytes suggest that marker sequences in plasmid and bacterial chromosomal DNA can be successfully detected.     

Genomic DNA fingerprinting of clinical isolates of Helicobacter pylori using short oligonucleotide probes containing repetitive sequences

The ability of oligonucleotide probes containing short repetitive sequence motifs to differentiate between isolates of Helicobacter pylori was investigated. Genomic DNA preparations from H. pylori were digested with the restriction enzyme Hind III, electrophoresed in agarose gels and transferred to nylon filters. Five separate oligonucleotide probes were tested for hybridization sequentially to fingerprint the digested DNA from a panel of 29 clinical isolates and one type strain of H. pylori , and their relative discriminatory abilities were assessed. Four probes, (GACA)₄, (GT)₈, (GTG)₅ and (GGAT)₄, were each shown to yield highly informative hybridization band profiles allowing differentiation of H. pylori isolates. The DNA fingerprints of individual isolates obtained with each probe were distinct and reproducible. Direct comparison with ribotyping revealed that oligonucleotide fingerprinting had far superior discriminatory power. Computer-assisted similarity analysis of (GGAT)₄-generated hybridization profiles of pairwise combinations of H. pylori isolates revealed that there was no correlation between ribotype and oligonucleotide fingerprint patterns. The results of this study demonstrate that oligonucleotide probes containing microsatellite sequences provide a new and powerful tool for isolate discrimination of H. pylori .     

Affinity capture of specific DNA fragments with short synthetic sequences

The ability of short peptide nucleic acid (PNA) oligomers and oligonucleotides containing modified residues of 5-methylcitidine, 2-aminoadenosine, and 5-propynyl-2′-deoxyuridine (strong binding oligonucleotides, SBO) to affinity capture the target double-stranded DNA fragment from mixture by means of the end invasion was compared. Both types of probes were highly effective at the conditions used. The SBO-based probes may represent a handy and easily prepared alternative to PNA for selection of target DNA fragments in mixtures.     

Rapid determination of short DNA sequences by the use of MALDI-MS

We have developed a protocol for rapid sequencing of short DNA stretches (15–20 nt) using MALDI-TOF-MS. The protocol is based on the Sanger concept with the modification that double-stranded template DNA is used and all four sequencing reactions are performed in one reaction vial. The sequencing products are separated and detected by MALDI-TOF-MS and the sequence is determined by comparing measured molecular mass differences to expected values. The protocol is optimized for low costs and broad applicability. One reaction typically includes 300 fmol template, 10 pmol primer and 200 pmol each nucleotide monomer. Neither the primer nor any of the nucleotide monomers are labeled. Solid phase purification, concentration and mass spectrometric sample preparation of the sequencing products are accomplished in a few minutes and parallel processing of 96 samples is possible. The mass spectrometric analyses and subsequent sequence read-out require only a few seconds per template.     

Detection byin-situ fluorescence of short,single-copy sequences of chromosomal DNA

Short single-copy probes have been widely used in plant molecular biology. However, they have rarely been effective in plant research usingin situ hybridization techniques, possibly due to limitations imposed by the cell wall. We recently developed two fluorescencein situ hybridization protocols for the single-copy sequence detection in soybean. By enzymatically removing the cell wall, single-copy sequences as short as 1 kb were detected by probes using standard fluorescencein situ hybridization or PCR-primedin situ hybridization (PCR-PRINS). Such technology is useful for genome analysis, in plant molecular, cellular, and biotechnological research.     

Identifying fern gametophytes using DNA sequences

Identification of fern gametophytes is generally hampered by low morphological complexity. Here we explore an alternative: DNA‐based identification. We obtained a plastid rbcL sequence from a sterile gametophyte of unknown origin (cultivated for more than 30 years) and employed blast to determine its affinities. Using this approach, we identified the gametophyte as Osmunda regalis. To evaluate the robustness of this determination, and the usefulness of rbcL in differentiating among species, we conducted a phylogenetic analysis of osmundaceous fern sequences. Based on our results, it is evident that DNA‐based identification has considerable potential in exploring the ecology of fern gametophytes.     

Stretching short sequences of DNA with constant force axial optical tweezers

Single-molecule techniques for stretching DNA of contour lengths less than a kilobase are fraught with experimental difficulties. However, many interesting biological events such as histone binding and protein-mediated looping of DNA, occur on this length scale. In recent years, the mechanical properties of DNA have been shown to play a significant role in fundamental cellular processes like the packaging of DNA into compact nucleosomes and chromatin fibers. Clearly, it is then important to understand the mechanical properties of short stretches of DNA. In this paper, we provide a practical guide to a single-molecule optical tweezing technique that we have developed to study the mechanical behavior of DNA with contour lengths as short as a few hundred basepairs. The major hurdle in stretching short segments of DNA is that conventional optical tweezers are generally designed to apply force in a direction lateral to the stage (see Fig. 1). In this geometry, the angle between the bead and the coverslip, to which the DNA is tethered, becomes very steep for submicron length DNA. The axial position must now be accounted for, which can be a challenge, and, since the extension drags the microsphere closer to the coverslip, steric effects are enhanced. Furthermore, as a result of the asymmetry of the microspheres, lateral extensions will generate varying levels of torque due to rotation of the microsphere within the optical trap since the direction of the reactive force changes during the extension. Alternate methods for stretching submicron DNA run up against their own unique hurdles. For instance, a dual-beam optical trap is limited to stretching DNA of around a wavelength, at which point interference effects between the two traps and from light scattering between the microspheres begin to pose a significant problem. Replacing one of the traps with a micropipette would most likely suffer from similar challenges. While one could directly use the axial potential to stretch the DNA, an active feedback scheme would be needed to apply a constant force and the bandwidth of this will be quite limited, especially at low forces. We circumvent these fundamental problems by directly pulling the DNA away from the coverslip by using a constant force axial optical tweezers. This is achieved by trapping the bead in a linear region of the optical potential, where the optical force is constant-the strength of which can be tuned by adjusting the laser power. Trapping within the linear region also serves as an all optical force-clamp on the DNA that extends for nearly 350 nm in the axial direction. We simultaneously compensate for thermal and mechanical drift by finely adjusting the position of the stage so that a reference microsphere stuck to the coverslip remains at the same position and focus, allowing for a virtually limitless observation period.     

Comparison of different melting temperature calculation methods for short DNA sequences

MOTIVATION: The overall performance of several molecular biology techniques involving DNA/DNA hybridization depends on the accurate prediction of the experimental value of a critical parameter: the melting temperature Tm. Till date, many computer software programs based on different methods and/or parameterizations are available for the theoretical estimation of the experimental Tm value of any given short oligonucleotide sequence. However, in most cases, large and significant differences in the estimations of Tm were obtained while using different methods. Thus, it is difficult to decide which Tm value is the accurate one. In addition, it seems that most people who use these methods are unaware about the limitations, which are well described in the literature but not stated properly or restricted the inputs of most of the web servers and standalone software programs that implement them. RESULTS: A quantitative comparison on the similarities and differences among some of the published DNA/DNA Tm calculation methods is reported. The comparison was carried out for a large set of short oligonucleotide sequences ranging from 16 to 30 nt long, which span the whole range of CG-content. The results showed that significant differences were observed in all the methods, which in some cases depend on the oligonucleotide length and CG-content in a non-trivial manner. Based on these results, the regions of consensus and disagreement for the methods in the oligonucleotide feature space were reported. Owing to the lack of sufficient experimental data, a fair and complete assessment of accuracy for the different methods is not yet possible. Inspite of this limitation, a consensus Tm with minimal error probability was calculated by averaging the values obtained from two or more methods that exhibit similar behavior to each particular combination of oligonucleotide length and CG-content class. Using a total of 348 DNA sequences in the size range between 16mer and 30mer, for which the experimental Tm data are available, we demonstrated that the consensus Tm is a robust and accurate measure. It is expected that the results of this work would be constituted as a useful set of guidelines to be followed for the successful experimental implementation of various molecular biology techniques, such as quantitative PCR, multiplex PCR and the design of optimal DNA microarrays.     

Features in short guanine-rich sequences that stimulate DNA polymerization in vitro

We have discovered that short guanine-rich oligonucleotides are able to self-associate into higher order structures that stimulate DNA synthesis in vitro without the addition of a conventional template [Ying, J., Bradley, R. K., Jones, L. B., Reddy, M. S., Colbert, D. T., Smalley, R. E., and Hardin, S. H. (1999) Biochemistry 38, 16461-16468]. Our initial analysis indicated the importance of the presence of three contiguous guanines (G) in an oligonucleotide that stimulates DNA polymerization. To gain insight into and to refine sequence requirements for the unexpected DNA synthesis, we analyzed a 231-member guanine-rich octamer library in a fluorescent nucleotide polymerization assay. We observe that, in addition to three contiguous Gs, the presence of a secondary G cluster within the octamer is essential. Furthermore, the location of the primary G cluster in the center of the molecule is most stimulatory. The majority of the octamers that form extended DNA products have a single non-G base separating the primary and secondary G clusters, the identity of which is predominantly thymine (T). Further, a T 5' or 3' of the primary G cluster positively influences the stimulatory function of the oligonucleotide. Overall, the occurrence of bases in the octamer is in the descending order of G > T > A > C. Our studies demonstrate that structures stabilized by noncanonical base pairings are recognized by a DNA polymerase in vitro, and these findings may have relevance within the cell. In particular, the features of these G-rich stimulatory sequences show striking similarities to telomeric sequences that form diverse G-quartet structures in vitro.     

Site2genome: locating short DNA sequences in whole genomes

SUMMARY: Many biological papers describe short, functional DNA sites without specifying their exact positions in the genome. We have developed a Web server that automates the tedious task of locating such sites in eukaryotic genomes, thus giving access to the context of rich annotations that are increasingly available for genome sequences. AVAILABILITY: http://zlab.bu.edu/site2genome/     

Ultrafast and memory-efficient alignment of short DNA sequences to the human genome

Bowtie is an ultrafast, memory-efficient alignment program for aligning short DNA sequence reads to large genomes. For the human genome, Burrows-Wheeler indexing allows Bowtie to align more than 25 million reads per CPU hour with a memory footprint of approximately 1.3 gigabytes. Bowtie extends previous Burrows-Wheeler techniques with a novel quality-aware backtracking algorithm that permits mismatches. Multiple processor cores can be used simultaneously to achieve even greater alignment speeds. Bowtie is open source .     

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by ±10° in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3⁷) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150–325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CX_nG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5′ end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10⁶ conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar–phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques. © 1993 John Wiley & Sons, Inc.     

Simulations using random-generated DNA and RNA sequences

Using a very simple computer program written in BASIC, a very large number of random-generated DNA or RNA sequences are obtained. Students individually and as a group use these sequences to perform a wide variety of simulated events including predicting complementary sequences and trans-lational products, evaluating base compositions, determining frequencies of particular triplet codons, and suggesting possible secondary structures.     

Statistical assignment of DNA sequences using Bayesian phylogenetics

We provide a new automated statistical method for DNA barcoding based on a Bayesian phylogenetic analysis. The method is based on automated database sequence retrieval, alignment, and phylogenetic analysis using a custom-built program for Bayesian phylogenetic analysis. We show on real data that the method outperforms Blast searches as a measure of confidence and can help eliminate 80% of all false assignment based on best Blast hit. However, the most important advance of the method is that it provides statistically meaningful measures of confidence. We apply the method to a re-analysis of previously published ancient DNA data and show that, with high statistical confidence, most of the published sequences are in fact of Neanderthal origin. However, there are several cases of chimeric sequences that are comprised of a combination of both Neanderthal and modern human DNA.     

Detection of short repeated genomic sequences on metaphase chromosomes using padlock probes and target primed rolling circle DNA synthesis

Background  

In situ detection of short sequence elements in genomic DNA requires short probes with high molecular resolution and powerful specific signal amplification. Padlock probes can differentiate single base variations. Ligated padlock probes can be amplified in situ by rolling circle DNA synthesis and detected by fluorescence microscopy, thus enhancing PRINS type reactions, where localized DNA synthesis reports on the position of hybridization targets, to potentially reveal the binding of single oligonucleotide-size probe molecules. Such a system has been presented for the detection of mitochondrial DNA in fixed cells, whereas attempts to apply rolling circle detection to metaphase chromosomes have previously failed, according to the literature.     

Detection of specific DNA sequences using antibodies recognizing UV-labelled DNA.

This non-isotopic method for detection of nucleic acids is based on the in situ labelling of the nucleic acid by exposure to UV-irradiation. The different UV-induced photoproducts, mainly of the thymidine dimer type, are recognized by purified rabbit antibodies specific to the lesions introduced. The UV-labelled nucleic acids can then be visualized by conventional immunostaining procedures. A major advantage of the technique is the low cost and the ease by which the DNA is specifically labelled. The purified rabbit antibodies were shown to be specific for UV-irradiated DNA, and the method was applied for detection of specific DNA sequences hybridized to homologous target DNA on membrane support. We believe that the sensitivity of the method can be improved, and the significance of using different UV-doses, immunostaining methods and membrane types is discussed.     

Enforcement of late replication origin firing by clusters of short G-rich DNA sequences

Previous studies in budding yeast suggested that the default firing time of most DNA replication origins is early in S phase and that origins can be forced to fire later by proximity to certain cis-acting sequences. However, these cis-acting sequences were not well defined. We have attempted to characterize cis-acting sequences that affect replication timing in the fission yeast. We identified a stretch of 200 bp that was sufficient to compel nearby origins to fire late. The 200-bp stretch was able to force an origin to fire late whether adjacent to the origin or approximately 800 bp away in opposite orientation. The stretch contains a cluster of three close matches to a G-rich, 10-bp late consensus sequence (LCS). The three LCS elements cooperate with each other and with other sequences within the 200-bp stretch to enforce late replication. Although only a few origins that fire in very late S phase have been identified in fission yeast, all of them are located close to a cluster of LCS elements.