Genome Mapping by Fluorescent Fingerprinting

The construction of sequence-ready maps of overlapping genomic clones is central to large-scale genome sequencing. We have implemented a method for fluorescent fingerprinting of bacterial clones to assemble contig maps. The method utilizes three spectrally distinct fluorescently tagged dideoxy ATPs to specifically label the HindIII termini in HindIII and Sau3AI restriction digests of clones that are multiplexed prior to electrophoresis and data collection. There is excellent reproducibility of raw data, improved resolution of large fragments, and concordance between the results obtained using this and the equivalent radioactive protocol. This method also allows detection of smaller overlaps between clones when compared to the analysis of restriction digests on nondenaturing agarose gels.     

Internet Contig Explorer (iCE)—A Tool for Visualizing Clone Fingerprint Maps

Fingerprinted clone physical maps have proven useful in various applications, supporting both whole-genome and region-specific DNA sequencing as well as gene cloning studies. Fingerprint maps have been generated for several genomes, including those of human, mouse, rat, the nematodes Caenorhabditis elegans and Caenorhabditis briggsae, Arabidopsis thaliana and rice. Fingerprint maps of other genomes, including those of fungi, bacteria, poplar, and the cow, are being generated. The increasing use of fingerprint maps in genomic research has spawned a need in the research community for intuitive computer tools that facilitate viewing of the maps and the underlying fingerprint data. In this report we describe a new Java-based application called iCE (Internet Contig Explorer) that has been designed to provide views of fingerprint maps and associated data. Users can search for and display individual clones, contigs, clone fingerprints, clone insert sizes and markers. Users can also load into the software lists of particular clones of interest and view their fingerprints. iCE is being used at our Genome Centre to offer up to the research community views of the mouse, rat, bovine, C. briggsae, and several fungal genome bacterial artificial chromosome (BAC) fingerprint maps we have either completed or are currently constructing. We are also using iCE as part of the Rat Genome Sequencing Project to manage our provision of rat BAC clones for sequencing at the Human Genome Sequencing Center at the Baylor College of Medicine.     

A 7.5 Mb sequence-ready PAC contig and gene expression map of human chromosome 11p13-p14.1

The region p13 of the short arm of human chromosome 11 has been studied intensely during the search for genes involved in the etiology of the Wilms' tumor, aniridia, genitourinary abnormalities, mental retardation (WAGR) syndrome, and related conditions. The gene map for this region is far from being complete, however, strengthening the need for additional gene identification efforts. We describe the extension of an existing contig map with P1-derived artificial chromosomes (PACs) to cover 7.5 Mb of 11p13-14.1. The extended sequence-ready contig was established by end probe walking and fingerprinting and consists of 201 PAC clones. Utilizing bins defined by overlapping PACs, we generated a detailed gene map containing 20 genes as well as 22 anonymous ESTs which have been identified by searching the RH databases. RH maps and our established gene map show global correlation, but the limits of resolution of the current RH panels are evident at this scale. Initial expression studies on the novel genes have been performed by Northern blot analyses. To extend these expression profiles, corresponding mouse cDNA clones were identified by database search and employed for Northern blot analyses and RNA in situ hybridizations to mouse embryo sections. Genomic sequencing of clones along a minimal tiling path through the contig is currently under way and will facilitate these expression studies by in silico gene identification approaches.     

A contiguous 3-Mb sequence-ready map in the S3-MX region on 21q22.2 based on high- throughput nonisotopic library screenings

