Sequence-tagged site (STS) content mapping of human chromosomes: theoretical considerations and early experiences

Human cells involved in atherosclerosis have a sex

The influence of sex has been largely described in cardiovascular diseases. Atherosclerosis is a complex process that involves many cell types such as vessel cells, immune cells and endothelial progenitor cells; however, many, if not all, studies do not report the sex of the cells. This review focuses on sex differences in human cells involved in the atherosclerotic process, emphasizing the role of sex hormones. Furthermore, we report sex differences and issues related to the processes that determine the fate of the cells such as apoptotic and autophagic mechanisms. The analysis of the data reveals that there are still many gaps in our knowledge regarding sex influences in atherosclerosis, largely for the cell types that have not been well studied, stressing the urgent need for a clear definition of experimental conditions and the inclusion of both sexes in preclinical studies.

Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps

The genome of potato, a major global food crop, was recently sequenced. The work presented here details the integration of the potato reference genome (DM) with a new sequence-tagged site marker−based linkage map and other physical and genetic maps of potato and the closely related species tomato. Primary anchoring of the DM genome assembly was accomplished by the use of a diploid segregating population, which was genotyped with several types of molecular genetic markers to construct a new ~936 cM linkage map comprising 2469 marker loci. In silico anchoring approaches used genetic and physical maps from the diploid potato genotype RH89-039-16 (RH) and tomato. This combined approach has allowed 951 superscaffolds to be ordered into pseudomolecules corresponding to the 12 potato chromosomes. These pseudomolecules represent 674 Mb (~93%) of the 723 Mb genome assembly and 37,482 (~96%) of the 39,031 predicted genes. The superscaffold order and orientation within the pseudomolecules are closely collinear with independently constructed high density linkage maps. Comparisons between marker distribution and physical location reveal regions of greater and lesser recombination, as well as regions exhibiting significant segregation distortion. The work presented here has led to a greatly improved ordering of the potato reference genome superscaffolds into chromosomal “pseudomolecules”.

Advances in BAC-based physical mapping and map integration strategies in plants

In the advent of next-generation sequencing (NGS) platforms, map-based sequencing strategy has been recently suppressed being too expensive and laborious. The detailed studies on NGS drafts alone indicated these assemblies remain far from gold standard reference quality, especially when applied on complex genomes. In this context the conventional BAC-based physical mapping has been identified as an important intermediate layer in current hybrid sequencing strategy. BAC-based physical map construction and its integration with high-density genetic maps have benefited from NGS and high-throughput array platforms. This paper addresses the current advancements of BAC-based physical mapping and high-throughput map integration strategies to obtain densely anchored well-ordered physical maps. The resulted maps are of immediate utility while providing a template to harness the maximum benefits of the current NGS platforms.

Whole Genome Profiling provides a robust framework for physical mapping and sequencing in the highly complex and repetitive wheat genome

Background

Sequencing projects using a clone-by-clone approach require the availability of a robust physical map. The SNaPshot technology, based on pair-wise comparisons of restriction fragments sizes, has been used recently to build the first physical map of a wheat chromosome and to complete the maize physical map. However, restriction fragments sizes shared randomly between two non-overlapping BACs often lead to chimerical contigs and mis-assembled BACs in such large and repetitive genomes. Whole Genome Profiling (WGP™) was developed recently as a new sequence-based physical mapping technology and has the potential to limit this problem.

Results

A subset of the wheat 3B chromosome BAC library covering 230 Mb was used to establish a WGP physical map and to compare it to a map obtained with the SNaPshot technology. We first adapted the WGP-based assembly methodology to cope with the complexity of the wheat genome. Then, the results showed that the WGP map covers the same length than the SNaPshot map but with 30% less contigs and, more importantly with 3.5 times less mis-assembled BACs. Finally, we evaluated the benefit of integrating WGP tags in different sequence assemblies obtained after Roche/454 sequencing of BAC pools. We showed that while WGP tag integration improves assemblies performed with unpaired reads and with paired-end reads at low coverage, it does not significantly improve sequence assemblies performed at high coverage (25x) with paired-end reads.

Conclusions

Our results demonstrate that, with a suitable assembly methodology, WGP builds more robust physical maps than the SNaPshot technology in wheat and that WGP can be adapted to any genome. Moreover, WGP tag integration in sequence assemblies improves low quality assembly. However, to achieve a high quality draft sequence assembly, a sequencing depth of 25x paired-end reads is required, at which point WGP tag integration does not provide additional scaffolding value. Finally, we suggest that WGP tags can support the efficient sequencing of BAC pools by enabling reliable assignment of sequence scaffolds to their BAC of origin, a feature that is of great interest when using BAC pooling strategies to reduce the cost of sequencing large genomes.

Cerebral cavernous malformations: from genes to proteins to disease

Over the past half century molecular biology has led to great advances in our understanding of angio- and vasculogenesis and in the treatment of malformations resulting from these processes gone awry. Given their sporadic and familial distribution, their developmental and pathological link to capillary telangiectasias, and their observed chromosomal abnormalities, cerebral cavernous malformations (CCMs) are regarded as akin to cancerous growths. Although the exact pathological mechanisms involved in the formation of CCMs are still not well understood, the identification of 3 genetic loci has begun to shed light on key developmental pathways involved in CCM pathogenesis. Cavernous malformations can occur sporadically or in an autosomal dominant fashion. Familial forms of CCMs have been attributed to mutations at 3 different loci implicated in regulating important processes such as proliferation and differentiation of angiogenic precursors and members of the apoptotic machinery. These processes are important for the generation, maintenance, and pruning of every vessel in the body. In this review the authors highlight the latest discoveries pertaining to the molecular genetics of CCMs, highlighting potential new therapeutic targets for the treatment of these lesions.

