Free energy of DNA duplex formation on short oligonucleotide microarrays

DNA/DNA duplex formation is the basic mechanism that is used in genome tiling arrays and SNP arrays manufactured by Affymetrix. However, detailed knowledge of the physical process is still lacking. In this study, we show a free energy analysis of DNA/DNA duplex formation these arrays based on the positional-dependent nearest-neighbor (PDNN) model, which was developed previously for describing DNA/RNA duplex formation on expression microarrays. Our results showed that the two ends of a probe contribute less to the stability of the duplexes and that there is a microarray surface effect on binding affinities. We also showed that free energy cost of a single mismatch depends on the bases adjacent to the mismatch site and obtained a comprehensive table of the cost of a single mismatch under all possible combination of adjacent bases. The mismatch costs were found to be correlated with those determined in aqueous solution. We further demonstrate that the DNA copy number estimated from the SNP array correlates negatively with the target length; this is presumably caused by inefficient PCR amplification for long fragments. These results provide important insights into the molecular mechanisms of microarray technology and have implications for microarray design and the interpretation of observed data.     

Molecular evolution and phylogenetic implications of internal transcribed spacer sequences of Ribosomal DNA in Winteraceae

The internal transcribed spacers and the 5.8S coding region of nuclear ribosomal DNA were sequenced and analyzed to address questions of generic relationships in Winteraceae. The molecular data generated a single tree that is congruent with one based on morphological data. The sequences of ITS 1 in the family range from 235 to 252 bases in size and of ITS 2 from 213 to 226 bases. The size of the 5.8S coding region is 164 bases. The range of ITS 1 and ITS 2 sequence divergence between pairs of genera within Winteraceae is relatively low in comparison to other plant families. Two types of ITS 1 and ITS 2 sequences were observed in the same individual for some taxa. Sequence variations between the two arrays are 4.7%–6.3% for ITS 1 and 5.1%–7.0% for ITS 2. Both arrays of sequences, however, generate the same phylogenetic relationships. Rates of nucleotide substitutions for the internal transcribed spacers are 3.2–5.2 × 10^-10 substitution per site per year estimated in ITS 1 and 3.6–5.7 × 10^-10 in ITS 2.     

Multiple copies of a retroposon interrupt spliced leader RNA genes in the African trypanosome, Trypanosoma gambiense.

The 140-nucleotide spliced leader (SL) RNA, involved in mRNA maturation in the African trypanosomes and in other kinetoplastida, is encoded by a tandem array of spliced leader genes. We show that the 1.4-kb SL gene repeat unit in Trypanosoma gambiense is organized in tandem arrays confined to two large (minimum size 350-450 kb) restriction fragments. SL genes in both arrays are interrupted by a total of eight conserved insertion elements. Cleavage of genomic DNA at restriction sites present within the insertion element but not in the SL gene repeat, releases variable numbers of SL genes from the tandem array. Since the insertion element contains a terminal poly(A) track of 36 bases and because a 49-bp duplication of target DNA has occurred at the integration site, we conclude that it is a retroposon. This retropson is uniquely associated with the SL gene clusters. These retroposons presumably originated from a single insertion event after which their copy number increased, possibly through unequal sister chromatid exchange.     

A measure of DNA periodicity

A Fourier transform g(n) of a sequence of bases along a given stretch of DNA is defined. The transform is invariant to the labelling of the bases and can therefore be used as a measure of periodicity for segments of DNA with differing base content. It can also be conveniently used to search for base periodicities within large DNA data bases.     

Microevolution of the mitochondrial DNA control region in the Japanese brown bear (Ursus arctos) population.

We investigated nucleotide sequences of the mitochondrial DNA control region to describe natural genetic variations and to assess the relationships between subpopulations of the brown bear Ursus arctos on Hokkaido Island, Japan. Using the polymerase chain reaction product-direct sequencing technique, partial sequences (about 930 bases) of the control region were determined for 56 brown bears sampled throughout Hokkaido Island. A sequence alignment revealed that the brown bear control region included a variable sequence on the 5' side and a repetitive region on the 3' side. Phylogenetic trees reconstructed from the 5' variable region (696-702 bases) exhibited 17 haplotypes, which were clustered into three groups (Clusters A, B, and C). The distribution of each group did not overlap with those of the others, and the three different areas were located in separate mountainous forests of Hokkaido Island. Furthermore, most of the phylogenetically close haplotypes within each group were distributed geographically close to each other. In addition, the 3' repetitive region (arrays of 10 bases) exhibited a much faster mutation rate than the 5' variable region, resulting in heteroplasmy. Such mitochondrial DNA divergence in each group could have occurred after the brown bears migrated from the continent to Hokkaido and became fixed in the different areas.     

