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.     

SPC-P1: a pathogenicity-associated prophage of Salmonella paratyphi C

Background

The presence of closely related genomes in polyploid species makes the assembly of total genomic sequence from shotgun sequence reads produced by the current sequencing platforms exceedingly difficult, if not impossible. Genomes of polyploid species could be sequenced following the ordered-clone sequencing approach employing contigs of bacterial artificial chromosome (BAC) clones and BAC-based physical maps. Although BAC contigs can currently be constructed for virtually any diploid organism with the SNaPshot high-information-content-fingerprinting (HICF) technology, it is currently unknown if this is also true for polyploid species. It is possible that BAC clones from orthologous regions of homoeologous chromosomes would share numerous restriction fragments and be therefore included into common contigs. Because of this and other concerns, physical mapping utilizing the SNaPshot HICF of BAC libraries of polyploid species has not been pursued and the possibility of doing so has not been assessed. The sole exception has been in common wheat, an allohexaploid in which it is possible to construct single-chromosome or single-chromosome-arm BAC libraries from DNA of flow-sorted chromosomes and bypass the obstacles created by polyploidy.

Results

The potential of the SNaPshot HICF technology for physical mapping of polyploid plants utilizing global BAC libraries was evaluated by assembling contigs of fingerprinted clones in an in silico merged BAC library composed of single-chromosome libraries of two wheat homoeologous chromosome arms, 3AS and 3DS, and complete chromosome 3B. Because the chromosome arm origin of each clone was known, it was possible to estimate the fidelity of contig assembly. On average 97.78% or more clones, depending on the library, were from a single chromosome arm. A large portion of the remaining clones was shown to be library contamination from other chromosomes, a feature that is unavoidable during the construction of single-chromosome BAC libraries.

Conclusions

The negligibly low level of incorporation of clones from homoeologous chromosome arms into a contig during contig assembly suggested that it is feasible to construct contigs and physical maps using global BAC libraries of wheat and almost certainly also of other plant polyploid species with genome sizes comparable to that of wheat. Because of the high purity of the resulting assembled contigs, they can be directly used for genome sequencing. It is currently unknown but possible that equally good BAC contigs can be also constructed for polyploid species containing smaller, more gene-rich genomes.     

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.     

Expression analysis and physical mapping of low-molecular-weight glutenin loci in hexaploid wheat (Triticum aestivum L.).

Gliadins and glutenins are storage proteins important in determining the bread-, noodle-, and pasta-making quality of wheat. Glutenins consist of HMW and LMW subunits. The Glu-A3, Glu-B3, and Glu-D3 loci on the short arms of chromosomes 1A, 1B, and 1D, respectively, are the major loci for LMW glutenins. To construct physical maps of the Glu-3 loci, a set of 24 high-density filters representing a 3.1x genome coverage hexaploid wheat BAC library was screened by hybridization using a probe made of 3 LMW glutenin sequences. After 2 rounds of hybridization, a subset of 536 BAC clones were selected and fingerprinted. Three developing seed cDNA libraries were also constructed. A total of 5000-6000 ESTs were generated from each library, assembled into contigs and searched by homology for LMW glutenin sequences. In total, 90 full-length LMW glutenin sequences were found to cluster into 8 distinct groups representing at least 21 different LMW glutenin subunits. A set of 24 pairs of PCR primers was designed from these groups and used as markers on the BAC clones. The combined fingerprinting and marker data were used to build the physical maps using FPC software. A total of 91 contigs comprising 254 clones were obtained and 282 clones remained singletons.     

The generation of ordered sets of cosmid DNA clones from human chromosome region 11p.

We describe progress in a continuing project aimed at the generation of an overlapping cosmid DNA clone map of the short arm of human chromosome 11. The automated procedures used to prepare DNA samples and the computerized data collection and recording systems are described. We also demonstrate the use of the clones as reagents for the rapid isolation of genomic DNAs containing smaller probed regions. We have isolated approximately 4700 human cosmid DNA clones from mouse/human hybrid cell lines that contain predominantly human chromosomal region 11p. Of the DNA in the cell lines, 60% is derived from this chromosomal region, and the remaining 40% is derived from regions of chromosomes 3, 19, and 20. A total of 4159 clones have been fingerprinted to identify potential overlaps, and we have developed 535 sets ("contigs"). Using random modeling, it is estimated that 65% of 11p must be contained in the analyzed cosmids. The database of clones has been used to identify single or overlapping clones from noncosmid DNA probes. Examples are presented. It is proposed that cosmid reference filters be distributed to requesting laboratories.     

