Studies on locus content mapping

Systematic evaluation of map quality: human chromosome 22

Marker positions on nine genetic linkage, radiation hybrid, and integrated maps of human chromosome 22 were compared with their corresponding positions in the completed DNA sequence. The proportion of markers whose map position is <250 kb from their respective sequence positions ranges from 100% to 35%. Several discordant markers were identified, as well as four regions that show common inconsistencies across multiple maps. These shared discordant regions surround duplicated DNA segments and may indicate mapping or assembly errors due to sequence homology. Recombination-rate distributions along the chromosome were also evaluated, with male and female meioses showing significantly different patterns of recombination, including an 8-Mb male recombination desert. The distributions of radiation-induced chromosome breakage for the GB4 and the G3 radiation hybrid panels were also evaluated. Both panels show fluctuations in breakage intensity, with different regions of significantly elevated rates of breakage. These results provide support for the common assumption that radiation-induced breaks are generally randomly distributed. The present studies detail the limitations of these important map resources and should prove useful for clarifying potential problems in the human maps and sequence assemblies, as well as for mapping and sequencing projects in and across other species.

The first linkage disequilibrium (LD) maps: delineation of hot and cold blocks by diplotype analysis

Linkage disequilibrium (LD) provides information about positional cloning, linkage, and evolution that cannot be inferred from other evidence, even when a correct sequence and a linkage map based on more than a handful of families become available. We present theory to construct an LD map for which distances are additive and population-specific maps are expected to be approximately proportional. For this purpose, there is only a modest difference in relative efficiency of haplotypes and diplotypes: resolving the latter into 2-locus haplotypes has significant cost or error and increases information by about 50%. LD maps for a cold spot in 19p13.3 and a more typical region in 3q21 are optimized by interval estimates. For a random sample and trustworthy map the value of LD at large distance can be predicted reliably from information over a small distance and does not depend on the evolutionary variance unless the sample size approaches the population size. Values of the association probability that can be distinguished from the value at large distance are determined not by population size but by time since a critical bottleneck. In these examples, omission of markers with significant Hardy-Weinberg disequilibrium does not improve the map, and widely discrepant draft sequences have similar estimates of the genetic parameters. The LD cold spot in 19p13.3 gives an unusually high estimate of time, supporting an argument that this relationship is general. As predicted for a region with ancient haplotypes or uniformly high recombination, there is no clear evidence of LD clustering. On the contrary, the 3q21 region is resolved into alternating blocks of stable and decreasing LD, as expected from crossover clustering. Construction of a genomewide LD map requires data not yet available, which may be complemented but not replaced by a catalog of haplotypes.

A sequence-based integrated map of chromosome 22

The near-completion of the sequence for chromosome 22q revolutionizes map integration. We describe a sequence-based integrated map containing 968 loci including 516 known or predicted gene sequences, 317 STSs not included in these sequences, and 135 nonexpressed multinucleotide polymorphisms. The published sequence spans 34.6 Mb, inclusive of gaps estimated to total 1.1 Mb, compared with a top-down estimate of 43 Mb. This discrepancy is discussed, but will not be resolved until more of the genome is analyzed. The radiation hybrid map has 5% error in order and 34% error in location exceeding 1 Mb. The utility of a composite location based on evidence other than sequence is limited to regions not yet sequenced. A genetic map conditional on sequence order was constructed from pairwise lods. Its length of 74.8 cM in males and 80.2 cM in females is slightly less than the previous estimate not constrained by sequence order. Five recombination hot spots are detected, with differences in location between the sexes. Male recombination correlates with repetitive DNA, whereas female recombination does not. It remains to be seen whether this is true for other human chromosomes. An algorithm to improve the fit of cytogenetic bands sequence location reduces the discrepancies in cytogenetic assignment from 61 to 38. This sequence-based integrated map is represented in the genetic location database (LDB2000), which is available at http://cedar.genetics.soton.ac.uk/public_html/LDB2000.html.

Microsatellite variation and recombination rate in the human genome

Background (purifying) selection on deleterious mutations is expected to remove linked neutral mutations from a population, resulting in a positive correlation between recombination rate and levels of neutral genetic variation, even for markers with high mutation rates. We tested this prediction of the background selection model by comparing recombination rate and levels of microsatellite polymorphism in humans. Published data for 28 unrelated Europeans were used to estimate microsatellite polymorphism (number of alleles, heterozygosity, and variance in allele size) for loci throughout the genome. Recombination rates were estimated from comparisons of genetic and physical maps. First, we analyzed 61 loci from chromosome 22, using the complete sequence of this chromosome to provide exact physical locations. These 61 microsatellites showed no correlation between levels of variation and recombination rate. We then used radiation-hybrid and cytogenetic maps to calculate recombination rates throughout the genome. Recombination rates varied by more than one order of magnitude, and most chromosomes showed significant suppression of recombination near the centromere. Genome-wide analyses provided no evidence for a strong positive correlation between recombination rate and polymorphism, although analyses of loci with at least 20 repeats suggested a weak positive correlation. Comparisons of microsatellites in lowest-recombination and highest-recombination regions also revealed no difference in levels of polymorphism. Together, these results indicate that background selection is not a major determinant of microsatellite variation in humans.

