Genome-wide linkage analysis for alcohol dependence: a comparison between single-nucleotide polymorphism and microsatellite marker assays

Both theoretical and applied studies have proven that the utility of single nucleotide polymorphism (SNP) markers in linkage analysis is more powerful and cost-effective than current microsatellite marker assays. Here we performed a whole-genome scan on 115 White, non-Hispanic families segregating for alcohol dependence, using one 10.3-cM microsatellite marker set and two SNP data sets (0.33-cM, 0.78-cM spacing). Two definitions of alcohol dependence (ALDX1 and ALDX2) were used. Our multipoint nonparametric linkage analysis found alcoholism was nominal linked to 12 genomic regions. The linkage peaks obtained by using the microsatellite marker set and the two SNP sets had a high degree of correspondence in general, but the microsatellite marker set was insufficient to detect some nominal linkage peaks. The presence of linkage disequilibrium between markers did not significantly affect the results. Across the entire genome, SNP datasets had a much higher average linkage information content (0.33 cM: 0.93, 0.78 cM: 0.91) than did microsatellite marker set (0.57). The linkage peaks obtained through two SNP datasets were very similar with some minor differences. We conclude that genome-wide linkage analysis by using approximately 5,000 SNP markers evenly distributed across the human genome is sufficient and might be more powerful than current 10-cM microsatellite marker assays.     

Investigation of altering single-nucleotide polymorphism density on the power to detect trait loci and frequency of false positive in nonparametric linkage analyses of qualitative traits

Genome-wide linkage analysis using microsatellite markers has been successful in the identification of numerous Mendelian and complex disease loci. The recent availability of high-density single-nucleotide polymorphism (SNP) maps provides a potentially more powerful option. Using the simulated and Collaborative Study on the Genetics of Alcoholism (COGA) datasets from the Genetics Analysis Workshop 14 (GAW14), we examined how altering the density of SNP marker sets impacted the overall information content, the power to detect trait loci, and the number of false positive results. For the simulated data we used SNP maps with density of 0.3 cM, 1 cM, 2 cM, and 3 cM. For the COGA data we combined the marker sets from Illumina and Affymetrix to create a map with average density of 0.25 cM and then, using a sub-sample of these markers, created maps with density of 0.3 cM, 0.6 cM, 1 cM, 2 cM, and 3 cM. For each marker set, multipoint linkage analysis using MERLIN was performed for both dominant and recessive traits derived from marker loci. Our results showed that information content increased with increased map density. For the homogeneous, completely penetrant traits we created, there was only a modest difference in ability to detect trait loci. Additionally, as map density increased there was only a slight increase in the number of false positive results when there was linkage disequilibrium (LD) between markers. The presence of LD between markers may have led to an increased number of false positive regions but no clear relationship between regions of high LD and locations of false positive linkage signals was observed.     

Handling marker-marker linkage disequilibrium: pedigree analysis with clustered markers

Single-nucleotide polymorphisms (SNPs) are rapidly replacing microsatellites as the markers of choice for genetic linkage studies and many other studies of human pedigrees. Here, we describe an efficient approach for modeling linkage disequilibrium (LD) between markers during multipoint analysis of human pedigrees. Using a gene-counting algorithm suitable for pedigree data, our approach enables rapid estimation of allele and haplotype frequencies within clusters of tightly linked markers. In addition, with the use of a hidden Markov model, our approach allows for multipoint pedigree analysis with large numbers of SNP markers organized into clusters of markers in LD. Simulation results show that our approach resolves previously described biases in multipoint linkage analysis with SNPs that are in LD. An updated version of the freely available Merlin software package uses the approach described here to perform many common pedigree analyses, including haplotyping and haplotype frequency estimation, parametric and nonparametric multipoint linkage analysis of discrete traits, variance-components and regression-based analysis of quantitative traits, calculation of identity-by-descent or kinship coefficients, and case selection for follow-up association studies. To illustrate the possibilities, we examine a data set that provides evidence of linkage of psoriasis to chromosome 17.     

Comparing single-nucleotide polymorphism marker-based and microsatellite marker-based linkage analyses

