Extensive Linkage Disequilibrium Research Articles

Linkage disequilibrium and demographic history of the isolated population of the Faroe Islands.

The isolated population of the Faroe Islands has a history of recent expansion after being limited to a small size for centuries. Such an isolated population may be ideal for linkage disequilibrium mapping of disease genes if linkage disequilibrium (LD) extends over large regions. Analyses of 18 markers on 12q24.3, spanning a region of 4.3 Mb (16 cM), revealed extensive LD in the Faroese population. Maximum LD was found between marker pairs separated by more than 3.8 Mb. The same region had a maximum LD of only 1.2 and 1.4 Mb respectively in two outbred Danish and British populations analysed here for comparison. The analyses of gene diversity excess at 15 unlinked microsatellite markers did not reveal any sign of a severe bottleneck to have occurred within approximately 1200 years' history of the Faroese population. The extensive LD in this population may, therefore, have arisen primarily by random genetic drift. The implications for future gene mapping studies are discussed.

Extensive Linkage Disequilibrium in Small Human Populations in Eurasia

The extent of linkage disequilibrium (LD) was studied in two small food-gathering populations-Evenki and Saami-and two larger food-producing populations-Finns and Swedes-in northern Eurasia. In total, 50 single-nucleotide polymorphisms (SNPs) from five genes were genotyped using real-time pyrophosphate DNA sequencing, whereas 14 microsatellites were genotyped in two X-chromosomal regions. In addition, hypervariable region I of the mtDNA was sequenced to shed light on the demographic history of the populations. The SNP data, as well as the microsatellite data, reveal extensive levels of LD in Evenki and Saami when compared to Finns and Swedes. mtDNA-sequence variation is compatible with constant population size over time in Evenki and Saami but indicates population expansion in Finns and Swedes. Furthermore, the similarity between Finns and Swedes in SNP allele- and haplotype-frequency distributions indicate that these two populations may share a recent common origin. These findings suggest that populations such as the Evenki and the Saami, rather than the Finns, may be particularly suited for the initial coarse mapping of common complex diseases.

The extent of linkage disequilibrium in Arabidopsis thaliana.

Linkage disequilibrium (LD), the nonrandom occurrence of alleles in haplotypes, has long been of interest to population geneticists. Recently, the rapidly increasing availability of genomic polymorphism data has fueled interest in LD as a tool for fine-scale mapping, in particular for human disease loci. The chromosomal extent of LD is crucial in this context, because it determines how dense a map must be for associations to be detected and, conversely, limits how finely loci may be mapped. Arabidopsis thaliana is expected to harbor unusually extensive LD because of its high degree of selfing. Several polymorphism studies have found very strong LD within individual loci, but also evidence of some recombination. Here we investigate the pattern of LD on a genomic scale and show that in global samples, LD decays within approximately 1 cM, or 250 kb. We also show that LD in local populations may be much stronger than that of global populations, presumably as a result of founder events. The combination of a relatively high level of polymorphism and extensive haplotype structure bodes well for developing a genome-wide LD map in A. thaliana.

Sporadic—but Not Variant—Creutzfeldt-Jakob Disease Is Associated with Polymorphisms Upstream of PRNP Exon 1

