Mapping using gene encyclopaedias

Identification of phosphate-regulated genes by differential expression in the UV-irradiated host system

The UV-irradiated host system has been used for identifying protein products of genes cloned in a phage vector. By starving the host cells for phosphate immediately before UV-irradiation, we demonstrate that phosphate-regulated genes can be easily identified. By employing this new technique, we also provide evidence showing that the gpsA gene might be a new member of the phosphate starvation-inducible (psi) genes of E. coli.

Genome encyclopedias and their use for comparative analysis of Rhodobacter capsulatus strains

This paper consists of two components: the use of gene encyclopedias in genomic studies and Rhodobacter capsulatus genome project. A survey of vectors used for encyclopedia construction includes a brief discussion of their relative advantages and limitations. Projects employing various methods of encyclopedia assembly including the comparison of restriction patterns, restriction maps, linking by hybridization, oligonucleotide fingerprinting, sequence tagged site (STS) fingerprinting and encyclopedias derived from genetic maps are listed and briefly described. The R. capsulatus SB 1003 genome project started with the construction of its cosmid encyclopedia, which comprises 192 cosmids covering the chromosome and the 134 kbp plasmid in strain SB 1003, with the exact map coordinates of each cosmid. In a pilot sequencing study, several cosmids were individually subcloned using the vector M13mp18 and merged into one 189 kbp contig. About 160 open reading frames (ORFs) identified by the CodonUse program were subjected to similarity searches. The biological functions of eighty ORFs could be assigned reliably using the WIT (what is there) genome investigation environment. Eighty percent of these recognizable ORFs were organized in functional clusters, which simplified assignment decisions and increased the strength of the predictions. A set of 26 genes for cobalamin biosynthesis, genes for polyhydroxyalkanoic acid metabolism, DNA replication and recombination, and DNA gyrase were among those identified. Recently, another 1.2 Mbp genome fragment of the Rhodobacter genome was sequenced using a slightly modified approach. These results together with some genome investigation tools, have been placed at our web site (http://capsulapedia.uchicago.edu). The sequence of R. capsulatus is expected to be completed by summer 1998. A project to construct a systematic set of deletion strains of R. capsulatus in order to assign functions to unknown ORFs has been started. Preliminary data demonstrate the extreme convenience of the unique gene transfer agent (GTA) system to perform such work.

Design and analysis of chromosome physical mapping experiments

Mathematical and statistical aspects of constructing ordered-clone physical maps of chromosomes are reviewed. Three broad problems are addressed: analysis of fingerprint data to identify configurations of overlapping clones, prediction of the rate of progress of a mapping strategy and optimal design of pooling schemes for screening large clone libraries.

Bacterial genomics

During the last decade, great advances have been made in the study of bacterial genomes which is perhaps better described by the term bacterial genomics. The application of powerful techniques, such as pulsed-field gel electrophoresis of macro-restriction fragments of genomic DNA, has freed the characterisation of the chromosomes of many bacteria from the constraints imposed by classical genetic analysis. It is now possible to analyse the genome of virtually every microorganism by direct molecular methods and to construct detailed physical and gene maps. In this review, the various practical approaches are compared and contrasted, and some of the emerging themes of bacterial genomics, such as the size, shape, number and organisation of chromosomes are discussed.

A large-conductance mechanosensitive channel in E. coli encoded by mscL alone

All cellular organisms respond to vibration, touch, gravity or changes in osmolarity, although the molecules on which such mechanosensations depend are unknown. Candidates include certain channels that gate in response to membrane stretch. Patch-clamp experiments with Escherichia coli envelope have revealed a mechanosensitive channel with very large conductance (MscL) and one with a smaller conductance (MscS) which may be important in osmoregulation. Here we have solubilized and fractionated the envelope, reconstituted the MscL activity in vitro, and traced it to a small protein, whose gene, mscL, we then cloned. Insertional disruption of mscL removes the channel activity, whereas re-expression of mscL borne on an expression plasmid restores it. MscL-channel activities were observed in material from a cell-free expression system with mscL as the only template. The mscL nucleotide sequence predicts a unique protein of only 136 amino acids, with a highly hydrophobic core and very different from porins or other known proteins.

Sequence analysis of four new heat-shock genes constituting the hslTS/ibpAB and hslVU operons in Escherichia coli

Sequences of four new heat-shock (HS) genes of Escherichia coli organized into two operons were determined. The operon at 83 min specifies two proteins of 15.8 kDa (HslT) and 16.1 kDa (HslS), which are identical to IbpA and IbpB, respectively. Expression of mRNA from a sigma 32-dependent promoter of the hslTS/ibpAB operon is stimulated 30-75-fold upon temperature upshift. The transcription start point (tsp) is located at a G, 96 bp upstream from the AUG start codon of hslT/ibpA. The deduced amino acid sequences of HslT/IbpA and HslS/IbpB are 48% identical to each other and were found to be remotely related to the chloroplast low-molecular-weight HS protein, which is highly conserved among plants. The second hs operon is much less actively stimulated by temperature upshift, although it has a hs promoter that perfectly matches the consensus of promoters recognized by sigma 32. Located at 88.9 min, the hslVU operon specifies proteins of 19.1 kDa (HslV) and 49.6 kDa (HslU). Multiple tsp were found in this operon. HslV is remotely related to the eukaryotic proteasome proteins, and HslU is very similar to a Pasteurella haemolytica protein of unknown function. Both HslU and the P. haemolytica protein share a ATP/GTP-binding motif near their N-termini. The two operons described here are transcribed counterclockwise on the standard genetic map.

Global regulation of gene expression in Escherichia coli

Global transcription responses of Escherichia coli to various stimuli or genetic defects were studied by measuring mRNA levels in about 400 segments of the genome. Measuring mRNA levels was done by analyzing hybridization to DNA dot blots made with overlapping lambda clones spanning the genome of E. coli K-12. Conditions examined included isopropyl-beta-D-thiogalactopyranoside (IPTG) induction, heat shock, osmotic shock, starvation for various nutrients, entrance of cells into the stationary phase of growth, anaerobic growth in a tube, growth in the gnotobiotic mouse gut, and effects of pleiotropic mutations rpoH, himA, topA, and crp. Most mapped genes known to be regulated by a particular situation were successfully detected. In addition, many chromosomal regions containing no previously known regulated genes were discovered that responded to various stimuli. This new method for studying globally regulated genetic systems in E. coli combines detection, cloning, and physical mapping of a battery of coregulated genes in one step.

Analysis of the Escherichia coli genome: DNA sequence of the region from 84.5 to 86.5 minutes

The DNA sequence of 91.4 kilobases of the Escherichia coli K-12 genome, spanning the region between rrnC at 84.5 minutes and rrnA at 86.5 minutes on the genetic map (85 to 87 percent on the physical map), is described. Analysis of this sequence identified 82 potential coding regions (open reading frames) covering 84 percent of the sequenced interval. The arrangement of these open reading frames, together with the consensus promoter sequences and terminator-like sequences found by computer searches, made it possible to assign them to proposed transcriptional units. More than half the open reading frames correlated with known genes or functions suggested by similarity to other sequences. Those remaining encode still unidentified proteins. The sequenced region also contains several RNA genes and two types of repeated sequence elements were found. Intergenic regions include three "gray holes," 0.6 to 0.8 kilobases, with no recognizable functions.

The application of Markov chain analysis to oligonucleotide frequency prediction and physical mapping of Drosophila melanogaster

Here we compare several methods for predicting oligonucleotide frequencies in 691 kb of Drosophila melanogaster DNA. As in previous work on Escherichia coli and Saccharomyces cerevisiae, a relatively simple equation based on tetranucleotide frequencies can be used in predicting frequencies of higher order oligonucleotides. For example, the mean of observed/expected abundances of 4,096 hexamers was 1.07 with a sample standard deviation of .55. This simple predictor arises by considering each base on the sense strand of D. melanogaster to depend only on the three bases 5' to it (a 3rd order Markov chain) and is more accurate than the random predictor. This equation is useful in predicting restriction enzyme fragment sizes, selecting restriction enzymes that cut preferentially in coding vs noncoding regions, and in selecting probes to fingerprint clones in contig mapping. Once again, this equation well predicts the occurrence of higher order oligonucleotides, supporting our hypothesis that this predictor holds in evolutionarily diverse organisms. When ranked from highest to lowest abundance, the observed frequencies of oligomers of a given length are closely tracked by the predicted abundances of a 3rd order Markov chain. Through use of the dependence of oligomer frequencies on base composition, we report a list of oligomers that will be useful for the completion of a cosmid physical map of D. melanogaster. Presently, the library is such that it will be possible to construct large contigs using only 30 oligonucleotide probes to fingerprint cosmids.

SfiI genomic cleavage map of Escherichia coli K-12 strain MG1655

An SfiI restriction map of Escherichia coli K-12 strain MG1655 is presented. The map contains thirty-one cleavage sites separating fragments ranging in size from 407 kb to 3.7 kb. Several techniques were used in the construction of this map, including CHEF pulsed field gel electrophoresis; physical analysis of a set of twenty-six auxotrophic transposon insertions; correlation with the restriction map of Kohara and coworkers using the commercially available E. coli Gene Mapping Membranes; analysis of publicly available sequence information; and correlation of the above data with the combined genetic and physical map developed by Rudd, et al. The combination of these techniques has yielded a map in which all but one site can be localized within a range of +/- 2 kb, and over half the sites can be localized precisely by sequence data. Two sites present in the EcoSeq5 sequence database are not cleaved in MG1655 and four sites are noted to be sensitive to methylation by the dcm methylase. This map, combined with the NotI physical map of MG1655, can aid in the rapid, precise mapping of several different types of genetic alterations, including transposon mediated mutations and other insertions, inversions, deletions and duplications.

Distribution of moderately repetitive sequences pTR5 and LF1 in Xq24-q28 human DNA and their use in assembling YAC contigs

Xq24-q28 DNA, from a hamster/human hybrid cell containing only that portion of the human X chromosome, was found to contain 56 TaqI restriction fragments that hybridized to the moderately repetitive sequence pTR5. Using the pTR5 sequence as a probe in colony hybridization, 136 cognate yeast artificial chromosome (YAC) clones were detected among a collection of 820 containing about three genomic equivalents of the Xq24-q28 DNA. The YACs were then grouped into 48 contigs and single clones containing one or more of the TaqI fragments. Overlaps were confirmed both by fingerprinting YACs with AluI and L1 probes and by additional information. A less complete analysis was also carried out with a second moderately repetitive sequence, LF1, and some smaller contigs were merged into larger ones. Moderately repetitive sequences can thus be used as probes for multiple loci in single hybridization experiments and can help to organize and confirm YAC overlaps during the development of maps with long-range contiguity.

Statistical analysis of DNA fingerprint data for ordered clone physical mapping of human chromosomes

A statistical framework is proposed for analysing DNA fingerprint data from experiments aimed at constructing ordered clone physical maps of chromosomes. The fingerprint data consists of the lengths and hybridization states of restriction digest fragments and the paper develops a solution to the fundamental problem of deciding whether or not two randomly selected clones overlap. Overlap probabilities are calculated using Bayes' rule together with appropriate statistical descriptions of the chromosome and experimental procedure. The analysis is flexible, allowing a variety of assumptions to account for experimental errors and difficulties, such as unobserved fragments. The approach described here provides a basis for predicting the rate of progress of an experimental protocol and hence for comparing alternate protocols. It is readily generalized to related problems with a wide range of possible data. Results are presented for the clone mapping protocol currently being employed at Los Alamos National Laboratory on human chromosome 16 (Stallings et al., 1990, Proc. natl. Acad. Sci. U.S.A., 87, 6218-6222).

Mapping sequenced E.coli genes by computer: software, strategies and examples

Methods are presented for organizing and integrating DNA sequence data, restriction maps, and genetic maps for the same organism but from a variety of sources (databases, publications, personal communications). Proper software tools are essential for successful organization of such diverse data into an ordered, cohesive body of information, and a suite of novel software to support this endeavor is described. Though these tools automate much of the task, a variety of strategies is needed to cope with recalcitrant cases. We describe such strategies and illustrate their application with numerous examples. These strategies have allowed us to order, analyze, and display over one megabase of E. coli DNA sequence information. The integration task often exposes inconsistencies in the available data, perhaps caused by strain polymorphisms or human oversight, necessitating the application of sound biological judgment. The examples illustrate both the level of expertise required of the database curator and the knowledge gained as apparent inconsistencies are resolved. The software and mapping methods are applicable to the study of any genome for which a high resolution restriction map is available. They were developed to support a weakly coordinated sequencing effort involving many laboratories, but would also be useful for highly orchestrated sequencing projects.

Physical mapping of human chromosomes by repetitive sequence fingerprinting

We have developed an approach for identifying overlapping cosmid clones by exploiting the high density of repetitive sequences in complex genomes. Individual clones are fingerprinted, using a combination of restriction enzyme digestions followed by hybridization with selected classes of repetitive sequences. This "repeat fingerprinting" technique allows small regions of clone overlap (10-20%) to be unambiguously assigned. We demonstrate the utility of this approach, using the fingerprinting of 3145 cosmid clones (1.25 x coverage), containing one or more (GT)n repeats, from human chromosome 16. A statistical analysis was used to link these clones into 460 contiguous sequences (contigs), averaging 106 kilobases (kb) in length and representing approximately 54% (48.7 Mb) of the euchromatic arms of this chromosome. These values are consistent with theoretical calculations and indicate that 150- to 200-kb contigs can be generated with 1.5 x coverage. This strategy requires the fingerprinting of approximately one-fourth as many cosmids as random strategies requiring 50% minimum overlap for overlap detection. By "nucleating" at specific regions in the human genome, and exploiting the high density of interspersed sequences, this approach allows (i) the rapid generation of large (greater than 100-kb) contigs in the early stages of contig mapping and (ii) the production of a contig map with useful landmarks for rapid integration of the genetic and physical maps.

The complete AvrII restriction map of the Escherichia coli genome and comparisons of several laboratory strains

The complete 13 site AvrII restriction map of the genome of E coli strain MG1655 is presented and compared with several other E. coli strains. The map was determined primarily by isolating individual AvrII fragments from pulsed-field gels, and hybridizing these large probes to a battery of mapped E. coli clones in lambda vectors. AvrII restriction patterns for eight other laboratory strains were determined and maps for seven of them deduced from the gel and comparisons between the strain genotypes, the MG1655 map, and AvrII sites in E. coli sequences taken from Genbank.

The expected fraction of clonable genomic DNA

Random clone mapping of genomic DNA is a subject of great interest in molecular biology. E. coli has just been mapped and work is progressing on some human chromosomes. In this paper we give estimates of the fraction of genomic DNA which is not clonable by partial digest with a restriction enzyme.

DNA fingerprinting by sampled sequencing

We describe a method for characterizing DNA segments that combines limited sequencing with size separation of restriction fragments. As part of a multistep procedure, 5' overhangs of unknown sequence are generated by cleavage with a class IIS restriction enzyme. After labeling of these ends by using dideoxynucleotides tagged with distinctive fluorescent dyes, the restriction fragments are analyzed by polyacrylamide gel electrophoresis and detection of fluorescent emissions using a commercially available DNA sequencer. The nucleotide-specific fluorescent signatures permit determination of the terminal sequence for each labeled end. The set of labeled fragments, characterized by both size and terminal sequence, constitutes a fingerprint that can be used to compare DNA segments for overlap or relatedness. The inclusion of terminal sequence data dramatically increases the information content of the fingerprint, making comparisons more reliable and efficient than those based upon size alone.

Completion of the detailed restriction map of the E. coli genome by the isolation of overlapping cosmid clones

Ordered sets of cosmids derived from E. coli K-12 803 overlap the 6 remaining gaps left in the physical map of strain W3110. We present detailed restriction maps of the gaps and surrounding regions, thus providing a comparison of about 30% of the genome of the two E. coli strains. Our analysis shows that there is a high degree of homology between the strains, with only occasional restriction fragment differences. However, the large inversion occurring between rrnD (72.1') and rrnE (90.4') in strain W3110 is absent in strain 803. Instead, a new inversion and adjacent deletion near argF is present in strain 803. The distribution of cosmid clones at, and adjacent to, the gaps shows that all gaps except one were difficult to clone in both lambda and cosmid clones. A low copy number cosmid vector, pOU61cos, developed previously, was essential for cloning 3 of the 8 gaps.

Regulation of enterobacterial cephalosporinase production: the role of a membrane-bound sensory transducer

In clinical isolates of Enterobacter cloacae, resistance to the newer beta-lactam antibiotics often results from overproduction of a cephalosporinase encoded by the beta-lactam-inducible ampC gene. Regulation of ampC is controlled by the divergently expressed activator gene, ampR, and a second unlinked locus. In this presentation we show that although Escherichia coli has lost its ampR gene it has retained the second regulatory locus and that this comprises the bicistronic ampDE operon. Genetic and biochemical studies define the ampD gene as encoding a repressor for ampC transcription whereas the ampE gene product is a cytoplasmic membrane protein. Inactivation of the AmpD protein by mutation causes massive overproduction of cephalosporinase which, in E. cloacae, can terminate in therapeutic failure. In contrast, loss of AmpE results in a total block in induction, despite the presence of the activator, AmpR.

Cosmid-derived map of E. coli strain BHB2600 in comparison to the map of strain W3110

A physical map for the genome of E. coli K12 strain BHB2600 was constructed by use of 570 cloned DNA elements (CDEs) withdrawn from a cosmid library. Dot blot hybridisation was applied to establish contig interrelations with subsequent fine mapping achieved by analysis of EcoR1 restriction patterns on Southern blots. The derived map covers nearly 95% of the E. coli genome resulting in 12 minor gaps. It may be compared to the almost complete map for strain W3110 of Kohara et al. (1). Except for one tiny gap (lpp,36.5') remaining gaps in BHB2600 do not coincide with those in W3110 so that both maps complement each other establishing an essentially complete clone represented map. Besides numerous minute differences (site and fragment gains and losses) both strains harbour at differing positions extended rearrangements flanked by mutually inverted repetitive elements, in our case insertion elements (IS1 and IS5).

Transverse alternating field electrophoresis and applications to mammalian genome mapping

The transverse alternating field electrophoresis system is a pulsed field gel apparatus that has been used to separate DNA molecules that range in size from a few thousand to approximately 7 million base pairs. This apparatus uses a vertical gel and a simple electrode arrangement to produce electric fields that are uniform across all lanes of the gel. The velocity of identical molecules does not vary from lane to lane, and hence there is no distortion in the paths of the DNA. The performance of this system is illustrated here using the chromosomes from S. pombe and S. cerevisiae, and restriction enzyme digested mammalian DNA. The mobility of molecules up to 1100 kilobase pairs is linear with size and can be accomplished in overnight runs. Resolution of very large molecules requires electrophoresis for several days, but molecules from 200 to 7000 kilobase pairs can be separated on a single gel. This electrophoresis system has been used extensively in the construction of a physical map of human chromosome 21, and examples of this application are discussed.

Construction of an ordered cosmid collection of the Escherichia coli K-12 W3110 chromosome

A cosmid library of the Escherichia coli K-12 W3110 chromosome was constructed in which clones were assigned to locations on the chromosome map by hybridization and genetic marker complementation tests. Approximately 70% of the genome was represented by this library. The identified clones can be maintained in the homologous system and would facilitate genetic studies of E. coli.

Mono- through hexanucleotide composition of the sense strand of yeast DNA: a Markov chain analysis

Here we compare several methods for predicting oligonucleotide frequencies in 392 kb of yeast DNA. As in previous work on E. coli, a relatively simple equation based on tetranucleotide frequencies can be used in predicting the frequencies of longer oligonucleotides. For example, the mean of observed/expected abundances of 4,096 hexamers was 1.00 with a sample standard deviation of .18. This simple predictor arises by considering each base on the sense strand of yeast to depend only on the three bases 5' to it (a 3rd order Markov chain) and is more accurate in estimating oligonucleotide frequencies than other statistical methods examined. This equation is useful in predicting restriction enzyme fragment sizes, selecting restriction enzymes that cut preferentially in coding vs noncoding regions, and in constructing detailed physical maps of whole genomes. When ranked highest to lowest abundance, the observed frequencies of oligomers of a given length (up to 6 bases) are closely tracked by the predicted abundances of a 3rd or 4th order Markov chain. These ordered abundance curves have a power curve shape with a broad linear range with a sharp break at the top end of the curve. There is also a strong disparity between the most and least abundant oligomer with for example a 79-fold variation between the most and least abundant hexamer. The curves reveal a strong dependence of oligomer frequencies on base composition. Unlike E. Coli, there is no sharp downturn at the low end of the curves and hence, no class of oligomers rare relative to other oligomers of the same length.

Genomic mapping by fingerprinting random clones: a mathematical analysis

Results from physical mapping projects have recently been reported for the genomes of Escherichia coli, Saccharomyces cerevisiae, and Caenorhabditis elegans, and similar projects are currently being planned for other organisms. In such projects, the physical map is assembled by first "fingerprinting" a large number of clones chosen at random from a recombinant library and then inferring overlaps between clones with sufficiently similar fingerprints. Although the basic approach is the same, there are many possible choices for the fingerprint used to characterize the clones and the rules for declaring overlap. In this paper, we derive simple formulas showing how the progress of a physical mapping project is affected by the nature of the fingerprinting scheme. Using these formulas, we discuss the analytic considerations involved in selecting an appropriate fingerprinting scheme for a particular project.

Molecular cloning and expression of the Escherichia coli dimethyl sulfoxide reductase operon

The dimethyl sulfoxide (DMSO) reductase operon coding for a membrane-bound iron-sulfur, molybdoenzyme, which functions as a terminal reductase in Escherichia coli, has been isolated and cloned from an E. coli gene bank. Two clones, MV12(pLC19-36) and MV12(pLC43-43), overexpressed both DMSO and trimethylamine N-oxide (TMAO) reductase activities 13- to 15-fold compared with wild-type cells. Amplification was highest in cells grown anaerobically on fumarate, while cells grown on DMSO or TMAO displayed reduced levels of enzyme amplification. Growth on nitrate or aerobic growth repressed expression of the enzyme. A 6.5-kilobase-pair DNA restriction endonuclease fragment was subcloned from pLC19-36 into the vector pBR322, yielding a recombinant DMSO reductase plasmid, pDMS159. Two polypeptides were amplified and identified on sodium dodecyl sulfate-polyacrylamide gels of proteins from E. coli HB101 harboring pDMS159: a membrane-bound protein with molecular weight 82,600 and a soluble polypeptide with molecular weight 23,600. Three plasmid-encoded polypeptides with molecular weights of 87,500, 23,300, and 22,600 were detected by in vivo transcription/translation studies. The smallest subunit was poorly defined and not detectable by Coomassie blue staining. The DMSO reductase operon was localized to the 20.0-min position on the E. coli linkage map.

Construction of long-range restriction maps in human DNA using pulsed field gel electrophoresis

Pulsed field gel electrophoresis (PFGE) is a powerful new tool for genetic analysis that can be applied to a variety of problems concerning genome structure and organization. This technique uses an agarose gel matrix to separate DNA molecules in a size range from 40 kb to 2,000 kb, molecules far larger than the maximum separable using standard agarose gel electrophoresis. The PFGE method can be used to separate the intact chromosomes from lower eukaryotes or to separate very large DNA fragments from higher eukaryotes generated by digestion with restriction endonucleases whose cleavage sites are rare. This paper describes the use of PFGE for construction of long-range restriction maps in the human genome and includes detailed methods for all steps. A pulsed field gel device that utilizes a rotating platform for altering the applied electric field is also described. Map construction is illustrated using a cloned DNA fragment (D3S2) from human chromosome 3. Several technical problems specific for mammalian genomes are discussed.

The physical map of the whole E. coli chromosome: application of a new strategy for rapid analysis and sorting of a large genomic library

Thirty-four hundred lambda phage clones containing segments of the E. coli chromosome were isolated and used to construct a 4700 kb long integrated restriction map for eight six-base-recognizing enzymes by a rapid mass-analysis method. Our strategy was to measure the sizes of partial restriction enzyme digests by hybridization with a vector probe in a manner analogous to nucleotide sequencing. The data were sorted into groups by a computer program and the boundary clones were further correlated with each other using a mass hybridization method. These clones can be exploited for the isolation of any desired E. coli genes if their map positions are known. Also, the strategy is applicable to analyses of the genomes of other organisms.