We compared linkage analysis results for an alcoholism trait, ALDX1 (DSM-III-R and Feigner criteria) using a nonparametric linkage analysis method, which takes into account allele sharing among several affected persons, for both microsatellite and single-nucleotide polymorphism (SNP) markers (Affymetrix and Illumina) in the Collaborative Study on the Genetics of Alcoholism (COGA) dataset provided to participants at the Genetic Analysis Workshop 14 (GAW14). The two sets of linkage results from the dense Affymetrix SNP markers and less densely spaced Illumina SNP markers are very similar. The linkage analysis results from microsatellite and SNP markers are generally similar, but the match is not perfect. Strong linkage peaks were found on chromosome 7 in three sets of linkage analyses using both SNP and microsatellite marker data. We also observed that for SNP markers, using the given genetic map and using the map by converting 1 megabase pair (1 Mb) to 1 centimorgan (cM), did not change the linkage results. We recommend the use of the 1 Mb-to-1 cM converted map in a first round of linkage analysis with SNP markers in which map integration is an issue.     

Comparative informativeness for linkage of multiple SNPs and single microsatellites

Single nucleotide polymorphisms (SNPs), or biallelic markers, are popular in genetic linkage studies due to their abundance in the genome, stability, and ease of scoring. We determined the 'information ratio' (IR) of closely spaced SNPs in simulated nuclear families and affected sib pairs (ASPs). (The IR is the ratio of actual average maximum lod score to the maximum lod score attainable if the marker were fully informative.) The nuclear families included parental information, whereas the ASPs did not. We analyzed these SNPs in two ways: (1) using multipoint analysis, and (2) treating the SNPs as 'composite markers' (i.e., haplotypes, as assigned by GENEHUNTER). (3) We also calculated the IR of a single microsatellite marker with multiple alleles and compared with the IR from the SNPs. For each set of input conditions, we simulated 1000 nuclear families, of 2, 3, 4, or 5 children each, as well as 1000 ASPs. We generated SNP marker data for strings of k = 1, 2, 3, 5, 7, and 10 SNP loci, with no recombination (theta = 0) and no linkage disequilibrium among the SNPs. The MAF (minor allele frequency) was either 0.5 or 0.25, and allele frequencies were the same for all k loci in any analysis. We also generated marker data for one single-locus microsatellite marker, with m = 3, 4, 5, 6, 7, and 9 equally frequent alleles. In all simulations, the disease was fully penetrant dominant, and there was no recombination or linkage disequilibrium among markers or between marker and disease. When multipoint analysis was used, we found that 5-7 closely spaced SNPs were usually enough to yield an IR of approximately 100%, for nuclear families of any size. However, for the ASPs, even 7-10 SNPs yielded an IR of only 70-80%. A microsatellite with 9 equally frequent alleles yielded about the same IR (86-88%) as a string of 4-5 SNPs, in nuclear families. SNPs analyzed as 'composite markers' analyses performed worse, due to the inherent ambiguity of SNP haplotyping.     

Comparison of the power between microsatellite and single-nucleotide polymorphism markers for linkage and linkage disequilibrium mapping of an electrophysiological phenotype

We performed linkage and linkage disequilibrium (LD) mapping analyses to compare the power between microsatellite and single nucleotide polymorphism (SNP) markers. Chromosome-wide analyses were performed for a quantitative electrophysiological phenotype, ttth1, on chromosome 7. Multipoint analysis of microsatellite markers using the variance component (VC) method showed the highest LOD score of 4.20 at 162 cM, near D7S509 (163.7 cM). Two-point analysis of SNPs using the VC method yielded the highest LOD score of 3.98 in the Illumina SNP data and 3.45 in the Affymetrix SNP data around 152-153 cM. In family-based single SNP and SNP haplotype LD analysis, we identified seven SNPs associated with ttth1. We searched for any potential candidate genes in the location of the seven SNPs. The SNPs rs1476640 and rs768055 are located in the FLJ40852 gene (a hypothetical protein), and SNP rs1859646 is located in the TAS2R5 gene (a taste receptor). The other four SNPs are not located in any known or annotated genes. We found the high density SNP scan to be superior to microsatellites because it is effective in downstream fine mapping due to a better defined linkage region. Our study proves the utility of high density SNP in genome-wide mapping studies.     

