The Use of Simulated Annealing in Chromosome Reconstruction Experiments Based on Binary Scoring

We present a method of combinatorial optimization, simulated annealing, to order clones in a library with respect to their position along a chromosome. This ordering method relies on scoring each clone for the presence or absence of specific target sequences, thereby assigning a digital signature to each clone. Specifically, we consider the hybridization of oligonucleotide probes to a clone to constitute the signature. In that the degree of clonal overlap is reflected in the similarity of their signatures, it is possible to construct maps based on the minimization of the differences in signatures across a reconstructed chromosome. Our simulations show that with as few as 30 probes and a clonal density of 4.5 genome equivalents, it is possible to assemble a small eukaryotic chromosome into 33 contiguous blocks of clones (contigs). With higher clonal densities and more probes, this number can be reduced to less than 5 contigs per chromosome.     

On the design of genome mapping experiments using short synthetic oligonucleotides.

The DNA of an organism can be digested into smaller fragments, stored individually as clones in phage, for example, to create a clone library, and retrieved later, when needed. The original ordering of fragments is lost in the process of creating the library. Hence, it is important to be able to place clones in order according to their position along chromosome(s), and this process is referred to as "in vitro reconstruction" or "contig mapping" of an organismal genome. Clones in the phage library can be assigned binary call numbers by scoring each clone for hybridization (0 or 1) with a battery of short manufactured DNA sequences called synthetic oligonucleotides or with restriction enzyme digests of each clone. Those clones with similar call numbers are placed close together in the ordered library. We address the design question of how many clones and probes to use to carry out in vitro reconstruction of an organism's chromosomes. This physical mapping problem is placed in the context of coverage problems in geometrical probability. Various statistics are developed to summarize how an ordered library covers a chromosome, the extent of clone overlap, and the similarity between clone call numbers. Several tests for whether clones overlap are given, together with their power properties. A simulation study is used to determine how robust some of the tests for clone overlap are to model violations. Tables are presented for researchers to choose the number of clones and probes on the basis of both power and technical considerations surrounding the hybridization experiments.     

A combinatorial approach for fast,high-resolution mapping

High-resolution physical maps can be used as a scaffold for several subsequent studies, such as sequencing projects and positional cloning of disease genes and genetic elements that regulate gene expression. Here we describe a method for fast, high-resolution physical mapping on stretched DNA molecules, based on a combinatorial multi-FISH approach. Fluorescent labels are assigned to a binary code and probes are identified by a binary tag according to their labeling. To validate the approach, we have mapped eight probes covering a region of about 300 kb on human chromosome 11 with three hybridization assays. This approach enables one to determine the structural organization of a large region by means of the order of its clones, without ambiguities. The structure established in a control cell constitutes a reference for further studies, to detect rearrangements displayed by disease cells and to find differences shown by different cell types and organisms.     

Development of one set of chromosome-specific microsatellite-containing BACs and their physical mapping in Gossypium hirsutum L.

Fluorescence in situ hybridization (FISH), using bacterial artificial chromosome (BAC) clone as probe, is a reliable cytological technique for chromosome identification. It has been used in many plants, especially in those containing numerous small chromosomes. We previously developed eight chromosome-specific BAC clones from tetraploid cotton, which were used as excellent cytological markers for chromosomes identification. Here, we isolated the other chromosome-specific BAC clones to make a complete set for the identification of all 26 chromosome-pairs by this technology in tetraploid cotton (Gossypium hirsutum L.). This set of BAC markers was demonstrated to be useful to assign each chromosome to a genetic linkage group unambiguously. In addition, these BAC clones also served as convenient and reliable landmarks for establishing physical linkage with unknown targeted sequences. Moreover, one BAC containing an EST, with high sequence similarity to a G. hirsutum ethylene-responsive element-binding factor was located physically on the long arm of chromosome A7 with the help of a chromosome-A7-specific BAC FISH marker. Comparative analysis of physical marker positions in the chromosomes by BAC-FISH and genetic linkage maps demonstrated that most of the 26 BAC clones were localized close to or at the ends of their respective chromosomes, and indicated that the recombination active regions of cotton chromosomes are primarily located in the distal regions. This technology also enables us to make associations between chromosomes and their genetic linkage groups and re-assign each chromosome according to the corresponding genetic linkage group. This BAC clones and BAC-FISH technology will be useful for us to evaluate grossly the degree to which a linkage map provides adequate coverage for developing a saturated genetic map, and provides a powerful resource for cotton genomic researches.     

