Towards fully automated genotyping: use of an X linked recessive spastic paraplegia family to test alternative analysis methods

Advances in dinucleotide-based genetic maps open possibilities for large scale genotyping at high resolution. The current rate-limiting steps in use of these dense maps is data interpretation (allele definition), data entry, and statistical calculations. We have recently reported automated allele identification methods. Here we show that a 10-cM framework map of the human X chromosome can be analyzed on two lanes of an automated sequencer per individual (10–12 loci per lane). We use this map and analysis strategy to generate allele data for an X-linked recessive spastic paraplegia family with a known PLP mutation. We analyzed 198 genotypes in a single gel and used the data to test three methods of data analysis: manual meiotic breakpoint mapping, automated concordance analysis, and whole chromosome multipoint linkage analysis. All methods pinpointed the correct location of the gene. We propose that multipoint exclusion mapping may permit valid inflation of LOD scores using the equation max LOD — (next best LOD).     

Integration of the PiGMaP and USD A maps for porcine chromosome 14

In order to align two previously published genetic linkage maps, a set of four of the United States Department of Agriculture (USDA) microsatellite linkage markers was mapped in the International Pig Gene Mapping Project (PiGMaP) reference families. Two-point linkage analysis was used between these USDA markers and the set of genes and markers previously mapped on the PiGMaP chromosome 14 map-Markers with threshold lod scores of three or greater were used for multipoint map construction. The USDA and PigGMaP linkage maps of chromosome 14 were aligned using the four USDA microsatellite markers along with three markers that are common to both maps. The PiGMaP genetic linkage map order for chromosome 14 was confirmed and the map was expanded to 193 cM with addition of the new markers.     

A detailed multipoint map of human chromosome 4 provides evidence for linkage heterogeneity and position-specific recombination rates

Utilizing the CEPH reference panel and genotypic data for 53 markers, we have constructed a 20-locus multipoint genetic map of human chromosome 4. New RFLPs are reported for four loci. The map integrates a high-resolution genetic map of 4p16 into a continuous map extending to 4q31 and an unlinked cluster of three loci at 4q35. The 20 linked markers form a continuous linkage group of 152 cM in males and 202 cM in females. Likely genetic locations are provided for 25 polymorphic anonymous sequences and 28 gene-specific RFLPs. The map was constructed employing the LINKAGE and CRIMAP computational methodologies to build the multipoint map via a stepwise algorithm. A detailed 10-point map of the 4p16 region constructed from the CEPH panel provides evidence for heterogeneity in the linkage maps constructed from families segregating for Huntington disease (HD). It additionally provides evidence for position-specific recombination frequencies in the telomeric region of 4p.     

Comparison of microsatellites, single-nucleotide polymorphisms (SNPs) and composite markers derived from SNPs in linkage analysis

There is growing evidence that a map of dense single-nucleotide polymorphisms (SNPs) can outperform a map of sparse microsatellites for linkage analysis. There is also argument as to whether a clustered SNP map can outperform an evenly spaced SNP map. Using Genetic Analysis Workshop 14 simulated data, we compared for linkage analysis microsatellites, SNPs, and composite markers derived from SNPs. We encoded the composite markers in a two-step approach, in which the maximum identity length contrast method was employed to allow for recombination between loci. A SNP map 2.3 times as dense as a microsatellite map (approximately 2.9 cM compared to approximately 6.7 cM apart) provided slightly less information content (approximately 0.83 compared to approximately 0.89). Most inheritance information could be extracted when the SNPs were spaced < 1 cM apart. Comparing the linkage results on using SNPs or composite markers derived from them based on both 3 cM and 0.3 cM resolution maps, we showed that the inter-SNP distance should be kept small (< 1 cM), and that for multipoint linkage analysis the original markers and the derived composite markers had similar power; but for single point linkage analysis the resulting composite markers lead to more power. Considering all factors, such as information content, flexibility of analysis method, map errors, and genotyping errors, a map of clustered SNPs can be an efficient design for a genome-wide linkage scan.     

MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations

With the advent of RFLPs, genetic linkage maps are now being assembled for a number of organisms including both inbred experimental populations such as maize and outbred natural populations such as humans. Accurate construction of such genetic maps requires multipoint linkage analysis of particular types of pedigrees. We describe here a computer package, called MAPMAKER, designed specifically for this purpose. The program uses an efficient algorithm that allows simultaneous multipoint analysis of any number of loci. MAPMAKER also includes an interactive command language that makes it easy for a geneticist to explore linkage data. MAPMAKER has been applied to the construction of linkage maps in a number of organisms, including the human and several plants, and we outline the mapping strategies that have been used.     

