A BAC-based physical map of the channel catfish genome

Catfish is the major aquaculture species in the United States. To enhance its genome studies involving genetic linkage and comparative mapping, a bacterial artificial chromosome (BAC) contig-based physical map of the channel catfish (Ictalurus punctatus) genome was generated using four-color fluorescence-based fingerprints. Fingerprints of 34,580 BAC clones (5.6x genome coverage) were generated for the FPC assembly of the BAC contigs. A total of 3307 contigs were assembled using a cutoff value of 1x10(-20). Each contig contains an average of 9.25 clones with an average size of 292 kb. The combined contig size for all contigs was 0.965 Gb, approximately the genome size of the channel catfish. The reliability of the contig assembly was assessed by both hybridization of gene probes to BAC clones contained in the fingerprinted assembly and validation of randomly selected contigs using overgo probes designed from BAC end sequences. The presented physical map should greatly enhance genome research in the catfish, particularly aiding in the identification of genomic regions containing genes underlying important performance traits.     

Characterization of three maize bacterial artificial chromosome libraries toward anchoring of the physical map to the genetic map using high-density bacterial artificial chromosome filter hybridization

Three maize (Zea mays) bacterial artificial chromosome (BAC) libraries were constructed from inbred line B73. High-density filter sets from all three libraries, made using different restriction enzymes (HindIII, EcoRI, and MboI, respectively), were evaluated with a set of complex probes including the 185-bp knob repeat, ribosomal DNA, two telomere-associated repeat sequences, four centromere repeats, the mitochondrial genome, a multifragment chloroplast DNA probe, and bacteriophage lambda. The results indicate that the libraries are of high quality with low contamination by organellar and lambda-sequences. The use of libraries from multiple enzymes increased the chance of recovering each region of the genome. Ninety maize restriction fragment-length polymorphism core markers were hybridized to filters of the HindIII library, representing 6x coverage of the genome, to initiate development of a framework for anchoring BAC contigs to the intermated B73 x Mo17 genetic map and to mark the bin boundaries on the physical map. All of the clones used as hybridization probes detected at least three BACs. Twenty-two single-copy number core markers identified an average of 7.4 +/- 3.3 positive clones, consistent with the expectation of six clones. This information is integrated into fingerprinting data generated by the Arizona Genomics Institute to assemble the BAC contigs using fingerprint contig and contributed to the process of physical map construction.     

A Physical Map of the Short Arm of Wheat Chromosome 1A

Bread wheat (Triticum aestivum) has a large and highly repetitive genome which poses major technical challenges for its study. To aid map-based cloning and future genome sequencing projects, we constructed a BAC-based physical map of the short arm of wheat chromosome 1A (1AS). From the assembly of 25,918 high information content (HICF) fingerprints from a 1AS-specific BAC library, 715 physical contigs were produced that cover almost 99% of the estimated size of the chromosome arm. The 3,414 BAC clones constituting the minimum tiling path were end-sequenced. Using a gene microarray containing ∼40 K NCBI UniGene EST clusters, PCR marker screening and BAC end sequences, we arranged 160 physical contigs (97 Mb or 35.3% of the chromosome arm) in a virtual order based on synteny with Brachypodium, rice and sorghum. BAC end sequences and information from microarray hybridisation was used to anchor 3.8 Mbp of Illumina sequences from flow-sorted chromosome 1AS to BAC contigs. Comparison of genetic and synteny-based physical maps indicated that ∼50% of all genetic recombination is confined to 14% of the physical length of the chromosome arm in the distal region. The 1AS physical map provides a framework for future genetic mapping projects as well as the basis for complete sequencing of chromosome arm 1AS.     

Evaluation of a cosmid contig physical map of human chromosome 16.