Gene expression profiling in the rhesus macaque: methodology, annotation and data interpretation

Gene microarray analyses represent potentially effective means for high-throughput gene expression profiling in non-human primates. In the companion article, we emphasize effective experimental design based on the in vivo physiology of the rhesus macaque, whereas this article emphasizes considerations for gene annotation and data interpretation using gene microarray platforms from Affymetrix. Initial annotation of the rhesus genome array was based on Affymetrix human GeneChips. However, annotation revisions improve the precision with which rhesus transcripts are identified. Annotation of the rhesus GeneChip is under continuous revision with large percentages of probesets under multiple annotation systems having undergone multiple reassignments between March 2007 and November 2008. It is also important to consider that quantitation and comparison of gene expression levels across multiple chips requires appropriate normalization. External corroboration of microarray results using PCR-based methodology also requires validation of appropriate internal reference genes for normalization of expression values. Many tools are now freely available to aid investigators with microarray normalization and selection of internal reference genes to be used for independent corroboration of microarray results.

Uprobe: a genome-wide universal probe resource for comparative physical mapping in vertebrates

Interspecies comparisons are important for deciphering the functional content and evolution of genomes. The expansive array of >70 public vertebrate genomic bacterial artificial chromosome (BAC) libraries can provide a means of comparative mapping, sequencing, and functional analysis of targeted chromosomal segments that is independent and complementary to whole-genome sequencing. However, at the present time, no complementary resource exists for the efficient targeted physical mapping of the majority of these BAC libraries. Universal overgo-hybridization probes, designed from regions of sequenced genomes that are highly conserved between species, have been demonstrated to be an effective resource for the isolation of orthologous regions from multiple BAC libraries in parallel. Here we report the application of the universal probe design principal across entire genomes, and the subsequent creation of a complementary probe resource, Uprobe, for screening vertebrate BAC libraries. Uprobe currently consists of whole-genome sets of universal overgo-hybridization probes designed for screening mammalian or avian/reptilian libraries. Retrospective analysis, experimental validation of the probe design process on a panel of representative BAC libraries, and estimates of probe coverage across the genome indicate that the majority of all eutherian and avian/reptilian genes or regions of interest can be isolated using Uprobe. Future implementation of the universal probe design strategy will be used to create an expanded number of whole-genome probe sets that will encompass all vertebrate genomes.

Genetic and physical maps of jerker (Espn(je)) on mouse chromosome 4

The jerker mutation causes degeneration of cochlea and vestibular sensory hair cells in mice. A frame-shift mutation in the actin bundling gene Espin (Espn) leads to hair bundle defects by disrupting the actin filament assembly in stereocilia. Previously, jerker was mapped to distal mouse chromosome 4. Here, analyzing 2536 informative meioses derived from two intersubspecific intercrosses, we localize jerker to a 0.51+/-0.14cM interval on chromosome 4. The following order and distances of genes and markers were determined: D4Mit180-0.44+/-0.13cM-Hes2, Espn(je)-0.08+/-0.06cM-D4Mit356-0.28+/-0.1cM-D4Mit208. A 300kb physical bacterial artificial chromosome (BAC) contig was generated containing the Espn(je) locus. The human homologous region maps to 1p36.31. We present a detailed high-resolution genetic and physical map of markers located at distal chromosome 4 and demonstrate concordance of Espn with jerker.

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.

Positional cloning of the mouse saccharin preference (Sac) locus

Differences in sweetener intake among inbred strains of mice are partially determined by allelic variation of the saccharin preference (Sac) locus. Genetic and physical mapping limited a critical genomic interval containing Sac to a 194 kb DNA fragment. Sequencing and annotation of this region identified a gene (Tas1r3) encoding the third member of the T1R family of putative taste receptors, T1R3. Introgression by serial backcrossing of the 194 kb chromosomal fragment containing the Tas1r3 allele from the high-sweetener-preferring C57BL/6ByJ strain onto the genetic background of the low-sweetener-preferring 129P3/J strain rescued its low-sweetener-preference phenotype. Polymorphisms of Tas1r3 that are likely to have functional significance were identified using analysis of genomic sequences and sweetener-preference phenotypes of genealogically distant mouse strains. Tas1r3 has two common haplotypes, consisting of six single nucleotide polymorphisms: one haplotype was found in mouse strains with elevated sweetener preference and the other in strains relatively indifferent to sweeteners. This study provides compelling evidence that Tas1r3 is equivalent to the Sac locus and that the T1R3 receptor responds to sweeteners.

Physical and radiation hybrid mapping of canine chromosome 12, in a region corresponding to human chromosome 6p12-q12

The positional cloning of the hypocretin receptor 2, the gene for autosomal recessive canine narcolepsy, has led to the development of a physical map spanning a large portion of canine chromosome 12 (CFA12), in a region corresponding to human chromosome 6p12-q13. More than 40 expressed sequence tags (ESTs) were used in homology search experiments, together with chromosome walking, to build both physical and radiation hybrid maps of the CFA12 13-21 region. The resulting map of bacterial artificial chromosome ends, ESTs, and microsatellite markers represents the longest continuous high-density map of the dog genome reported to date. These data further establish the dog as a system for studying disease genes of interest to human populations and highlight feasible approaches for positional cloning of disease genes in organisms where genomic resources are limited.

Contig selection in physical mapping