Linkage analysis using single nucleotide polymorphisms

We performed multipoint linkage analysis using 83 markers from the SNP Consortium (TSC) SNP linkage map in 3 regions covering 190 cM previously scanned with microsatellite markers and found to be linked to type 2 diabetes. Since the average linkage disequilibrium present in the TSC SNP marker clusters is relatively low, we assumed the intracluster genetic distances were a reasonable small nonzero distance (0.03 cM) and performed linkage analysis using GENEHUNTER PLUS and ASM linkage analysis software. We found that for the pedigree structures and missing data patterns in our samples the average information content in all three regions and the LOD score curves in two regions obtained from the TSC SNP markers were similar to results obtained from microsatellite marker maps with 10 cM average spacing. We also give an algorithm which extends the Lander-Green algorithm to permit multipoint linkage analysis of clusters of tightly linked markers with arbitrarily high levels of intracluster linkage disequilibrium.     

Development and applications of a molecular genetic linkage map of the mouse genome

Interspecific mouse backcrosses provide almost limitless genetic variation for gene mapping. We have used interspecific backcrosses to develop the first comprehensive molecular genetic linkage map of the mouse genome. More than 600 loci have been positioned on the map; the current average map resolution is less than 3 cM. Since all loci were mapped using a single backcross panel, gene order can be determined unambiguously. With this level of resolution, it is now possible to position any new locus on the linkage map with virtually 100% certainty. In this article, we review how interspecific linkage maps are constructed, the salient features of our linkage map, and some of the many applications of interspecific linkage maps, in general, for future research.     

Genetic Linkage Map of Fishes of the Genus Xiphophorus (Teleostei: Poeciliidae)

Analysis of genotypes of 76 polymorphic loci in more than 2600 backcross hybrid individuals derived from intra- and interspecific genetic crosses of fishes of the genus Xiphophorus (Poeciliidae) resulted in the identification of 17 multipoint linkage groups containing 55 protein-coding loci and one sex chromosome-linked pigment pattern gene. Multipoint linkage analyses identified highly probable gene orders for 10 linkage groups. The total genome length was estimated to be approximately 18 Morgans. Comparisons of the Xiphophorus linkage map with those of other fishes, amphibians and mammals suggested that fish gene maps are remarkably similar and probably retain many syntenic groups present in the ancestor of all vertebrates.     

Examination of a region showing linkage map discrepancies across sheep breeds

The availability of accurate linkage maps is an important step for the localization of genetic variants of interest. However, most studies in livestock assume the published map is applicable in their population despite the large differences between the breeds of a species. A region of sheep Chromosome 1 was previously identified as providing evidence for a marker order inconsistent with the published linkage map. In this study the identified region was investigated in more detail. Four microsatellite markers covering the central 5 cM of the inconsistent region and two flanking markers were genotyped in three sheep breeds, a commercial population (Charollais), an experimental population (Scottish Blackface), and a feral population (Soay). With the inclusion of the published linkage map, this provided evidence for three different marker orders across four sheep populations. Evidence for selection in this region was investigated using both a single-point allelic competition model and a multipoint allele-sharing statistic. Only the Charollais population provided evidence for selection, with significant transmission bias observed at marker BM7145. The implications of variation in linkage maps on the design and analysis of fine-mapping studies are discussed.     

A primary linkage map of the human chromosome 11q22-23 region

We have constructed a genetic map of the human chromosomal region 11q22-23 by multipoint linkage analysis of 13 DNA polymorphisms that we have condensed into eight loci. An analysis for linkage disequilibrium between tightly linked probe/enzyme systems allows us to make specific recommendations for future DNA typing at these loci. The resulting sex-averaged multipoint map spans approximately 80 cM and differs considerably from previously reported genetic maps of this region. Our mathematically derived "most likely order" of the markers is compatible with physical mapping data using somatic cell hybrids. The known localizations of at least 14 functional genes and several disease loci to 11q22-23, including ataxia telangiectasia, make the mapping of this region especially relevant to studies of disease pathogenesis.     

Genetic and physical maps of human chromosome 4 based on dinucleotide repeats.

