SimPed: a simulation program to generate haplotype and genotype data for pedigree structures

With the widespread availability of SNP genotype data, there is great interest in analyzing pedigree haplotype data. Intermarker linkage disequilibrium for microsatellite markers is usually low due to their physical distance; however, for dense maps of SNP markers, there can be strong linkage disequilibrium between marker loci. Linkage analysis (parametric and nonparametric) and family-based association studies are currently being carried out using dense maps of SNP marker loci. Monte Carlo methods are often used for both linkage and association studies; however, to date there are no programs available which can generate haplotype and/or genotype data consisting of a large number of loci for pedigree structures. SimPed is a program that quickly generates haplotype and/or genotype data for pedigrees of virtually any size and complexity. Marker data either in linkage disequilibrium or equilibrium can be generated for greater than 20,000 diallelic or multiallelic marker loci. Haplotypes and/or genotypes are generated for pedigree structures using specified genetic map distances and haplotype and/or allele frequencies. The simulated data generated by SimPed is useful for a variety of purposes, including evaluating methods that estimate haplotype frequencies for pedigree data, evaluating type I error due to intermarker linkage disequilibrium and estimating empirical p values for linkage and family-based association studies.     

Incorporating genotypes of relatives into a test of linkage disequilibrium.

Genetic data from autosomal loci in diploids generally consist of genotype data for which no phase information is available, making it difficult to implement a test of linkage disequilibrium. In this paper, we describe a test of linkage disequilibrium based on an empirical null distribution of the likelihood of a sample. Information on the genotypes of related individuals is explicitly used to help reconstruct the gametic phase of the independent individuals. Simulation studies show that the present approach improves on estimates of linkage disequilibrium gathered from samples of completely independent individuals but only if some offspring are sampled together with their parents. The failure to incorporate some parents sharply decreases the sensitivity and accuracy of the test. Simulations also show that for multiallelic data (more than two alleles) our testing procedure is not as powerful as an exact test based on known haplotype frequencies, owing to the interaction between departure from Hardy-Weinberg equilibrium and linkage disequilibrium.     

Joint modeling of linkage and association: identifying SNPs responsible for a linkage signal

Once genetic linkage has been identified for a complex disease, the next step is often association analysis, in which single-nucleotide polymorphisms (SNPs) within the linkage region are genotyped and tested for association with the disease. If a SNP shows evidence of association, it is useful to know whether the linkage result can be explained, in part or in full, by the candidate SNP. We propose a novel approach that quantifies the degree of linkage disequilibrium (LD) between the candidate SNP and the putative disease locus through joint modeling of linkage and association. We describe a simple likelihood of the marker data conditional on the trait data for a sample of affected sib pairs, with disease penetrances and disease-SNP haplotype frequencies as parameters. We estimate model parameters by maximum likelihood and propose two likelihood-ratio tests to characterize the relationship of the candidate SNP and the disease locus. The first test assesses whether the candidate SNP and the disease locus are in linkage equilibrium so that the SNP plays no causal role in the linkage signal. The second test assesses whether the candidate SNP and the disease locus are in complete LD so that the SNP or a marker in complete LD with it may account fully for the linkage signal. Our method also yields a genetic model that includes parameter estimates for disease-SNP haplotype frequencies and the degree of disease-SNP LD. Our method provides a new tool for detecting linkage and association and can be extended to study designs that include unaffected family members.     

Allele and Haplotype Diversity of 26 X-STR Loci in Four Nationality Populations from China

Background

Haplotype analysis of closely associated markers has proven to be a powerful tool in kinship analysis, especially when short tandem repeats (STR) fail to resolve uncertainty in relationship analysis. STR located on the X chromosome show stronger linkage disequilibrium compared with autosomal STR. So, it is necessary to estimate the haplotype frequencies directly from population studies as linkage disequilibrium is population-specific.

Methodology and Findings