In physical mapping, one orders a set of genetic landmarks or a library of cloned fragments of DNA according to their position in the genome. Our approach to physical mapping divides the problem into smaller and easier subproblems by partitioning the probe set into independent parts (probe contigs). For this purpose we introduce a new distance function between probes, the averaged rank distance (ARD) derived from bootstrap resampling of the raw data. The ARD measures the pairwise distances of probes within a contig and smoothes the distances of probes across different contigs. It shows distinct jumps at contig borders. This makes it appropriate for contig selection by clustering. We have designed a physical mapping algorithm that makes use of these observations and seems to be particularly well suited to the delineation of reliable contigs. We evaluated our method on data sets from two physical mapping projects. On data from the recently sequenced bacterium Xylella fastidiosa, the probe contig set produced by the new method was evaluated using the probe order derived from the sequence information. Our approach yielded a basically correct contig set. On this data we also compared our method to an approach which uses the number of supporting clones to determine contigs. Our map is much more accurate. In comparison to a physical map of Pasteurella haemolytica that was computed using simulated annealing, the newly computed map is considerably cleaner. The results of our method have already proven helpful for the design of experiments aimed at further improving the quality of a map.

Integration of the Aedes aegypti mosquito genetic linkage and physical maps

Two approaches were used to correlate the Aedes aegypti genetic linkage map to the physical map. STS markers were developed for previously mapped RFLP-based genetic markers so that large genomic clones from cosmid libraries could be found and placed to the metaphase chromosome physical maps using standard FISH methods. Eight cosmids were identified that contained eight RFLP marker sequences, and these cosmids were located on the metaphase chromosomes. Twenty-one cDNAs were mapped directly to metaphase chromosomes using a FISH amplification procedure. The chromosome numbering schemes of the genetic linkage and physical maps corresponded directly and the orientations of the genetic linkage maps for chromosomes 2 and 3 were inverted relative to the physical maps. While the chromosome 2 linkage map represented essentially 100% of chromosome 2, approximately 65% of the chromosome 1 linkage map mapped to only 36% of the short p-arm and 83% of the chromosome 3 physical map contained the complete genetic linkage map. Since the genetic linkage map is a RFLP cDNA-based map, these data also provide a minimal estimate for the size of the euchromatic regions. The implications of these findings on positional cloning in A. aegypti 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 bacterial artificial chromosome library for sequencing the complete human genome

A 30-fold redundant human bacterial artificial chromosome (BAC) library with a large average insert size (178 kb) has been constructed to provide the intermediate substrate for the international genome sequencing effort. The DNA was obtained from a single anonymous volunteer, whose identity was protected through a double-blind donor selection protocol. DNA fragments were generated by partial digestion with EcoRI (library segments 1--4: 24-fold) and MboI (segment 5: sixfold) and cloned into the pBACe3.6 and pTARBAC1 vectors, respectively. The quality of the library was assessed by extensive analysis of 169 clones for rearrangements and artifacts. Eighteen BACs (11%) revealed minor insert rearrangements, and none was chimeric. This BAC library, designated as "RPCI-11," has been used widely as the central resource for insert-end sequencing, clone fingerprinting, high-throughput sequence analysis and as a source of mapped clones for diagnostic and functional studies.

Genetic and physical mapping in mosquitoes: molecular approaches

The genetic background of individual mosquito species and populations within those species influences the transmission of mosquito-borne pathogens to humans. Technical advances in contemporary genomics are contributing significantly to the detailed genetic analysis of this mosquito-pathogen interaction as well as all other aspects of mosquito biology, ecology, and evolution. A variety of DNA-based marker types are being used to develop genetic maps for a number of mosquito species. Complex phenotypic traits such as vector competence are being dissected into their discrete genetic components, with the intention of eventually using this information to develop new methods to prevent disease transmission. Both genetic- and physical-mapping techniques are being used to define and compare genome architecture among and within mosquito species. The integration of genetic- and physical-map information is providing a sound framework for map-based positional cloning of target genes of interest. This review focuses on advances in genome-based analysis and their specific applications to mosquitoes.

Application of bootstrap techniques to physical mapping

Ordering genetic markers or clones from a genomic library into a physical map is a central problem in genetics. In the presence of errors, there is no efficient algorithm known that solves this problem. Based on a standard heuristic algorithm for it, we present a method to construct a confidence neighborhood for a computed solution. We compute a confidence value for putative local solutions derived from bootstrap replicates of the original solution. In the reliable parts, the confidence neighborhood and the computed solution tend to coincide. In regions that are ill-defined by the data, the neighborhood contains additional reasonable alternatives. This offers the possibility of designing further experiments for the badly defined regions to improve the quality of the physical map. We analyze our approach by a simulation study and by application to a dataset of the genome of the bacterium Xylella fastidiosa.

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.

A BAC-based physical map of the major autosomes of Drosophila melanogaster

We constructed a bacterial artificial chromosome (BAC)-based physical map of chromosomes 2 and 3 of Drosophila melanogaster, which constitute 81% of the genome. Sequence tagged site (STS) content, restriction fingerprinting, and polytene chromosome in situ hybridization approaches were integrated to produce a map spanning the euchromatin. Three of five remaining gaps are in repeat-rich regions near the centromeres. A tiling path of clones spanning this map and STS maps of chromosomes X and 4 was sequenced to low coverage; the maps and tiling path sequence were used to support and verify the whole-genome sequence assembly, and tiling path BACs were used as templates in sequence finishing.

Comparative genomic sequence analysis of the human and mouse cystic fibrosis transmembrane conductance regulator genes