Progress in complete genomic sequencing of human chromosome 21 relies on the construction of high-quality bacterial clone maps spanning large chromosomal regions. To achieve this goal, we have applied a strategy based on nonradioactive hybridizations to contig building. A contiguous sequence-ready map was constructed in the Down syndrome congenital heart disease (DS-CHD) region in 21q22.2, as a framework for large-scale genomic sequencing and positional candidate gene approach. Contig assembly was performed essentially by high throughput nonisotopic screenings of genomic libraries, prior to clone validation by (1) restriction digest fingerprinting, (2) STS analysis, (3) Southern hybridizations, and (4) FISH analysis. The contig contains a total of 50 STSs, of which 13 were newly isolated. A minimum tiling path (MTP) was subsequently defined that consists of 20 PACs, 2 BACs, and 5 cosmids covering 3 Mb between D21S3 and MX1. Gene distribution in the region includes 9 known genes (c21-LRP, WRB, SH3BGR, HMG14, PCP4, DSCAM, MX2, MX1, and TMPRSS2) and 14 new additional gene signatures consisting of cDNA selection products and ESTs. Forthcoming genomic sequence information will unravel the structural organization of potential candidate genes involved in specific features of Down syndrome pathogenesis.     

Comparative genome mapping in the sequence-based era: early experience with human chromosome 7

The success of the ongoing Human Genome Project has resulted in accelerated plans for completing the human genome sequence and the earlier-than-anticipated initiation of efforts to sequence the mouse genome. As a complement to these efforts, we are utilizing the available human sequence to refine human-mouse comparative maps and to assemble sequence-ready mouse physical maps. Here we describe how the first glimpses of genomic sequence from human chromosome 7 are directly facilitating these activities. Specifically, we are actively enhancing the available human-mouse comparative map by analyzing human chromosome 7 sequence for the presence of orthologs of mapped mouse genes. Such orthologs can then be precisely positioned relative to mapped human STSs and other genes. The chromosome 7 sequence generated to date has allowed us to more than double the number of genes that can be placed on the comparative map. The latter effort reveals that human chromosome 7 is represented by at least 20 orthologous segments of DNA in the mouse genome. A second component of our program involves systematically analyzing the evolving human chromosome 7 sequence for the presence of matching mouse genes and expressed-sequence tags (ESTs). Mouse-specific hybridization probes are designed from such sequences and used to screen a mouse bacterial artificial chromosome (BAC) library, with the resulting data used to assemble BAC contigs based on probe-content data. Nascent contigs are then expanded using probes derived from newly generated BAC-end sequences. This approach produces BAC-based sequence-ready maps that are known to contain a gene(s) and are homologous to segments of the human genome for which sequence is already available. Our ongoing efforts have thus far resulted in the isolation and mapping of >3,800 mouse BACs, which have been assembled into >100 contigs. These contigs include >250 genes and represent approximately 40% of the mouse genome that is homologous to human chromosome 7. Together, these approaches illustrate how the availability of genomic sequence directly facilitates studies in comparative genomics and genome evolution.     

Parallel construction of orthologous sequence-ready clone contig maps in multiple species

Comparison is a fundamental tool for analyzing DNA sequence. Interspecies sequence comparison is particularly powerful for inferring genome function and is based on the simple premise that conserved sequences are likely to be important. Thus, the comparison of a genomic sequence with its orthologous counterpart from another species is increasingly becoming an integral component of genome analysis. In ideal situations, such comparisons are performed with orthologous sequences from multiple species. To facilitate multispecies comparative sequence analysis, a robust and scalable strategy for simultaneously constructing sequence-ready bacterial artificial chromosome (BAC) contig maps from targeted genomic regions has been developed. Central to this approach is the generation and utilization of "universal" oligonucleotide-based hybridization probes ("overgo" probes), which are designed from sequences that are highly conserved between distantly related species. Large collections of these probes are used en masse to screen BAC libraries from multiple species in parallel, with the isolated clones assembled into physical contig maps. To validate the effectiveness of this strategy, efforts were focused on the construction of BAC-based physical maps from multiple mammalian species (chimpanzee, baboon, cat, dog, cow, and pig). Using available human and mouse genomic sequence and a newly developed computer program to design the requisite probes, sequence-ready maps were constructed in all species for a series of targeted regions totaling approximately 16 Mb in the human genome. The described approach can be used to facilitate the multispecies comparative sequencing of targeted genomic regions and can be adapted for constructing BAC contig maps in other vertebrates.     