Twenty-six X-STR loci including six clusters of linked markers DXS6807-DXS8378-DXS9902(Xp22), DXS7132-DXS10079-DXS10074-DXS10075-DXS981 (Xq12), DXS6801-DXS6809-DXS6789-DXS6799(Xq21), DXS7424-DXS101-DXS7133(Xq22), DXS6804-GATA172D05(Xq23), DXS8377-DXS7423 (Xq28) and the loci DXS6800, DXS6803, DXS9898, GATA165B12, DXS6854, HPRTB and GATA31E08 were typed in four nationality (Han, Uigur, Kazakh and Mongol) samples from China (n = 1522, 876 males and 646 females). Allele and haplotype frequency as well as linkage disequilibrium data for kinship calculation were observed. The allele frequency distribution among different populations was compared. A total of 5–20 alleles for each locus were observed and altogether 289 alleles for all the selected loci were found. Allele frequency distribution for most X-STR loci is different in different populations. A total of 876 male samples were investigated by haplotype analysis and for linkage disequilibrium. A total of 89, 703, 335, 147, 39 and 63 haplotypes were observed. Haplotype diversity was 0.9584, 0.9994, 0.9935, 0.9736, 0.9427 and 0.9571 for cluster I, II, III, IV, V and VI, respectively. Eighty-two percent of the haplotype of cluster IIwas found only once. And 94% of the haplotype of cluster III show a frequency of <1%.

Conclusions

These results indicate that allele frequency distribution for most X-STR loci is population-specific and haplotypes of six clusters provide a powerful tool for kinship testing and relationship investigation. So it is necessary to obtain allele frequency and haplotypes data of the linked loci for forensic application.     

Toward localization of the Werner syndrome gene by linkage disequilibrium and ancestral haplotyping: lessons learned from analysis of 35 chromosome 8p11.1-21.1 markers.

Werner syndrome (WS) is an autosomal recessive disorder characterized by premature onset of a number of age-related diseases. The gene for WS, WRN, has been mapped to the 8p 11.1-21.1 region with further localization through linkage disequilibrium mapping. Here we present the results of linkage disequilibrium and ancestral haplotype analyses of 35 markers to further refine the location of WRN. We identified an interval in this region in which 14 of 18 markers tested show significant evidence of linkage disequilibrium in at least one of the two populations tested. Analysis of extended and partial haplotypes covering 21 of the markers studied supports the existence of both obligate and probable ancestral recombinant events which localize WRN almost certainly to the interval between D8S2196 and D8S2186, and most likely to the narrower interval between D8S2168 and D8S2186. These haplotype analyses also suggest that there are multiple WRN mutations in each of the two populations under study. We also present a comparison of approaches to performing disequilibrium tests with multiallelic markers, and show that some commonly used approximations for such tests perform poorly in comparison to exact probability tests. Finally, we discuss some of the difficulties introduced by the high mutation rate at microsatellite markers which influence our ability to use ancestral haplotype analysis to localize disease genes.     

Linkage disequilibrium detected between myotonic dystrophy and the anonymous marker D19S63 in the Spanish population

Summary We used the following polymorphic markers: APOC2 (BanI, BglI, TaqI), CKMM (NcoI, TaqI), and D19S63 (PstI) to haplotype 33 Spanish myotonic dystrophy (DM) families. We analysed the allele and haplotype frequencies of our sample, and the possible association of alleles or haplotypes with the disease. We found a slight linkage disequilibrium between APOC2 (BanI) and DM, but no disequilibrium when using all other APOC2 and CKMM RFLPs; this agrees with data previously reported. In addition, we found a very strong linkage disequilibrium when using D19S63 (PstI), the + allele being associated with the DM locus. This disequilibrium in the Spanish population indicates that D19S63 is very close to the DM locus.     

Linkage disequilibrium at the Angelman syndrome gene UBE3A in autism families.