The identification of the cystic fibrosis transmembrane conductance regulator gene (CFTR) in 1989 represents a landmark accomplishment in human genetics. Since that time, there have been numerous advances in elucidating the function of the encoded protein and the physiological basis of cystic fibrosis. However, numerous areas of cystic fibrosis biology require additional investigation, some of which would be facilitated by information about the long-range sequence context of the CFTR gene. For example, the latter might provide clues about the sequence elements responsible for the temporal and spatial regulation of CFTR expression. We thus sought to establish the sequence of the chromosomal segments encompassing the human CFTR and mouse Cftr genes, with the hope of identifying conserved regions of biologic interest by sequence comparison. Bacterial clone-based physical maps of the relevant human and mouse genomic regions were constructed, and minimally overlapping sets of clones were selected and sequenced, eventually yielding approximately 1.6 Mb and approximately 358 kb of contiguous human and mouse sequence, respectively. These efforts have produced the complete sequence of the approximately 189-kb and approximately 152-kb segments containing the human CFTR and mouse Cftr genes, respectively, as well as significant amounts of flanking DNA. Analyses of the resulting data provide insights about the organization of the CFTR/Cftr genes and potential sequence elements regulating their expression. Furthermore, the generated sequence reveals the precise architecture of genes residing near CFTR/Cftr, including one known gene (WNT2/Wnt2) and two previously unknown genes that immediately flank CFTR/Cftr.

An integrated physical map for the short arm of human chromosome 5

The short arm of human chromosome 5 contains approximately 48 Mb of DNA and comprises 1.5% of the genome. We have constructed a mega-YAC/ STS map of this region that includes 436 YACs anchored by 216 STSs. By combining and integrating our map with the 5p maps of other groups using the same recombinant DNA library, a comprehensive map was constructed that includes 552 YACs and 504 markers. The YAC map covers >94% of 5p in four YAC contigs, bridges the centromere, and includes an additional 5 Mb of 5q DNA. The average marker density is 95 kb. This integrated 5p map will serve as a resource for the continuing localization of genes on the short arm of human chromosome 5 and as a framework for both generating and aligning the DNA sequence of this region.

Sequence-tagged connectors: a sequence approach to mapping and scanning the human genome

The sequence-tagged connector (STC) strategy proposes to generate sequence tags densely scattered (every 3.3 kilobases) across the human genome by arraying 450,000 bacterial artificial chromosomes (BACs) with randomly cleaved inserts, sequencing both ends of each, and preparing a restriction enzyme fingerprint of each. The STC resource, containing end sequences, fingerprints, and arrayed BACs, creates a map where the interrelationships of the individual BAC clones are resolved through their STCs as overlapping BAC clones are sequenced. Once a seed or initiation BAC clone is sequenced, the minimum overlapping 5' and 3' BAC clones can be identified computationally and sequenced. By reiterating this "sequence-then-map by computer analysis against the STC database" strategy, a minimum tiling path of clones can be sequenced at a rate that is primarily limited by the sequencing throughput of individual genome centers. As of February 1999, we had deposited, together with The Institute for Genomic Research (TIGR), into GenBank 314,000 STCs ( approximately 135 megabases), or 4.5% of human genomic DNA. This genome survey reveals numerous genes, genome-wide repeats, simple sequence repeats (potential genetic markers), and CpG islands (potential gene initiation sites). It also illustrates the power of the STC strategy for creating minimum tiling paths of BAC clones for large-scale genomic sequencing. Because the STC resource permits the easy integration of genetic, physical, gene, and sequence maps for chromosomes, it will be a powerful tool for the initial analysis of the human genome and other complex genomes.

The sleep disorder canine narcolepsy is caused by a mutation in the hypocretin (orexin) receptor 2 gene

Narcolepsy is a disabling sleep disorder affecting humans and animals. It is characterized by daytime sleepiness, cataplexy, and striking transitions from wakefulness into rapid eye movement (REM) sleep. In this study, we used positional cloning to identify an autosomal recessive mutation responsible for this sleep disorder in a well-established canine model. We have determined that canine narcolepsy is caused by disruption of the hypocretin (orexin) receptor 2 gene (Hcrtr2). This result identifies hypocretins as major sleep-modulating neurotransmitters and opens novel potential therapeutic approaches for narcoleptic patients.

An algorithm for the assembly of robust physical maps based on a combination of multi-level hybridization data and fingerprinting data

We have developed an algorithm which combines data obtained from restriction digestion experiments and hybridization experiments to construct robust physical maps of whole chromosomes. The algorithm has been incorporated into a program which accepts hybridization data consisting of an unordered hybridization matrix and fingerprinting data containing band coordinates for each clone. The combined data is used to produce a non-redundant, ordered matrix which can be further reduced to represent a minimum tile coverage of the chromosome. In addition, the method also takes into account multi-level hybridization events which allows for an improved treatment of the hybridization data. The program is evaluated against several other contig building programs using simulated and real data sets. Finally, it is applied to construct a physical map of the 4.1 Mb genome of Ochrobactrum anthropi based on 1387 clones and 70 probes, as well as 624 fingerprints.

A 1-Mb BAC/PAC-based physical map of the autosomal recessive polycystic kidney disease gene (PKHD1) region on chromosome 6

