Strength Of Linkage Disequilibrium Research Articles

Human haplotype block sizes are negatively correlated with recombination rates.

The International Haplotype Map ("HapMap") Project is motivated, in part, by the belief that the organization of the human genome, the mechanics of recombination, and the population-level behavior of alleles at adjacent loci should allow researchers to parse the genome into small segments, or "blocks," that show strong linkage disequilibrium (LD) between alleles at loci within those segments. The discovery and evidence for these blocks is to be based solely on the observed LD strength and patterns between alleles at adjacent loci throughout the genome. Although there are many factors that contribute to LD strength, we assessed the correlation between block structure, in terms of length and percentage of the genome assembled into blocks within a region, and recombination rate obtained from two independent sources. We found evidence of a striking negative correlation between the average recombination rate and average block length, suggesting that recombination rate is a strong contributor to haplotype block structure within the genome. We discuss the potential implications of this negative correlation in the context of the organization, properties, and potential ubiquity of a block-like structure in the human genome.

TNXB locus may be a candidate gene predisposing to schizophrenia.

We report here on the detection of nine single nucleotide polymorphisms (SNPs) near to the NOTCH4 locus in the search for schizophrenia susceptibility genes in the class III region of the human major histocompatibility complex (MHC). We totally analyzed 122 family trios recruited in the UK. The TDT analysis demonstrated that of the nine SNPs, three were associated with schizophrenia, including rs1009382 (P = 0.00047), rs204887 (P = 0.007), and rs8283 (P = 0.015). Both rs1009382 and rs204887 are present in the TNXB locus. The rs1009382 is a non-synonymous SNP located in exon 23 of the gene and its A to G base change causes a Glu2578Gly substitution. The goodness-of-fit test showed that genotypic distribution of rs1009382 was deviated from Hardy-Weinberg equilibrium due to homozygote excess in the patient group (P = 0.01), suggesting that a double dose of a genetic risk may be involved. Possibly, rs1009382 is a candidate SNP predisposing to a schizophrenic illness. Moreover, the test for linkage disequilibrium (LD) between paired SNPs showed that the nine SNPs studied may be in the same LD block with an unexpected pattern as the strength of LD was not correlated with the distance between paired SNPs. The haplotype analysis suggested that there might be more than one disease-related allele located in the class III region of the MHC, and that these alleles possibly confer either susceptibility or resistance to schizophrenia.

Dopamine D2 receptor TaqIA and TaqIB polymorphisms in Parkinson's disease.

In a case control study, we examined the association of DRD2 Taq1A and Taq1B polymorphisms and risk of PD, and evaluated the strength of linkage disequilibrium of the polymorphisms. The Taq1A and Taq1B polymorphisms were in strong linkage disequilibrium. There was, however, no significant association of the two polymorphisms with PD.

Designing an optimum genetic association study using dense SNP markers and family-based sample.

Genetic association analysis using thousands of single nucleotide polymorphism (SNP) markers has become a promising alternative to genome-wide linkage scan. Analysis based on linkage-disequilibrium (LD) is more efficient because meiotic information of past generations is utilized. However, in addition to the physical distance between the disease locus and a marker locus, numerous other factors such as admixture, genetic drift, and multiple mutations can affect the observed value of LD. The effect of these factors in a genomic LD association study must be carefully analyzed to obtain an efficient study design. In the following review, we consider studies using family-based data and carefully study the effects of some of these important design factors, including the sample size, frequency of SNP markers, and marker density. For example, we conclude that (1) for reasonably frequent SNP markers, a moderately large sample of 500 families is appropriate for a moderately stringent significance level (alpha = 0.00009); (2) to maintain a power of 80%, maximal difference in allele frequencies between the disease gene and a SNP marker varies between 0.1 (under additive model) and 0.5 (multiplicative); (3) a map density of 10 cM is appropriate only under idea scenario (moderately large sample size, equal trait/marker allele frequencies, maximum LD strength etc.). Results shown here should have practical implications to designing efficient LD association studies using dense SNP markers.

Evaluation of linkage disequilibrium between chromosome 22q11 single nucleotide polymorphisms in a large outbred population.