A cosmid contig physical map of human chromosome 16 has been developed by repetitive sequence finger-printing of approximately 4000 cosmid clones obtained from a chromosome 16-specific cosmid library. The arrangement of clones in contigs is determined by (1) estimating cosmid length and determining the likelihoods for all possible pairwise clone overlaps, using the fingerprint data, and (2) using an optimization technique to fit contig maps to these estimates. Two important questions concerning this contig map are how much of chromosome 16 is covered and how accurate are the assembled contigs. Both questions can be addressed by hybridization of single-copy sequence probes to gridded arrays of the cosmids. All of the fingerprinted clones have been arrayed on nylon membranes so that any region of interest can be identified by hybridization. The hybridization experiments indicate that approximately 84% of the euchromatic arms of chromosome 16 are covered by contigs and singleton cosmids. Both grid hybridization (26 contigs) and pulsed-field gel electrophoresis experiments (11 contigs) confirmed the assembled contigs, indicating that false positive overlaps occur infrequently in the present map. Furthermore, regional localization of 93 contigs and singleton cosmids to a somatic cell hybrid mapping panel indicates that there is no bias in the coverage of the euchromatic arms.     

Integrating sequence with FPC fingerprint maps

Recent advances in both clone fingerprinting and draft sequencing technology have made it increasingly common for species to have a bacterial artificial clone (BAC) fingerprint map, BAC end sequences (BESs) and draft genomic sequence. The FPC (fingerprinted contigs) software package contains three modules that maximize the value of these resources. The BSS (blast some sequence) module provides a way to easily view the results of aligning draft sequence to the BESs, and integrates the results with the following two modules. The MTP (minimal tiling path) module uses sequence and fingerprints to determine a minimal tiling path of clones. The DSI (draft sequence integration) module aligns draft sequences to FPC contigs, displays them alongside the contigs and identifies potential discrepancies; the alignment can be based on either individual BES alignments to the draft, or on the locations of BESs that have been assembled into the draft. FPC also supports high-throughput fingerprint map generation as its time-intensive functions have been parallelized for Unix-based desktops or servers with multiple CPUs. Simulation results are provided for the MTP, DSI and parallelization. These features are in the FPC V9.3 software package, which is freely available.     

Construction of BAC library for the amphibian Xenopus tropicalis

Xenopus tropicalis has become an alternative model to the amphibian Xenopus laevis because it is better suited for genetic and genomic studies. We have constructed a genomic BAC library consisting of over 100,000 clones from sperm of Xenopus tropicalis. Analysis by pulsed field gel electrophoresis of representative BAC clones indicated the average size of insert DNA to be 100 kb, and we estimated the library covers 6 times the Xenopus tropicalis genome of 1.7 x 10(9) base pairs. To evaluate the BAC library, we attempted to isolate BAC clones which contain a protocadherin gamma (Pcdh gamma) gene and found that the isolated BAC clones are assembled as two separate contigs. This result suggests the presence of at least two clusters for the Pcdh gamma gene in the genome of X. tropicalis.     

Towards a physical map of the Drosophila melanogaster genome: mapping of cosmid clones within defined genomic divisions.

A physical map of the D. melanogaster genome is being constructed, in the form of overlapping cosmid clones that are assigned to specific polytene chromosome sites. A master library of ca. 20,000 cosmids is screened with probes that correspond to numbered chromosomal divisions (ca. 1% of the genome); these probes are prepared by microdissection and PCR-amplification of individual chromosomes. The 120 to 250 cosmids selected by each probe are fingerprinted by Hinfl digestion and gel electrophoresis, and overlaps are detected by computer analysis of the fingerprints, permitting us to assemble sets of contiguous clones (contigs). Selected cosmids, both from contigs and unattached, are then localized by in situ hybridization to polytene chromosomes. Crosshybridization analysis using end probes links some contigs, and hybridization to previously cloned genes relates the physical to the genetic map. This approach has been used to construct a physical map of the 3.8 megabase DNA in the three distal divisions of the x chromosome. The map is represented by 181 canonical cosmids, of which 108 clones in contigs and 32 unattached clones have been mapped individually by in situ hybridization to chromosomes. Our current database of in situ hybridization results also includes the beginning of a physical map for the rest of the genome: 162 cosmids have been assigned by in situ hybridization to 129 chromosomal subdivisions elsewhere in the genome, representing 5 to 6 megabases of additional mapped DNA.     

