Linkage disequilibrium between mutation and RFLP haplotype at the phenylalanine hydroxylase locus in the German population

Summary Restriction fragment length polymorphism (RFLP) haplotypes at the phenylalanine hydroxylase (PAH) locus have been determined in 60 German families with PAH deficiency. Similar to the Danish population, about 90% of the mutant alleles are confined to four distinct haplotypes. There are however, differences in the frequency distributiion of these haplotypes among the mutant alleles between the two populations. Using an oligonucleotide probe for the splicing mutation associated with mutant haplotype 3 in the Danish population, a tight association between the mutation and the RFLP haplotype has also been observed in Germany. The results provide strong evidence that the splicing mutation occurred on a haplotype 3 chromosome and that the mutant allele has spread into different populations smong Caucasians.     

Linkage disequilibrium and haplotype analysis in German Friedreich ataxia families

The main mutation causing Friedreich ataxia (FRDA) is the expansion of a GAA repeat localized within the intron between exon 1 and exon 2 of the gene X25. This expansion has been observed in 98% of FRDA chromosomes. To analyze frequencies of markers tightly linked to the Friedreich ataxia gene and to investigate wheter a limited number of ancestral chromosomes are shared by German FRDA families, a detailed analysis employing nine polymorphic markers was performed. We found strong linkage disequilibria and association of FRDA expansions with a few haplotypes. FRDA haplotypes differ significantly from control haplotypes. Our results confirm that GAA repeat expansions in intron 1 of the frataxin gene are limited to a few chromosomes and indicate an obvious founder effect in German patients. Based on these analyses, we estimate a minimum age of the mutation of 107 generations.     

Linkage disequilibrium decay and haplotype block structure in the pig

Linkage disequilibrium (LD) may reveal much about domestication and breed history. An investigation was conducted, to analyze the extent of LD, haploblock partitioning, and haplotype diversity within haploblocks across several pig breeds from China and Europe and in European wild boar. In total, 371 single-nucleotide-polymorphisms located in three genomic regions were genotyped. The extent of LD differed significantly between European and Chinese breeds, extending up to 2 cM in Europe and up to 0.05 cM in China. In European breeds, LD extended over large haploblocks up to 400 kb, whereas in Chinese breeds the extent of LD was smaller and generally did not exceed 10 kb. The European wild boar showed an intermediate level of LD between Chinese and European breeds. In Europe, the extent of LD also differed according to genomic region. Chinese breeds showed a higher level of haplotype diversity and shared high levels of frequent haplotypes with Large White, Landrace, and Duroc. The extent of LD differs between both centers of pig domestication, being higher in Europe. Two hypotheses can explain these findings. First, the European ancestral stock had a higher level of LD. Second, modern breeding programs increased the extent of LD in Europe and caused differences of LD between genomic regions. Large White, Landrace, and Duroc showed evidence of past introgression from Chinese breeds.     

Linkage disequilibrium between phenylketonuria and RFLP haplotype 1 at the phenylalanine hydroxylase locus in Portugal

Summary RFLPs of 36 normal and 41 mutant alleles at the phenylalanine hydroxylase locus were determined in 31 Portuguese kindreds. A total of 14 haplotypes including 10 normal and 7 mutant alleles were observed. Almost 75% of all mutant alleles were confined within only two haplotypes, namely haplotype 9 (17.1%) and haplotype 1 (56.1%). This frequency of mutant haplotype 1 in Portugal is, to our knowledge, the highest for this mutant haplotype in all studies reported to date. Other mutant haplotypes were either rare (haplotype 2, 9.7%) or totally absent (haplotype 3, 0%). Only 24.5% of all mutant alleles were found to consistently carry identified mutations, particularly R261Q (9.8%), R252W (3.3%), R408W (1.6%) and I94 (3.3%). A new mutation, L249F, located in the seventh exon of the gene, accounted for 6.5% of all mutant alleles in our series. Interestingly, this mutant genotype was consistently associated with mutant haplotype 1 (P<0.01), as also observed for the R261Q mutation. It appears, therefore, that mutant haplotype 1 is genotypically heterogeneous in Portugal and that more than two mutations account for its prevalence in this country.     

Linkage disequilibrium in a population undergoing periodic fragmentation and admixture