Human prion diseases have inherited, sporadic, and acquired etiologies. The appearance of the novel acquired prion disease, variant Creutzfeldt-Jakob disease (vCJD), and the demonstration that it is caused by the same prion strain as that causing bovine spongiform encephalopathy, has led to fears of a major human epidemic. The etiology of classical (sporadic) CJD, which has a worldwide incidence, remains obscure. A common human prion-protein-gene (PRNP) polymorphism (encoding either methionine or valine at codon 129) is a strong susceptibility factor for sporadic and acquired prion disease. However, a quantitative-trait-locus study of prion incubation periods in mice has demonstrated an important factor that is close to Prnp but is independent of its coding sequence or that of the nearby prion-like doppel gene (Prnd). We have analyzed the PRNP locus for such tightly linked susceptibility factors. Fifty-six polymorphic sites have been identified within 25 kb of the PRNP open reading frame, including sites within the PRNP promoter and the PRNP 3' untranslated region. These have been characterized in 61 Centre d'Etude du Polymorphisme Humain (CEPH) families, demonstrating extensive linkage disequilibrium around PRNP and the existence of 11 major European PRNP haplotypes. Haplotype frequencies estimated in healthy U.K. control individuals were very similar to those deduced in the CEPH families. A common haplotype was overrepresented in patients with sporadic CJD (sCJD). Through use of a log-linear modeling approach to simultaneously model Hardy-Weinberg and linkage disequilibria, a significant independent association was found between sCJD and a polymorphism upstream of PRNP exon 1 (P=.005), in addition to the strong susceptibility conferred by codon 129 (P=2x10(-8)). However, although our sample size was necessarily small, no association was found between these polymorphisms and vCJD or iatrogenic CJD, in keeping with their having distinct disease mechanisms. In addition, there was no evidence of a PRNP founder effect in the first reported geographical cluster of vCJD.

The Genetic Dissection of Complex Traits in a Founder Population

We estimated broad heritabilities (H(2)) and narrow heritabilities (h(2)) and conducted genomewide screens, using a novel association-based mapping approach for 20 quantitative trait loci (QTLs) among the Hutterites, a founder population that practices a communal lifestyle. Heritability estimates ranged from.21 for diastolic blood pressure (DBP) to.99 for whole-blood serotonin levels. Using a multipoint method to detect association under a recessive model we found evidence of major QTLs for six traits: low-density lipoprotein (LDL), triglycerides, lipoprotein (a) (Lp[a]), systolic blood pressure (SBP), serum cortisol, and whole-blood serotonin. Second major QTLs for Lp(a) and for cortisol were identified using a single-point method to detect association under a general two-allele model. The heritabilities for these six traits ranged from.37 for triglycerides to.99 for serotonin, and three traits (LDL, SBP, and serotonin) had significant dominance variances (i.e., H(2) > h(2)). Surprisingly, there was little correlation between measures of heritability and the strength of association on a genomewide screen (P>.50), suggesting that heritability estimates per se do not identify phenotypes that are influenced by genes with major effects. The present study demonstrates the feasibility of genomewide association studies for QTL mapping. However, even in this young founder population that has extensive linkage disequilibrium, map densities <<5 cM may be required to detect all major QTLs.

Global Analysis of ATM Polymorphism Reveals Significant Functional Constraint

ATM, the gene that is mutated in ataxia-telangiectasia, is associated with cerebellar degeneration, abnormal proliferation of small blood vessels, and cancer. These clinically important manifestations have stimulated interest in defining the sequence variation in the ATM gene. Therefore, we undertook a comprehensive survey of sequence variation in ATM in diverse human populations. The protein-encoding exons of the gene (9,168 bp) and the adjacent intron and untranslated sequences (14,661 bp) were analyzed in 93 individuals from seven major human populations. In addition, the coding sequence was analyzed in one chimpanzee, one gorilla, one orangutan, and one Old World monkey. In human ATM, 88 variant sites were discovered by denaturing high-performance liquid chromatography, which is 96%-100% sensitive for detection of DNA sequence variation. ATM was compared to 14 other autosomal genes for nucleotide diversity. The noncoding regions of ATM had diversity values comparable to other genes, but the coding regions had very low diversity, especially in the last 29% of the protein sequence. A test of the neutral evolution hypothesis, through use of the Hudson/Kreitman/Aguadé statistic, revealed that this region of the human ATM gene was significantly constrained relative to that of the orangutan, the Old World monkey, and the mouse, but not relative to that of the chimpanzee or the gorilla. ATM displayed extensive linkage disequilibrium, consistent with suppression of meiotic recombination at this locus. Seven haplotypes were defined. Two haplotypes accounted for 82% of all chromosomes analyzed in all major populations; two others carrying the same D126E missense polymorphism accounted for 33% of chromosomes in Africa but were never observed outside of Africa. The high frequency of this polymorphism may be due either to a population expansion within Africa or to selective pressure.