The PKHD1 (polycystic kidney and hepatic disease 1) gene responsible for autosomal recessive polycystic kidney disease has been mapped to 6p21.1-p12 to an approximately 1-cM interval flanked by the markers D6S1714/D6S243 and D6S1024. We have developed a sequence-ready BAC/PAC-based contig map of this region as the next step for the positional cloning of PKHD1. This contig comprising 52 clones spanning approximately 1 Mb was established by content mapping of 44 BAC/PAC-end-derived STSs, 3 known genetic markers, 5 YAC-end-derived STSs, 3 random STSs, 1 previously mapped gene, and 1 EST. The average depth per marker is 6.3 clones, and the average STS density is 20 kb. The genomic clone overlaps were confirmed by restriction fragment fingerprint analysis. A high-resolution BAC/PAC-based contig map is essential to the ultimate goal of identifying the PKHD1 gene.

Error Checking and Graphical Representation of Multiple–Complete–Digest (MCD) Restriction-Fragment Maps

Genetic and physical maps display the relative positions of objects or markers occurring within a target DNA molecule. In constructing maps, the primary objective is to determine the ordering of these objects. A further objective is to assign a coordinate to each object, indicating its distance from a reference end of the target molecule. This paper describes a computational method and a body of software for assigning coordinates to map objects, given a solution or partial solution to the ordering problem. We describe our method in the context of multiple–complete–digest (MCD) mapping, but it should be applicable to a variety of other mapping problems. Because of errors in the data or insufficient clone coverage to uniquely identify the true ordering of the map objects, a partial ordering is typically the best one can hope for. Once a partial ordering has been established, one often seeks to overlay a metric along the map to assess the distances between the map objects. This problem often proves intractable because of data errors such as erroneous local length measurements (e.g., large clone lengths on low-resolution physical maps). We present a solution to the coordinate assignment problem for MCD restriction-fragment mapping, in which a coordinated set of single-enzyme restriction maps are simultaneously constructed. We show that the coordinate assignment problem can be expressed as the solution of a system of linear constraints. If the linear system is free of inconsistencies, it can be solved using the standard Bellman–Ford algorithm. In the more typical case where the system is inconsistent, our program perturbs it to find a new consistent system of linear constraints, close to those of the given inconsistent system, using a modified Bellman–Ford algorithm. Examples are provided of simple map inconsistencies and the methods by which our program detects candidate data errors and directs the user to potential suspect regions of the map.

Construction and validation of yeast artificial chromosome contig maps by RecA-assisted restriction endonuclease cleavage

RecA-assisted restriction endonuclease (RARE) cleavage is an "Achilles' heel" approach to restriction mapping whereby a RecA-protein-oligodeoxynucleotide complex protects an individual restriction site from methylation, thus limiting subsequent digestion to a single, predetermined site. We have used RARE cleavage to cut yeast artificial chromosomes (YACs) at specific EcoRI sites located within or adjacent to sequence-tagged sites (STSs). Each cleavage reaction produces two YAC fragments whose sizes are a direct measure of the position of the STS in the YAC. In this fashion, we have positioned 45 STSs within a contig of 19 independent YACs and constructed a detailed RARE-cleavage map that represents 8.4 Mbp of human chromosome 6p21.3-22. By comparing maps of overlapping YACs, we were able to detect seven internal deletions that ranged from approximately 75 kbp to approximately 1 Mbp in size. Thirteen pairs of EcoRI sites were targeted for double RARE cleavage in uncloned total human DNA. The excised fragments, up to 2 Mbp in size, were resolved by pulsed-field gel electrophoresis and were detected by hybridization. In general, the genomic RARE-cleavage results support the YAC-based map. In one case, the distance in uncloned DNA between the two terminal EcoRI sites of a YAC insert was approximately 1 Mbp larger than the YAC itself, indicating a major deletion. The general concept of RARE-cleavage mapping as well as its applications and limitations are discussed.

Evidence for a novel osteosarcoma tumor-suppressor gene in the chromosome 18 region genetically linked with Paget disease of bone

Paget disease of bone, or "osteitis deformans," is a bone disorder characterized by rapid bone remodeling resulting in abnormal bone formation. It is the second most common metabolic bone disease after osteoporosis, affecting 3%-5% of subjects aged >40 years. Recent evidence suggests that predisposition to Paget disease may have a genetic component. Genetic linkage analysis of families with multigenerational Paget disease shows linkage to a region of chromosome 18q near the polymorphic locus D18S42. Approximately 1% of Paget patients develop osteosarcoma, which represents an increase in risk that is several thousandfold over that of the general population. Osteosarcoma in Paget patients is the underlying basis for a significant fraction of osteosarcomas occurring after age 60 years. Our analysis of tumor-specific loss of constitutional heterozygosity (LOH) in 96 sporadic osteosarcomas has identified a putative tumor-suppressor locus that maps to chromosome 18q. We have localized this tumor-suppressor locus between D18S60 and D18S42, a region tightly linked to familial Paget disease. Analysis of osteosarcomas from patients with Paget disease revealed that these tumors also undergo LOH in this region. These findings suggest that the association between Paget disease and osteosarcoma is the result of a single gene or two tightly linked genes on chromosome 18.

Impact and implications of yeast and human artificial chromosomes

Artificial chromosomes have been developed in the last 10 years to sustain genome mapping and, more recently, to begin initiating functional studies and some approaches to gene therapy. The use of yeast artificial chromosomes (YACs) in mapping the human X chromosome is reported as an example. The requirements which have postponed the development of human artificial chromosomes have now been relatively met, and some prospects are previewed here.

Contig maps and genomic sequencing identify candidate genes in the usher 1C locus