Fine mapping and SNP analysis of positional candidates at the preeclampsia susceptibility locus (PREG1) on chromosome 2

Genome scans in Icelandic, Australian and New Zealand, and Finnish families have localized putative susceptibility loci for preeclampsia/ eclampsia to chromosome 2. The locus mapped in the Australian and New Zealand study (designated PREG1) was thought to be the same locus as that identified in the Icelandic study. In both these studies, two distinct quantitative trait locus (QTL) regions were evident on chromosome 2. Here, we describe our fine mapping of the PREG1 locus and a genetic analysis of two positional candidate genes. Twenty-five additional microsatellite markers were genotyped within the 74-cM linkage region defined by the combined Icelandic and Australian and New Zealand genome scans. The overall position and shape of the localization evidence obtained using nonparametric multipoint analysis did not change from that seen previously in our 10-cM resolution genome scan; two peaks were displayed, one on chromosome 2p at marker D2S388 (107.46 cM) and the other on chromosome 2q at 151.5 cM at marker D2S2313. Using the robust two-point linkage analysis implemented in the Analyze program, all 25 markers gave positive LOD scores with significant evidence of linkage being seen at marker D2S2313 (151.5 cM), achieving a LOD score of 3.37 under a strict diagnostic model. Suggestive evidence of linkage was seen at marker D2S388 (107.46 cM) with a LOD score of 2.22 under the general diagnostic model. Two candidate genes beneath the peak on chromosome 2p were selected for further analysis using public single nucleotide polymorphisms (SNPs) within these genes. Maximum LOD scores were obtained for an SNP in TACR1 (LOD = 3.5) and for an SNP in TCF7L1 (LOD = 3.33), both achieving genome-wide significance. However, no evidence of association was seen with any of the markers tested. These data strongly support the presence of a susceptibility gene on chromosome 2p11-12 and substantiate the possibility of a second locus on chromosome 2q23.     

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.     

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.     

Comparison of marker types and map assumptions using Markov chain Monte Carlo-based linkage analysis of COGA data

We performed multipoint linkage analysis of the electrophysiological trait ECB21 on chromosome 4 in the full pedigrees provided by the Collaborative Study on the Genetics of Alcoholism (COGA). Three Markov chain Monte Carlo (MCMC)-based approaches were applied to the provided and re-estimated genetic maps and to five different marker panels consisting of microsatellite (STRP) and/or SNP markers at various densities. We found evidence of linkage near the GABRB1 STRP using all methods, maps, and marker panels. Difficulties encountered with SNP panels included convergence problems and demanding computations.     

The extent and distribution of linkage disequilibrium in a multi-hierarchic outbred canine pedigree

A canine integrated linkage-radiation map has been recently constructed by using microsatellite markers. This map, with a good coverage of the canine genome, allows for a genome-wide search for the extent and distribution of linkage disequilibrium derived from linkage and evolutionary forces. In this study, we genotyped an outbred pedigree between Labrador retriever and Greyhound breeds with a set of microsatellite markers (240) from the canine linkage map. Linkage disequilibrium was measured between all syntenic and nonsyntenic marker pairs. Analysis of syntenic pairs revealed a significant correlation (–0.229, P < 0.001) between linkage disequilibrium and genetic distance (log transformed). Significant linkage disequilibria were observed more frequently between syntenic pairs spaced <40 cM than those paced >40 cM. There is a clear trend for linkage disequilibrium to decline with marker distance. From our results, a genome-wide screen with markers at low to moderate density (1–2 per 10 cM) should take full advantage of linkage disequilibrium for quantitative trait locus mapping in dogs. This study supports the appropriateness of linkage disequilibrium analysis to detect and map quantitative trait loci underlying complex traits in dogs.     

A comparison of univariate, bivariate, and trivariate whole-genome linkage screens of genetically correlated electrophysiological endophenotypes

