D5S2500 is an ambiguously characterized STR: Identification and description of forensic microsatellites in the genomics age.

In the process of establishing short tandem repeat (STR) sequence variant nomenclature guidelines in anticipation of expanded forensic multiplexes for massively parallel sequencing (MPS), it was discovered that the STR D5S2500 has multiple positions and genomic characteristics reported. This ambiguity is because the marker named D5S2500 consists of two different microsatellites forming separate components in the capillary electrophoresis multiplexes of Qiagen's HDplex (Hilden, Germany) and AGCU ScienTech's non-CODIS STR 21plex (Wuxi, Jiangsu, China). This study outlines the genomic details used to identify each microsatellite and reveals the D5S2500 marker in HDplex has the correctly assigned STR name, while the D5S2500 marker in the AGCU 21plex, closely positioned a further 1643 nucleotides in the human reference sequence, is an unnamed microsatellite. The fact that the D5S2500 marker has existed as two distinct STR loci undetected for almost ten years, even with reported discordant genotypes for the standard control DNA, underlines the need for careful scrutiny of the genomic properties of forensic STRs, as they become adapted for sequence analysis with MPS systems. We make the recommendation that precise chromosome location data must be reported for any forensic marker under development but not in common use, so that the genomic characteristics of the locus are validated to the same level of accuracy as its allelic variation and forensic performance. To clearly differentiate each microsatellite, we propose the name D5S2800 be used to identify the Chromosome-5 STR in the AGCU 21plex.

Forensic ancestry analysis with two capillary electrophoresis ancestry informative marker (AIM) panels: Results of a collaborative EDNAP exercise.

