Map error reduction: using genetic and sequence-based physical maps to order closely linked markers

The Marshfield comprehensive genetic maps are frequently used for linkage and association studies, however, for some regions of these maps the marker order has low level of likelihood ratio support. In order to investigate the level of statistical support and the accuracy of the genetic maps compared to sequence-based physical maps, two approximately 30 cM autosomal regions were selected. The first region was selected from chromosome 3 and consisted predominately of draft sequence. The second region was selected from chromosome 21 and consisted of finished sequence data. The physical order of these markers was based upon their position on Celera (CEL) and Human Genome Project-Santa Cruz (HGP-sc) sequence-based physical maps. The chromosome 3 and 21 regions contained 100 and 61 markers, respectively, on the Marshfield genetic map. The genetic and physical map order was consistent for 88.9 and 89.2% of the markers in the region on chromosome 3 and 21, respectively. Using a novel scoring criterion to assess inconsistent marker order between genetic and physical maps, it was determined that the physical order was likely the correct order for 3.3 and 7.1% of the markers in the chromosome 3 and 21 regions, respectively. To increase the accuracy of the order of markers selected for fine mapping a method is presented which combines information from genetic and sequence-based physical maps.     

Localization of a novel autosomal recessive non-syndromic hearing impairment locus (DFNB38) to 6q26-q27 in a consanguineous kindred from Pakistan

For autosomal recessive nonsyndromic hearing impairment over 30 loci have been mapped and 19 genes have been identified. DFNB38, a novel locus for autosomal recessive nonsyndromic hearing impairment, was localized in a consanguineous Pakistani kindred to 6q26-q27. The affected family members present with profound prelingual sensorineural hearing impairment and use sign language for communications. Linkage was established to microsatellite markers located on chromosome 6q26-q27 (Multipoint lod score 3.6). The genetic region for DFNB38 spans 10.1 cM according to the Marshfield genetic map and is bounded by markers D6S980 and D6S1719. This genetic region corresponds to 3.4 MB on the sequence-based physical map.     

MAP-O-MAT: internet-based linkage mapping

SUMMARY: MAP-O-MAT is a web-based server for automated linkage mapping of human polymorphic DNA markers. MAP-O-MAT facilitates the verification of order and map distances for custom mapping sets using genotype data from the CEPH database, and from the Marshfield, SNP Consortium and Rutgers linkage maps (exclusive to the deCODE genotyping data). The CRI-MAP program is used for likelihood calculations and some mapping algorithms, and physical map positions are provided from the human genome assembly. AVAILABILITY: MAP-O-MAT is located at http://compgen.rutgers.edu/mapomat/ CONTACT: matise@biology.rutgers.edu.     

Integration of the Physical and Genetic Linkage Map for Human Chromosome 13

We have used a panel of somatic cell hybrids containing different rearrangements of human chromosome 13 to integrate genetic and physical maps of this chromosome. The positions of 17 translocation/deletion breakpoints on human chromosome 13 have been determined relative to the microsatellite markers on the genetic linkage map compiled by Généthon. Because markers on maps from several other Consortium groups have also been analyzed using many of the same hybrids, it was possible to relate these with the Généthon map. The position of all of the chromosome breakpoints have been placed, wherever possible, between two adjacent markers on the genetic linkage maps using PCR analysis for the presence/absence of the markers in the somatic cell hybrids. The positions of the breakpoints have already been determined cytogenetically, and some of these breakpoints are located at landmark positions on the chromosome. The relative density of markers along the chromosome differs between independently derived maps, and, based on the known locations of certain breakpoints in the physical map, inconsistencies in the genetic maps have been identified.     

A Radiation Hybrid Map of 95 STSs Spanning Human Chromosome 13q

We have constructed a high-resolution physical map of the long arm of human chromosome 13 using a panel of 94 radiation hybrids. A comprehensive map of 95 chromosome 13-specific sequence tagged sites (STSs) spanning 13q from the presumed centromere at D13Z1 to the known telomere was obtained by multipoint maximum likelihood statistical methods. The 95 markers have an average retention frequency of 10%, with markers closer to the centromere having much greater retention frequencies (22-49%) than distal 13q markers (2-12%) The most likely radiation hybrid map localized the 95 STSs into 54 unique map positions, 34 with odds of 1000:1 or greater; the comprehensive map localized all but 17 STSs with odds exceeding 10:1. The total map length of 13q was 1302 cR₉₀₀₀ (range 6.4-94.4 cR₉₀₀₀) and a physical distance of 98 Mb, so that 1% breakage in the RH panel corresponds to 75 kb. A comparison of the comprehensive RH map to genetic maps of chromosome 13q shows identical locus orders for the common markers, with two exceptions over 1-cM distances. We discuss the possible relationships between the genetic and the radiation hybrid maps.     