DNA packaging and organization in mammalian spermatozoa: comparison with somatic cells

Mammalian sperm DNA is the most tightly compacted eukaryotic DNA, being at least sixfold more highly condensed than the DNA in mitotic chromosomes. To achieve this high degree of packaging, sperm DNA interacts with protamines to form linear, side-by-side arrays of chromatin. This differs markedly from the bulkier DNA packaging of somatic cell nuclei and mitotic chromosomes, in which the DNA is coiled around histone octamers to form nucleosomes. The overall organization of mammalian sperm DNA, however, resembles that of somatic cells in that both the linear arrays of sperm chromatin and the 30-nm solenoid filaments of somatic cell chromatin are organized into loop domains attached at their bases to a nuclear matrix. In addition to the sperm nuclear matrix, sperm nuclei contain a unique structure termed the sperm nuclear annulus to which the entire complement of DNA appears to be anchored when the nuclear matrix is disrupted during decondensation. In somatic cells, proper function of DNA is dependent upon the structural organization of the DNA by the nuclear matrix, and the structural organization of sperm DNA is likely to be just as vital to the proper functioning of the spermatozoa.     

Hybridization of DNA targets to glass-tethered oligonucleotide probes

Hybridization of nucleic acids to surface-tethered oligonucleotide probes has numerous potential applications in genome mapping and DNA sequence analysis. In this article, we describe a simple standard protocol for routine preparation of terminal amine-derivatized 9-mer oligonucleotide arrays on ordinary microscope slides and hybridization conditions with DNA target strands of up to several hundred bases in length with good discrimination against mismatches. Additional linker arms separating the glass surface from the probe sequence are not necessary. The technique described here offers a powerful tool for the detection of specific genetic mutations.     

Assessment of oxidative base damage to isolated and cellular DNA by HPLC-MS/MS measurement

Oxidation reactions that involve several oxygen and nitrogen reactive species together with nucleobase radical cations give rise among various classes of lesions to modified bases. About 70 of oxidized nucleosides that include diastereomeric forms have been characterized in mechanistic studies involving isolated DNA and related model compounds. However, only eight modified bases have been accurately measured within cellular DNA upon exposure to either gamma or UVA radiations. Emphasis is placed in this survey on recent developments of HPLC associated with tandem mass spectrometry (MS/MS) operating in the mild electrospray ionization mode. Interestingly, the HPLC-MS/MS assay in the multiple reaction monitoring mode appears to be the more sensitive and accurate method currently available for singling out several oxidized nucleosides including 8-oxo-7,8-dihydro-2'-deoxyguanosine, 8-oxo-7,8-dihydro-2'-deoxyadenosine, 5-formyl-2'-deoxyuridine, 5-(hydroxymethyl-2'-deoxyuridine, 5-hydroxy-2'-deoxyuridine, and the four diastereomers of 5,6-dihydroxy-5,6-dihydrothymidine within isolated and cellular DNA. However, one limitation of the assay that also applied to all chromatographic methods is the slight side-oxidation of normal bases during DNA extraction and subsequent work-up. This explains why the combined use of DNA repair glycosylases with either the comet assay or the alkaline elution technique is a better alternative to monitor the formation of low levels of oxidized bases within cellular DNA.     

An analysis of the relationship between hydration and protein-DNA interactions.

Eleven protein-DNA crystal structures were analyzed to test the hypothesis that hydration sites predicted in the first hydration shell of DNA mark the positions where protein residues hydrogen-bond to DNA. For nine of those structures, protein atoms, which form hydrogen bonds to DNA bases, were found within 1.5 A of the predicted hydration positions in 86% of the interactions. The correspondence of the predicted hydration sites with the hydrogen-bonded protein side chains was significantly higher for bases inside the conserved DNA recognition sequences than outside those regions. In two CAP-DNA complexes, predicted base hydration sites correctly marked 71% (within 1.5 A) of protein atoms, which form hydrogen bonds to DNA bases. Phosphate hydration was compared to actual protein binding sites in one CAP-DNA complex with 78% marked contacts within 2.0 A. These data suggest that hydration sites mark the binding sites at protein-DNA interfaces.     