There is increasing interest in forensic ancestry tests, which are part of a growing number of DNA analyses that can enhance routine profiling by obtaining additional genetic information about unidentified DNA donors. Nearly all ancestry tests use single nucleotide polymorphisms (SNPs), but these currently rely on SNaPshot single base extension chemistry that can fail to detect mixed DNA. Insertion-deletion polymorphism (Indel) tests have been developed using dye-labeled primers that allow direct capillary electrophoresis detection of PCR products (PCR-to-CE). PCR-to-CE maintains the direct relationship between input DNA and signal strength as each marker is detected with a single dye, so mixed DNA is more reliably detected. We report the results of a collaborative inter-laboratory exercise of 19 participants (15 from the EDNAP European DNA Profiling group) that assessed a 34-plex SNP test using SNaPshot and a 46-plex Indel test using PCR-to-CE. Laboratories were asked to type five samples with different ancestries and detect an additional mixed DNA sample. Statistical inference of ancestry was made by participants using the Snipper online Bayes analysis portal plus an optional PCA module that analyzes the genotype data alongside calculation of Bayes likelihood ratios. Exercise results indicated consistent genotyping performance from both tests, reaching a particularly high level of reliability for the Indel test. SNP genotyping gave 93.5% concordance (compared to the organizing laboratory's data) that rose to 97.3% excluding one laboratory with a large number of miscalled genotypes. Indel genotyping gave a higher concordance rate of 99.8% and a reduced no-call rate compared to SNP analysis. All participants detected the mixture from their Indel peak height data and successfully assigned the correct ancestry to the other samples using Snipper, with the exception of one laboratory with SNP miscalls that incorrectly assigned ancestry of two samples and did not obtain informative likelihood ratios for a third. Therefore, successful ancestry assignments were achieved by participants in 92 of 95 Snipper analyses. This exercise demonstrates that ancestry inference tests based on binary marker sets can be readily adopted by laboratories that already have well-established CE regimes in place. The Indel test proved to be easy to use and allowed all exercise participants to detect the DNA mixture as well as achieving complete and concordant profiles in nearly all cases. Lastly, two participants successfully ran parallel next-generation sequencing analyses (each using different systems) and achieved high levels of genotyping concordance using the exercise PCR primer mixes unmodified.

Revision of the SNPforID 34-plex forensic ancestry test: Assay enhancements, standard reference sample genotypes and extended population studies.

A revision of an established 34 SNP forensic ancestry test has been made by swapping the under-performing rs727811 component SNP with the highly informative rs3827760 that shows a near-fixed East Asian specific allele. We collated SNP variability data for the revised SNP set in 66 reference populations from 1000 Genomes and HGDP-CEPH panels and used this as reference data to analyse four U.S. populations showing a range of admixture patterns. The U.S. Hispanics sample in particular displayed heterogeneous values of co-ancestry between European, Native American and African contributors, likely to reflect in part, the way this disparate group is defined using cultural as well as population genetic parameters. The genotyping of over 700 U.S. population samples also provided the opportunity to thoroughly gauge peak mobility variation and peak height ratios observed from routine use of the single base extension chemistry of the 34-plex test. Finally, the genotyping of the widely used DNA profiling Standard Reference Material samples plus other control DNAs completes the audit of the 34-plex assay to allow forensic practitioners to apply this test more readily in their own laboratories.

Forensic performance of two insertion-deletion marker assays.

Improving the amplification and analysis of highly degraded DNA extracts has been a longstanding area of research in forensic genetics. One of the most promising recent developments in analysis of degraded DNA is the availability of short, biallelic insertion-deletion length polymorphisms (InDels) in highly multiplexed assays. InDels share many of the favourable characteristics of single-nucleotide polymorphisms (SNPs) that make them ideal markers for analysis of degraded DNA, including: analysis in short amplicon size ranges, high multiplexing capability and low mutation rates. In addition, as length-based polymorphisms, InDels can be analysed with the same simple dye-labelled PCR primer methods as standard forensic short tandem repeats. Separation and detection of fluorescently dye-labelled PCR products by capillary electrophoresis eliminate the multiple step protocols required by SNP typing with single-base extension assays and provide a closer relationship between the input DNA and the profile peak height ratios. Therefore InDel genotyping represents an effective new approach for human identification that adds informative new loci to the existing battery of forensic markers. To assess the utility of InDels for forensic analysis, we characterised population variation with two InDel identification assays: the 30-plex Qiagen DIPplex panel and a 38-plex panel developed by Pereira et al. in 2009. Allele frequencies were generated for the 68 markers in US African American, Caucasian, East Asian and Hispanic samples. We made a thorough assessment of the individual and combined performance of the InDel sets, as well as characterising profile artifacts and other issues related to the routine use of these newly developed forensic assays based on artificially degraded DNA and mixed source samples.

Autosomal SNP typing of forensic samples with the GenPlex™ HID System: results of a collaborative study.

The GenPlex™ HID System (Applied Biosystems - AB) offers typing of 48 of the 52 SNPforID SNPs and amelogenin. Previous studies have shown a high reproducibility of the GenPlex™ HID System using 250-500pg DNA of good quality. An international exercise was performed by 14 laboratories (9 in Europe and 5 in the US) in order to test the robustness and reliability of the GenPlex™ HID System on forensic samples. Three samples with partly degraded DNA and 10 samples with low amounts of DNA were analyzed in duplicates using various amounts of DNA. In order to compare the performance of the GenPlex™ HID System with the most commonly used STR kits, 500pg of partly degraded DNA from three samples was typed by the laboratories using one or more STR kits. The median SNP typing success rate was 92.3% with 500pg of partly degraded DNA. Three of the fourteen laboratories counted for more than two thirds of the locus dropouts. The median percentage of discrepant results was 0.2% with 500pg degraded DNA. An increasing percentage of locus dropouts and discrepant results were observed when lower amounts of DNA were used. Different success rates were observed for the various SNPs. The rs763869 SNP was the least successful. With the exception of the MiniFiler™ kit (AB), GenPlex™ HID performed better than five other tested STR kits. When partly degraded DNA was analyzed, GenPlex™ HID showed a very low mean mach probability, while all STR kits except MiniFiler™ had very limited discriminatory power.

Forensic typing of autosomal SNPs with a 29 SNP-multiplex--results of a collaborative EDNAP exercise.

We report the results of an inter-laboratory exercise on typing of autosomal single nucleotide polymorphisms (SNP) for forensic genetic investigations in crime cases. The European DNA Profiling Group (EDNAP), a working group under the International Society for Forensic Genetics (ISFG), organised the exercise. A total of 11 European and one US forensic genetic laboratories tested a subset of a 52 SNP-multiplex PCR kit developed by the SNPforID consortium. The 52 SNP-multiplex kit amplifies 52 DNA fragments with 52 autosomal SNP loci in one multiplex PCR. The 52 SNPs are detected in two separate single base extension (SBE) multiplex reactions with 29 and 23 SNPs, respectively, using SNaPshot kit, capillary electrophoresis and multicolour fluorescence detection. For practical reasons, only the 29 SBE multiplex reaction was carried out by the participating laboratories. A total of 11 bloodstains on FTA cards including a sample of poor quality and a negative control were sent to the laboratories together with the essential reagents for the initial multiplex PCR and the multiplex SBE reaction. The total SNP locus dropout rate was 2.8% and more than 50% of the dropouts were observed with the poor quality sample. The overall rate of discrepant SNP allele assignments was 2.0%. Two laboratories reported 60% of all the discrepancies. Two laboratories reported all 29 SNP alleles in all 10 positive samples correctly. The results of the collaborative exercise were surprisingly good and demonstrate that SNP typing with SBE, capillary electrophoresis and multicolour detection methods can be developed for forensic genetics.

Forensic application of the Affymetrix human mitochondrial resequencing array.

In the field of forensic DNA testing, sequencing regions of the mitochondrial genome is performed when insufficient genomic DNA is present for traditional autosomal short tandem repeat (STR) testing. Sequencing coding region polymorphisms in the mitochondrial genome can be useful for resolving individuals who have the identical HV1 and HV2 control region sequence. Various methods and strategies have been established to interrogate coding region polymorphisms. These range from SNP assays probing sites most likely to differentiate individuals based on their HV1/HV2 sequence to the use of mass spectrometry to pyrosequencing. Here we evaluate the potential of the Affymetrix GeneChip Mitochondrial Resequencing Array (version 2.0) for forensic applications.

The impact of additional Y-STR loci on resolving common haplotypes and closely related individuals.

Commercial Y-STR kits have permitted laboratories to go beyond the original nine minimal haplotype loci (MHL) and to discover the advantage of additional Y-STR loci in resolving common haplotypes. In an effort to examine the impact of Y-STR markers beyond the 17 loci now available in commercial kit form, new Y-STR loci are being investigated on a common set of samples representative of the major U.S. population groups. Additional Y-STRs can also increase the power of discrimination between closely related male individuals, which is important not only in forensics but also in the paternity and genetic genealogy communities.

Analysis of artificially degraded DNA using STRs and SNPs--results of a collaborative European (EDNAP) exercise.

Recently, there has been much debate about what kinds of genetic markers should be implemented as new core loci that constitute national DNA databases. The choices lie between conventional STRs, ranging in size from 100 to 450 bp; mini-STRs, with amplicon sizes less than 200 bp; and single nucleotide polymorphisms (SNPs). There is general agreement by the European DNA Profiling Group (EDNAP) and the European Network of Forensic Science Institutes (ENFSI) that the reason to implement new markers is to increase the chance of amplifying highly degraded DNA rather than to increase the discriminating power of the current techniques. A collaborative study between nine European and US laboratories was organised under the auspices of EDNAP. Each laboratory was supplied with a SNP multiplex kit (Foren-SNPs) provided by the Forensic Science Service, two mini-STR kits provided by the National Institute of Standards and Technology (NIST) and a set of degraded DNA stains (blood and saliva). Laboratories tested all three multiplex kits, along with their own existing DNA profiling technique, on the same sets of degraded samples. Results were collated and analysed and, in general, mini-STR systems were shown to be the most effective. Accordingly, the EDNAP and ENFSI working groups have recommended that existing STR loci are reengineered to provide smaller amplicons, and the adoption of three new European core loci has been agreed.

Calculating sequence-dependent melting stability of duplex DNA oligomers and multiplex sequence analysis by graphs.

The analytical methods for characterizing DNA sequence-dependent thermodynamic stability have been reviewed. A set of n-n sequence stability parameters is presented. Examples in which these values are used to calculate the thermodynamic stability of short duplex DNA oligomers are presented. The problem of determining sets of isothermal sequences is addressed by representing DNA sequences as graphs. Representing DNA sequences by a graph descriptor with special mathematical properties minimizes the computational difficulty of determining the number of DNA sequences with identical predicted thermodynamic stability. This is achieved by replacement of a whole set of sequences by a single representative. Applications of this concept were demonstrated for sequences assembled from individual bases and sequences assembled from oligomeric blocks.

Quality control of PCR primers used in multiplex STR amplification reactions.

Reliable amplification of short tandem repeat (STR) DNA markers with the polymerase chain reaction (PCR) is dependent on high quality PCR primers. The particular primer combinations and concentrations are especially important with multiplex amplification reactions where multiple STR loci are simultaneously copied. Commercially available kits are now widely used for STR amplification and subsequent DNA typing. We present here the use of high performance liquid chromatography (HPLC) and time-of-flight mass spectrometry (TOF-MS) methods for characterization of commercially available STR kits.

Capillary electrophoresis as a tool for optimization of multiplex PCR reactions.

Copying multiple regions of a DNA molecule is routinely performed today using the polymerase chain reaction (PCR) in a process commonly referred to as multiplex PCR. The development of a multiplex PCR reaction involves designing primer sets and examining various combinations of those primer sets and different reaction components and/or thermal cycling conditions. The process of optimizing a multiplex PCR reaction in order to obtain a well-balanced set of amplicons can be time-consuming and labor-intensive. The rapid separation and quantitation capabilities of capillary electrophoresis make it an efficient technique to help in the multiplex PCR optimization process.

Genotyping of two mutations in the HFE gene using single-base extension and high-performance liquid chromatography.

Currently, a major focus of human genetics is the utilization of single-nucleotide polymorphisms for clinical diagnostics, whole-genome linkage disequilibrium screens to identify common disease genes such as Alzheimer disease, determination of the recent evolutionary history of a species, and the process of speciation. We have examined single-nucleotide extension coupled with high-performance liquid chromatography as a method to simultaneously genotype two SNPs occurring in the coding region of the HFE gene that produce clinical effects. This assay allows concurrent genotyping of the C282Y and H63D mutations in 11 min and is 100% concordant with current testing methods for both of these mutations.

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in four 40 base pair deoxyoligonucleotides.

Effects of different end sequences on melting, circular dichroism spectra (CD), and enzyme binding properties were investigated for four 40 base pair, non-self-complementary duplex DNA oligomers. The center sequences of these oligoduplexes have either of two 22 base pair modules flanked on both sides by sequences differing in AT content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from a van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming that the melting transition is two-state. Melting free energies (20 degrees C) evaluated from DSC melting experiments on the four duplex DNAs ranged from -52.2 to -77.5 kcal/mol. Free energies based on the van't Hoff analysis were -37.9 to -58.8 kcal/mol. Although the values are different, trends in the melting free energies of the four duplex DNAs as a function of sequence were identical in both DSC and optical analyses. Subject to several assumptions, values for the initiation free energy were estimated for each duplex, defined as DeltaG(int) = DeltaG(cal) - DeltaG(pred), where DeltaG(cal) is the experimental free energy at 20 degrees C determined from the experimentially measured values of the transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal). The predicted free energy of the sequence, DeltaG(pred)(20 degrees C), is obtained using published nearest-neighbor sequence stability values. For three of the four duplexes, values of DeltaG(int) are essentially nil. In contrast, the duplex with 81.8% GC has a considerably higher estimate of DeltaG(int) = 7.1 kcal/mol. The CD spectra for the six duplexes collected over the wavelength range from 200 to 320 nm are also sequence-dependent. Factor analysis of the CD spectra by singular value decomposition revealed that the experimental CD spectra could be reconstructed from linear combinations of two minor and one major subspectra. Changes in the coefficients of the major subspectrum for different sequences reflect incremental sequence-dependent variations of the CD spectra. Equilibrium binding by BamHI restriction endonuclease to the 40 base pair DNAs whose central eight base pairs contain the recognition sequence for BamHI restriction enzyme bounded by A.T base pairs, 5'-A-GGATCC-A-3' was investigated. Binding assays were performed by titering BamHI against a constant concentration of each of the duplex DNA substrates, in the absence of Mg(2+), followed by analysis by gel retardation. Under the conditions employed, the enzyme binds but does not cleave the DNAs. Results of the assays revealed two binding modes with retarded gel mobilities. Binding isotherms for the fraction of bound DNA species versus enzyme concentration for each binding mode were constructed and analyzed with a simple two-step equilibrium binding model. This analysis provided semiquantitative estimates on the equilibrium binding constants for BamHI to the four DNAs. Values obtained for the binding constants varied only 7-fold and ranged from 6 x 10(-)(8) to 42 x 10(-)(8) M, with binding free energies from -8.6 to -9.7 (+/- 0.2) kcal/mol depending on the sequence that flanks the enzyme binding site. Unlike what was found earlier in binding studies of the 22 base pair duplexes that constitute the core modules of the present 40-mers [Riccelli, P. V., Vallone, P. M., Kashin, I., Faldasz, B. D., Lane, M. J., and Benight, A. S. (1999) Biochemistry 38, 11197-11208], no obvious relationship between binding and stability was found for these longer DNAs. Apparently, effects of flanking sequence stability on restriction enzyme binding may only be measurable in very short duplex deoxyoligonucl

The role of the loop in binding of an actinomycin D analog to hairpins formed by single-stranded DNA.

Our recent work has indicated that the potent antibiotic and antitumor agent actinomycin D has the ability to selectively bind and stabilize single-stranded DNA that is capable of adopting a hairpin conformation. This mechanism of DNA binding has been implicated in the drug's ability to inhibit transcription by HIV reverse transcriptase from single-stranded DNA templates. In this report, we studied the importance of the hairpin loop on the ability of the 7-amino analog of actinomycin D to selectively bind DNA hairpins. Binding dissociation constant (Kd) values were determined to be 0.22 +/- 0.11 microM for the hairpin formed from the single-stranded DNA 5'-AAAAAAATAGTTTTAAATATTTTTTT-3' (dubbed HP1). The hairpin stem without the loop resulted in binding with Kd = 2.6 +/- 0.9 microM. The drug showed low affinity for the HP1 strand fully duplexed to its complementary sequence (estimated to be at least Kd > 21 microM). Evaluation of 7-aminoactinomycin D binding to a library of thermodynamically characterized DNA hairpins revealed an affinity for the hairpin-forming sequence 5'-GGATACCCCCGTATCC-3' (dubbed ACC4) of Kd = 6.8 +/- 2.2 microM. Replacement of the terminal guanines of this sequence to generate 5'-ATATACCCCCGTATAT-3' resulted in a 10-fold increase in affinity for this hairpin compared to ACC4, to Kd = 0.74 +/- 0.06 microM. A molecular model of the ACC4actinomycin D complex reveals that significant interactions between the hairpin loop and the pentapeptide rings of the drug must occur during drug binding. Taken together, our data indicate that the composition of the stem-loop interface is critical for the selectivity of actinomycin D and its 7-amino analog for DNA hairpins and suggests that novel drugs may be designed based on selection for the desired hairpin composition.

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in six 22-base pair deoxyoligonucleotides.

Effects of different end sequences on stability, circular dichroism spectra (CD), and enzyme binding properties were investigated for six 22-base pair, non-self-complementary duplex DNA oligomers. The center sequences of these deoxyoligonucleotides have 8-14 base pairs in common and are flanked on both sides by sequences differing in context and A-T content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming a two-state melting transition. Melting free energies (20 degrees C) of the six DNAs evaluated from DSC experiments ranged from -18.7 to -32.7 kcal/mol. van't Hoff estimates of the free energies ranged from -18.5 to -48.0 kcal/mol. With either method, the trends in free energy as a function of sequence were identical. Equilibrium binding by BamHI restriction endonuclease to the 22-base pair DNAs was also investigated. The central eight base pairs of all six molecules, 5'-A-GGATCC-A-3', contained a BamHI recognition sequence bounded by A-T base pairs. Magnesium free binding assays were performed by titering BamHI against a constant concentration of each of the deoxyoligonucleotide substrates and analyzing reaction products by gel retardation. Binding isotherms of the total amount of bound DNA versus protein concentration were constructed which provided semiquantitative estimates of the equilibrium dissociation constants for dissociation of BamHI from the six DNA oligomers. Dissociation constants ranged from 0.5 x 10(-)(9) to 12.0 x 10(-)(9) M with corresponding binding free energies of -12.5 to -10.6 (+/-0. 1) kcal/mol. An inverse relationship is found when binding and stability are compared.

Melting studies of short DNA hairpins containing the universal base 5-nitroindole.

Effects of the universal base 5-nitroindole on the thermodynamic stability of DNA hairpins having a 6 bp stem and four base loops were investigated by optical absorbance and differential scanning calorimetry techniques. Melting studies were conducted in buffer containing 115 mM Na(+). Five different modified versions of DNA hairpins containing a 5-nitroindole base or bases substituted at different positions in the stem and loop regions were examined. Thermo-dynamic parameters of the melting transitions estimated from a two-state analysis of optical melting curves and measured directly by calorimetry revealed that the presence of 5-nitroindole bases in the duplex stem or loop regions of short DNA hairpins significantly affects both their enthalpic and entropic melting components in a compensating manner, while the transition free energy varies linearly with the transition temperature. The calorimetrically determined enthalpy and entropy values of the modified hairpins were considerably smaller (43-53%) than the two-state optical parameters, suggesting that solvent effects may be significant in the melting processes of these hairpins. Results of circular dichroism measurements also revealed slight differences between the modified hairpins and the control in both the duplex and melted states, suggesting subtle structural differences between the control and DNA hairpins containing a 5-nitroindole base or bases.

Melting studies of short DNA hairpins: influence of loop sequence and adjoining base pair identity on hairpin thermodynamic stability.

Spectroscopic and calorimetric melting studies of 28 DNA hairpins were performed. These hairpins form by intramolecular folding of 16 base self-complementary DNA oligomer sequences. Sequence design dictated that the hairpin structures have a six base pair duplex linked by a four base loop and that the first five base pairs in the stem are the same in every molecule. Only loop sequence and identity of the duplex base pair closing the loop vary for the set of hairpins. For these DNA samples, melting studies were carried out to investigate effects of the variables on hairpin stability. Stability of the 28 oligomers was ascertained from their temperature-induced melting transitions in buffered 115 mM Na(+) solvent, monitored by ultraviolet absorbance and differential scanning calorimetry (DSC). Experiments revealed the melting temperatures of these molecules range from 32.4 to 60.5 degrees C and are concentration independent over strand concentrations of 0.5 to 260 microM; thus, as expected for hairpins, the melting transitions are apparently unimolecular. Model independent thermodynamic transition parameters, DeltaH(cal), DeltaS(cal), and DeltaG(cal), were determined from DSC measurements. Model dependent transition parameters, DeltaH(vH), DeltaS(vH), and DeltaG(vH) were estimated from a van't Hoff (two-state) analysis of optical melting transitions. Results of these studies reveal a significant sequence dependence to DNA hairpin stability. Thermodynamic parameters evaluated by either procedure reveal the transition enthalpy, DeltaH(cal) (DeltaH(vH)) can differ by as much as 20 kcal/mol depending on sequence. Similarly, values of the transition entropy DeltaS(cal) (DeltaS(vH)) can differ by as much as 60 cal/Kmol (eu) for different molecules. Differences in free energies DeltaG(cal) (DeltaG(vH)) are as large as 4 kcal/mol for hairpins with different sequences. Comparisons between the model independent calorimetric values and the thermodynamic parameters evaluated assuming a two-state model reveal that 10 of the 28 hairpins display non-two-state melting behavior. The database of sequence-dependent melting free energies obtained for the hairpins was employed to extract a set of n-n (nearest-neighbor) sequence dependent loop parameters that were able to reproduce the input data within error (with only two exceptions). Surprisingly, this suggests that the thermodynamic stability of the DNA hairpins can in large part be reasonably represented in terms of sums of appropriate nearest-neighbor loop sequence parameters.

Studies of DNA dumbbells VII: evaluation of the next-nearest-neighbor sequence-dependent interactions in duplex DNA.

Melting experiments were conducted on 22 DNA dumbbells as a function of solvent ionic strength from 25-115 mM Na(+). The dumbbell molecules have short duplex regions comprised of 16-20 base pairs linked on both ends by T(4) single-strand loops. Only the 4-8 central base pairs of the dumbbell stems differ for different molecules, and the six base pairs on both sides of the central sequence and adjoining loops on both ends are the same in every molecule. Results of melting analysis on the 22 new DNA dumbbells are combined with our previous results on 17 other DNA dumbbells, with stem lengths containing from 14-18 base pairs, reported in the first article of this series (Doktycz, Goldstein, Paner, Gallo, and Benight, Biopoly 32, 1992, 849-864). The combination of results comprises a database of optical melting parameters for 39 DNA dumbbells in ionic strengths from 25-115 mM Na(+). This database is employed to evaluate the thermodynamics of singlet, doublet, and triplet sequence-dependent interactions in duplex DNA. Analysis of the 25 mM Na(+) data reveals the existence of significant sequence-dependent triplet or next-nearest-neighbor interactions. The enthalpy of these interactions is evaluated for all possible triplets. Some of the triplet enthalpy values are less than the uncertainty in their evaluation, indicating no measurable interaction for that particular sequence. This finding suggests that the thermodynamic stability of duplex DNA depends on solvent ionic strength in a sequence-dependent manner. As a part of the analysis, the nearest-neighbor (base pair doublet) interactions in 55, 85, and 115 mM Na(+) are also reevaluated from the larger database.

Predicting sequence-dependent melting stability of short duplex DNA oligomers.

Many important applications of DNA sequence-dependent hybridization reactions have recently emerged. This has sparked a renewed interest in analytical calculations of sequence-dependent melting stability of duplex DNA. In particular, for many applications it is often desirable to accurately predict the transition temperature, or tm of short duplex DNA oligomers (approximately 20 base pairs or less) from their sequence and concentration. The thermodynamic analytical method underlying these predictive calculations is based on the nearest-neighbor model. At least 11 sets of nearest-neighbor sequence-dependent thermodynamic parameters for DNA have been published. These sets are compared. Use of the nearest-neighbor sets in predicting tm from the DNA sequence is demonstrated, and the ability of the nearest-neighbor parameters to provide accurate predictions of experimental tm's of short duplex DNA oligomers is assessed.