The Physical and Genetic Framework of the Maize B73 Genome

Maize is a major cereal crop and an important model system for basic biological research. Knowledge gained from maize research can also be used to genetically improve its grass relatives such as sorghum, wheat, and rice. The primary objective of the Maize Genome Sequencing Consortium (MGSC) was to generate a reference genome sequence that was integrated with both the physical and genetic maps. Using a previously published integrated genetic and physical map, combined with in-coming maize genomic sequence, new sequence-based genetic markers, and an optical map, we dynamically picked a minimum tiling path (MTP) of 16,910 bacterial artificial chromosome (BAC) and fosmid clones that were used by the MGSC to sequence the maize genome. The final MTP resulted in a significantly improved physical map that reduced the number of contigs from 721 to 435, incorporated a total of 8,315 mapped markers, and ordered and oriented the majority of FPC contigs. The new integrated physical and genetic map covered 2,120 Mb (93%) of the 2,300-Mb genome, of which 405 contigs were anchored to the genetic map, totaling 2,103.4 Mb (99.2% of the 2,120 Mb physical map). More importantly, 336 contigs, comprising 94.0% of the physical map (∼1,993 Mb), were ordered and oriented. Finally we used all available physical, sequence, genetic, and optical data to generate a golden path (AGP) of chromosome-based pseudomolecules, herein referred to as the B73 Reference Genome Sequence version 1 (B73 RefGen_v1).     

HumanMSD and MouseMSD: generating genetic maps for human and murine microsatellite markers

We present two Web interfaces that generate genetic maps for given sets of human or mouse microsatellite markers. The genetic maps are generated from available databases using linear interpolation of physical map distances to infer genetic map positions for missing markers in these databases.     

DFNB44, a novel autosomal recessive non-syndromic hearing impairment locus, maps to chromosome 7p14.1-q11.22

The genetic etiology for many forms of hearing impairment (HI) is very diverse. Non-syndromic HI (NSHI) is one of the most heterogeneous traits known. Autosomal recessive forms of prelingual HI account for approximately 75% of hereditary cases. A novel autosomal recessive NSHI locus, DFNB44, was mapped to a 20.9 cM genetic interval on chromosome 7p14.1-q11.22, according to the Marshfield genetic map, in a consanguineous Pakistani family. Multipoint linkage analysis resulted in a maximum LOD score of 5.0 at marker D7S1818. The 3-unit support interval ranged from marker D7S2209 to marker D7S2435, spanning a 30.1 Mb region on the sequence-based physical map.     

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.     

A 3.9-centimorgan-resolution human single-nucleotide polymorphism linkage map and screening set

Recent advances in technologies for high-throughout single-nucleotide polymorphism (SNP)-based genotyping have improved efficiency and cost so that it is now becoming reasonable to consider the use of SNPs for genomewide linkage analysis. However, a suitable screening set of SNPs and a corresponding linkage map have yet to be described. The SNP maps described here fill this void and provide a resource for fast genome scanning for disease genes. We have evaluated 6,297 SNPs in a diversity panel composed of European Americans, African Americans, and Asians. The markers were assessed for assay robustness, suitable allele frequencies, and informativeness of multi-SNP clusters. Individuals from 56 Centre d'Etude du Polymorphisme Humain pedigrees, with >770 potentially informative meioses altogether, were genotyped with a subset of 2,988 SNPs, for map construction. Extensive genotyping-error analysis was performed, and the resulting SNP linkage map has an average map resolution of 3.9 cM, with map positions containing either a single SNP or several tightly linked SNPs. The order of markers on this map compares favorably with several other linkage and physical maps. We compared map distances between the SNP linkage map and the interpolated SNP linkage map constructed by the deCode Genetics group. We also evaluated cM/Mb distance ratios in females and males, along each chromosome, showing broadly defined regions of increased and decreased rates of recombination. Evaluations indicate that this SNP screening set is more informative than the Marshfield Clinic's commonly used microsatellite-based screening set.     

A physical map of the Chinese chestnut (Castanea mollissima) genome and its integration with the genetic map

Three Chinese chestnut bacterial artificial chromosome (BAC) libraries were developed and used for physical map construction. Specifically, high information content fingerprinting was used to assemble 126,445 BAC clones into 1,377 contigs and 12,919 singletons. Integration of the dense Chinese chestnut genetic map with the physical map was achieved via high-throughput hybridization using overgo probes derived from sequence-based genetic markers. A total of 1,026 probes were anchored to the physical map including 831 probes corresponding to 878 expressed sequence tag-based markers. Within the physical map, three BAC contigs were anchored to the three major fungal blight-resistant quantitative trait loci on chestnut linkage groups B, F, and G. A subset of probes corresponding to orthologous genes in poplar showed only a limited amount of conserved gene order between the poplar and chestnut genomes. The integrated genetic and physical map of Chinese chestnut is available at www.fagaceae.org/physical_maps.     