Usher syndrome 1C (USH1C) is a congenital condition manifesting profound hearing loss, the absence of vestibular function, and eventual retinal degeneration. The USH1C locus has been mapped genetically to a 2- to 3-cM interval in 11p14-15.1 between D11S899 and D11S861. In an effort to identify the USH1C disease gene we have isolated the region between these markers in yeast artificial chromosomes (YACs) using a combination of STS content mapping and Alu-PCR hybridization. The YAC contig is approximately 3.5 Mb and has located several other loci within this interval, resulting in the order CEN-LDHA-SAA1-TPH-D11S1310-(D11S1888/KCNC1 )-MYOD1-D11S902D11S921-D11S 1890-TEL. Subsequent haplotyping and homozygosity analysis refined the location of the disease gene to a 400-kb interval between D11S902 and D11S1890 with all affected individuals being homozygous for the internal marker D11S921. To facilitate gene identification, the critical region has been converted into P1 artificial chromosome (PAC) clones using sequence-tagged sites (STSs) mapped to the YAC contig, Alu-PCR products generated from the YACs, and PAC end probes. A contig of >50 PAC clones has been assembled between D11S1310 and D11S1890, confirming the order of markers used in haplotyping. Three PAC clones representing nearly two-thirds of the USH1C critical region have been sequenced. PowerBLAST analysis identified six clusters of expressed sequence tags (ESTs), two known genes (BIR, SUR1) mapped previously to this region, and a previously characterized but unmapped gene NEFA (DNA binding/EF hand/acidic amino-acid-rich). GRAIL analysis identified 11 CpG islands and 73 exons of excellent quality. These data allowed the construction of a transcription map for the USH1C critical region, consisting of three known genes and six or more novel transcripts. Based on their map location, these loci represent candidate disease loci for USH1C. The NEFA gene was assessed as the USH1C locus by the sequencing of an amplified NEFA cDNA from an USH1C patient; however, no mutations were detected.

Application of interspersed repetitive sequence polymerase chain reaction for construction of yeast artificial chromosome contigs

Construction of physical maps across candidate regions is one of the rate-limiting steps of positional cloning projects. To date, most physical maps have been constructed by polymerase chain reaction (PCR)-based sequence-tagged site (STS) content mapping. While effective, this technique has a number of disadvantages including problems with yeast artificial chromosome (YAC) chimerism, the time and effort required to generate new STSs from YAC ends, the cost of primer synthesis for large contiging projects, and the time, effort, and expense necessary for screening each STS in the two-tiered hierarchical YAC library screening format. An alternative strategy, interspersed repetitive sequence (IRS) PCR genomics, alleviates many of these constraints. Clonal overlap is detected by hybridization of individual IRS-PCR products to IRS-PCR product pools of the three-dimensional coordinate pools of YAC libraries in dot-blot format. Entire libraries can be screened in a single step, and multiple libraries can be screened simultaneously. Cloning YAC fragments, sequencing, and primer generation are eliminated, increasing the efficiency of contig construction and reducing the expense. In addition, the genomic location of the individual IRS-PCR products can also be simultaneously determined by screening either interspecific backcrosses or radiation hybrid panels, in dot-blot format, confirming contig extension in the region of interest.

High-resolution physical map of the immunoglobulin lambda variant gene cluster assembled by quantitative DNA fiber mapping

Quantitative DNA fiber mapping (QDFM) allows rapid construction of near-kilobase-resolution physical maps by hybridizing specific probes to individual stretched DNA molecules. We evaluated the utility of QDFM for the large-scale physical mapping of a rather unstable, repeat-rich 850-kb region encompassing the immunoglobulin lambda variant (IGLV) gene segments. We mapped a minimal tiling path composed of 32 cosmid clones to three partially overlapping yeast artificial chromosome (YAC) clones and determined the physical size of each clone, the extent of overlap between clones, and contig orientation, as well as the sizes of gaps between adjacent contigs. Regions of germline DNA for which we had no YAC coverage were characterized by cosmid to cosmid hybridizations. Compared to other methods commonly used for physical map assembly, QDFM is a rapid, versatile technique delivering unambiguous data necessary for map closure and preparation of sequence-ready minimal tiling paths.

A physical map of human chromosome 7: an integrated YAC contig map with average STS spacing of 79 kb

The construction of highly integrated and annotated physical maps of human chromosomes represents a critical goal of the ongoing Human Genome Project. Our laboratory has focused on developing a physical map of human chromosome 7, a approximately 170-Mb segment of DNA that corresponds to an estimated 5% of the human genome. Using a yeast artificial chromosome (YAC)-based sequence-tagged site (STS)-content mapping strategy, 2150 chromosome 7-specific STSs have been established and mapped to a collection of YACs highly enriched for chromosome 7 DNA. The STSs correspond to sequences generated from a variety of DNA sources, with particular emphasis placed on YAC insert ends, genetic markers, and genes. The YACs include a set of relatively nonchimeric clones from a human-hamster hybrid cell line as well as clones isolated from total genomic libraries. For map integration, we have localized 260 STSs corresponding to Genethon genetic markers and 259 STSs corresponding to markers orders by radiation hybrid (RH) mapping on our YAC contigs. Analysis of the data with the program SEGMAP results in the assembly of 22 contigs that are "anchored" on the Genethon genetic map, the RH map, and/or the cytogenetic map. These 22 contigs are ordered relative to one another, are (in all but 3 cases) oriented relative to the centromere and telomeres, and contain > 98% of the mapped STSs. The largest anchored YAC contig, accounting for most of 7p, contains 634 STSs and 1260 YACs. An additional 14 contigs, accounting for approximately 1.5% of the mapped STSs, are assembled but remain unanchored on either the genetic or RH map. Therefore, these 14 "orphan" contigs are not ordered relative to other contigs. In our contig maps, adjacent STSs are connected by two or more YACs in > 95% of cases. With 2150 mapped STSs, our map provides an average STS spacing of approximately 79 kb. The physical map we report here exceeds the goal of 100-kb average STS spacing and should provide an excellent framework for systematic sequencing of the chromosome.