A zinc finger protein array for the visual detection of specific DNA sequences for diagnostic applications

The visual detection of specific double-stranded DNA sequences possesses great potential for the development of diagnostics. Zinc finger domains provide a powerful scaffold for creating custom DNA-binding proteins that recognize specific DNA sequences. We previously demonstrated sequence-enabled reassembly of TEM-1 β-lactamase (SEER–LAC), a system consisting of two inactive fragments of β-lactamase each linked to engineered zinc finger proteins (ZFPs). Here the SEER–LAC system was applied to develop ZFP arrays that function as simple devices to identify bacterial double-stranded DNA sequences. The ZFP arrays provided a quantitative assay with a detection limit of 50 fmol of target DNA. The method could distinguish target DNA from non-target DNA within 5 min. The ZFP arrays provided sufficient sensitivity and high specificity to recognize specific DNA sequences. These results suggest that ZFP arrays have the potential to be developed into a simple and rapid point-of-care (POC) diagnostic for the multiplexed detection of pathogens.     

Analysis of difference spectra of protonated DNA: determination of degree of protonation of nitrogen bases and the fractions of disordered nucleotide pairs.

The titration curves of nitrogen bases and fractions of disordered nucleotide pairs are obtained during DNA protonation. It is shown that purine bases are the first sites of the DNA double helix protonation. The cytosine protonation is due to proton-induced conformational transition within GC pairs with the sequence proton transfer from (N-7) of guanine to (N-3) of cytosine. Within DNA with unwound regions the bases are protonated in the following order: cytosine, adenine, guanine. It is shown that GC pairs are the primary centres in which the unwinding of protonated DNAs occurs.     

Human rDNA: evolutionary patterns within the genes and tandem arrays derived from multiple chromosomes

Human rDNA forms arrays on five chromosome pairs and is homogenized by concerted evolution through recombination and gene conversion (loci RNR1, RNR2, RNR3, RNR4, RNR5, OMIM: 180450). Homogenization is not perfect, however, so that it becomes possible to study its efficiency within genes, within arrays, and between arrays by measuring and comparing DNA sequence variation. Previous studies with randomly cloned genomic DNA fragments showed that different parts of the gene evolve at different rates but did not allow comparison of rDNA sequences derived from specific chromosomes. We have now cloned and sequenced rDNA fragments from specific acrocentric chromosomes to (1) study homogenization along the rDNA and (2) compare homogenization within chromosomes and between homologous and nonhomologous chromosomes. Our results show high homogeneity among regulatory and coding regions of rDNA on all chromosomes, a surprising homogeneity among adjacent distal non-rDNA sequences, and the existence of one to three very divergent intergenic spacer classes within each array.     

Handling long targets and errors in sequencing by hybridization.

Sequencing by hybridization (SBH) is a DNA sequencing technique, in which the sequence is reconstructed using its k-mer content. This content, which is called the spectrum of the sequence, is obtained by hybridization to a universal DNA array. Standard universal arrays contain all k-mers for some fixed k, typically 8 to 10. Currently, in spite of its promise and elegance, SBH is not competitive with standard gel-based sequencing methods. This is due to two main reasons: lack of tools to handle realistic levels of hybridization errors and an inherent limitation on the length of uniquely reconstructible sequence by standard universal arrays. In this paper, we deal with both problems. We introduce a simple polynomial reconstruction algorithm which can be applied to spectra from standard arrays and has provable performance in the presence of both false negative and false positive errors. We also propose a novel design of chips containing universal bases that differs from the one proposed by Preparata et al. (1999). We give a simple algorithm that uses spectra from such chips to reconstruct with high probability random sequences of length lower only by a squared log factor compared to the information theoretic bound. Our algorithm is very robust to errors and has a provable performance even if there are both false negative and false positive errors. Simulations indicate that its sensitivity to errors is also very small in practice.