A high-resolution consensus linkage map of the rat, integrating radiation hybrid and genetic maps

We have constructed a high-resolution consensus genetic map of the rat in a single large intercross, which integrates 747 framework markers and 687 positions of our whole-genome radiation hybrid (RH) map of the rat. We selected 136 new gene markers from the GenBank database and assigned them either genetically or physically to rat chromosomes to evaluate the accuracy of the integrated linkage-RH maps in the localization of new markers and to enrich existing comparative mapping data. These markers and 631 D-Got- markers, which are physically mapped but still uncharacterized for evidence of polymorphism, were tested for allele variations in a panel of 16 rat strains commonly used in genetic studies. The consensus linkage map constructed in the GK x BN cross now comprises 1620 markers of various origins, defining 840 resolved genetic positions with an average spacing of 2.2 cM between adjacent loci, and includes 407 gene markers. This whole-genome genetic map will contribute to the advancement of genetic studies in the rat by incorporating gene/EST maps, physical mapping information, and sequence data generated in rat and other mammalian species into genetic intervals harboring disease susceptibility loci identified in rat models of human genetic disorders.     

In silico mapping of important genes and markers available in the public domain for efficient sorghum breeding

Crop genome sequencing projects generate massive amounts of genomic sequence information, and the utilization of this information in applied crop improvement programs has been augmented by the availability of sophisticated bioinformatics tools. Here, we present the possible direct utilization of sequence data from a sorghum genome sequencing project in applied crop breeding programs. Based on sequence homology, we aligned all publicly available simple sequence repeat markers on a sequence-based physical map for sorghum. Linking this physical map with already existing linkage map(s) provides better options for applied molecular breeding programs. When a new set of markers is made available, the new markers can be first aligned on a sequence-based physical map, and those located near the quantitative trait locus (QTL) can be identified from this map, thereby reducing the number of markers to be tested in order to identify polymorphic flanking markers for the QTL for any given donor × recurrent parent combination. Polymorphic markers that are expected (on the basis of their position on the sequence-based physical map) to be closely linked to the target can be used for foreground selection in marker-assisted breeding. This map facilitates the identification of a set of markers representing the entire genome, which would provide better resolution in diversity analyses and further linkage disequilibrium mapping. Filling the gaps in existing linkage maps and fine mapping can be achieved more efficiently by targeting the specific genomic regions of interest. It also opens up new exciting opportunities for comparative mapping and for the development of new genomic resources in related crops, both of which are lagging behind in the current genomic revolution. This paper also presents a number of examples of potential applications of sequence-based physical map for sorghum.     

The challenge of pharmacogenetics and pharmacogenomics

The Human Genome Project along with a number of rival commercial and academic efforts such as that of Celera Genomics looks set to provide a comprehensive map and sequence of the human genome within the next two to three years. The development of techniques in genetic engineering has given us the ability to identify and isolate many genes that are important to the understanding of disease. The pharmaceutical industry is being flooded with information as it enters a period both of prodigious complexity and unprecedented opportunity. The challenge for the industry will be in selecting the best therapeutic targets among the huge number that present themselves. The future of drug development may well lie in the efficient identification of these targets that are most receptive to therapeutic attack.     

Genetic mapping of 14 short tandem repeat polymorphisms on human chromosome 22

We have constructed a linkage map of 14 short tandem repeat polymorphisms (11 with heterozygosity > 70%) on the long arm of human chromosome 22 using 23 non-CEPH pedigrees. Twelve of the markers could be positioned uniquely with a likelihood of at least 1,000:1, and distributed at an average distance of 6.62 cM (range 1.5–16.1 cM). The sex-combined map covers a total of 79.6 cM, the female map 93.2 cM and the male map 64.6 cM. Based on comparisons between physical maps and other genetic maps, we estimate that our map covers 70%–80% of the chromosome. The map integrates markers from previous genetic maps and uniquely positions one marker (D22S307). Data from physical mapping on the location of four genetic markers correlates well with our linkage map, and provides information on an additional marker (D22S315). This map will facilitate high resolution mapping of additional polymorphic loci and disease genes on chromosome 22, and act as a reference for building and verifying physical maps.     

A first-generation porcine whole-genome radiation hybrid map