Sequence diversity and linkage disequilibrium within the merozoite surface protein-1 (Msp-1) locus of Plasmodium falciparum: a longitudinal study in Brazil.

The merozoite surface protein-1 (MSP-1) is a major vaccine candidate for the asexual blood stage of malaria. We examined both the extent of sequence diversity in block 17, the 3' end of Msp-1 gene coding for a 19-kDa polypeptide (MSP-1(19)) putatively involved in red blood cell binding, and the patterns of linkage disequilibrium between polymorphic sites throughout the Msp-1 locus. The parasite population sample consisted of Plasmodium falciparum isolates collected between 1985 and 1998 in Rondĵnia, an area of hypoendemic malaria transmission in the southwestern Brazilian Amazon. Results were summarized as follows. (1) Seven block-17 sequence variants or haplotypes were found among 130 isolates, including two new haplotypes (novel combinations of previously reported amino acid replacements), here named Brazil-1 (E-TSR-F) and Brazil-2 (Q-TSR-F). (2) As previously shown for other Msp-1 polymorphisms, frequencies of block-17 haplotypes displayed significant temporal variation. (3) Extensive linkage disequilibrium was demonstrated between neighboring dimorphic sites within block 17, as well as between polymorphisms at the 5' and 3' ends of Msp-1 (map distance range: 3.83-4.99 kb). (4) The overall patterns of linkage disequilibrium within Msp-1 remained stable over a period of nearly one decade, and examples of possible 'epidemic' expansion of parasites carrying particular Msp-1 alleles were found in the 1980s and 1990s. These results are discussed in relation to the population biology of P. falciparum and the development of malaria vaccines based on MSP-1.

Allelic association with SNPs: metrics, populations, and the linkage disequilibrium map.

Comparison of different metrics, using three large samples of haplotypes from different populations, demonstrates that rho is the most efficient measure of association between pairs of single nucleotide polymorphisms (SNPs). Pairwise data can be modeled, using composite likelihood, to describe the decline in linkage disequilibrium with distance (the Malecot model). The evidence from more isolated populations (Finland, Sardinia) suggests that linkage disequilibrium extends to 427-893 kb but, even in samples representative of large heterogeneous populations, such as CEPH, the extent is 385 kb or greater. This suggests that isolated populations are not essential for linkage disequilibrium mapping of common diseases with SNPs. The in parameter of the Malecot model (recombination and time), evaluated at each SNP, indicates regions of the genome with extensive and less extensive disequilibrium (low and high values of in respectively). When plotted against the physical map, the regions with extensive and less extensive linkage disequilibrium may correspond to recombination cold and hot spots. This is discussed in relation to the Xq25 cytogenetic band and the HFE gene region.

Haplotypes at ATM Identify Coding-Sequence Variation and Indicate a Region of Extensive Linkage Disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.

Genetic studies: look no further than your own backyard

Eaves, I.A. et al. (2000) The genetically isolated populations of Finland and Sardinia may not be a panacea for linkage disequilibrium mapping of common disease genes. Nat. Genet. 25, 320–323 Taillon-Miller, P. et al. (2000) Juxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28. Nat. Genet. 25, 324–328

Genetic studies: look no further than your own backyard