Insight into Trichoderma reesei's genome content, organization and evolution revealed through BAC library characterization

Trichoderma reesei is an important industrial fungus known for its ability to efficiently secrete large quantities of protein as well as its wide variety of biomass degrading enzymes. Past research on this fungus has primarily focused on extending its protein production capabilities, leaving the structure of its 33 Mb genome essentially a mystery. To begin to address these deficiencies and further our knowledge of T. reesei's secretion and cellulolytic potential, we have created a genomic framework for this fungus. We constructed a BAC library containing 9216 clones with an average insert size of 125 kb which provides a coverage of 28 genome equivalents. BAC ends were sequenced and annotated using publicly available software which identified a number of genes not seen in previously sequenced EST datasets. Little evidence was found for repetitive sequence in T. reesei with the exception of several copies of an element with similarity to the Podospora anserina transposon, PAT. Hybridization of 34 genes involved in biomass degradation revealed five groups of co-located genes in the genome. BAC clones were fingerprinted and analyzed using fingerprinted contigs (FPC) software resulting in 334 contigs covering 28 megabases of the genome. The assembly of these FPC contigs was verified by congruence with hybridization results.     

Algorithms and software tools for ordering clone libraries: application to the mapping of the genome of Schizosaccharomyces pombe.

A complete set of software tools to aid the physical mapping of a genome has been developed and successfully applied to the genomic mapping of the fission yeast Schizosaccharomyces pombe. Two approaches were used for ordering single-copy hybridisation probes: one was based on the simulated annealing algorithm to order all probes, and another on inferring the minimum-spanning subset of the probes using a heuristic filtering procedure. Both algorithms produced almost identical maps, with minor differences in the order of repetitive probes and those having identical hybridisation patterns. A separate algorithm fitted the clones to the established probe order. Approaches for handling experimental noise and repetitive elements are discussed. In addition to these programs and the database management software, tools for visualizing and editing the data are described. The issues of combining the information from different libraries are addressed. Also, ways of handling multiple-copy probes and non-hybridisation data are discussed.     

ODS2: a multiplatform software application for creating integrated physical and genetic maps

A contig map is a physical map that shows the native order of a library of overlapping genomic clones. One common method for creating such maps involves using hybridization to detect clone overlaps. False- positive and false-negative hybridization errors, the presence of chimeric clones, and gaps in library coverage lead to ambiguity and error in the clone order. Genomes with good genetic maps, such as Neurospora crassa, provide a means for reducing ambiguities and errors when constructing contig maps if clones can be anchored with genetic markers to the genetic map. A software application called ODS2 for creating contig maps based on clone-clone hybridization data is presented. This application is also designed to exploit partial ordering information provided by anchorage of clones to a genetic map. This information, along with clone-clone hybridization data, is used by a clone ordering algorithm and is represented graphically, allowing users to interactively align physical and genetic maps. ODS2 has a graphical user interface and is implemented entirely in Java, so it runs on multiple platforms. Other features include the flexibility of storing data in a local file or relational database and the ability to create full or minimum tiling contig maps.     

A comprehensive BAC resource