To assess the utility of linkage disequilibrium (LD) as a tool for fine-mapping disease genes in non-isolated populations, we have assessed the linkage disequilibrium strength among a series of single nucleotide polymorphisms (SNPs) in an approximate 1 Mb region of human chromosome 22q11. Nineteen random SNPs were discovered and tested across this region with an average spacing of 57 kb (range=1.4-289 kb). These 19 SNPs were genotyped in a population consisting of 444 unrelated pedigrees that were largely collected in the U.S. and U.K. Haplotypes for all pedigrees were derived from pedigree data and over 1,400 haplotypes from unrelated individuals were evaluated for linkage disequilibrium between marker alleles. In addition, linkage disequilibrium between marker alleles was also evaluated using estimated haplotypes without genealogical information (i.e., without parental genotype information). Every marker pair combination was tested for a total of 171 tests and 2x2 contingency tables were constructed to measure LD strength. In general the haplotypes derived from pedigree data provided a more conservative estimate of LD strength. Using genealogical information for estimates of D', 59% (10/17) of marker pairs less than 50 kb apart had D' values >0.30. Finally, we observed a 60 kb region with non-significant LD, which could reflect increased recombination in this region.

Accuracy of Haplotype Frequency Estimation for Biallelic Loci, via the Expectation-Maximization Algorithm for Unphased Diploid Genotype Data

Haplotype analyses have become increasingly common in genetic studies of human disease because of their ability to identify unique chromosomal segments likely to harbor disease-predisposing genes. The study of haplotypes is also used to investigate many population processes, such as migration and immigration rates, linkage-disequilibrium strength, and the relatedness of populations. Unfortunately, many haplotype-analysis methods require phase information that can be difficult to obtain from samples of nonhaploid species. There are, however, strategies for estimating haplotype frequencies from unphased diploid genotype data collected on a sample of individuals that make use of the expectation-maximization (EM) algorithm to overcome the missing phase information. The accuracy of such strategies, compared with other phase-determination methods, must be assessed before their use can be advocated. In this study, we consider and explore sources of error between EM-derived haplotype frequency estimates and their population parameters, noting that much of this error is due to sampling error, which is inherent in all studies, even when phase can be determined. In light of this, we focus on the additional error between haplotype frequencies within a sample data set and EM-derived haplotype frequency estimates incurred by the estimation procedure. We assess the accuracy of haplotype frequency estimation as a function of a number of factors, including sample size, number of loci studied, allele frequencies, and locus-specific allelic departures from Hardy-Weinberg and linkage equilibrium. We point out the relative impacts of sampling error and estimation error, calling attention to the pronounced accuracy of EM estimates once sampling error has been accounted for. We also suggest that many factors that may influence accuracy can be assessed empirically within a data set-a fact that can be used to create "diagnostics" that a user can turn to for assessing potential inaccuracies in estimation.

The effect of epistasis on the structure of hybrid zones.

Within hybrid zones that are maintained by a balance between selection and dispersal, linkage disequilibrium is generated by the mixing of divergent populations. This linkage disequilibrium causes selection on each locus to act on all other loci, thereby steepening clines, and generating a barrier to gene flow. Diffusion models predict simple relations between the strength of linkage disequilibrium and the dispersal rate, sigma, and between the barrier to gene flow, B, and the reduction in mean fitness, W. The aim of this paper is to test the accuracy of these predictions by comparison with an exact deterministic model of unlinked loci (r = 0.5). Disruptive selection acts on the proportion of alleles from the parental populations (p,q): W = exp[-S(4pq)beta], such that the least fit genotype has fitness e-s. Where beta << 1, fitness is reduced for a wide range of intermediate genotypes; where beta >> 1, fitness is only reduced for those genotypes close to p = 0.5. Even with strong epistasis, linkage disequilibria are close to sigma 2p'ip'j/rij, where p'i, p'j are the gradients in allele frequency at loci i, j. The barrier to gene flow, which is reflected in the steepening of neutral clines, is given by [formula: see text] where r, the harmonic mean recombination rate between the neural and selected loci, is here 0.5. This is a close approximation for weak selection, but underestimates B for strong selection. The barrier is stronger for small beta, because hybrid fitness is then reduced over a wider range of p. The widths of the selected clines are harder to predict: though simple approximations are accurate for beta = 1, they become inaccurate for extreme beta because, then, fitness changes sharply with p. Estimates of gene number, made from neutral clines on the assumption that selection acts against heterozygotes, are accurate for weak selection when beta = 1; however, for strong selection, gene number is overestimated. For beta > 1, gene number is systematically overestimated and, conversely, when beta < 1, it is underestimated.

Linkage and association study of the CTLA-4 region in coeliac disease for Italian and Tunisian populations.