Linkage mapping and phenotypic analysis of autosomal dominant Pallister-Hall syndrome

Pallister-Hall syndrome is a human developmental disorder that is inherited in an autosomal dominant pattern. The phenotypic features of the syndrome include hypothalamic hamartoma, polydactyly, imperforate anus, laryngeal clefting, and other anomalies. Here we describe the clinical characterisation of a family with 22 affected members and the genetic mapping of the corresponding locus. Clinical, radiographic, and endoscopic evaluations showed that this disorder is a fully penetrant trait with variable expressivity and low morbidity. By analysing 60 subjects in two families using anonymous STRP markers, we have established linkage to 7p13 by two point analysis with D7S691 resulting in a lod score of 7.0 at theta = 0, near the GLI3 locus. Deletions and translocations in GLI3 are associated with the Greig cephalopolysyndactyly syndrome. Although Greig cephalopolysyndactyly syndrome has some phenotypic overlap with Pallister-Hall syndrome, these two disorders are clinically distinct. The colocalisation of loci for these distinct phenotypes led us to analyse GLI3 for mutations in patients with Pallister-Hall syndrome. We have previously shown GLI3 mutations in two other small, moderately affected families with Pallister-Hall syndrome. The linkage data reported here suggest that these larger, mildly affected families may also have mutations in GLI3.

Binning clones by hybridization with complex probes: statistical refinement of an inner product mapping method

Molecular methods that use long-range information to solve genomics problems (i.e., top-down strategies) efficiently have become increasingly prominent in the genomics literature. One such method, an implementation of inner product mapping (IPM), uses noisy, long-range radiation hybrid (RH)/YAC overlap data and relatively noise-free RH/STS overlap data to localize clones to specific chromosomal regions. Because the molecular data are rarely noise-free, statistical models tailored to the top-down molecular methods make the methods far more effective. We develop two statistical models for IPM (or any other top-down strategy of similar form), a parametric logit model and a nonparametric order-restricted model, and show how these models can be implemented within a hierarchical Bayes framework. Using these models, we refine the chromosome 11 map reported in M. Perlin et al. (1995, Genomics 28: 315-327). Our analyses improve the IPM map, both in terms of successful localization of clones and in terms of the confidence with which they are localized.

A sequence-ready high-resolution physical map of the best macular dystrophy gene region in 11q12-q13

Best disease, an autosomal dominant inherited macular degenerative disorder, was previously localized between D11S1765 and UGB (uteroglobin) in 11q13 by genetic linkage analysis. Since this region was found to be refractory to cloning in YAC (yeast artificial chromosome)-based vectors, a P1 artificial chromosome (PAC) contig was assembled. Gridded PAC libraries representing a 16-fold genome equivalent were screened by hybridization using PCR products representing STSs derived from YAC end sequences, markers binned to 11q13, and PAC-derived insert ends. A highly marker dense approximately 1.7-Mb PAC contig that encompassed the disease gene region was constructed, allowing us to order accurately the markers throughout the region and to provide the most precise estimate of its physical size. Using this contig, thus far we have mapped seven anonymous ESTs and five known genes into this region. This high-resolution physical map will facilitate the isolation of polymorphic markers for refinement of the disease gene region, as well as the identification of candidate genes by exon trapping, cDNA selection, and gene prediction from PAC-derived genomic sequence.

X chromosome map at 75-kb STS resolution, revealing extremes of recombination and GC content

A YAC/STS map of the X chromosome has reached an inter-STS resolution of 75 kb. The map density is sufficient to provide YACs or other large-insert clones that are cross-validated as sequencing substrates across the chromosome. Marker density also permits estimates of regional gene content and a detailed comparison of genetic and physical map distances. Five regions are detected with relatively high G + C, correlated with gene richness; and a 17-Mb region with very low recombination is revealed between the Xq13.3 [XIST] and Xq21.3 XY homology loci.

A collection of 1814 human chromosome 7-specific STSs

An established goal of the ongoing Human Genome Project is the development and mapping of sequence-tagged sites (STSs) every 100 kb, on average, across all human chromosomes. En route to constructing such a physical map of human chromosome 7, we have generated 1814 chromosome 7-specific STSs. The corresponding PCR assays were designed by the use of DNA sequence determined in our laboratory (79%) or generated elsewhere (21%) and were demonstrated to be suitable for screening yeast artificial chromosome (YAC) libraries. This collection provides the requisite landmarks for constructing a physical map of chromosome 7 at < 100-kb average spacing of STSs.

A high-resolution physical map of human chromosome 11

The development of a highly reliable physical map with landmark sites spaced an average of 100 kbp apart has been a central goal of the Human Genome Project. We have approached the physical mapping of human chromosome 11 with this goal as a primary target. We have focused on strategies that would utilize yeast artificial chromosome (YAC) technology, thus permitting long-range coverage of hundreds of kilobases of genomic DNA, yet we sought to minimize the ambiguities inherent in the use of this technology, particularly the occurrence of chimeric genomic DNA clones. This was achieved through the development of a chromosome 11-specific YAC library from a human somatic cell hybrid line that has retained chromosome 11 as its sole human component. To maximize the efficiency of YAC contig assembly and extension, we have employed an Alu-PCR-based hybridization screening system. This system eliminates many of the more costly and time-consuming steps associated with sequence tagged site content mapping such as sequencing, primer production, and hierarchical screening, resulting in greater efficiency with increased throughput and reduced cost. Using these approaches, we have achieved YAC coverage for >90% of human chromosome 11, with an average intermarker distance of <100 kbp. Cytogenetic localization has been determined for each contig by fluorescent in situ hybridization and/or sequence tagged site content. The YAC contigs that we have generated should provide a robust framework to move forward to sequence-ready templates for the sequencing efforts of the Human Genome Project as well as more focused positional cloning on chromosome 11.