The feasibility of using single nucleotide polymorphisms and other types of genetic variation for disease-association and pharmacogenomic studies could be greatly enhanced if favorable levels of linkage disequilibrium (LD) existed in the population under study. LD, or correlation between nearby genetic markers, allows the survey of random genetic variants in the genome instead of having to find and evaluate the causative variant for a disease or drug response directly. Genetically isolated populations have long been coveted by geneticists because they are thought to be more genetically homogeneous, and because they have LD that extends over longer regions of the genome than mixed general populations, like those in the US and Europe. However, recent evidence suggests that these isolated populations might have fewer advantages than expected.Much of the early success in human genetics was the result of using large families from homogeneous isolated populations to uncover genes for monogenic disorders. As the field turns to the study of complex diseases, a number of challenges must be overcome including the access to large patient populations for association studies. Whereas isolated populations are thought to offer biological and logistical advantages, they rarely have the large numbers of patients needed to carry out large-scale disease-association studies and often do not participate in early clinical trials for pharmacogenomics. Furthermore, specific genetic variants underlying disease in an isolated population might only partially contribute to the disease or drug response in larger, more-outbred populations.Two recent papers in Nature Genetics demonstrate that genetically isolated populations might not be as advantageous as previously thought for studying complex traits1xThe genetically isolated populations of Finland and Sardinia may not be a panacea for linkage disequilibrium mapping of common disease genes. Eaves, I.A. et al. Nat. Genet. 2000; 25: 320–323Crossref | PubMed | Scopus (134)See all References, 2xJuxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28. Taillon-Miller, P. et al. Nat. Genet. 2000; 25: 324–328Crossref | PubMed | Scopus (208)See all References. Eaves et al. examined the extent of LD on chromosome 18q21 in two isolated populations (Finland and Sardinia). LD between 21 different polymorphic microsatellite markers was evaluated and compared with LD found in two mixed populations (US and UK). Neither isolate displayed a greater level of LD compared with the mixed populations. Taillon-Miller et al. observed no differences in the extent or general location of LD between 39 SNPs on the X chromosome when comparing Finnish and Sardinian populations with an outbred European population. Both of these studies provide the first empirical support of similar previous claims based on simulation studies.While these findings are intriguing, caution should be observed in generalizing the results to other isolated populations and other chromosomal regions. Nonetheless, they give pause to re-evaluate the need to conduct genetic association studies of complex traits in specialized populations and open the door for large-scale studies in more-readily accessible populations. No doubt, the challenges of finding associations are still significant, especially for those genes that lie in regions of low LD. However, now we might not need to look any further than our own backyards.

Natural selection mapping of the warfarin-resistance gene.

In theory, genes under natural selection can be revealed by unique patterns of linkage disequilibrium (LD) and polymorphism at physically linked loci. However, given the effects of recombination and mutation, the physical extent and persistence of LD patterns in natural populations is uncertain. To assess the LD signature of selection, we survey variation in 26 microsatellite loci spanning an approximately 32-cM region that includes the warfarin-resistance gene (Rw) in five wild rat populations having resistance levels between 0 and 95%. We find a high frequency of heterozygote deficiency at microsatellite loci in resistant populations, and a negative association between gene diversity (H) and resistance. Contrary to previous studies, these data suggest that directional rather than overdominant selection may predominate during periods of intense anticoagulant treatment. In highly resistant populations, extensive LD was observed over a chromosome segment spanning approximately 14% of rat chromosome 1. In contrast, LD in a moderately resistant population was more localized and, in conjunction with likelihood ratios, allowed assignment of Rw to a 2. 2-cM interval. Within this genomic window, a diagnostic marker, D1Rat219, assigned 91% of rats to the correct resistance category. These results further demonstrate that "natural selection mapping" in field populations can detect and map major fitness-related genes, and question overdominance as the predominant mode of selection in anticoagulant-resistant rat populations.

Juxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28