Coeliac disease (CD) is a multifactorial disease for which there is an intensive search for genetic risk factors. Some authors found an association between the CTLA-4 region and CD. In the present work, we investigate the possible implication of the CTLA-4 region as a genetic risk factor for CD, through two statistical approaches: the maximum likelihood score (MLS) test in a large Italian sample of affected sib-pairs using polymorphic genetic markers on chromosome 2, and the transmission disequilibrium test (TDT) in continental Italian and Tunisian families using the CTLA-4 exon 1 49 A/G polymorphism. None of these approaches provides evidence for linkage or association between the CTLA-4 region and CD. This might result from a difference in the CTLA-4 region from population to population, either in its involvement as a risk factor or in the strength of linkage disequilibrium.

Sequence Diversity in 36 Candidate Genes for Cardiovascular Disorders

Two strategies involving whole-genome association studies have been proposed for the identification of genes involved in complex diseases. The first one seeks to characterize all common variants of human genes and to test their association with disease. The second one seeks to develop dense maps of single-nucleotide polymorphisms (SNPs) and to detect susceptibility genes through linkage disequilibrium. We performed a molecular screening of the coding and/or flanking regions of 36 candidate genes for cardiovascular diseases. All polymorphisms identified by this screening were further genotyped in 750 subjects of European descent. In the whole set of genes, the lengths explored spanned 53.8 kb in the 5' regions, 68.4 kb in exonic regions, and 13 kb in the 3' regions. The strength of linkage disequilibrium within candidate regions suggests that genomewide maps of SNPs might be efficient ways to identify new disease-susceptibility genes, provided that the maps are sufficiently dense. However, the relatively large number of polymorphisms within coding and regulatory regions of candidate genes raises the possibility that several of them might be functional and that the pattern of genotype-phenotype association might be more complex than initially envisaged, as actually has been observed in some well-characterized genes. These results argue in favor of both genomewide association studies and detailed studies of the overall sequence variation of candidate genes, as complementary approaches.

Identification of a C/G polymorphism in the promoter region of the BRCA1 gene and its use as a marker for rapid detection of promoter deletions.

Reduced expression of BRCA1 has been implicated in sporadic breast cancer, although the mechanisms underlying this phenomenon remain unclear. To determine whether regulatory mutations could account for the reduced expression, we screened the promoter region by sequencing in 20 patients with sporadic disease. No mutations were detected; however, a new polymorphism consisting of a C-to-G base change within the β-promoter was identified, with the frequency of the G allele being 0.34. Close to complete linkage disequilibrium was found between this marker and the Pro871Leu polymorphism, situated in exon 11, which has previously been shown not to be associated with breast or ovarian cancer. This indicates that the C/G polymorphism is also unlikely to play a role in either disease. However, the strength of linkage disequilibrium between these markers permitted their use for rapid screening for genomic deletions within BRCA1. A series of 214 cases with familial breast cancer were analysed using this approach; 88/214 were heterozygous for the promoter polymorphism, thereby excluding a deletion in this region. Among the remaining patients, one hemizygous case reflecting a promoter deletion was successfully identified. Therefore, this study indicates that deletions within the β-promoter region of BRCA1 are an uncommon event in familial breast cancer. Furthermore, it suggests that mutations within the BRCA1 promoter are unlikely to account for the reported decreased expression of BRCA1 in sporadic disease. 1999 Cancer Research Campaign

Plasmodium falciparum antigenic diversity: evidence of clonal population structure.

Plasmodium falciparum, the agent of malignant malaria, is one of mankind's most severe scourges. Efforts to develop preventive vaccines or remedial drugs are handicapped by the parasite's rapid evolution of drug resistance and protective antigens. We examine 25 DNA sequences of the gene coding for the highly polymorphic antigenic circumsporozoite protein. We observe total absence of silent nucleotide variation in the two nonrepeated regions of the gene. We propose that this absence reflects a recent origin (within several thousand years) of the world populations of P. falciparum from a single individual; the amino acid polymorphisms observed in these nonrepeat regions would result from strong natural selection. Analysis of these polymorphisms indicates that: (i) the incidence of recombination events does not increase with nucleotide distance; (ii) the strength of linkage disequilibrium between nucleotides is also independent of distance; and (iii) haplotypes in the two nonrepeat regions are correlated with one another, but not with the central repeat region they span. We propose two hypotheses: (i) variation in the highly polymorphic central repeat region arises by mitotic intragenic recombination, and (ii) the population structure of P. falciparum is clonal--a state of affairs that persists in spite of the necessary stage of physiological sexuality that the parasite must sustain in the mosquito vector to complete its life cycle.