The Human Genome Project has generated extensive map and sequence data for a large number of Bacterial Artificial Chromosome (BAC) clones. In order to maximize the efficient use of the data and to minimize the redundant work for the research community, The Institute for Genomic Research (TIGR) comprehensive BAC resource (cBACr) (http://www.tigr.org/tdb/BacResource/BAC_resourc e_intro. html) was built as an expansion of the TIGR human BAC ends database. This resource collects, integrates and reports the information on library, maps, sequence, annotation and functions for each human and mouse BAC. The current database contains 635 016 human BACs and 265 617 mouse BACs that were characterized by various approaches, among which 22 705 human clones and 1000 mouse clones have sequence and annotation data.     

Construction and characterization of a soybean bacterial artificial chromosome library and use of multiple complementary libraries for genome physical mapping

Two plant-transformation-competent large-insert binary clone bacterial artificial chromosome (hereafter BIBAC) libraries were previously constructed for soybean cv. Forrest, using BamHI or HindIII. However, they are not well suited for clone-based genomic sequencing due to their larger ratio of vector to insert size (27.6 kbp:125 kbp). Therefore, we developed a larger-insert bacterial artificial chromosome (BAC) library for the genotype in a smaller vector (pECBAC1), using EcoRI. The BAC library contains 38,400 clones; about 99.1% of the clones have inserts; the average insert size is 157 kbp; and the ratio of vector to insert size is much smaller (7.5 kbp:157 kbp). Colony hybridization with probes derived from several chloroplast and mitochondrial genes showed that 0.89% and 0.45% of the clones were derived from the chloroplast and mitochondrial genomes, respectively. Considering these data, the library represents 5.4 haploid genomes of soybean. The library was hybridized with six RFLP marker probes, 5S rDNA and 18S-5.8S-25S rDNA, respectively. Each RFLP marker hybridized to about six clones, and the 5S and 18S-5.8S-25S rDNA probes collectively hybridized to 402 BACs—about 1.05% of the clones in the library. The BAC library complements the existing soybean Forrest BIBAC libraries by using different restriction enzymes and vector systems. Together, the BAC and BIBAC libraries encompass 13.2 haploid genomes, providing the most comprehensive clone resource for a single soybean genotype for public genome research. We show that the BAC library has enhanced the development of the soybean whole-genome physical map and use of three complementary BAC libraries improves genome physical mapping by fingerprint analysis of most of the clones of the library. The rDNA-containing clones were also fingerprinted to evaluate the feasibility of constructing contig maps of the rDNA regions. It was found that physical maps for the rDNA regions could not be readily constructed by fingerprint analysis, using one or two restriction enzymes. Additional data to fingerprints and/or different fingerprinting methods are needed to build contig maps for such highly tandem repetitive regions and thus, the physical map of the entire soybean genome.     

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.     

Construction and utility of a human chromosome 22-specific Fosmid library

We have previously demonstrated the capability of the Fosmid vector based on Escherichia coli F-factor replicon to stably propagate cosmid-sized human genomic DNA fragments. Using the Fosmid vector, we have constructed and arrayed a 10 × human chromosome 22-specific library, partly by picking human positive clones from a total Fosmid library constructed using DNA from human-hamster hybrid cell line containing human chromosome 22, and partly by using flow-sorted chromosomal DNA. The clones and physical contig maps of the clones in the library will serve as a valuable resource for detailed analysis of the chromosome by providing reliable materials for high resolution mapping and sequencing. In order to efficiently built physical maps for the chromosomal regions of interest spanning several hundred kilobases to a megabase, it is necessary to rapidly identify subsets of the Fosmid clones from the library that cover such regions. In this report, we describe a method of using random amplification products derived from YAC clones to rapidly identify a subset of Fosmid clones that cover a specific genomic subregion.     

Physical map and gene survey of the Ochrobactrum anthropi genome using bacterial artificial chromosome contigs

Bacterial artificial chromosome (BAC) clones are effective mapping and sequencing reagents for use with a wide variety of small and large genomes. This report describes research aimed at determining the genome structure of Ochrobactrum anthropi, an opportunistic human pathogen that has potential applications in biodegradation of hazardous organic compounds. A BAC library for O. anthropi was constructed that provides a 70-fold genome coverage based on an estimated genome size of 4.8 Mb. The library contains 3072 clones with an average insert size of 112 kb. High-density colony filters of the library were made, and a physical map of the genome was constructed using a hybridization without replacement strategy. In addition, 1536 BAC clones were fingerprinted with HindIII and analyzed using IMAGE and Fingerprint Contig software (FPC, Sanger Centre, U.K.). The FPC results supported the hybridization data, resulting in the formation of two major contigs representing the two major replicons of the O. anthropi genome. After determining a reduced tiling path, 138 BAC ends from the reduced tile were sequenced for a preliminary gene survey. A search of the public databases with the BLASTX algorithm resulted in 77 strong hits (E-value < 0.001), of which 89% showed similarity to a wide variety of prokaryotic genes. These results provide a contig-based physical map to assist the cloning of important genomic regions and the potential sequencing of the O. anthropi genome.     

Carbon-fate maps for metabolic reactions

MOTIVATION: Stable isotope labeling of small-molecule metabolites (e.g. (13)C-labeling of glucose) is a powerful tool for characterizing pathways and reaction fluxes in a metabolic network. Analysis of isotope labeling patterns requires knowledge of the fates of individual atoms and moieties in reactions, which can be difficult to collect in a useful form when considering a large number of enzymatic reactions. RESULTS: We report carbon-fate maps for 4605 enzyme-catalyzed reactions documented in the KEGG database. Every fate map has been manually checked for consistency with known reaction mechanisms. A map includes a standardized structure-based identifier for each reactant (namely, an InChI string); indices for carbon atoms that are uniquely derived from the metabolite identifiers; structural data, including an identification of homotopic and prochiral carbon atoms; and a bijective map relating the corresponding carbon atoms in substrates and products. Fate maps are defined using the BioNetGen language (BNGL), a formal model-specification language, which allows a set of maps to be automatically translated into isotopomer mass-balance equations. AVAILABILITY: The carbon-fate maps and software for visualizing the maps are freely available (http://cellsignaling.lanl.gov/FateMaps/).     

Interleukin-1 receptor cluster: gene organization of IL1R2, IL1R1, IL1RL2 (IL-1Rrp2), IL1RL1 (T1/ST2), and IL18R1 (IL-1Rrp) on human chromosome 2q

The family of interleukin-1 receptor-like genes currently has six known members. We have constructed a contig of 10 overlapping human PAC clones that covers 530 kb and includes five of the six family members. The termini of the contig were mapped to the interval between D2S373 and D2S176 (chromosome 2q12) by radiation hybrid mapping. The contig contains the genes (cen --> tel), in the order given, for the type II interleukin-1 (IL-1) receptor (IL1R2), the type I IL-1 receptor (IL1R1), the IL-1 receptor-related protein 2 (IL1RL2), T1/ST2/fit-1 (IL1RL1), and the IL-1 receptor-related protein 1, which has recently been shown to be a component of the IL-18 receptor (IL18R1). We show that all the genes are transcribed in the same direction, with IL1R2 being transcribed toward the cluster. The only known family member that is absent from the human contig is the IL-1 receptor accessory protein gene (IL1RAP), which maps to 3q28.     

Animal QTLdb: beyond a repository

Over the past ten years there have been a large number of publications that have described hundreds of quantitative trait loci (QTL) in livestock species. To facilitate the comparison of QTL results across experiments, the Animal QTL database (QTLdb) was developed to house all published QTL information as a public repository. The QTLdb was originally developed to serve the porcine genomics community (previously known as PigQTLdb). We have further developed the QTLdb to house QTL data from multiple species, including but not limited to cattle, chickens, and pigs. In addition, tools have been developed to allow QTL map alignments against consensus linkage maps, radiation hybrid (RH) maps, BAC fingerprinted contig (FPC) maps, single nucleotide polymorphism (SNP) location maps, and human maps. In addition, we have expanded the capabilities of the database such that research tools were developed where “private” preliminary QTL data could be entered and compared against all public data. This allows researchers to visualize data before publication and compare it with published results to aid in data interpretation. To serve this purpose, the database curator/editor tools also include functions that allow registered users to enter their own QTL data, make use of the QTLdb tools for data analysis, and use the QTLdb as a publishing tool (URL: ).     

Comparative mapping of human Chromosome 2 identifies segments of conserved synteny near the bovine mh locus

The ``double-muscling' (mh) locus has been localized to an interval between the centromere and the microsatellite marker TGLA44 on bovine Chromosome (Chr) 2 (BTA2). We identified segments of conserved synteny that correspond to this region of BTA2 by assigning large genomic clones containing bovine homologs of seven genes from the long arm of human Chr 2 (HSA2q). Polymorphic markers developed from these clones integrated the physical and linkage maps of BTA2 from 2q12 to 2q44 and extended genetic coverage towards the centromere. This comparative analysis suggests the mh locus resides on HSA2q near both the protein C and collagen type III alpha-1 genes. Overall, our data reveal a complex rearrangement of gene order between BTA2q12-44 and HSA2q14-37 that underscores the need to establish boundaries of conserved synteny when applying comparative mapping information to identify genes or traits of interest. Received: 3 March 1997 / Accepted: 12 May 1997     

Genome physical mapping from large-insert clones by fingerprint analysis with capillary electrophoresis: a robust physical map of Penicillium chrysogenum

Physical mapping with large-insert clones is becoming an active area of genomics research, and capillary electrophoresis (CE) promises to revolutionize the physical mapping technology. Here, we demonstrate the utility of the CE technology for genome physical mapping with large-insert clones by constructing a robust, binary bacterial artificial chromosome (BIBAC)-based physical map of Penicillium chrysogenum. We fingerprinted 23.1× coverage BIBAC clones with five restriction enzymes and the SNaPshot kit containing four fluorescent-ddNTPs using the CE technology, and explored various strategies to construct quality physical maps. It was shown that the fingerprints labeled with one or two colors, resulting in 40–70 bands per clone, were assembled into much better quality maps than those labeled with three or four colors. The selection of fingerprinting enzymes was crucial to quality map construction. From the dataset labeled with ddTTP–dROX, we assembled a physical map for P.chrysogenum, with 2–3 contigs per chromosome and anchored the map to its chromosomes. This map represents the first physical map constructed using the CE technology, thus providing not only a platform for genomic studies of the penicillin-producing species, but also strategies for efficient use of the CE technology for genome physical mapping of plants, animals and microbes.     

A Composite Linkage Map from Three Crosses Between Commercial Clones of Cacao, Theobroma cacao L.

In this paper, we report the construction of the first composite map of cacao from linkage data of one F₂ and two F₁ mapping populations with a high number of codominant markers in common. The combination of linkage information from all three maps results in the currently most precise estimates of marker locations and distances between markers, especially in densely marked areas. JoinMap®V4 software was used for all marker quality assessment and mapping. Individual (sub-composite) maps and the composite map contained 10 major linkage groups, corresponding to the number of cacao chromosomes. Homogeneity of marker placement was very high among sub-composite maps, the composite map, and the designated “reference” map. Care was exercised in the re-creation of sub-composite maps and the composite map to include only markers with acceptable mapping quality parameters. The composite map places more markers with higher precision than any individual map. This research clearly demonstrates for the first time a very high level of marker homogeneity among commercial cacao clones compared to other species. The observed homogeneity between different maps, including the composite one, is probably due to a narrow genetic base of commercial cacao clones. Markers linked to identified quantitative trait loci (QTLs) are more likely to retain linkage in other commercial clones, rendering the QTLs in cacao potentially more stable than in other species.     

Differences between the radiation hydrid and genetic linkage maps of bovine Chromosome 5 resolved with a quasi-phylogenetic method of analysis

Two major differences were detected in gene order between the radiation hybrid map and the genetic linkage map of bovine Chromosome (Chr) 5, and these were resolved by analyzing the raw radiation hybrid data by a quasi-phylogenetic method. Seventeen loci were typed on the new cattle whole genome radiation hybrid panel. Most of these loci are framework loci and include AGLA293, BM315, BM6026, BP1, BZRP, CD9, CSSM22, CSSM34, CYP2D@, ETH2, ETH10, ETH152, IGF1, LALBA, SLC2A3, SYT1, and TPI1. BP1 was found to be closer to the centromere than either BM6026 or SYT1 with two standard computer software packages for analyzing radiation hybrid panel data. This is inconsistent with any of the genetic linkage maps as well as their consensus. CYP2D@ was placed between ETH2 and BZRP, and this is also inconsistent with the genetic linkage maps, since CYP2D@ should be the most telomeric of the loci tested in this study. Resolution was reached by analyzing the raw radiation hybrid data for clones that bind some but not all of the loci, and the binding pattern was more consistent with the linkage maps and less consistent with the software-generated radiation hybrid map. The comparative mapping data confirm the relative inversion of gene order of SYT1 compared with humans and mice. A non-polymorphic fragment for CD9 indicates the conservation of gene order for three loci located on human Chr 12p. The genes of bovine Chr 5 conserved on human Chr 12p are located separately from the genes conserved on human Chr 12q. It is recommended that the raw data for radiation hybrid maps be made publicly available so that conflicts in gene order can be evaluated explicity. Received: 19 February 1999 / Accepted: 3 January 2000