Glacial and interglacial cycles are considered to have caused the fragmentation and admixture of populations in many organisms. A simple model incorporating such periodic changes of the population structure is analysed in order to investigate the behaviour of neutral genetic variation at one and two loci. The equilibrium is reached very quickly in terms of cycles if the length of a cycle is long, as would be expected of the glaciation cycles. Heterozygosity and linkage disequilibrium are shown to depend on the length of time of the fragmented and admixed phases, population sizes, and number (n) of subpopulations in the fragmented phase. If the population size is small in the fragmented phase and its duration is long, the squared correlation coefficient of two loci (a measure of linkage disequilibrium) just after the admixture is approximated by 1/(n-1) for n > 1. After admixture, the correlation decays at a rate of approximately twice the recombination rate. Therefore, if post-glaciation admixture created linkage disequilibrium, we expect to observe linkage disequilibrium even between moderately linked loci, and its decay pattern along the chromosome is very different from that in a random mating population at equilibrium. This is especially true in organisms with long generation times such as trees.     

Linkage disequilibrium and haplotype mapping of a skin cancer susceptibility locus in outbred mice

Car-R (carcinogenesis-resistant) and Car-S (carcinogenesis-susceptible) outbred mice, obtained by phenotypic selection from an initial intercross of eight inbred strains, show a >100-fold difference in their susceptibility to two-stage skin tumorigenesis. We found that the lines carry a high degree of genetic polymorphism, with an average heterozygosity of 0.39. This polymorphism allowed the use of linkage disequilibrium (LD) and haplotype analysis for the mapping of a skin cancer modifier locus on Chr 7, in a short region of 6 cM, around the Tyr gene. Car-S mice inherited the susceptibility allele at this locus from the A/J, BALB/c, SJL/J, and SWR/J strains. Our results demonstrate the feasibility and usefulness of mapping disease genes by LD in phenotypically selected, genetically heterogeneous animals. Received: 16 March 2000 / Accepted: 9 June 2000     

Linkage disequilibrium in a finite population that is partially selfing