Construction of a BAC contig map of chromosome 16q by two-dimensional overgo hybridization

We have used sequence-based markers from an integrated YAC STS-content/somatic cell hybrid breakpoint physical map and radiation hybrid maps of human chromosome 16 to construct a new sequence-ready BAC map of the long arm of this chromosome. The integrated physical map was generated previously in our laboratory and contains 1150 STSs, providing a marker on average every 78 kb on the euchromatic arms of chromosome 16. The other two maps used for this effort were the radiation hybrid maps of chromosome 16 from Whitehead Institute and Stanford University. To create large sequenceable targets of this chromosome, we used a systematic approach to screen high-density BAC filters with probes generated from overlapping oligonucleotides (overgos). We first identified all available sequences in the three maps. These include sequences from genes, ESTs, STSs, and cosmid end sequences. We then used BLASTto identify 36-bp unique fragments of DNA for overgo probes. A total of 906 overgos were selected from the long arm of chromosome 16. Hybridizations occurred in three stages: (1) superpool hybridizations against the 12x coverage human BAC library (RPCI-11); (2) two-dimensional hybridizations against rearrayed positive BACs identified in the superpool hybridizations; and (3) pooled tertiary hybridizations for those overgos that had ambiguous positives remaining after the two-dimensional hybridization. For the superpool hybridizations, up to 236 overgos have been pooled in a single hybridization against the 12x BAC library. A total of 5187 positive BACs from chromosome 16q were identified as a result of five superpool hybridizations. These positive clones were rearrayed on membranes and hybridized with 161 two-dimensional subpools of overgos to determine which BAC clones were positive for individual overgos. An additional 46 tertiary hybridizations were required to resolve ambiguous overgo-BAC relationships. Thus, after a total of 212 hybridizations, we have constructed an initial probe-content BAC map of chromosome 16q consisting of 828 overgo markers and 3363 BACs providing >85% coverage of the long arm of this chromosome. The map has been confirmed by the fingerprinting data and BAC end PCR screening.     

A clone-array pooled shotgun strategy for sequencing large genomes

A simplified strategy for sequencing large genomes is proposed. Clone-Array Pooled Shotgun Sequencing (CAPSS) is based on pooling rows and columns of arrayed genomic clones, for shotgun library construction. Random sequences are accumulated, and the data are processed by sequential comparison of rows and columns to assemble the sequence of clones at points of intersection. Compared with either a clone-by-clone approach or whole-genome shotgun sequencing, CAPSS requires relatively few library constructions and only minimal computational power for a complete genome assembly. The strategy is suitable for sequencing large genomes for which there are no sequence-ready maps, but for which relatively high resolution STS maps and highly redundant BAC libraries are available. It is immediately applicable to the sequencing of mouse, rat, zebrafish, and other important genomes, and can be managed in a cooperative fashion to take advantage of a distributed international DNA sequencing capacity.     

High-Resolution Landmark Framework for the Sequence-Ready Mapping of Xq23–q26.1

We have established a landmark framework map over 20-25 Mb of the long arm of the human X chromosome using yeast artificial chromosome (YAC) clones. The map has approximately one landmark per 45 kb of DNA and stretches from DXS7531 in proximal Xq23 to DXS895 in proximal Xq26, connecting to published framework maps on its proximal and distal sides. There are three gaps in the framework map resulting from the failure to obtain clone coverage from the YAC resources available. Estimates of the maximum sizes of these gaps have been obtained. The four YAC contigs have been positioned and oriented using somatic-cell hybrids and fluorescence in situ hybridization, and the largest is estimated to cover approximately 15 Mb of DNA. The framework map is being used to assemble a sequence-ready map in large-insert bacterial clones, as part of an international effort to complete the sequence of the X chromosome. PAC and BAC contigs currently cover 18 Mb of the region, and from these, 12 Mb of finished sequence is available.     