Construction of a binary BAC library for an apomictic monosomic addition line of<Emphasis Type="Italic"> Beta corolliflora</Emphasis> in sugar beet and identification of the clones derived from the alien chromosome

A plant-transformation-competent binary BAC library was constructed from the genomic DNA of the chromosome 9 monosomic addition line of Beta corolliflora Zoss. in sugar beet (B. vulgaris. L). This monosomic addition line (designated M14) is characterized by diplosporic reproduction caused by the alien chromosome carrying the gene(s) responsible for diplospory. The library consists of 49,920 clones with an average insert size of 127 kb, representing approximately 7.5 haploid genome equivalents and providing a greater than 99% probability of isolating a single-copy DNA sequence from the library. To develop the scaffold of a physical map for the alien chromosome, B. corolliflora genome-specific dispersed repetitive DNA sequences were used as probes to isolate BAC clones derived from the alien chromosome in the library. A total of 2,365 positive clones were obtained and arrayed into a sublibrary specific for B. corolliflora chromosome 9 (designated bcBAC-IX). The bcBAC-IX sublibrary was further screened with a subtractive cDNA pool generated from the ovules of M14 and the floral buds of B. vulgaris by the suppression subtractive hybridization method. One hundred and three positive binary BACs were obtained, which potentially contain the genes of the alien chromosome specifically expressed during the ovule and embryo development of M14, and may be associated with apomictic reproduction. Thus, these binary BAC clones will be useful for identification of the genes for apomixis by genetic transformation.Communicated by H. C. Becker     

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.     

Systematic generation of sequence-tagged sites for physical mapping of human chromosomes: application to the mapping of human chromosome 7 using yeast artificial chromosomes.

Basic to the development of long-range physical maps of DNA are the detection and localization of landmarks within recombinant clones. Sequence-tagged sites (STSs), which are short stretches of DNA that can be specifically detected by the polymerase chain reaction (PCR), can be used as such landmarks. Our interest is to construct physical maps of whole human chromosomes by localizing STSs within yeast artificial chromosome (YAC) clones. Here we describe a generalized strategy for the systematic generation of large numbers of STSs specific for human chromosome 7. These STSs can be detected by PCR assays developed following the sequencing of anonymous pieces of chromosome 7 DNA, which was derived from flow-sorted chromosomes or from lambda clones made from DNA of a human-hamster hybrid cell line. Our approach for STS generation is tailored for the development of PCR assays capable of screening a large YAC library. In this study, we report the generation of 100 new STSs specific to human chromosome 7.     

Identification of genetically linked RGAs by BAC screening in maize and implications for gene cloning,mapping and MAS

The resistance gene analogue (RGA) pic19 in maize, a candidate for sugarcane mosaic virus (SCMV) resistance gene (R gene) Scmv1, was used to screen a maize BAC library to identify homologous sequences in the maize genome and to investigate their genomic organisation. Fifteen positive BAC clones were identified and could be classified into five physically independent contigs consisting of overlapping clones. Genetic mapping clustered three contigs into the same genomic region as Scmv1 on chromosome 6S. The two remaining contigs mapped to the same region as a QTL for SCMV resistance on chromosome 1. Thus, RGAs mapping to a target region can be successfully used to identify further-linked candidate sequences. The pic19 homologous sequences of these clones revealed a sequence similarity of 94-98% on the nucleotide level. The high sequence similarity reveals potential problems for the use of RGAs as molecular markers. Their application in marker-assisted selection (MAS) and the construction of high-density genetic maps is complicated by the existence of closely linked homologues resulting in 'ghost' marker loci analogous to 'ghost' QTLs. Therefore, implementation of genomic library screening, including genetic mapping of potential homologues, seems necessary for the safe application of RGA markers in MAS and gene isolation.     

