Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

Automated parallel DNA sequencing on multiple channel microchips

We report automated DNA sequencing in 16-channel microchips. A microchip prefilled with sieving matrix is aligned on a heating plate affixed to a movable platform. Samples are loaded into sample reservoirs by using an eight-tip pipetting device, and the chip is docked with an array of electrodes in the focal plane of a four-color scanning detection system. Under computer control, high voltage is applied to the appropriate reservoirs in a programmed sequence that injects and separates the DNA samples. An integrated four-color confocal fluorescent detector automatically scans all 16 channels. The system routinely yields more than 450 bases in 15 min in all 16 channels. In the best case using an automated base-calling program, 543 bases have been called at an accuracy of >99%. Separations, including automated chip loading and sample injection, normally are completed in less than 18 min. The advantages of DNA sequencing on capillary electrophoresis chips include uniform signal intensity and tolerance of high DNA template concentration. To understand the fundamentals of these unique features we developed a theoretical treatment of cross-channel chip injection that we call the differential concentration effect. We present experimental evidence consistent with the predictions of the theory.

DNA sequencing and genotyping by transcriptional synthesis of chain-terminated RNA ladders and MALDI-TOF mass spectrometry

Sets of RNA ladders can be synthesized by transcription of a bacteriophage-encoded RNA polymerase using 3′-deoxynucleotides as chain terminators. These ladders can be used for sequencing of DNA. Using a nicked form of phage SP6 RNA polymerase in this study substantially enhanced yields of transcriptional sequencing ladders. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) of chain-terminated RNA ladders allowed DNA sequence determination of up to 56 nt. It is also demonstrated that A→G and C→T variations in heterozygous and homozygous samples can be unambiguously identified by the mass spectrometric analysis. As a step towards single-tube sequencing reactions, α-thiotriphosphate nucleotide analogs were used to overcome problems caused by chain terminator-independent, premature termination and by the small mass difference between natural pyrimidine nucleotides.

DNA sequencing by delayed extraction-matrix-assisted laser desorption/ionization time of flight mass spectrometry.

Matrix-assisted laser desorption/ionization (MALDI) time of flight mass spectrometry was used to detect and order DNA fragments generated by Sanger dideoxy cycle sequencing. This was accomplished by improving the sensitivity and resolution of the MALDI method using a delayed ion extraction technique (DE-MALDI). The cycle sequencing chemistry was optimized to produce as much as 100 fmol of each specific dideoxy terminated fragment, generated from extension of a 13-base primer annealed on 40- and 50-base templates. Analysis of the resultant sequencing mixture by DE-MALDI identified the appropriate termination products. The technique provides a new non-gel-based method to sequence DNA which may ultimately have considerable speed advantages over traditional methodologies.

A DNA sequencing strategy that requires only five bases of known terminal sequence for priming.

We have previously reported an enhanced version of sequencing by hybridization (SBH), termed positional SBH (PSBH). PSBH uses partially duplex probes containing single-stranded 3' overhangs, instead of simple single-stranded probes. Stacking interactions between the duplex probe and a single-stranded target allow us to reduce the probe sizes required to 5-base single-stranded overhangs. Here we demonstrate the use of PSBH to capture relatively long single-stranded DNA targets and perform standard solid-state Sanger sequencing on these primer-template complexes without ligation. Our results indicate that only 5 bases of known terminal sequence are required for priming. In addition, the partially duplex probes have the ability to capture their specific target from a mixture of five single-stranded targets with different 3'-terminal sequences. This indicates the potential utility of the PSBH approach to sequence mixtures of DNA targets without prior purification.

Fluorescence energy transfer dye-labeled primers for DNA sequencing and analysis.

Fluorescent dye-labeled DNA primers have been developed that exploit fluorescence energy transfer (ET) to optimize the absorption and emission properties of the label. These primers carry a fluorescein derivative at the 5' end as a common donor and other fluorescein and rhodamine derivatives attached to a modified thymidine residue within the primer sequence as acceptors. Adjustment of the donor-acceptor spacing through the placement of the modified thymidine in the primer sequence allowed generation of four primers, all having strong absorption at a common excitation wavelength (488 nm) and fluorescence emission maxima of 525, 555, 580, and 605 nm. The ET efficiency of these primers ranges from 65% to 97%, and they exhibit similar electrophoretic mobilities by gel electrophoresis. With argon-ion laser excitation, the fluorescence of the ET primers and of the DNA sequencing fragments generated with ET primers is 2- to 6-fold greater than that of the corresponding primers or fragments labeled with single dyes. The higher fluorescence intensity of the ET primers allows DNA sequencing with one-fourth of the DNA template typically required when using T7 DNA polymerase. With single-stranded M13mp18 DNA as the template, a typical sequencing reaction with ET primers on a commercial sequencer provided DNA sequences with 99.8% accuracy in the first 500 bases. ET primers should be generally useful in the development of other multiplex DNA sequencing and analysis methods.