A high-throughput AFLP-based method for constructing integrated genetic and physical maps: progress toward a sorghum genome map

Sorghum is an important target for plant genomic mapping because of its adaptation to harsh environments, diverse germplasm collection, and value for comparing the genomes of grass species such as corn and rice. The construction of an integrated genetic and physical map of the sorghum genome (750 Mbp) is a primary goal of our sorghum genome project. To help accomplish this task, we have developed a new high-throughput PCR-based method for building BAC contigs and locating BAC clones on the sorghum genetic map. This task involved pooling 24,576 sorghum BAC clones ( approximately 4x genome equivalents) in six different matrices to create 184 pools of BAC DNA. DNA fragments from each pool were amplified using amplified fragment length polymorphism (AFLP) technology, resolved on a LI-COR dual-dye DNA sequencing system, and analyzed using Bionumerics software. On average, each set of AFLP primers amplified 28 single-copy DNA markers that were useful for identifying overlapping BAC clones. Data from 32 different AFLP primer combinations identified approximately 2400 BACs and ordered approximately 700 BAC contigs. Analysis of a sorghum RIL mapping population using the same primer pairs located approximately 200 of the BAC contigs on the sorghum genetic map. Restriction endonuclease fingerprinting of the entire collection of sorghum BAC clones was applied to test and extend the contigs constructed using this PCR-based methodology. Analysis of the fingerprint data allowed for the identification of 3366 contigs each containing an average of 5 BACs. BACs in approximately 65% of the contigs aligned by AFLP analysis had sufficient overlap to be confirmed by DNA fingerprint analysis. In addition, 30% of the overlapping BACs aligned by AFLP analysis provided information for merging contigs and singletons that could not be joined using fingerprint data alone. Thus, the combination of fingerprinting and AFLP-based contig assembly and mapping provides a reliable, high-throughput method for building an integrated genetic and physical map of the sorghum genome.     

Software for automated analysis of DNA fingerprinting gels

Here we describe software tools for the automated detection of DNA restriction fragments resolved on agarose fingerprinting gels. We present a mathematical model for the location and shape of the restriction fragments as a function of fragment size, with model parameters determined empirically from "marker" lanes containing molecular size standards. Automated identification of restriction fragments involves several steps, including: image preprocessing, to put the data in a form consistent with a linear model; marker lane analysis, for determination of the model parameters; and data lane analysis, a procedure for detecting restriction fragment multiplets while simultaneously determining the amplitude curve that describes restriction fragment amplitude as a function of mobility. In validation experiments conducted on fingerprinted and sequenced Bacterial Artificial Chromosome (BAC) clones, sensitivity and specificity of restriction fragment identification exceeded 96% on restriction fragments ranging in size from 600 base pairs (bp) to 30,000 bp. The integrated suite of software tools, written in MATLAB and collectively called BandLeader, is in use at the BC Cancer Agency Genome Sciences Centre (GSC) and the Washington University Genome Sequencing Center, and has been provided to the Wellcome Trust Sanger Institute and the Whitehead Institute. Employed in a production mode at the GSC, BandLeader has been used to perform automated restriction fragment identification for more than 850,000 BAC clones for mouse, rat, bovine, and poplar fingerprint mapping projects.     

A 3-Mb High-Resolution BAC/PAC Contig of 12q22 Encompassing the 830-kb Consensus Minimal Deletion in Male Germ Cell Tumors