Characterization of inherited variations within tandem arrays of dinucleotide repeats has substantially advanced the construction of genetic maps using linkage approaches over the last several years. Using a backbone of 10 newly identified microsatellite repeats on human chromosome 4 and 6 previously identified short tandem repeat element polymorphisms, we have constructed several genetic maps and a physical map of human chromosome 4. The genetic and physical maps are in complete concordance with each other. The genetic maps include a 15-locus microsatellite-based linkage map, a framework map of high support incorporating a total of 39 independent loci, a 25-locus high-heterozygosity, easily used index map, and a gene-based comprehensive map that provides the best genetic location for 35 genes mapped to chromosome 4. The 16 microsatellite markers are each localized to one of nine regions of chromosome 4, delineated by a panel of somatic cell hybrids. These results demonstrate the utility of PCR-based repeat elements for the construction of genetic maps and provide a valuable resource for continued high-resolution mapping of chromosome 4 and of genetic disorders to this chromosome.     

Modeling segregation distortion for viability selection I. Reconstruction of linkage maps with distorted markers

Molecular markers have been widely used to map quantitative trait loci (QTL). The QTL mapping partly relies on accurate linkage maps. The non-Mendelian segregation of markers, which affects not only the estimation of genetic distance between two markers but also the order of markers on a same linkage group, is usually observed in QTL analysis. However, these distorted markers are often ignored in the real data analysis of QTL mapping so that some important information may be lost. In this paper, we developed a multipoint approach via Hidden Markov chain model to reconstruct the linkage maps given a specified gene order while simultaneously making use of distorted, dominant and missing markers in an F₂ population. The new method was compared with the methods in the MapManager and Mapmaker programs, respectively, and verified by a series of Monte Carlo simulation experiments along with a working example. Results showed that the adjusted linkage maps can be used for further QTL or segregation distortion locus (SDL) analysis unless there are strong evidences to prove that all markers show normal Mendelian segregation.     

Linkage analysis using sex-specific recombination fractions with GENEHUNTER-MODSCORE

MOTIVATION: Sex-specific marker maps have become increasingly available. We have implemented the usage of sex-specific recombination frequencies in the GENEHUNTER-MODSCORE program that performs multipoint linkage analysis. Furthermore, we have devised a consistent method to choose the combinations of male and female genetic positions at which linkage scores should be calculated. Marker coordinates can be read automatically from publicly available genetic maps. RESULTS: In a MOD-score analysis of the COGA dataset provided for Genetic Analysis Workshop 14, the highest linkage peak on chromosome 1 further increases when using sex-specific maps, while some smaller peaks are decreased. Simulations confirm that the MOD score can be biased when a sex-averaged instead of the correct sex-specific map is employed. This shows that an adequate modeling of the female:male ratio of genetic distances is important, especially for complex traits. AVAILABILITY: The new version of GENEHUNTER-MODSCORE can be downloaded from the following website: http://www.staff.uni-marburg.de/~strauchk/software.html     

A comparison of RFLP maps based on anther culture derived, selfed, and hybrid progenies of Solanum chacoense.

Comparative RFLP linkage maps were constructed using five segregating populations derived from two self-incompatible lines (termed PI 230582 and PI 458314) of diploid tuber-bearing Solanum chacoense Bitt. The analysis was based on 84 RFLP loci identified by 73 different cDNA clones. Distortion of expected Mendelian segregation ratios was observed; less than 10% of the markers showed a skewed segregation in the gametes forming the F1, hybrid population compared with 30% in the selfed population and 46 and 70% in the two populations produced by anther culture. For the anther culture derived populations, most of the skewed loci were scattered throughout the genome, whereas in the populations derived from selfing, they were found primarily in linkage group 1, around the S locus. In this study, we also found that the rate of meiotic recombination could differ between the male and female gametes produced by our parental lines. Thus, male gametes of line PI 458314 showed significantly less recombination as assessed by the total length of the map (206 cM for male gametes vs. 375 cM for female gametes) and the phenomenon was genome-wide. In contrast, the maps from the gametes of PI 230582 had about the same length, but some linkage groups were longer in the female gametes, while others were longer in the male gametes. Key words : Solanum chacoense, RFLP, anther culture, skewed segregation, self-incompatibility, sex differences in recombination.     

A 12,000-rad porcine radiation hybrid (IMNpRH2) panel refines the conserved synteny between SSC12 and HSA17