Linkage disequilibrium (LD), or the non-random association of alleles, is poorly understood in the human genome. Population genetic theory suggests that LD is determined by the age of the markers, population history, recombination rate, selection and genetic drift. Despite the uncertainties in determining the relative contributions of these factors, some groups have argued that LD is a simple function of distance between markers. Disease-gene mapping studies and a simulation study gave differing predictions on the degree of LD in isolated and general populations. In view of the discrepancies between theory and experimental observations, we constructed a high-density SNP map of the Xq25-Xq28 region and analysed the male genotypes and haplotypes across this region for LD in three populations. The populations included an outbred European sample (CEPH males) and isolated population samples from Finland and Sardinia. We found two extended regions of strong LD bracketed by regions with no evidence for LD in all three samples. Haplotype analysis showed a paucity of haplotypes in regions of strong LD. Our results suggest that, in this region of the X chromosome, LD is not a monotonic function of the distance between markers, but is more a property of the particular location in the human genome.

Extensive genome-wide linkage disequilibrium in cattle.

A genome-wide linkage disequilibrium (LD) map was generated using microsatellite genotypes (284 autosomal microsatellite loci) of 581 gametes sampled from the dutch black-and-white dairy cattle population. LD was measured between all marker pairs, both syntenic and nonsyntenic. Analysis of syntenic pairs revealed surprisingly high levels of LD that, although more pronounced for closely linked marker pairs, extended over several tens of centimorgan. In addition, significant gametic associations were also shown to be very common between nonsyntenic loci. Simulations using the known genealogies of the studied sample indicate that random drift alone is likely to account for most of the observed disequilibrium. No clear evidence was obtained for a direct effect of selection ("Bulmer effect"). The observation of long range disequilibrium between syntenic loci using low-density marker maps indicates that LD mapping has the potential to be very effective in livestock populations. The frequent occurrence of gametic associations between nonsyntenic loci, however, encourages the combined use of linkage and linkage disequilibrium methods to avoid false positive results when mapping genes in livestock.

Genetic Structure of theBovicola ovis(Mallophaga: Trichodectidae) in Southwestern Australia

The sheep biting louse, Bovicola ovis (Schrank), is an economically important, worldwide ectoparasite of sheep. In Australia, up to 30% of sheep flocks are infested with lice. The usual method of control has been synthetic pyrethroids applied as pour-on along the back of the sheep, but treatment failures have become widespread since 1985 because of the development of resistance. We used allozyme markers to study the distribution of genetic variation within and among populations of lice on different farms in Western Australia. Genetic variation within populations was similar to previously reported values for other ectoparasitic arthropods. Heterozygote deficiencies were found at 1 locus in a number of population and another 2 loci in 1 other population. However, another variable locus conformed to Hardy-Weinberg equilibrium and there was little evidence of extensive linkage disequilibrium between loci. Further studies are necessary to establish the breeding system. Genetic differences among populations were not related to geographic separation, which is consistent with an island model of population structure. A small but significant proportion (2.8%) of the total genetic variation was distributed among populations' equivalent under the island model to a gene flow of 8.7 individuals exchanged per generation. The implications of this result are discussed in terms of controlling and managing synthetic pyrethroid resistance in sheep biting lice.

HLA typing in the United Kingdom multiple sclerosis genome screen.

The United Kingdom multiple sclerosis genome screen demonstrated a peak maximum lod score of 2.8 in the HLA region, together with statistically significant excess transmission of the 121-base pair (bp) allele of the tumour necrosis factor-a marker. In order to determine whether this association is independent of the established HLA association, or simply a consequence of the 121-bp allele being part of the same haplotype, we HLA-DR and -DQ typed the 227 sibling-pair families used in the original screen. The expected associations of multiple sclerosis with the DR15 (p=8.7E-18), DQ6 (p=2.0E-09) and DR51 (p=2.8E-16) phenotypes were confirmed, and excess transmission of the DRB1*1501 and DQB1*0602 alleles was demonstrated. Combining HLA typing with the original microsatellite data demonstrated extensive linkage disequilibrium between the 121-bp allele and the 1501-0602 haplotype. Outside this extended haplotype (121-1501-0602), none of the alleles demonstrated significant transmission distortion. Having established the importance of this extended haplotype, we reanalysed the entire genome screen data after excluding those sibling pairs sharing the extended haplotype (n=27). Conditioning the full genome screen data on the basis of identity by state sharing showed that some potential linkage regions identified in the original screen clustered in families, in which the extended haplotype was shared (1p, 2p and 17q), whereas others grouped with those in which it was not (5cen, 7p and Xq). This suggests complexity in the genetics of multiple sclerosis.