The linkage disequilibrium expected in a finite, partially selfing population is analyzed, assuming the infinite allele model. Formulas for the expected sum of squares of the linkage disequilibria and the squared standard linkage disequilibrium are derived from the equilibrium values of sixteen inbreeding coefficients required to describe the behavior of the system. These formulas are identical to those obtained with random mating if the effective population size N_e = (1-½S)N and the effective recombination value r_e = (1-S)r/(1-½S), where S is the proportion of selfing, are substituted for the population size and the recombination value. Therefore, the effect of partial selfing at equilibrium is to reduce the population size by a factor 1-½S and the recombination value by a factor (1-S)/(1-½S).     

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium ( r ²   >   0.8) in genomic regions of ∼50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26–113 kb by various definitions, which was larger than that observed in humans (∼10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.     

Linkage disequilibrium under skewed offspring distribution among individuals in a population

Correlations in coalescence times between two loci are derived under selectively neutral population models in which the offspring of an individual can number on the order of the population size. The correlations depend on the rates of recombination and random drift and are shown to be functions of the parameters controlling the size and frequency of these large reproduction events. Since a prediction of linkage disequilibrium can be written in terms of correlations in coalescence times, it follows that the prediction of linkage disequilibrium is a function not only of the rate of recombination but also of the reproduction parameters. Low linkage disequilibrium is predicted if the offspring of a single individual frequently replace almost the entire population. However, high linkage disequilibrium can be predicted if the offspring of a single individual replace an intermediate fraction of the population. In some cases the model reproduces the standard Wright-Fisher predictions. Contrary to common intuition, high linkage disequilibrium can be predicted despite frequent recombination, and low linkage disequilibrium under infrequent recombination. Simulations support the analytical results but show that the variance of linkage disequilibrium is very large.     

Linkage disequilibrium analyses of natriuretic peptide precursor B locus reveal risk haplotype conferring high plasma BNP levels

BACKGROUND: Brain natriuretic peptide (BNP) has been widely used for the diagnosis and prognostic evaluation of chronic heart failure (CHF). In the present study, we performed association study of single nucleotide polymorphisms (SNPs) surrounding the natriuretic peptide precursor B (NPPB) gene with plasma BNP levels in 2970 adult Japanese. METHODS AND RESULTS: Association analysis between SNPs of the NPPB gene and plasma BNP revealed significant associations of the 8 SNPs surrounding the entire NPPB gene with plasma BNP levels. For instance, as to SNP rs198389 (T-381C), plasma BNP levels among the three genotypic categories, i.e., 2189 homozygous T-allele carriers (BNP 26.4+/-0.6pg/ml), 697 heterozygous carriers (35.0+/-1.1pg/ml), and 52 homozygous C-allele carriers (46.0+/-4.1pg/ml) indicated a co-dominant effect of the minor C-allele on elevating plasma BNP levels (P<0.0001). Linkage disequilibrium (LD) analysis among the 8 SNPs revealed that the region consisted of two, 5' major and 3' minor, LD blocks. Haplotype-based association analysis demonstrated that plasma BNP levels were associated closely with the haplotypes-1 and -2 of the major LD block. CONCLUSION: These results suggest that genetic variation at the primary locus NPPB gene, represented by definition of risk haplotypes, may be an important determinant of plasma BNP levels.     

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 and haplotype analysis among Polish families with spinal muscular atrophy.

Spinal Muscular Atrophy (SMA) is an inherited degenerative disorder of anterior horn cells that results in progressive muscle weakness and atrophy. The autosomal recessive forms of childhood-onset SMA have been mapped to chromosome 5q11.2-13.3, in a number of studies examining different populations. A total of 9 simple sequence repeat markers were genotyped against 32 Polish families with SMA. The markers span an approximately 0.7 cM region defined by the SMA flanking markers D5S435 and MAP1B. Significant linkage disequilibrium (corrected P < .05) was detected at four of these markers, with D/Dmax values of < or = .89. Extended haplotype analysis revealed a predominant haplotype associated with SMA. The apparently high mutation rate of some of the markers has resulted in a number of haplotypes that vary slightly from this predominant haplotype. The predominant haplotype and these closely related patterns represent 25% of the disease chromosomes and none of the nontransmitted parental chromosomes. This predominant haplotype is present both in patients with acute (type I) and in chronic (types II and III) forms of SMA and occurs twice in a homozygous state, both times in children with chronic SMA.     

Linkage disequilibrium and historical effective population size in the Thoroughbred horse

Many genomic methodologies rely on the presence and extent of linkage disequilibrium (LD) between markers and genetic variants underlying traits of interest, but the extent of LD in the horse has yet to be comprehensively characterized. In this study, we evaluate the extent and decay of LD in a sample of 817 Thoroughbreds. Horses were genotyped for over 50,000 single nucleotide polymorphism (SNP) markers across the genome, with 34,848 autosomal SNPs used in the final analysis. Linkage disequilibrium, as measured by the squared correlation coefficient (r(2)), was found to be relatively high between closely linked markers (>0.6 at 5 kb) and to extend over long distances, with average r(2) maintained above non-syntenic levels for single nucleotide polymorphisms (SNPs) up to 20 Mb apart. Using formulae which relate expected LD to effective population size (N(e)), and assuming a constant actual population size, N(e) was estimated to be 100 in our population. Values of historical N(e), calculated assuming linear population growth, suggested a decrease in N(e) since the distant past, reaching a minimum twenty generations ago, followed by a subsequent increase until the present time. The qualitative trends observed in N(e) can be rationalized by current knowledge of the history of the Thoroughbred breed, and inbreeding statistics obtained from published pedigree analyses are in agreement with observed values of N(e). Given the high LD observed and the small estimated N(e), genomic methodologies such as genomic selection could feasibly be applied to this population using the existing SNP marker set.     

Linkage disequilibrium analysis in the genetically isolated Norfolk Island population

Norfolk Island is a human genetic isolate, possessing unique population characteristics that could be utilized for complex disease gene localization. Our intention was to evaluate the extent and strength of linkage disequilibrium (LD) in the Norfolk isolate by investigating markers within Xq13.3 and the NOS2A gene encoding the inducible nitric oxide synthase. A total of six microsatellite markers spanning approximately 11 Mb were assessed on chromosome Xq13.3 in a group of 56 men from Norfolk Island. Additionally, three single nucleotide polymorphisms (SNPs) localizing to the NOS2A gene were analyzed in a subset of the complex Norfolk pedigree. With the exception of two of the marker pairs, one of which is the most distantly spaced marker, all the Xq13.3 marker pairs were found to be in significant LD indicating that LD extends up to 9.5-11.5 Mb in the Norfolk Island population. Also, all SNPs studied showed significant LD in both Norfolk Islanders and Australian Caucasians, with two of the marker pairs in complete LD in the Norfolk population only. The Norfolk Island study population possesses a unique set of characteristics including founder effect, geographical isolation, exhaustive genealogical information and phenotypic data of use to cardiovascular disease risk traits. With LD extending up to 9.5-11 Mb, the Norfolk isolate should be a powerful resource for the localization of complex disease genes.     

Linkage disequilibrium, selection and recombination at three Loci

Limits to the relationship among linkage disequilibrium, selection and recombination at equilibrium in three-locus, two-allele, deterministic, discrete generation models are determined using linear programming techniques. These results show that the commonly used measures of linkage disequilibrium are not appropriate for a multilocus setting. Additionally, interactions among three loci are important in reducing the strength of selection necessary to maintain a given level of disequilibrium, relative to a two-locus model.     

Deep haplotype divergence and long-range linkage disequilibrium at xp21.1 provide evidence that humans descend from a structured ancestral population

Fossil evidence links human ancestry with populations that evolved from modern gracile morphology in Africa 130,000-160,000 years ago. Yet fossils alone do not provide clear answers to the question of whether the ancestors of all modern Homo sapiens comprised a single African population or an amalgamation of distinct archaic populations. DNA sequence data have consistently supported a single-origin model in which anatomically modern Africans expanded and completely replaced all other archaic hominin populations. Aided by a novel experimental design, we present the first genetic evidence that statistically rejects the null hypothesis that our species descends from a single, historically panmictic population. In a global sample of 42 X chromosomes, two African individuals carry a lineage of noncoding 17.5-kb sequence that has survived for >1 million years without any clear traces of ongoing recombination with other lineages at this locus. These patterns of deep haplotype divergence and long-range linkage disequilibrium are best explained by a prolonged period of ancestral population subdivision followed by relatively recent interbreeding. This inference supports human evolution models that incorporate admixture between divergent African branches of the genus Homo.     

SNP identification, linkage disequilibrium, and haplotype analysis for a 200-kb genomic region in a Korean population

Understanding patterns of linkage disequilibrium (LD) across genomes may facilitate association mapping studies to localize genetic variants influencing complex diseases, a recognition that led to the International Haplotype Mapping Project (HapMap). Divergent patterns of haplotype frequency and LD across global populations require that the HapMap database be supplemented with haplotype and LD data from additional populations. We conducted a pilot study of the LD and haplotype structure of a genomic region in a Korean population. A total of 165 SNPs were identified in a 200-kb region of 22q13.2 by direct sequencing. Unphased genotype data were generated for 76 SNPs in 90 unrelated Korean individuals. LD, haplotype diversity, and recombination rates were assessed in this region and compared with the HapMap database. The pattern of LD and haplotype frequencies of Korean samples showed a high degree of similarity with Japanese data. There was a strong correlation between high LD and low recombination frequency in this region. We found considerable similarities in local LD patterns between three Asian populations (Han Chinese, Japanese, and Korean) and the CEPH population. Haplotype frequencies were, however, significantly different between them. Our results should further the understanding of distinctive Korean genomic features and assist in designing appropriate association studies.     

Inferring coancestry in population samples in the presence of linkage disequilibrium

In both pedigree linkage studies and in population-based association studies there has been much interest in the use of modern dense genetic marker data to infer segments of gene identity by descent (ibd) among individuals not known to be related, to increase power and resolution in localizing genes affecting complex traits. In this article, we present a hidden Markov model (HMM) for ibd among a set of chromosomes and describe methods and software for inference of ibd among the four chromosomes of pairs of individuals, using either phased (haplotypic) or unphased (genotypic) data. The model allows for missing data and typing error, but does not model linkage disequilibrium (LD), because fitting an accurate LD model requires large samples from well-studied populations. However, LD remains a major confounding factor, since LD is itself a reflection of coancestry at the population level. To study the impact of LD, we have developed a novel simulation approach to generate realistic dense marker data for the same set of markers but at varying levels of LD. Using this approach, we present results of a study of the impact of LD on the sensitivity and specificity of our HMM model in estimating segments of ibd among sets of four chromosomes and between genotype pairs. We show that, despite not incorporating LD, our model has been quite successful in detecting segments as small as 10(6) bp (1 Mpb); we present also comparisons with fastIBD which uses an LD model in estimating ibd.