Bacterial artificial chromosome-based physical map of the rice genome constructed by restriction fingerprint analysis

Genome-wide physical mapping with bacteria-based large-insert clones (e.g., BACs, PACs, and PBCs) promises to revolutionize genomics of large, complex genomes. To accelerate rice and other grass species genome research, we developed a genome-wide BAC-based map of the rice genome. The map consists of 298 BAC contigs and covers 419 Mb of the 430-Mb rice genome. Subsequent analysis indicated that the contigs constituting the map are accurate and reliable. Particularly important to proficiency were (1) a high-resolution, high-throughput DNA sequencing gel-based electrophoretic method for BAC fingerprinting, (2) the use of several complementary large-insert BAC libraries, and (3) computer-aided contig assembly. It has been demonstrated that the fingerprinting method is not significantly influenced by repeated sequences, genome size, and genome complexity. Use of several complementary libraries developed with different restriction enzymes minimized the "gaps" in the physical map. In contrast to previous estimates, a clonal coverage of 6.0-8.0 genome equivalents seems to be sufficient for development of a genome-wide physical map of approximately 95% genome coverage. This study indicates that genome-wide BAC-based physical maps can be developed quickly and economically for a variety of plant and animal species by restriction fingerprint analysis via DNA sequencing gel-based electrophoresis.     

An integrated BAC and genome sequence physical map of Phytophthora sojae

Phytophthora spp. are serious pathogens that threaten numerous cultivated crops, trees, and natural vegetation worldwide. The soybean pathogen P. sojae has been developed as a model oomycete. Here, we report a bacterial artificial chromosome (BAC)-based, integrated physical map of the P. sojae genome. We constructed two BAC libraries, digested 8,681 BACs with seven restriction enzymes, end labeled the digested fragments with four dyes, and analyzed them with capillary electrophoresis. Fifteen data sets were constructed from the fingerprints, using individual dyes and all possible combinations, and were evaluated for contig assembly. In all, 257 contigs were assembled from the XhoI data set, collectively spanning approximately 132 Mb in physical length. The BAC contigs were integrated with the draft genome sequence of P. sojae by end sequencing a total of 1,440 BACs that formed a minimal tiling path. This enabled the 257 contigs of the BAC map to be merged with 207 sequence scaffolds to form an integrated map consisting of 79 superscaffolds. The map represents the first genome-wide physical map of a Phytophthora sp. and provides a valuable resource for genomics and molecular biology research in P. sojae and other Phytophthora spp. In one illustration of this value, we have placed the 350 members of a superfamily of putative pathogenicity effector genes onto the map, revealing extensive clustering of these genes.     

A BAC-based physical map of the apple genome

Genome-wide physical mapping is an essential step toward investigating the genetic basis of complex traits as well as pursuing genomics research of virtually all plant and animal species. We have constructed a physical map of the apple genome from a total of 74,281 BAC clones representing approximately 10.5x haploid genome equivalents. The physical map consists of 2702 contigs, and it is estimated to span approximately 927 Mb in physical length. The reliability of contig assembly was evaluated by several methods, including assembling contigs using variable stringencies, assembling contigs using fingerprints from individual libraries, checking consensus maps of contigs, and using DNA markers. Altogether, the results demonstrated that the contigs were properly assembled. The apple genome-wide BAC-based physical map represents the first draft genome sequence not only for any member of the large Rosaceae family, but also for all tree species. This map will play a critical role in advanced genomics research for apple and other tree species, including marker development in targeted chromosome regions, fine-mapping and isolation of genes/QTL, conducting comparative genomics analyses of plant chromosomes, and large-scale genomics sequencing.     