Autistic disorder is a neurodevelopmental disorder with a complex genetic etiology. Observations of maternal duplications affecting chromosome 15q11-q13 in patients with autism and evidence for linkage and linkage disequilibrium to markers in this region in chromosomally normal autism families indicate the existence of a susceptibility locus. We have screened the families of the Collaborative Linkage Study of Autism for several markers spanning a candidate region covering approximately 2 Mb and including the Angelman syndrome gene (UBE3A) and a cluster of gamma-aminobutyric acid (GABA(A)) receptor subunit genes (GABRB3, GABRA5, and GABRG3). We found significant evidence for linkage disequilibrium at marker D15S122, located at the 5' end of UBE3A. This is the first report, to our knowledge, of linkage disequilibrium at UBE3A in autism families. Characterization of null alleles detected at D15S822 in the course of genetic studies of this region showed a small (approximately 5-kb) genomic deletion, which was present at somewhat higher frequencies in autism families than in controls.     

Genetic linkage and transmission disequilibrium of marker haplotypes at chromosome 1q41 in human systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies to a wide range of self-antigens. Recent genome screens have implicated numerous chromosomal regions as potential SLE susceptibility loci. Among these, the 1q41 locus is of particular interest, because evidence for linkage has been found in several independent SLE family collections. Additionally, the 1q41 locus appears to be syntenic with a susceptibility interval identified in the NZM2410 mouse model for SLE. Here, we report the results of genotyping of 11 microsatellite markers within the 1q41 region in 210 SLE sibpair and 122 SLE trio families. These data confirm the modest evidence for linkage at 1q41 in our family collection (LOD = 1.21 at marker D1S2616). Evidence for significant linkage disequilibrium in this interval was also found. Multiple markers in the region exhibit transmission disequilibrium, with the peak single marker multiallelic linkage disequilibrium noted at D1S490 (pedigree disequilibrium test [PDT] global P value = 0.0091). Two- and three-marker haplotypes from the 1q41 region similarly showed strong transmission distortion in the collection of 332 SLE families. The finding of linkage together with significant transmission disequilibrium provides strong evidence for a susceptibility locus at 1q41 in human SLE.     

Disequilibrium Pattern Analysis. I. Theory

We have developed a method, disequilibrium pattern analysis, for examining the disequilibrium distribution of the entire array of two locus multiallelic haplotypes in a population. It is shown that a selected haplotype will produce a distinct pattern of linkage disequilibrium values for all generations while the selection is acting. This pattern will also presumably be maintained for many generations after the selection event, until the disequilibrium pattern is eventually broken down by genetic drift and recombination. Related haplotypes, sharing an allele with a selected haplotype, assume a value of linkage disequilibrium proportional to the frequency of the unshared allele and have a single negative value of the normalized linkage disequilibrium. The analysis assumes zero linkage disequilibrium for all allelic combinations initially. The same basic results continue to apply if the selection involves a new mutant, the occurrence of which creates linkage disequilibrium for some haplotypes. The disequilibrium pattern predicted under selection is robust with respect to the influence of migration and random genetic drift. This method is applicable to population data having linked polymorphic loci including that determined from protein or DNA sequencing.     

The inhibitory effects of camptothecin, a topoisomerase I inhibitor, on collagen synthesis in fibroblasts from patients with systemic sclerosis

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies to a wide range of self-antigens. Recent genome screens have implicated numerous chromosomal regions as potential SLE susceptibility loci. Among these, the 1q41 locus is of particular interest, because evidence for linkage has been found in several independent SLE family collections. Additionally, the 1q41 locus appears to be syntenic with a susceptibility interval identified in the NZM2410 mouse model for SLE. Here, we report the results of genotyping of 11 microsatellite markers within the 1q41 region in 210 SLE sibpair and 122 SLE trio families. These data confirm the modest evidence for linkage at 1q41 in our family collection (LOD = 1.21 at marker D1S2616). Evidence for significant linkage disequilibrium in this interval was also found. Multiple markers in the region exhibit transmission disequilibrium, with the peak single marker multiallelic linkage disequilibrium noted at D1S490 (pedigree disequilibrium test [PDT] global P value = 0.0091). Two- and three-marker haplotypes from the 1q41 region similarly showed strong transmission distortion in the collection of 332 SLE families. The finding of linkage together with significant transmission disequilibrium provides strong evidence for a susceptibility locus at 1q41 in human SLE.