A high-resolution cytogenetic map of 168 cosmid DNA markers for human chromosome 11.

We have constructed a high-resolution cytogenetic map with 168 DNA markers, including 90 RFLP markers for human chromosome 11. The cosmid clones were mapped by fluorescence in situ suppression hybridization, in which discrete fluorescent signals can be detected directly on prometaphase R-banded chromosomes. Although these cosmid clones were distributed throughout the chromosome, they had some tendency to localize in the regions of R-positive band, such as 11p15, 11p11.2, 11q13, 11q23, and 11q25. Since these regions of chromosome 11 are considered to contain genes responsible for certain genetic diseases, cancer breakpoints involved in chromosome rearrangements, and tumor-suppressor genes, this high-resolution cytogenetic map will contribute to the molecular characterization of such genes. This map will also provide many landmarks essential for construction of the complete physical map with contigs of cosmid and YAC clones.     

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.     

Characterization and application of soybean YACs to molecular cytogenetics

Yeast artificial chromosomes (YACs) are widely used in the physical analysis of complex genomes. In addition to their value in chromosome walking for map-based cloning, YACs represent excellent probes for chromosome mapping using fluorescence in situ hybridization (FISH). We have screened such a library for low-copy-number clones by hybridization to total genomic DNA. Four clones were chosen for chromosome tagging based upon their low or moderate signal. By using degenerate oligonucleotide-primed PCR (DOP-PCR), we were able to use relatively small amounts of soybean YAC DNA, isolated directly by preparative pulsed-field gel electrophoresis, as FISH probes for both metaphase chromosome spreads and interphase nuclei. FISH chromosomal analysis using the three of the clones as probes resulted in relatively simple hybridization patterns consistent with a single homologous locus or two homoeologous loci. The fourth YAC probe resulted in a diffuse hybridization pattern with signal on all metaphase chromosomes. We conclude that YACs represent a valuable source of probes for chromosomal analysis in soybean.     

An integrated map of Oryza sativa L. chromosome 5

The developments of molecular marker-based genetic linkage maps are now routine. Physical maps based on contigs of large insert genomic clones have been established in several plant species. However, integration of genetic, physical, and cytological maps is still a challenge for most plant species. Here we present an integrated map of rice (Oryza sativa L.) chromosome 5, developed by fluorescence in situ hybridization mapping of 18 bacterial artificial chromosome (BAC) clones or PI-derived artificial chromosome (PAC) clones on meiotic pachytene chromosomes. Each BAC/PAC clone was anchored by a restriction fragment length polymorphism marker mapped to the rice genetic linkage map. This molecular cytogenetic map shows the genetic recombination and sequence information of a physical map, correlated to the cytological features of rice chromosome 5. Detailed comparisons of the distances between markers on genetic, cytological, and physical maps, revealed the distributions of recombination events and molecular organization of the chromosomal features of rice chromosome 5 at the pachytene stage. Discordance of distances between the markers was found among the different maps. Our results revealed that neither the recombination events nor the degree of chromatin condensation were evenly distributed along the entire length of chromosome 5. Detailed comparisons of the correlative positions of markers on the genetic, cytological, and physical maps of rice chromosome 5 provide insight into the molecular architecture of rice chromosome 5, in relation to its cytological features and recombination events on the genetic map. The prospective applications of such an integrated cytogenetic map are discussed.     

Comparative genomic analysis of sequences sampled from a small region on soybean (Glycine max) molecular linkage group G.