Reverse or bidirectional Zoo-FISH suggests that synteny between porcine chromosome 12 (SSC12) and human chromosome 17 (HSA17) is completely conserved. The construction of a high-resolution radiation hybrid (RH) map for SSC12 provides a unique opportunity to determine whether chromosomal synteny is reflected at the molecular level by comparative gene mapping of SSC12 and HSA17. We report an initial, high-resolution RH map of SSC12 on the 12,000-rad IMNpRH2 panel using CarthaGene software. This map contains a total of 320 markers, including 20 microsatellites and 300 ESTs/genes, covering approximately 4836.9 cR12,000. The markers were ordered in 16 linkage groups at LOD 6.0 using framework markers previously mapped on the IMpRH7000-rad SSC12 and porcine genetic maps. Ten linkage groups ordered more than 10 markers, with the largest containing 101 STSs. The resolution of the current RH map is approximately 15.3 kb/cR on SSC12, a significant improvement over the second-generation EST SSC12 RH7000-rad map of 103 ESTs and 15 framework markers covering approximately 2287.2 cR7000. Compared to HSA17, six distinct segments were identified, revealing macro-rearrangements within the apparently complete synteny between SSC12 and HSA17. Further analysis of the order of 245 genes (ESTs) on HSA17 and SSC12 also revealed several micro-rearrangements within a synteny segment. A high-resolution SSC12 RH12,000-rad map will be useful in fine-mapping QTL and as a scaffold for sequencing this chromosome.     

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.     

Construction of a high-density,high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.     

Rearranged gene order between pig and human in a quantitative trait loci region on SSC3

A quantitative trait locus (QTL) for ovulation rate on chromosome 3 that peaks at 36 cM has been identified in a Meishan-White composite resource population with an additive effect of 2.2 corpora lutea. As part of an effort to identify the responsible gene(s), typing of additional genes on the INRA-University of Minnesota porcine radiation hybrid (IMpRH) map of SSC3 and comparative analysis of gene order was conducted. We placed 52 known genes and expressed sequence tags, two BAC-end sequences and one microsatellite (SB42) on a framework map that fills gaps on previous RH maps. Data were analysed for two-point and multipoint linkage with the IMpRH mapping tool and were submitted to the IMpRH database (http://imprh.toulouse.inra.fr/). Gene order was confirmed for 42 loci residing in the QTL region (spanning c. 17 Mb of human sequence) by using the high-resolution IMpRH2 panel. Carthagène (http://www.inra.fr/internet/departments/MIA/T/CarthaGene) was used to estimate multipoint marker distance and order using all public markers on SSC3 in the IMpRH database and those typed in this study. For the high-resolution map, only data for markers typed in both panels were used. Comparative analysis of human and porcine maps identified conservation of gene order for SSC3q and multiple blocks of conserved segments for SSC3p, which included six distinct segments of HSA7 and two segments of HSA16. The results of this study allow significant refinement of the SSC3p region that contains an ovulation rate QTL.     

Improving estimates of genetic maps: a maximum likelihood approach

As a result of previous large, multipoint linkage studies there is a substantial amount of existing marker data. Due to the increased sample size, genetic maps estimated from these data could be more accurate than publicly available maps. However, current methods for map estimation are restricted to data sets containing pedigrees with a small number of individuals, or cannot make full use of marker data that are observed at several loci on members of large, extended pedigrees. In this article, a maximum likelihood (ML) method for map estimation that can make full use of the marker data in a large, multipoint linkage study is described. The method is applied to replicate sets of simulated marker data involving seven linked loci, and pedigree structures based on the real multipoint linkage study of Abkevich et al. (2003, American Journal of Human Genetics 73, 1271-1281). The variance of the ML estimate is accurately estimated, and tests of both simple and composite null hypotheses are performed. An efficient procedure for combining map estimates over data sets is also suggested.     

A detailed multipoint gene map of chromosome 1q

Utilizing genotyping data for 23 markers, we have constructed a 21-locus multipoint genetic map of the long arm of chromosome 1. Five new RFLPs are reported. The map integrates anonymous loci from previous primary linkage maps and incorporates markers for 10 coding sequences. These markers form a continuous linkage group of 85 cM in males and 141 cM in females. The map was constructed employing the LINKAGE and CRIMAP computational methodologies via a stepwise algorithm.