Direct sequencing of bacterial and P1 artificial chromosome-nested deletions for identifying position-specific single-nucleotide polymorphisms

A loxP-transposon retrofitting strategy for generating large nested deletions from one end of the insert DNA in bacterial artificial chromosomes and P1 artificial chromosomes was described recently [Chatterjee, P. K. & Coren, J. S. (1997) Nucleic Acids Res. 25, 2205–2212]. In this report, we combine this procedure with direct sequencing of nested-deletion templates by using primers located in the transposon end to illustrate its value for position-specific single-nucleotide polymorphism (SNP) discovery from chosen regions of large insert clones. A simple ampicillin sensitivity screen was developed to facilitate identification and recovery of deletion clones free of transduced transposon plasmid. This directed approach requires minimal DNA sequencing, and no in vitro subclone library generation; positionally oriented SNPs are a consequence of the method. The procedure is used to discover new SNPs as well as physically map those identified from random subcloned libraries or sequence databases. The deletion templates, positioned SNPs, and markers are also used to orient large insert clones into a contig. The deletion clone can serve as a ready resource for future functional genomic studies because each carries a mammalian cell-specific antibiotic resistance gene from the transposon. Furthermore, the technique should be especially applicable to the analysis of genomes for which a full genome sequence or radiation hybrid cell lines are unavailable.

Fluorescence-based directed termination PCR: direct mutation characterization without sequencing

We describe a fluorescence-based directed termination PCR (fluorescent DT–PCR) that allows accurate determination of actual sequence changes without dideoxy DNA sequencing. This is achieved using near infrared dye-labeled primers and performing two PCR reactions under low and unbalanced dNTP concentrations. Visualization of resulting termination fragments is accomplished with a dual dye Li-cor DNA sequencer. As each DT–PCR reaction generates two sets of terminating fragments, a pair of complementary reactions with limiting dATP and dCTP collectively provide information on the entire sequence of a target DNA, allowing an accurate determination of any base change. Blind analysis of 78 mutants of the supF reporter gene using fluorescent DT–PCR not only correctly determined the nature and position of all types of substitution mutations in the supF gene, but also allowed rapid scanning of the signature sequences among identical mutations. The method provides simplicity in the generation of terminating fragments and 100% accuracy in mutation characterization. Fluorescent DT–PCR was successfully used to generate a UV-induced spectrum of mutations in the supF gene following replication on a single plate of human DNA repair-deficient cells. We anticipate that the automated DT–PCR method will serve as a cost-effective alternative to dideoxy sequencing in studies involving large-scale analysis for nucleotide sequence changes.

Cloning and expression in Escherichia coli of the OGG1 gene of Saccharomyces cerevisiae, which codes for a DNA glycosylase that excises 7,8-dihydro-8-oxoguanine and 2,6-diamino-4-hydroxy-5-N-methylformamidopyrimidine.

A spontaneous mutator strain of Escherichia coli (fpg mutY) was used to clone the OGG1 gene of Saccharomyces cerevisiae, which encodes a DNA glycosylase activity that excises 7,8-dihydro-8-oxoguanine (8-OxoG). E. coli (fpg mutY) was transformed by a yeast DNA library, and clones that showed a reduced spontaneous mutagenesis were selected. The antimutator activity was associated with pYSB10, an 11-kbp recombinant plasmid. Cell-free extracts of E. coli (fpg mutY) harboring pYSB10 possess an enzymatic activity that cleaves a 34-mer oligonucleotide containing a single 8-oxoG opposite a cytosine (8-OxoG/C). The yeast DNA fragment of 1.7 kbp that suppresses spontaneous mutagenesis and overproduces the 8-OxoG/C cleavage activity was sequenced and mapped to chromosome XIII. DNA sequencing identified an open reading frame, designated OGG1, which encodes a protein of 376 amino acids with a molecular mass of 43 kDa. The OGG1 gene was inserted in plasmid pUC19, yielding pYSB110. E. coli (fpg) harboring pYSB110 was used to purify the Ogg1 protein of S. cerevisiae to apparent homogeneity. The Ogg1 protein possesses a DNA glycosylase activity that releases 8-OxoG and 2,6-diamino-4-hydroxy-5-N-methylformamidopyrimidine. The Ogg1 protein preferentially incises DNA that contains 8-OxoG opposite cytosine (8-OxoG/C) or thymine (8-OxoG/T). In contrast, Ogg1 protein does not incise the duplex where an adenine is placed opposite 8-OxoG (8-OxoG/A). The mechanism of strand cleavage by Ogg1 protein is probably due to the excision of 8-OxoG followed by a beta-elimination at the resulting apurinic/apyrimidinic site.