Eight DNA markers spanning an interval of approximately 10 centimorgans (cM) on soybean (Glycine max) molecular linkage group G (MLG-G) were used to identify bacterial artificial chromosome (BAC) clones. Twenty-eight BAC clones in eight distinct contiguous groups (contigs) were isolated from this genome region, along with 59 BAC clones on 17 contigs homoeologous to those on MLG-G. BAC clones in four of the MLG-G contigs were also digested to produce subclones and detailed physical maps. All of the BAC-ends were sequenced, as were the subclones, to estimate proportions in different sequence categories, compare similarities among homoeologs, and explore microsynteny with Arabidopsis. Homoeologous BAC contigs were enriched in repetitive sequences compared with those on MLG-G or the soybean genome as a whole. Fingerprint and cross-hybridization comparisons between MLG-G and homoeologous contigs revealed cases of highly similar physical organization between soybean duplicates, as did DNA sequence comparisons. Twenty-seven out of 78 total sequences on soybean MLG-G showed significant similarity to Arabidopsis. The homologs mapped to six compact genome segments in Arabidopsis, with the longest containing seven homologs spanning two million base pairs. These results extend previous observations of large-scale duplication and selective gene loss in Arabidopsis, suggesting that networks of conserved synteny between Arabidopsis and other angiosperm families can stretch over long physical distances.     

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.     

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.     

Genetic analysis of the chromosome of alkaliphilic Bacillus halodurans C-125

Seventeen Sse8387I linking clones isolated from the chromosome of Bacillus halodurans C-125 for the purpose of constructing a physical map were sequenced and analyzed by comparison with the BSORF database and the nonredundant protein databank. The orientations of Sse8387I or AscI linking clones serving to join adjacent fragments were determined by southern blot analysis using specific DNA probes. One-third of the open reading frames (ORFs) identified in the Sse8387I linking clones showed no significant similarity to any protein so far reported. The ORFs showing significant similarities to those of Bacillus subtilis were mapped in the chromosome of strain C-125, and the locations of the putative genes on the map were not well conserved between B. halodurans C-125 and B. subtilis. Received: March 26, 1999 / Accepted: April 27, 1999     

Identification of dynein heavy chain genes expressed in human and mouse testis: chromosomal localization of an axonemal dynein gene

Dynein heavy chains are involved in microtubule-dependent transport processes. While cytoplasmic dyneins are involved in chromosome or vesicle movement, axonemal dyneins are essential for motility of cilia and flagella. Here we report the isolation of dynein heavy chain (DHC)-like sequences in man and mouse. Using polymerase chain reaction and reverse-transcribed human and mouse testis RNA cDNA fragments encoding the conserved ATP binding region of dynein heavy chains were amplified. We identified 11 different mouse and eight human dynein-like sequences in testis which show high similarity to known dyneins of different species such as rat, sea urchin or green algae. Sequence similarities suggest that two of the mouse clones and one human clone encode putative cytoplasmic dynein heavy chains, whereas the other sequences show higher similarity to axonemal dyneins. Two of nine axonemal dynein isoforms identified in the mouse testis are more closely related to known outer arm dyneins, while seven clones seem to belong to the inner arm dynein group. Of the isolated human isoforms three clones were classified as outer arm and four clones as inner arm dynein heavy chains. Each of the DHC cDNAs corresponds to an individual gene as determined by Southern blot experiments. The alignment of the deduced protein sequences between human (HDHC) and mouse (MDHC) dynein fragments reveals higher similarity between single human and mouse sequences than between two sequences of the same species. Human and mouse cDNA fragments were used to isolate genomic clones. Two of these clones, gHDHC7 and gMDHC7, are homologous genes encoding axonemal inner arm dyneins. While the human clone is assigned to 3p21, the mouse gene maps to chromosome 14.     

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.     

A chromosome 5-specific repetitive DNA sequence in rice (Oryza sativa L)

Repetitive DNA sequences in the rice genome comprise more than half of the nuclear DNA. The isolation and characterization of these repetitive DNA sequences should lead to a better understanding of rice chromosome structure and genome organization. We report here the characterization and chromosome localization of a chromosome 5-specific repetitive DNA sequence. This repetitive DNA sequence was estimated to have at least 900 copies. DNA sequence analysis of three genomic clones which contain the repeat unit indicated that the DNA sequences have two sub-repeat units of 37 bp and 19 bp, connected by 30-to 90-bp short sequences with high similarity. RFLP mapping and physical mapping by fluorescence in situ hybridization (FISH) indicated that almost all copies of the repetitive DNA sequence are located in the centromeric heterochromatic region of the long arm of chromosome 5. The strategy for cloning such repetitive DNA sequences and their uses in rice genome research are discussed.