A whole-genome radiation hybrid (WG-RH) panel was used to generate a first-generation radiation map of the porcine (Sus scrofa) genome. Over 900 Type I and II markers were used to amplify the INRA-University of Minnesota porcine Radiation Hybrid panel (IMpRH) comprised of 118 hybrid clones. Average marker retention frequency of 29.3% was calculated with 757 scorable markers. The RHMAP program established 128 linkage groups covering each chromosome (n = 19) at a lod ≥ 4.8. Fewer than 10% of the markers (59) could not be placed within any linkage group at a lod score ≥4.8. Linkage group order for each chromosome was determined by incorporating linkage data from the swine genetic map as well as physical assignments. The current map has an estimated ratio of ∼70 kb/cR and a maximum theoretical resolution of 145 kb. This initial map forms a template for establishing accurate YAC and BAC contigs and eventual positional cloning of genes associated with complex traits. Received: 8 January 1999 / Accepted: 13 April 1999     

A novel autosomal recessive non-syndromic hearing impairment locus (DFNB47) maps to chromosome 2p25.1-p24.3

Hereditary hearing impairment (HI) displays extensive genetic heterogeneity. Autosomal recessive (AR) forms of prelingual HI account for ~75% of cases with a genetic etiology. A novel AR non-syndromic HI locus (DFNB47) was mapped to chromosome 2p25.1-p24.3, in two distantly related Pakistani kindreds. Genome scan and fine mapping were carried out using microsatellite markers. Multipoint linkage analysis resulted in a maximum LOD score of 4.7 at markers D2S1400 and D2S262. The three-unit support interval was bounded by D2S330 and D2S131. The region of homozygosity was found within the three-unit support interval and flanked by markers D2S2952 and D2S131, which corresponds to 13.2 cM according to the Rutgers combined linkage-physical map. This region contains 5.3 Mb according to the sequence-based physical map. Three candidate genes, KCNF1, ID2 and ATP6V1C2 were sequenced, and were found to be negative for functional sequence variants.     

A genetic linkage map of human chromosome 5 with 60 RFLP loci.

A genetic map of human chromosome 5 that contains 60 restriction fragment length polymorphism (RFLP) loci in one linkage group has been constructed. Segregation data using these markers and 40 large multigenerational families supplied by the Centre d'Etude du Polymorphisme Humain have been collected. Linkage analyses were performed with the program package CRI-MAP; using odds greater than 1000:1, 30 RFLP loci could be placed on the map. This genetic map spans 289 cM sex-equal, 353 cM in females, and 244 cM in males. While the relative rate of recombination for female meioses is nearly twice that of males over much of the chromosome, several instances of statistically significant excess male recombination were observed. The order of probes on the genetic map has been confirmed by their physical order as determined by somatic cell hybrid lines containing deletions of normal chromosome 5. There is concordance between the physical positions of markers and their genetic positions. Our most distal probes on the genetic map are cytologically localized to the most distal portions of the chromosome. This suggests that our genetic map spans most of chromosome 5.     

Evaluation of linkage disequilibrium and its effect on non-parametric multipoint linkage analysis using two high density single-nucleotide polymorphism mapping panels

Genotype data from the Illumina Linkage III SNP panel (n = 4,720 SNPs) and the Affymetrix 10 k mapping array (n = 11,120 SNPs) were used to test the effects of linkage disequilibrium (LD) between SNPs in a linkage analysis in the Collaborative Study on the Genetics of Alcoholism pedigree collection (143 pedigrees; 1,614 individuals). The average r2 between adjacent markers across the genetic map was 0.099 +/- 0.003 in the Illumina III panel and 0.17 +/- 0.003 in the Affymetrix 10 k array. In order to determine the effect of LD between marker loci in a nonparametric multipoint linkage analysis, markers in strong LD with another marker (r2 > 0.40) were removed (n = 471 loci in the Illumina panel; n = 1,804 loci in the Affymetrix panel) and the linkage analysis results were compared to the results using the entire marker sets. In all analyses using the ALDX1 phenotype, 8 linkage regions on 5 chromosomes (2, 7, 10, 11, X) were detected (peak markers p < 0.01), and the Illumina panel detected an additional region on chromosome 6. Analysis of the same pedigree set and ALDX1 phenotype using short tandem repeat markers (STRs) resulted in 3 linkage regions on 3 chromosomes (peak markers p < 0.01). These results suggest that in this pedigree set, LD between loci with spacing similar to the SNP panels tested may not significantly affect the overall detection of linkage regions in a genome scan. Moreover, since the data quality and information content are greatly improved in the SNP panels over STR genotyping methods, new linkage regions may be identified due to higher information content and data quality in a dense SNP linkage panel.