A radiation hybrid map of BTA23: identification of a chromosomal rearrangement leading to separation of the cattle MHC class II subregions

Bovine chromosome 23 (BTA23) contains the bovine major histocompatibility complex (MHC) and is thus of particular interest because of the role of MHC genes in immunity. Previous studies have shown cattle MHC class II genes to be subdivided into two distinct subregions separated by a variable genetic distance of 15-30 cM. To elucidate the genetic events that resulted in the present organization of the class II and other MHC genes, a framework radiation hybrid (RH) map of BTA23 was developed by testing DNA samples from a 5000 rad whole genome RH panel. Twenty-six markers were screened with an average retention frequency of 0.27, ranging from 0.14 to 0.42. Total length of the chromosome was 220 cR5000, with 4.1 cR5000/cM when compared to linkage data. Gene orders for the markers common to both the RH framework map and the consensus framework linkage map are identical. Large centiray intervals, D23S23-D23S7, DYA-D23S24 and CYP21-D23S31, were observed compared to linkage distances. These data may indicate a much larger physical distance or suppression of recombination in the interval separating the class II subregions and also within the class I region than previously estimated. Comparison of 13 Type I genes conserved between BTA23 and the human homolog HSA6p suggests the occurrence of an inversion encompassing the centromeric half of the bovine chromosome, thus explaining the large distance between the bovine class IIa and IIb clusters. These results exemplify the power of RH mapping in solving problems in comparative genomics and evolution. Furthermore, noncongruence of the genetic and physical RH map distances indicates that caution must be observed in using either resource alone in searching for candidate genes controlling traits of economic importance.

Estimation of distances and map construction using radiation hybrids

A method of estimating distances between pairs of genetic markers is described that directly uses their observed joint frequency distribution in a panel of radiation hybrids (RHs). The distance measure is based on the strength of association between marker pairs, which is high for close markers and decays with distance. These distances are then submitted to a previous method that generates linear coordinates for the markers directly from the intermarker distance matrix. This method of map building from RH data is simpler than others, because it uses only the observed joint frequency distributions of markers in the panel, and does not attempt to model unobserved quantities such as the retention of different sized fragments that contain the markers. It also incorporates directly the observed variation in retention of different markers, without needing a model for differential fragment retention dependent on chromosomal location, which is generally not known. Only small, precise distances are used in map construction, thereby reducing any effects of different fragment retention frequencies and local variations in X-ray sensitivity. The method is tested by simulation, and known marker distances and locations are successfully recovered from RH raw data. The method is also applied to publicly available data sets related to the recent transcript map of the human genome.

A high-resolution metric HAPPY map of human chromosome 14

We have mapped 1001 novel sequence-tagged sites on human chromosome 14. The mean spacing between markers is approximately 90 kb, most markers are mapped with a resolution of better than 100 kb, and physical distances are determined. The map was produced using HAPPY mapping, a simple and widely applicable in vitro approach that is analogous to linkage or to radiation hybrid mapping, but that circumvents many of the difficulties and potential artifacts associated with these methods. We show also that the map serves as a robust scaffold for building physical maps using large-insert clones.

A metric map of humans: 23,500 loci in 850 bands

High-resolution maps integrated with the enhanced location data base software (LDB+) give improved estimates of genetic parameters and reveal characteristics of cytogenetic bands. Chiasma interference is intermediate between Kosambi and Carter-Falconer levels, as in Drosophila and the mouse. The autosomal genetic map is 2832 and 4348 centimorgans in males and females, respectively. Telomeric T-bands are strikingly associated with male recombination and gene density. Position and centromeric heterochromatin have large effects, but nontelomeric R-bands are not significantly different from G-bands. Several possible reasons are discussed. These regularities validate the maps, despite their high resolution and inevitable local errors. No other approach has been demonstrated to integrate such a large number of loci, which are increasing at about 45% per year. The maps and the data and software from which they are constructed are available through the Internet (http:@cedar.genetics.soton.ac.uk/public_html). Successive versions of this location data base may also be accessed on CD-ROM.