Mapping the RP10 locus for autosomal dominant retinitis pigmentosa on 7q: refined genetic positioning and localization within a well-defined YAC contig

Retinitis pigmentosa is a genetically heterogeneous disease that has autosomal dominant, autosomal recessive and X-linked forms. Autosomal dominant retinitis pigmentosa (adRP) has thus far been associated with eight distinct loci, including the rhodopsin and peripherin/RDS genes as well as unidentified genes on chromosomes 7p, 7q, 8q, 17p, 17q, and 19q. The RP10 locus for adRP on chromosome 7q was first mapped in a Spanish family; later, an unrelated American family was identified that also showed linkage to 7q. By combining the linkage results from both families, we are able to assign the disease gene to a 5-cM interval on 7q. Based on extensive physical mapping of this region, the genetic interval is now fully contained within a approximately 5-Mb segment on a well-defined YAC contig. These studies significantly reduce the size of the RP10 critical region, exclude a number of possible candidate genes, and provide the necessary cloned DNA for the positional cloning of the RP10 gene.

Generation of expressed sequence tags and sequence-tagged sites as physical landmarks in the mosquito, Aedes aegypti, genome

Twenty-three clones were developed as expressed sequence tags (ESTs) and sequence-tagged sites (STSs) that provide broad coverage of the mosquito, Aedes aegypti, genome at an average spacing of about 7.2 cM. Each of these clones had been mapped previously as a restriction fragment length polymorphism (RFLP) marker. Nineteen of these clones represent anonymous cDNAs and 4 clones represent known genes. Forty-seven percent of the anonymous cDNAs showed significant deduced amino acid sequence similarity to previously described sequences from a wide variety of organisms. STSs developed from RFLPs will provide effective tools for rapid screening of YAC and cosmid libraries, and will therefore, provide anchor points for contig construction for any region in the A. aegypti genome. Also, they facilitate the simultaneous integration of physical and genetic mapping data for A. aegypti.

The polymerase chain reaction: applications for the detection of foodborne pathogens

Faster methods for the detection of foodborne microbial pathogens are needed. The polymerase chain reaction (PCR) can amplify specific segments of DNA and is used to detect and identify bacterial genes responsible for causing diseases in humans. The major features and requirements for the PCR are described along with a number of important variations. A considerable number of PCR-based assays have been developed, but they have been applied most often to clinical and environmental samples and more rarely for the detection of foodborne microorganisms. Much of the difficulty in implementing PCR for the analysis of food samples lies in the problems encountered during the preparation of template DNAs from food matrices; a variety of approaches and considerations are examined. PCR methods developed for the detection and identification of particular bacteria, viruses, and parasites found in foods are described and discussed, and the major features of these reactions are summarized.

Refined localization of the cerebral cavernous malformation gene (CCM1) to a 4-cM interval of chromosome 7q contained in a well-defined YAC contig

Cerebral cavernous malformations (CCM) are vascular lesions present in some 20 million people worldwide that are responsible for seizures, migraine, hemorrhage, and other neurologic problems. Familial cases ofCCM can be inherited as an autosomal dominant disorder with variable expression. A gene for CCM (CCM/)was recently mapped to a 33-cM segment of chromosome 7q in a large Hispanic family (Dubovsky et al.1995). Here, the collection of several new short tandem repeat polymorphisms (STRPs) within the region of interest on 7q and the refinement of the marker order in this region using both linkage analysis in CEPH families and especially YAC-based STS content mapping are described. Affected members of three Hispanic families share allele haplotypes indicating a common ancestral mutation within these families. Using the shared haplotype information along with analysis of crossovers in affected individuals from both the Hispanic and Caucasian families, the region likely to contain the CCMI gene has been reduced to a 4-cM segment of 7q between D7S2410 and D7S689. All markers within the refined chromosomal segment were located on a single YAC contig estimated to be approximately 2 Mb in size. Four potential candidate genes have been mapped to this region.

The human obese (OB) gene: RNA expression pattern and mapping on the physical, cytogenetic, and genetic maps of chromosome 7

The recently identified mouse obese (ob) gene apparently encodes a secreted protein that may function in the signaling pathway of adipose tissue. Mutations in the mouse ob gene are associated with the early development of gross obesity. A detailed knowledge concerning the RNA expression pattern and precise genomic location of the human homolog, the OB gene, would facilitate examination of the role of this gene in the inheritance of human obesity. Northern blot analysis revealed that OB RNA is present at a high level in adipose tissue but at much lower levels in placenta and heart. OB RNA is undetectable in a wide range of other tissues. Comparative mapping of mouse and human DNA indicated that the ob gene is located within a region of mouse chromosome 6 that is homologous to a portion of human chromosome 7q. We mapped the human OB gene on a yeast artificial chromosome (YAC) contig from chromosome 7q31.3 that contains 43 clones and 19 sequence-tagged sites (STSs). Among the 19 STSs are eight corresponding to microsatellite-type genetic markers, including seven (CA)n repeat-type Genethon markers. Because of their close physical proximity to the human OB gene, these eight genetic markers represent valuable tools for analyzing families with evidence of hereditary obesity and for investigating the possible association between OB mutations and human obesity.