Construction and characterization of a papaya BAC library as a foundation for molecular dissection of a tree-fruit genome

A bacterial artificial chromosome (BAC) library was constructed from high-molecular-weight DNA isolated from young leaves of papaya (Carica papaya L.). This BAC library consists of 39168 clones from two separate ligation reactions. The average insert size of the library is 132 kb; 96.5% of the 18700 clones from the first ligation contained inserts that averaged 86 kb in size, 95.7% of the 20468 clones from the second ligation contained inserts that averaged 174 kb in size. Two sorghum chloroplast probes hybridized separately to the library and revealed a total of 504 chloroplast clones or 1.4% of the library. The entire BAC library was estimated to provide 13.7× papaya-genome equivalents, excluding the false-positive and chloroplast clones. High-density filters were made containing 94% or 36864 clones of the library with 12.7× papaya-genome equivalents. Eleven papaya-cDNA and ten Arabidopsis-cDNA probes detected an average of 22.8 BACs per probe in the library. Because of its relatively small genome (372 Mbp/1 C) and its ability to produce ripe fruit 9 to 15 months after planting, papaya shows promise as a model plant for studying genes that affect fruiting characters. A rapid approach to locating fruit-controlling genes will be to assemble a physical map based on BAC contigs to which ESTs have hybridized. A physical map of the papaya genome will significantly enhance our capacity to clone and manipulate genes of economic importance. Received: 11 April 2000 / Accepted: 28 July 2000     

Construction and characterization of a gridded cattle BAC library

