Search

We describe the analysis of errors and failure modes in the base-calling function in automated DNA sequencing, on instruments in which fluorescently-labeled Sanger dideoxy-sequencing ladders are detected via their times of migration past a fixed detector. A general approach entails the joint use of: (i) well-defined control samples such as M13mp18, and (ii) mathematical simulation of sequencing electropherograms, with the deliberate introduction of different types of distortion and noise. An algorithm, the electrophoretic trace simulator (ETS), is used to calculate electrophoresis traces corresponding to the output data stream of an automated fluorescent DNA sequencer. The ETS accepts a user-defined sequence of nucleotide bases (A, C, G, T) as input, and employs user-adjustable functions to compute the following critical parameters of an electropherogram: peak intensity, peak spacing, peak shape as a function of base number; background, noise, and spectral cross-talk correction (for a sequencer using multiple dyes). We use a combination of M13mp18 controls and simulated electropherograms to analyze two problems of considerable practical importance: (i) variation in electrophoretic migration rates between different lanes of a gel, and (ii) variation in signal intensity due to user-dependent loading artifacts. The issue of base-calling errors and failure modes, for electropherograms that contain noise and distortion, is addressed.

We present a new method for the linearization and alignment of data traces generated by multilane automated DNA sequencing instruments. Application of this method to data generated with the Visible Genetics Open Gene DNA sequencing system (using MicroCel 700 gel cassettes, with a 25 cm separation distance) allows read lengths of > 1,000 nucleotides to be routinely obtained with high confidence and > 97% accuracy. This represents an increase of 10-15% in average read length, relative to data from this system that have not been processed in the fashion described herein. Most importantly, the linearization and alignment method allows usable sequence to be obtained from a fraction of 10-15% of data sets which, because of original trace misalignment problems, would otherwise have to be discarded. Our method involves adding electrophoretic calibration standards to the DNA sequencing fragments. The calibration standards are labeled with a dye that differs spectrally from the dye attached to the sequencing fragments. The calibration standards are identical in all the lanes. Analysis of the mobilities of the calibration standards allows correction for both systematic and random variation of electrophoretic properties between gel lanes. We have successfully used this method with two-dye and three-dye DNA sequencing instruments.

The Visible Genetics Clipper sequencer is a new platform for automated DNA sequencing which employs disposable MicroCel cassettes and 50 microm thick polyacrylamide gels. Two DNA ladders can be analyzed simultaneously in each of 16 lanes on a gel, after labeling with far-red absorbing dyes such as Cy5 and Cy5.5. This allows a simultaneous bidirectional sequencing of four templates. We have evaluated the Clipper sequencer, by cycle-sequencing of an M13 single-stranded DNA standard, and by coupled amplification and sequencing (CLIP) of reverse-transcribed human immunodeficiency virus (HIV-1) RNA standards and clinical patient samples. (i) Limitations of instrument. We have examined basic instrument parameters such as detector stability, background, digital sampling rate, and gain. With proper usage, the optical and electronic subsystems of the Clipper sequencer do not limit the data collection or sequence-determination processes. (ii) Limitations of gel performance. We have also examined the physics of DNA band separation on 50 microm thick MicroCel gels. We routinely obtain well-resolved sequence which can be base-called with 98.5% accuracy to position approximately 450 on an 11 cm gel, and to position approximately 900 on a 25 cm gel. Resolution on 5 and 11 cm gels ultimately is limited by a sharp decrease in spacing between adjacent bands, in the biased reptation separation regime. Fick's (thermal) diffusion appears to be of minor importance on 6 cm or 11 cm gels, but becomes an additional resolution-limiting factor on 25 cm gels. (iii) Limitations of enzymology. Template quality, primer nesting, choice of DNA polymerase, and choice between dye primers and dye terminators are key determinants of the ability to detect mutations and polymorphisms on the Clipper sequencer, as on other DNA sequencers. When CLIP is used with dye-labeled primers and a DNA polymerase of the F667Y, delta(5'--> 3' exo) class, we can routinely detect single-nucleotide mutations and polymorphisms over the 0.35-0.65 heterozygosity range. We present an example of detecting therapeutically relevant mutations in a clinical HIV-1 RNA isolate.