We used a maximum-likelihood based multipoint linkage approach implemented in SOLAR to examine simultaneously linkage for three electrophysiological endophenotypes from the Collaborative Study of the Genetics of Alcoholism: TTTH1, TTTH2, and TTTH3. These endophenotypes have been identified as markers of alcohol dependence susceptibility. Data were from 905 individuals in 143 families. Measured covariates considered included sex, age at electrophysiology data collection, habitual smoking status, and the maximum number of drinks consumed in a 24-hour period. Comparisons were made among genome-wide univariate, bivariate, and trivariate linkage analyses using genotypes based on microsatellite markers supplied by the Center for Inherited Disease Research, and genotypes based on single-nucleotide polymorphism markers provided by Illumina. All LODs were corrected to a standard equivalent to 1 degree of freedom. Using the trivariate approach and the microsatellite-based genotypes, we estimated a maximum multipoint linkage signal of LOD = 2.66 on chromosome 7q at 157 cM. Analyses using the Illumina SNP genotypes produced similar results, yielding a maximum multipoint LOD of 2.95 on 7q at 174 cM. These regions of interest correspond to those identified in the univariate and bivariate linkage screens. Our results suggest that trivariate multipoint linkage analyses have utility in the further characterization of chromosomal regions potentially containing genes influencing the phenotypes being examined. Based on a comparison of the number of LOD scores achieving statistical significance, our results suggest that the microsatellite- and Illumina SNP-based genotypes have similar utility for detecting genomic regions of interest.     

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.     

Guidelines for genotyping in genomewide linkage studies: single-nucleotide-polymorphism maps versus microsatellite maps

Genomewide linkage scans have traditionally employed panels of microsatellite markers spaced at intervals of approximately 10 cM across the genome. However, there is a growing realization that a map of closely spaced single-nucleotide polymorphisms (SNPs) may offer equal or superior power to detect linkage, compared with low-density microsatellite maps. We performed a series of simulations to calculate the information content associated with microsatellite and SNP maps across a range of different marker densities and heterozygosities for sib pairs (with and without parental genotypes), sib trios, and sib quads. In the case of microsatellite markers, we varied density across 11 levels (1 marker every 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9, or 10 cM) and marker heterozygosity across 6 levels (2, 3, 4, 5, 10, or 20 equally frequent alleles), whereas, in the case of SNPs, we varied marker density across 4 levels (1 marker every 0.1, 0.2, 0.5, or 1 cM) and minor-allele frequency across 7 levels (0.5, 0.4, 0.3, 0.2, 0.1, 0.05, and 0.01). When parental genotypes were available, a map consisting of microsatellites spaced every 2 cM or a relatively sparse map of SNPs (i.e., at least 1 SNP/cM) was sufficient to extract most of the inheritance information from the map (>95% in most cases). However, when parental genotypes were unavailable, it was important to use as dense a map of markers as possible to extract the greatest amount of inheritance information. It is important to note that the information content associated with a traditional map of microsatellite markers (i.e., 1 marker every ~10 cM) was significantly lower than the information content associated with a dense map of SNPs or microsatellites. These results strongly suggest that previous linkage studies that employed sparse microsatellite maps could benefit substantially from reanalysis by use of a denser map of markers.     

Linkage disequilibrium inflates type I error rates in multipoint linkage analysis when parental genotypes are missing

OBJECTIVES: Describe the inflation in nonparametric multipoint LOD scores due to inter-marker linkage disequilibrium (LD) across many markers with varied allele frequencies. METHOD: Using simulated two-generation families with and without parents, we conducted nonparametric multipoint linkage analysis with 2 to 10 markers with minor allele frequencies (MAF) of 0.5 and 0.1. RESULTS: Misspecification of population haplotype frequencies by assuming linkage equilibrium caused inflated multipoint LOD scores due to inter-marker LD when parental genotypes were not included. Inflation increased as more markers in LD were included and decreased as markers in equilibrium were added. When marker allele frequencies were unequal, the r2 measure of LD was a better predictor of inflation than D'. CONCLUSION: This observation strongly supports the evaluation of LD in multipoint linkage analyses, and further suggests that unaccounted for LD may be suspected when two-point and multipoint linkage analyses show a marked disparity in regions with elevated r2 measures of LD. Given the increasing popularity of high-density genome-wide SNP screens, inter-marker LD should be a concern in future linkage studies.     