Haplotype Structure and Population Genetic Inferences from Nucleotide-Sequence Variation in Human Lipoprotein Lipase

Allelic variation in 9.7 kb of genomic DNA sequence from the human lipoprotein lipase gene (LPL) was scored in 71 healthy individuals (142 chromosomes) from three populations: African Americans (24) from Jackson, MS; Finns (24) from North Karelia, Finland; and non-Hispanic Whites (23) from Rochester, MN. The sequences had a total of 88 variable sites, with a nucleotide diversity (site-specific heterozygosity) of .002+/-.001 across this 9.7-kb region. The frequency spectrum of nucleotide variation exhibited a slight excess of heterozygosity, but, in general, the data fit expectations of the infinite-sites model of mutation and genetic drift. Allele-specific PCR helped resolve linkage phases, and a total of 88 distinct haplotypes were identified. For 1,410 (64%) of the 2,211 site pairs, all four possible gametes were present in these haplotypes, reflecting a rich history of past recombination. Despite the strong evidence for recombination, extensive linkage disequilibrium was observed. The number of haplotypes generally is much greater than the number expected under the infinite-sites model, but there was sufficient multisite linkage disequilibrium to reveal two major clades, which appear to be very old. Variation in this region of LPL may depart from the variation expected under a simple, neutral model, owing to complex historical patterns of population founding, drift, selection, and recombination. These data suggest that the design and interpretation of disease-association studies may not be as straightforward as often is assumed.

Genetic Structure of IDDM1: Two Separate Regions in the Major Histocompatibility Complex Contribute to Susceptibility or Protection

We analyzed 11 markers in the IDDM1 region in 120 IDDM patients and 83 healthy control subjects who were fully matched for the highest risk HLADQA1*O3O1-DQB1*O3O2/DQA1*O5O1-DQB1*O2O1 genotype. Our study provides strong evidence that two regions in the major histocompatibility complex contribute to IDDM susceptibility or protection. First, despite selection for highest IDDM-associated risk DQ genotypes, this region displays extensive linkage disequilibrium (LD) differences between IDDM patients and control subjects. A second critical region was mapped around the microsatellite locus D6S273 centromeric of TNF, and it is ∼200 kb in size. LD analysis shows that “diabetogenic haplotypes” may have resulted from a recombination telomeric of D6S1014 in the region of D6S273 and TNFa. Haplotype analysis using HLA and microsatellite loci refines IDDM risk assessment in carriers of the HLA-DQ highest risk genotype.

Genetic structure of IDDM1: two separate regions in the major histocompatibility complex contribute to susceptibility or protection. Belgian Diabetes Registry

We analyzed 11 markers in the IDDM1 region in 120 IDDM patients and 83 healthy control subjects who were fully matched for the highest risk HLA-DQA1*0301-DQB1 *0302/DQA1*0501-DQB1*0201 genotype. Our study provides strong evidence that two regions in the major histocompatibility complex contribute to IDDM susceptibility or protection. First, despite selection for highest IDDM-associated risk DQ genotypes, this region displays extensive linkage disequilibrium (LD) differences between IDDM patients and control subjects. A second critical region was mapped around the microsatellite locus D6S273 centromeric of TNF, and it is approximately 200 kb in size. LD analysis shows that "diabetogenic haplotypes" may have resulted from a recombination telomeric of D6S1014 in the region of D6S273 and TNFa. Haplotype analysis using HLA and microsatellite loci refines IDDM risk assessment in carriers of the HLA-DQ highest risk genotype.

Hereditary haemochromatosis as an immunological disease.

Hereditary haemochromatosis as an immunological disease.