Transposon tools for recombinant DNA manipulation: characterization of transcriptional regulators from yeast, Xenopus, and mouse.

Transposon Tn1000 has been adapted to deliver novel DNA sequences for manipulating recombinant DNA. The transposition procedure for these "tagged" Tn1000s is simple and applicable to most plasmids in current use. For yeast molecular biology, tagged Tn1000s introduce a variety of yeast selective markers and replication origins into plasmids and cosmids. In addition, the beta-globin minimal promoter and lacZ gene of Tn(beta)lac serve as a mobile reporter of eukaryotic enhancer activity. In this paper, Tn(beta)lac was used to localize a mouse HoxB-complex enhancer in transgenic mice. Other tagged transposons create Gal4 DNA-binding-domain fusions, in either Escherichia coli or yeast plasmids, for use in one- and two-hybrid tests of transcriptional activation and protein-protein interaction, respectively. With such fusions, the Saccharomyces cerevisiae Swi6 G1/S-phase transcription factor and the Xenopus laevis Pintallavis developmental regulator are shown to activate transcription. Furthermore, the same transposon insertions also facilitated mapping of the Swi6 and Pintallavis domains responsible for transcriptional activation. Thus, as well as introducing novel sequences, tagged transposons share the numerous other applications of transposition such as producing insertional mutations, creating deletion series, or serving as mobile primer sites for DNA sequencing.

Identification of a DNA transformation gene required for com101A+ expression and supertransformer phenotype in Haemophilus influenzae.

DNA sequencing, RNA mapping, and protein expression experiments revealed the presence of a gene, tfoX+, encoding a 24.9-kDa polypeptide, that is transcribed divergently from a common promoter region with the Haemophilus influenzae rec-1+ gene. H. influenzae strains mutant for tfoX failed to bind transforming DNA and were transformation deficient. Primer extension experiments utilizing in vivo total RNA from precompetent and competent H. influenzae cells demonstrated that transcription of tfoX+ increased immediately upon competence induction, suggesting that tfoX+ is an early competence gene. Similar experiments showed that the expression of the late competence-specific gene, com101A+, was tfoX+ dependent. Moreover, expression of plasmid-borne tfoX+ in H. influenzae resulted in constitutive competence. The addition of cyclic adenosine monophosphate (cAMP) to strains carrying a tfoX::lacZ operon fusion resulted in an immediate increase in beta-galactosidase activity that correlated with an increase in genetic transformability. Collectively, our results suggest that TfoX may play a key role in the development of genetic competence by regulating the expression of late competence-specific genes.

A mutation in the transmembrane/luminal domain of the ryanodine receptor is associated with abnormal Ca2+ release channel function and severe central core disease

Central core disease is a rare, nonprogressive myopathy that is characterized by hypotonia and proximal muscle weakness. In a large Mexican kindred with an unusually severe and highly penetrant form of the disorder, DNA sequencing identified an I4898T mutation in the C-terminal transmembrane/luminal region of the RyR1 protein that constitutes the skeletal muscle ryanodine receptor. All previously reported RYR1 mutations are located either in the cytoplasmic N terminus or in a central cytoplasmic region of the 5,038-aa protein. The I4898T mutation was introduced into a rabbit RYR1 cDNA and expressed in HEK-293 cells. The response of the mutant RyR1 Ca2+ channel to the agonists halothane and caffeine in a Ca2+ photometry assay was completely abolished. Coexpression of normal and mutant RYR1 cDNAs in a 1:1 ratio, however, produced RyR1 channels with normal halothane and caffeine sensitivities, but maximal levels of Ca2+ release were reduced by 67%. [3H]Ryanodine binding indicated that the heterozygous channel is activated by Ca2+ concentrations 4-fold lower than normal. Single-cell analysis of cotransfected cells showed a significantly increased resting cytoplasmic Ca2+ level and a significantly reduced luminal Ca2+ level. These data are indicative of a leaky channel, possibly caused by a reduction in the Ca2+ concentration required for channel activation. Comparison with two other coexpressed mutant/normal channels suggests that the I4898T mutation produces one of the most abnormal RyR1 channels yet investigated, and this level of abnormality is reflected in the severe and penetrant phenotype of affected central core disease individuals.