A bovine genomic large-insert bacterial artificial chromosome (BAC) library has been constructed from leukocytes of a Holstein-Friesian male. Size fractionated DpnII-digested genomic DNA was ligated to the dephosphorylated BamH1 ends of a pBACe3.6 vector. Approximately 8.3 x 10(4) individual BAC clones were picked into 384-well plates. Two-hundred and sixty-seven randomly chosen clones were characterized by pulsed-field gel electrophoresis (PFGE). The average insert size was 104 kb with a frequency of clones without inserts of 5.5%. Thirty-four BAC clones were mapped by fluorescence in situ hybridization (FISH) to cattle chromosomes. Three showed signals at more than one location, one of them on the centromeric regions of all autosomes, indicating that the clone contains centromeric repeats. A subset of these BAC clones was used for the development of sequence tagged sites. Both subcloning and direct sequencing of the BACs were used for generating sequence tagged site information. The clones from the library were gridded onto high-density membranes, and PCR superpools were produced from the same set of clones. Membranes and superpools are available through the Resource Centre of the German Human Genome Project in Berlin (http:// www.rzpd.de).     

The development of BAC-end sequence-based microsatellite markers and placement in the physical and genetic maps of soybean

The composite map of soybean shared among Soybase, LIS and SoyGD (March 2006) contained 3,073 DNA markers in the "Locus" class. Among the markers were 1,019 class I microsatellite markers with 2-3 bp simple sequence repeats (SSRs) of >10 iterations (BARC-SSR markers). However, there were few class II SSRs (2-5 bp repeats with <10 iterations; mostly SIUC-Satt markers). The aims here were to increase the number of classes I and II SSR markers and to integrate bacterial artificial chromosome (BAC) clones onto the soybean physical map using the markers. Used was 10 Mb of BAC-end sequence (BES) derived from 13,473 reads from 7,050 clones constituting minimum tile path 2 of the soybean physical map ( http://www.soybeangenome.siu.edu ; SoyGD). Identified were 1,053 1-6 bp motif, repeat sequences, 333 from class I (>10 repeats) and 720 from class II (<10 repeats). Potential markers were shown on the MTP_SSR track at Gbrowse. Primers were designed as 20-24 bp oligomers that had Tm of 55 +/- 1 C that would generate 100-500 bp amplicons. About 853 useful primer pairs were established. Motifs were not randomly distributed with biases toward AT rich motifs. Strong biases against the GC motif and all tetra-nucleotide repeats were found. The markers discovered were useful. Among the first 135 targeted for use in genetic map improvement about 60% of class II markers and 75% of class I markers were polymorphic among on the parents of four recombinant inbred line (RIL) populations. Many of the BES-based SSRs were located on the soybean genetic map in regions with few BARC-SSR markers. Therefore, BES-based SSRs represent useful tools for genetic map development in soybean. New members of a consortium to map the markers in additional populations are invited.     

Microsatellite discovery from BAC end sequences and genetic mapping to anchor the soybean physical and genetic maps.

Whole-genome sequencing of the soybean (Glycine max (L.) Merr. 'Williams 82') has made it important to integrate its physical and genetic maps. To facilitate this integration of maps, we screened 3290 microsatellites (SSRs) identified from BAC end sequences of clones comprising the 'Williams 82' physical map. SSRs were screened against 3 mapping populations. We found the AAT and ACT motifs produced the greatest frequency of length polymorphisms, ranging from 17.2% to 32.3% and from 11.8% to 33.3%, respectively. Other useful motifs include the dinucleotide repeats AG, AT, and AG, with frequency of length polymorphisms ranging from 11.2% to 18.4% (AT), 12.4% to 20.6% (AG), and 11.3% to 16.4% (GT). Repeat lengths less than 16 bp were generally less useful than repeat lengths of 40-60 bp. Two hundred and sixty-five SSRs were genetically mapped in at least one population. Of the 265 mapped SSRs, 60 came from BAC singletons not yet placed into contigs of the physical map. One hundred and ten originated in BACs located in contigs for which no genetic map location was previously known. Ninety-five SSRs came from BACs within contigs for which one or more other BACs had already been mapped. For these fingerprinted contigs (FPC) a high percentage of the mapped markers showed inconsistent map locations. A strategy is introduced by which physical and genetic map inconsistencies can be resolved using the preliminary 4x assembly of the whole genome sequence of soybean.     

A physical and genetic map of the Stigmatella aurantiaca DW4/3.1 chromosome

A physical map of the myxobacterium Stigmatella aurantiaca DW4/3.1 chromosome was constructed by pulsed-field gel (PFG) long-range mapping. One-and two-dimensional pulsed-field gel analyses were used together with reciprocal double-restriction, cross-hybridization and hybridization fingerprint analysis. These PFG results were confirmed by Smith-Birnstiel analysis, by Southern hybridization using linking clones and clones of a λ genomic library for the determination of adjacent restriction fragments and by transposon insertion mapping using defined genomic sequences for hybridization. It was thus possible to construct a circular restriction map of the single 9.35 Mbp chromosome of S. aurantiaca based on the endonucleases Asel and Spel. Genetic loci as well as the replication origin were located on the physical map by Southern hybridization using heterologous (derived from Myxococcus xanthus, Escherichia coli and Streptomyces lividans) and homologous probes that are mainly involved in development and ceil motility.     

An integrated map of Arabidopsis thaliana for functional analysis of its genome sequence.

The genome of the model plant species Arabidopsis thaliana has recently been sequenced. To accelerate its current genome research, we developed a whole-genome, BAC/BIBAC-based, integrated physical, genetic, and sequence map of the A. thaliana ecotype Columbia. This new map was constructed from the clones of a new plant-transformation-competent BIBAC library and is integrated with the existing sequence map. The clones were restriction fingerprinted by DNA sequencing gel-based electrophoresis, assembled into contigs, and anchored to an existing genetic map. The map consists of 194 BAC/BIBAC contigs, spanning 126 Mb of the 130-Mb Arabidopsis genome. A total of 120 contigs, spanning 114 Mb, were anchored to the chromosomes of Arabidopsis. Accuracy of the integrated map was verified using the existing physical and sequence maps and numerous DNA markers. Integration of the new map with the sequence map has enabled gap closure of the sequence map and will facilitate functional analysis of the genome sequence. The method used here has been demonstrated to be sufficient for whole-genome physical mapping from large-insert random bacterial clones and thus is applicable to rapid development of whole-genome physical maps for other species.     

A BAC-based physical map of Zhikong scallop (Chlamys farreri Jones et Preston)

Zhikong scallop (Chlamys farreri) is one of the most economically important aquaculture species in China. Physical maps are crucial tools for genome sequencing, gene mapping and cloning, genetic improvement and selective breeding. In this study, we have developed a genome-wide, BAC-based physical map for the species. A total of 81,408 clones from two BAC libraries of the scallop were fingerprinted using an ABI 3130xl Genetic Analyzer and a fingerprinting kit developed in our laboratory. After data processing, 63,641 (～5.8× genome coverage) fingerprints were validated and used in the physical map assembly. A total of 3,696 contigs were assembled for the physical map. Each contig contained an average of 10.0 clones, with an average physical size of 490 kb. The combined total physical size of all contigs was 1.81 Gb, equivalent to approximately 1.5 fold of the scallop haploid genome. A total of 10,587 BAC end sequences (BESs) and 167 markers were integrated into the physical map. We evaluated the physical map by overgo hybridization, BAC-FISH (fluorescence in situ hybridization), contig BAC pool screening and source BAC library screening. The results have provided evidence of the high reliability of the contig physical map. This is the first physical map in mollusc; therefore, it provides an important platform for advanced research of genomics and genetics, and mapping of genes and QTL of economical importance, thus facilitating the genetic improvement and selective breeding of the scallop and other marine molluscs.     

A 10-megabase physical map of human Xp21, including the Duchenne muscular dystrophy gene

Using pulsed-field gel electrophoresis and 12 Xp21-derived DNA probes, we have constructed a continuous restriction map spanning more than 4 million base pairs (4 Mbp), including the Duchenne muscular dystrophy gene of more than 2 Mbp. This detailed map is part of a less detailed map spanning 10 Mbp, also spanning the genes for glycerol kinase and congenital adrenal hypoplasia, constructed under electrophoresis conditions which separated DNA fragments in the range 200 to 4000 kbp. DNA from three different tissues was analyzed, and differential methylation was observed.     

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 ± 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa , version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.     

Physical mapping, BAC-end sequence analysis, and marker tagging of the soilborne nematicidal bacterium, Pseudomonas synxantha BG33R

A bacterial artificial chromosome (BAC) library was constructed for the genome of the rhizosphere-inhabiting fluorescent pseudomonad Pseudomonas synxantha BG33R. Three thousand BAC clones with an average insert size of 140 kbp and representing a 70-fold genomic coverage were generated and arrayed onto nylon membranes. EcoRI fingerprint analysis of 986 BAC clones generated 23 contigs and 75 singletons. Hybridization analysis allowed us to order the 23 contigs and condense them into a single contig, yielding an estimated genome size of 5.1 Mb for P. synxantha BG33R. A minimum-tile path of 47 BACs was generated and end-sequenced. The genetic loci involved in ring nematode egg-kill factor production in BG33R Tn5 mutants, 246 (vgrG homolog), 1122 (sensor kinase homolog), 1233 (UDP-galactose epimerase homolog), 1397 (ferrisiderophore receptor homolog), and 1917 (ribosomal subunit protein homolog), have been mapped onto the minimum-tile BAC library. Two of the genetic regions that flank Tn5 insertions in BG33R egg-kill-negative mutants 1233 and 1397 are separated by a single BAC clone. Fragments isolated by ligation-mediated PCR of the Tn5 mutagenized regions of 29 randomly selected, non-egg-kill-related, insertion mutants have been anchored onto the ordered physical map of P. synxantha.