Cytogenetic and molecular genetic analyses have shown that the 12q22 region is recurrently deleted in male germ cell tumors (GCTs), suggesting that this site may harbor a tumor suppressor gene (TSG). Previous loss of heterozygosity (LOH) analyses identified a consensus minimal deleted region between the markers D12S377 and D12S296, and a YAC clone contig covering the region was generated. Here, we describe a high-resolution sequence-ready physical map of this contig covering a 3-Mb region. The map comprised of 52 cosmids, 49 PACs, and 168 BACs that were anchored to the previous YAC contig; 99 polymorphic, nonpolymorphic, EST, and gene-based markers are now placed on this map in a unique order. Of these, 61 markers were isolated in the present study, including one that was polymorphic. In addition, we have narrowed the minimal deletion to approximately 830 kb between D12S1716 (proximal) and P382A8-AG (distal) by LOH analysis of 108 normal-tumor DNAs from GCT patients using 21 polymorphic STSs. These physical and deletion maps should prove useful for identification of the candidate TSG in GCTs, provide framework to generate complete DNA sequence, and ultimately generate a gene map of this segment of the chromosome 12. [The sequence data described in this paper have been submitted to the Genome Survey Sequence under accession nos. AQ254896-AQ254955 and AQ269251-AQ269266. Online supplementary material is available at http://www.genome.org]     

Characterization of Short Tandem Repeats from Thirty-One Human Telomeres

Completion of genetic and physical maps requires markers from the ends (telomeres) of every human chromosome. We have searched for short tandem repeats (microsatellites) in cosmid and P1 clones and generated 661 sequence-tagged sites (STS) from the terminal 300 kb of 31 human chromosome ends. PCR assays were successfully designed for 58 microsatellites and mapped both genetically and on radiation hybrids (RHs) to confirm their telomeric location. Sequence analysis revealed marked variation in sequence composition, consistent with the hypothesis that even very highly GC-rich chromosome bands (the T bands) are not homogenous. The STSs that we have generated will be a necessary resource for the construction of physical maps of these complex regions of the genome.     

Contigs built with fingerprints, markers, and FPC V4.7

Contigs have been assembled, and over 2800 clones selected for sequencing for human chromosomes 9, 10 and 13. Using the FPC (FingerPrinted Contig) software, the contigs are assembled with markers and complete digest fingerprints, and the contigs are ordered and localised by a global framework. Publicly available resources have been used, such as, the 1998 International Gene Map for the framework and the GSC Human BAC fingerprint database for the majority of the fingerprints. Additional markers and fingerprints are generated in-house to supplement this data. To support the scale up of building maps, FPC V4.7 has been extended to use markers with the fingerprints for assembly of contigs, new clones and markers can be automatically added to existing contigs, and poorly assembled contigs are marked accordingly. To test the automatic assembly, a simulated complete digest of 110 Mb of concatenated human sequence was used to create datasets with varying coverage, length of clones, and types of error. When no error was introduced and a tolerance of 7 was used in assembly, the largest contig with no false positive overlaps has 9534 clones with 37 out-of-order clones, that is, the starting coordinates of adjacent clones are in the wrong order. This paper describes the new features in FPC, the scenario for building the maps of chromosomes 9, 10 and 13, and the results from the simulation.     

Dynamic building of a BAC clone tiling path for the Rat Genome Sequencing Project

CLONEPICKER is a software pipeline that integrates sequence data with BAC clone fingerprints to dynamically select a minimal overlapping clone set covering the whole genome. In the Rat Genome Sequencing Project (RGSP), a hybrid strategy of "clone by clone" and "whole genome shotgun" approaches was used to maximize the merits of both approaches. Like the "clone by clone" method, one key challenge for this strategy was to select a low-redundancy clone set that covered the whole genome while the sequencing is in progress. The CLONEPICKER pipeline met this challenge using restriction enzyme fingerprint data, BAC end sequence data, and sequences generated from individual BAC clones as well as WGS reads. In the RGSP, an average of 7.5 clones was identified from each side of a seed clone, and the minimal overlapping clones were reliably selected. Combined with the assembled BAC fingerprint map, a set of BAC clones that covered >97% of the genome was identified and used in the RGSP.     

Characterization of Physical Gap Sizes at Human Telomeres

Genome-wide physical and genetic mapping efforts have not yet fully addressed the problem of closure at the telomeric ends of human chromosomes. Targeted efforts at cloning human and mouse telomeres have succeeded in identifying unique sequences at most telomeres, but gap sizes between these telomere clones and the distal markers on integrated genetic/physical maps remain largely unknown. As telomeric regions are known to be the most gene-rich regions of the human genome, filling these gaps should have a high priority in completion of the Human Genome Project. We reported previously a first generation set of unique sequence probes for human telomeric regions. Of 41 human telomere regions, 33 were represented by unique clones with a known distance (1 Mb, thus defining the physical mapping task for filling telomeric gaps.     

Mapping of the human P84 gene to the subtelomeric region of chromosome 20p

P84 is a novel neural adhesion molecule that may play an important role in synaptogenesis. We have recently cloned a murine cDNA encoding the P84 adhesion molecule. The human homologue of P84 has previously been isolated (by others) as a brain specific cDNA containing CCA repeats. We have mapped the human P84 gene to the subtelomeric region of chromosome 20p (20p13) by FISH. In addition, we have been albe to place P84 onto the high resolution physical map of the human genome by utilizing the Unigene database. P84 maps to several YAC clones, between STS markers IB255 and WI-9632, and very close to the polymorphic marker D20S199, in an interval of less than 1 Mb on 20p13. P84 is a strong candidate gene for neurological disorders which map into this region.     

The ins and outs of DNA fingerprinting the infectious fungi

DNA fingerprinting methods have evolved as major tools in fungal epidemiology. However, no single method has emerged as the method of choice, and some methods perform better than others at different levels of resolution. In this review, requirements for an effective DNA fingerprinting method are proposed and procedures are described for testing the efficacy of a method. In light of the proposed requirements, the most common methods now being used to DNA fingerprint the infectious fungi are described and assessed. These methods include restriction fragment length polymorphisms (RFLP), RFLP with hybridization probes, randomly amplified polymorphic DNA and other PCR-based methods, electrophoretic karyotyping, and sequencing-based methods. Procedures for computing similarity coefficients, generating phylogenetic trees, and testing the stability of clusters are then described. To facilitate the analysis of DNA fingerprinting data, computer-assisted methods are described. Finally, the problems inherent in the collection of test and control isolates are considered, and DNA fingerprinting studies of strain maintenance during persistent or recurrent infections, microevolution in infecting strains, and the origin of nosocomial infections are assessed in light of the preceding discussion of the ins and outs of DNA fingerprinting. The intent of this review is to generate an awareness of the need to verify the efficacy of each DNA fingerprinting method for the level of genetic relatedness necessary to answer the epidemiological question posed, to use quantitative methods to analyze DNA fingerprint data, to use computer-assisted DNA fingerprint analysis systems to analyze data, and to file data in a form that can be used in the future for retrospective and comparative studies.     

Integrated and sequence-ordered BAC- and YAC-based physical maps for the rat genome

As part of the effort to sequence the genome of Rattus norvegicus, we constructed a physical map comprised of fingerprinted bacterial artificial chromosome (BAC) clones from the CHORI-230 BAC library. These BAC clones provide ~13-fold redundant coverage of the genome and have been assembled into 376 fingerprint contigs. A yeast artificial chromosome (YAC) map was also constructed and aligned with the BAC map via fingerprinted BAC and P1 artificial chromosome clones (PACs) sharing interspersed repetitive sequence markers with the YAC-based physical map. We have annotated 95% of the fingerprint map clones in contigs with coordinates on the version 3.1 rat genome sequence assembly, using BAC-end sequences and in silico mapping methods. These coordinates have allowed anchoring 358 of the 376 fingerprint map contigs onto the sequence assembly. Of these, 324 contigs are anchored to rat genome sequences localized to chromosomes, and 34 contigs are anchored to unlocalized portions of the rat sequence assembly. The remaining 18 contigs, containing 54 clones, still require placement. The fingerprint map is a high-resolution integrative data resource that provides genome-ordered associations among BAC, YAC, and PAC clones and the assembled sequence of the rat genome.