Lethal paralysis of Caenorhabditis elegans by Pseudomonas aeruginosa

Identification of host factors that interact with pathogens is crucial to an understanding of infectious disease, but direct screening for host mutations to aid in this task is not feasible in mammals. The nematode Caenorhabditis elegans is a genetically tractable alternative for investigating the pathogenic bacterium Pseudomonas aeruginosa. A P. aeruginosa toxin, produced at high cell density under control of the quorum-sensing regulators LasR and RhlR, rapidly and lethally paralyzes C. elegans. Loss-of-function mutations in C. elegans egl-9, a gene required for normal egg laying, confer strong resistance to the paralysis. Thus, activation of EGL-9 or of a pathway that includes it may lead to the paralysis. The molecular identity of egl-9 was determined by transformation rescue and DNA sequencing. A mammalian homologue of EGL-9 is expressed in tissues in which exposure to P. aeruginosa could have clinical effects.

Genetic Separation of FK506 Susceptibility and Drug Transport in the Yeast Pdr5 ATP-binding Cassette Multidrug Resistance Transporter

Overexpression of the yeast Pdr5 ATP-binding cassette transporter leads to pleiotropic drug resistance to a variety of structurally unrelated cytotoxic compounds. To identify Pdr5 residues involved in substrate recognition and/or drug transport, we used a combination of random in vitro mutagenesis and phenotypic screening to isolate novel mutant Pdr5 transporters with altered substrate specificity. A plasmid library containing randomly mutagenized PDR5 genes was transformed into appropriate drug-sensitive yeast cells followed by phenotypic selection of Pdr5 mutants. Selected mutant Pdr5 transporters were analyzed with respect to their expression levels, subcellular localization, drug resistance profiles to cycloheximide, rhodamines, antifungal azoles, steroids, and sensitivity to the inhibitor FK506. DNA sequencing of six PDR5 mutant genes identified amino acids important for substrate recognition, drug transport, and specific inhibition of the Pdr5 transporter. Mutations were found in each nucleotide-binding domain, the transmembrane domain 10, and, most surprisingly, even in predicted extracellular hydrophilic loops. At least some point mutations identified appear to influence folding of Pdr5, suggesting that the folded structure is a major substrate specificity determinant. Surprisingly, a S1360F exchange in transmembrane domain 10 not only caused limited substrate specificity, but also abolished Pdr5 susceptibility to inhibition by the immunosuppressant FK506. This is the first report of a mutation in a yeast ATP-binding cassette transporter that allows for the functional separation of substrate transport and inhibitor susceptibility.

Subgenome chromosome walking in wheat: A 450-kb physical contig in Triticum monococcum L. spans the Lr10 resistance locus in hexaploid wheat (Triticum aestivum L.)

For many agronomically important plant genes, only their position on a genetic map is known. In the absence of an efficient transposon tagging system, such genes have to be isolated by map-based cloning. In bread wheat Triticum aestivum, the genome is hexaploid, has a size of 1.6 × 1010 bp, and contains more than 80% of repetitive sequences. So far, this genome complexity has not allowed chromosome walking and positional cloning. Here, we demonstrate that chromosome walking using bacterial artificial chromosome (BAC) clones is possible in the diploid wheat Triticum monococcum (Am genome). BAC end sequences were mostly repetitive and could not be used for the first walking step. New probes corresponding to rare low-copy sequences were efficiently identified by low-pass DNA sequencing of the BACs. Two walking steps resulted in a physical contig of 450 kb on chromosome 1AmS. Genetic mapping of the probes derived from the BAC contig demonstrated perfect colinearity between the physical map of T. monococcum and the genetic map of bread wheat on chromosome 1AS. The contig genetically spans the Lr10 leaf rust disease resistance locus in bread wheat, with 0.13 centimorgans corresponding to 300 kb between the closest flanking markers. Comparison of the genetic to physical distances has shown large variations within 350 kb of the contig. The physical contig can now be used for the isolation of the orthologous regions in bread wheat. Thus, subgenome chromosome walking in wheat can produce large physical contigs and saturate genomic regions to support positional cloning.