High-density single nucleotide polymorphism screen in a large multiplex neural tube defect family refines linkage to loci at 7p21.1-pter and 2q33.1-q35

BACKGROUND: Neural tube defects (NTDs) are considered complex, with both genetic and environmental factors implicated. To date, no major causative genes have been identified in humans despite several investigations. The first genomewide screen in NTDs demonstrated evidence of linkage to chromosomes 7 and 10. This screen included 44 multiplex families and consisted of 402 microsatellite markers spaced approximately 10 cM apart. Further investigation of the genomic screen data identified a single large multiplex family, pedigree 8776, as primarily driving the linkage results on chromosome 7. METHODS: To investigate this family more thoroughly, a high-density single nucleotide polymorphism (SNP) screen was performed. Two-point and multipoint linkage analyses were performed using both parametric and nonparametric methods. RESULTS: For both the microsatellite and SNP markers, linkage analysis suggested the involvement of a locus or loci proximal to the telomeric regions of chromosomes 2q and 7p, with both regions generating a LOD* score of 3.0 using a nonparametric identity by descent relative sharing method. CONCLUSIONS: The regions with the strongest evidence for linkage map proximal to the telomeres on these two chromosomes. In addition to mutations and/or variants in a major gene, these loci may harbor a microdeletion and/or translocation; potentially, polygenic factors may also be involved. This single family may be promising for narrowing the search for NTD susceptibility genes.     

Whole-genome scan, in a complex disease, using 11,245 single-nucleotide polymorphisms: comparison with microsatellites

Despite the theoretical evidence of the utility of single-nucleotide polymorphisms (SNPs) for linkage analysis, no whole-genome scans of a complex disease have yet been published to directly compare SNPs with microsatellites. Here, we describe a whole-genome screen of 157 families with multiple cases of rheumatoid arthritis (RA), performed using 11,245 genomewide SNPs. The results were compared with those from a 10-cM microsatellite scan in the same cohort. The SNP analysis detected HLA*DRB1, the major RA susceptibility locus (P=.00004), with a linkage interval of 31 cM, compared with a 50-cM linkage interval detected by the microsatellite scan. In addition, four loci were detected at a nominal significance level (P<.05) in the SNP linkage analysis; these were not observed in the microsatellite scan. We demonstrate that variation in information content was the main factor contributing to observed differences in the two scans, with the SNPs providing significantly higher information content than the microsatellites. Reducing the number of SNPs in the marker set to 3,300 (1-cM spacing) caused several loci to drop below nominal significance levels, suggesting that decreases in information content can have significant effects on linkage results. In contrast, differences in maps employed in the analysis, the low detectable rate of genotyping error, and the presence of moderate linkage disequilibrium between markers did not significantly affect the results. We have demonstrated the utility of a dense SNP map for performing linkage analysis in a late-age-at-onset disease, where DNA from parents is not always available. The high SNP density allows loci to be defined more precisely and provides a partial scaffold for association studies, substantially reducing the resource requirement for gene-mapping studies.     

A susceptibility gene for psoriatic arthritis maps to chromosome 16q: evidence for imprinting

Several genetic loci have been reported for psoriasis, but none has been specifically linked to psoriatic arthritis (PsA), a condition that affects >10% of patients with psoriasis. A genetic component for PsA is suggested by segregation within families and high concordance among identical twins. We performed a linkage scan to map genes contributing to PsA. We identified 178 patients with PsA out of 906 patients who were included in our genetic study of psoriasis. Using a comprehensive genealogy database, we were able to connect 100 of these into 39 families. We genotyped the patients using a framework marker set of 1,000 microsatellite markers, with an average density of 3 cM, and performed multipoint, affected-only, allele-sharing linkage analysis using the Allegro program. On the basis of the initial results, we genotyped more markers for the most prominent loci. A linkage with a LOD score of 2.17 was observed on chromosome 16q. The linkage analysis, conditioned on paternal transmission to affected individuals, gave a LOD score of 4.19, whereas a LOD score of only 1.03 was observed when conditioned for maternal transmission. A suggestive locus on chromosome 16q has previously been implicated in psoriasis. Our data indicate that a gene at this locus may be involved in paternal transmission of PsA.