Evaluation of linkage disequilibrium measures between multi-allelic markers as predictors of linkage disequilibrium between markers and QTL

Effectiveness of marker-assisted selection (MAS) and quantitative trait loci (QTL) mapping using population-wide linkage disequilibrium (LD) between markers and QTL depends on the extent of LD and how it declines with distance in a population. Because marker-QTL LD cannot be observed directly, the objective of this study was to evaluate alternative measures of observable LD between multi-allelic markers as predictors of usable LD of multi-allelic markers with presumed biallelic QTL. Observable LD between marker pairs was evaluated using eight existing measures and one new measure. These consisted of two pooled and standardized measures of LD between pairs of alleles at two markers based on Lewontin's LD measure, two pooled measures of squared correlations between alleles, one standardized measure using Hardy-Weinberg heterozygosities, and four measures based on the chi-square statistic for testing for association between alleles at two loci. In simulated populations with a range of LD generated by drift and a range of marker polymorphism, marker-marker LD measured by a standardized chi-square statistic (denoted chi(2')) was found to be the best predictor of useable marker-QTL LD for a group of multi-allelic markers. Estimates of the level and decline of marker-marker LD with distance obtained from chi(2') were linearly and highly correlated with usable LD of those markers with QTL across population structures and marker polymorphism. Corresponding relationships were poorer for the other marker-marker LD measures. Therefore, when LD is generated by drift, chi(2') is recommended to quantify the amount and extent of usable LD in a population for QTL mapping and MAS based on multi-allelic markers.     

A quantitative-trait analysis of human plasma-dopamine beta-hydroxylase activity: evidence for a major functional polymorphism at the DBH locus

Dopamine-beta-hydroxylase (D beta H) catalyzes the conversion of dopamine to norepinephrine and is released from sympathetic neurons into the circulation. Plasma-D beta H activity varies widely between individuals, and a subgroup of the population has very low activity levels. Mounting evidence suggests that the DBH structural gene is itself the major quantitative-trait locus (QTL) for plasma-D beta H activity, and a single unidentified polymorphism may account for a majority of the variation in activity levels. Through use of both sequencing-based mutational analysis of extreme phenotypes and genotype/phenotype correlations in samples from African American, European American (EA), and Japanese populations, we have identified a novel polymorphism (--1021C-->T), in the 5' flanking region of the DBH gene, that accounts for 35%--52% of the variation in plasma-D beta H activity in these populations. In EAs, homozygosity at the T allele predicted the very low D beta H-activity trait, and activity values in heterozygotes formed an intermediate distribution, indicating codominant inheritance. Our findings demonstrate that --1021C-->T is a major genetic marker for plasma-D beta H activity and provide new tools for investigation of the role of both D beta H and the DBH gene in human disease.     

Genetic diversity and the structure of linkage disequilibrium in the methylenetetrahydrofolate reductase locus

Investigation of linkage disequilibrium block architecture in human genome is modern, intensely investigated field of molecular genetics. In the present study, genetic differentiation and linkage disequilibrium pattern in the methylenetetrahydrofolate reductase (MTHFR) locus was examined in the populations of Russians, Tuvinians, and Northern and Southern Kyrgyzes. Methylenetetrahydrofolate reductase is the key enzyme of folate cycle, responsible for reduction of 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate. Decreased enzymatic activity of this protein often caused by certain associations of MTHFR alleles results in the increased plasma homocysteine levels. In the population groups examined, genotype and allele frequencies at five polymorphic MTHFR loci: rs17037397, rs4846052, rs1801133, rs1801131, and rs1537516 were evaluated. Statistically significant genetic differences between the population group of Southern Kyrgyzes and the other groups, as well as between Russians and Tuvinians, were demonstrated. In the MTHFR gene from the population of Southern Kyrgyzes one block was revealed; in the populations of Russians, Tuvinians, and Northern Kyrgyzes two blocks were detected. Thus, the structure of linkage disequilibrium in the MTHFR locus demonstrated population-specific pattern.     

Segregation and linkage studies of plasma dopamine-beta-hydroxylase (DBH), erythrocyte catechol-O-methyltransferase (COMT), and platelet monoamine oxidase (MAO): possible linkage between the ABO locus and a gene controlling DBH activity.

Measurements of dopamine-beta-hydroxylase (DBH), catechol-O-methyltransferase (COMT), and monoamine oxidase (MAO) along with 27 polymorphic marker phenotypes were available for 162 patients with major affective disorders and 1,125 of their relatives. Levels of enzymes were previously found not to be associated with illness. Pedigree analysis methods for quantitative traits are used to test single-gene hypotheses for segregation of DBH in 32 families with 411 individuals. COMT in 30 families with 351 individuals, and MAO in 50 families with 309 individuals. The familial distribution of both DBH and COMT are consistent with two codominant alleles at the same locus that account for 56% and 59% of the total variance, respectively. MAO activity cannot be shown to be segregating as a single major gene, but a purely nongenetic hypothesis is also rejected. A possible linkage of a locus for DBH to the ABO locus is indicated by a maximum lod score of 1.82 at 0% and 10% recombination fractions for males and females, respectively. A lod score of 0.61 at 0% recombination for a similar analysis in a single large pedigree was reported by Elston et al., making the combined lod score for the two studies equal to 2.32 at 0% recombination.     

The detection of linkage disequilibrium between closely linked markers: RFLPs at the AI-CIII apolipoprotein genes.

Study of very closely linked DNA variants at various loci has frequently shown linkage disequilibrium. We studied three closely linked RFLPs at the apolipoprotein AI-CIII locus. Two variants detected by MspI and SstI were in strong linkage disequilibrium; but when conventional statistical tests were used, a third variant (PstI), located between the MspI and SstI markers, appeared to be in linkage equilibrium with these two "outside" markers. Similar discrepancies from the expected monotone relationship between physical distance and linkage disequilibrium have been reported by others. To investigate these discrepancies, the power to detect linkage disequilibrium was calculated. It could be shown that, for the gene frequencies encountered, very large sample sizes would be required to demonstrate negative (i.e., repulsion-phase) linkage disequilibrium. Such numbers are usually very difficult to attain in human studies. Failure to demonstrate linkage disequilibrium by conventional methods therefore does not imply its absence. Appropriate nomograms and tables are provided.     

Haplotype distribution of and linkage disequilibrium between four polymorphic markers near the CFTR locus in Brazilian cystic fibrosis patients

To contribute to a better understanding of the origin and distribution of CFTR mutations in the Brazilian population, we have investigated the linkage between four polymorphic markers (XV2c, KM19, GATT, and TUB9) within or near the CFTR locus. The distribution of alleles for each polymorphism for both parental and cystic fibrosis (CF) chromosomes from Rio de Janeiro CF families were ascertained using a maximum-likelihood method. This same method was applied to study the distribution of the haplotypes defined by these markers. There was no significant association between the XV2c and KM19 loci on the parental and CF chromosomes. On the other hand, a strong association between GATT and TUB9 loci was observed on both CF and parental chromosomes, and striking linkage disequilibrium between the GATT-TUB9 pair and deltaF508 was observed (chi2 = 26.48, p < 0.0001). Remarkable linkage disequilibrium between the GATT-TUB9 marker pair and non-deltaF508 was also found (chi2 = 17.05, p < 0.0001). Our finding of a linkage disequilibrium between GATT-TUB9 and the CFTR locus could suggest that gene flow between different ethnic groups, mainly sub-Saharan and Mediterranean populations, with Brazilian populations could have resulted in some CF mutations originating on chromosomes that carried the GATT-TUB9 marker haplotype 7-2 (OR = 1.34 < 2.83 < 6.00; p = 0.0066).     

Selection of genetic markers for association analyses,using linkage disequilibrium and haplotypes

The genotyping of closely spaced single-nucleotide polymorphism (SNP) markers frequently yields highly correlated data, owing to extensive linkage disequilibrium (LD) between markers. The extent of LD varies widely across the genome and drives the number of frequent haplotypes observed in small regions. Several studies have illustrated the possibility that LD or haplotype data could be used to select a subset of SNPs that optimize the information retained in a genomic region while reducing the genotyping effort and simplifying the analysis. We propose a method based on the spectral decomposition of the matrices of pairwise LD between markers, and we select markers on the basis of their contributions to the total genetic variation. We also modify Clayton's "haplotype tagging SNP" selection method, which utilizes haplotype information. For both methods, we propose sliding window-based algorithms that allow the methods to be applied to large chromosomal regions. Our procedures require genotype information about a small number of individuals for an initial set of SNPs and selection of an optimum subset of SNPs that could be efficiently genotyped on larger numbers of samples while retaining most of the genetic variation in samples. We identify suitable parameter combinations for the procedures, and we show that a sample size of 50-100 individuals achieves consistent results in studies of simulated data sets in linkage equilibrium and LD. When applied to experimental data sets, both procedures were similarly effective at reducing the genotyping requirement while maintaining the genetic information content throughout the regions. We also show that haplotype-association results that Hosking et al. obtained near CYP2D6 were almost identical before and after marker selection.     

A global survey of haplotype frequencies and linkage disequilibrium at the DRD2 locus

A four-site haplotype system at the dopamine D2 receptor locus (DRD2) has been studied in a global sample of 28 distinct populations. The haplotype system spans about 25 kb, encompassing the coding region of the gene. The four individual markers include three TaqI restriction site polymorphisms (RSPs) – TaqI “A”, “B”, and “D” sites – and one dinucleotide short tandem repeat polymorphism (STRP). All four of the marker systems are polymorphic in all regions of the world and in most individual populations. The haplotype system shows the highest average heterozygosity in Africa, a slightly lower average heterozygosity in Europe, and the lowest average heterozygosities in East Asia and the Americas. Across all populations, 20 of the 48 possible haplotypes reached a frequency of at least 5% in at least one population sample. However, no single population had more than six haplotypes reaching that frequency. In general, African populations had more haplotypes present in each population and more haplotypes occurring at a frequency of at least 5% in that population. Permutation tests for significance of overall disequilibrium (all sites considered simultaneously) were highly significant (P<0.001) in all 28 populations. Except for three African samples, the pairwise disequilibrium between the outermost RSP markers, TaqI “B” and “A”, was highly significant with D’ values greater than 0.8; in two of those exceptions the RSP marker was not polymorphic. Except for those same two African populations, the 16-repeat allele at the STRP also showed highly significant disequilibrium with the TaqI “B” site in all populations, with D’ values usually greater than 0.7. Only four haplotypes account for more than 70% of all chromosomes in virtually all non-African populations, and two of those haplotypes account for more than 70% of all chromosomes in most East Asian and Amerindian populations. A new measure of the amount of overall disequilibrium shows least disequilibrium in African populations, somewhat more in European populations, and the greatest amount in East Asian and Amerindian populations. This pattern seems best explained by random genetic drift with low levels of recombination, a low mutation rate at the STRP, and essentially no recurrent mutation at the RSP sites, all in conjunction with an “Out of Africa” model for recent human evolution. Received: 14 January 1998 / Accepted 19 March 1998     

Detection of linkage disequilibrium between the myotonic dystrophy locus and a new polymorphic DNA marker.

We have examined the linkage of two new polymorphic DNA markers (D19S62 and D19S63) and a previously unreported polymorphism with an existing DNA marker (ERCC1) to the myotonic dystrophy (DM) locus. In addition, we have used pulsed-field gel electrophoresis to obtain a fine-structure map of this region. The detection of linkage disequilibrium between DM and one of these markers (D19S63) is the first demonstration of this phenomenon in a heterogeneous DM population. The results suggest that at least 58% of DM patients in the British population, as well as those in a French-Canadian subpopulation, are descended from the same ancestral DM mutation. We discuss the implications of this finding in terms of strategies for cloning the DM gene, for a possible role in modification of risk for prenatal and presymptomatic testing, and we speculate on the origin and number of existing mutations which may result in a DM phenotype.     

Extensive intragenic recombination and patterns of linkage disequilibrium at the CSN3 locus in European rabbit

Kappa-casein (CSN3) plays an important role in stabilising the Ca-sensitive caseins in the micelle. The European rabbit (Oryctolagus cuniculus) CSN3 has previously been shown to possess two alleles (A and B), which differ deeply in their intronic regions (indels of 100 and 1550 nucleotides in introns 1 and 4, respectively). Furthermore, a correlation between several reproductive performance traits and the different alleles was described. However, all these data were exclusively collected in rabbit domestic breeds, preventing a deeper understanding of the extensive polymorphism observed in the CSN3 gene. Additionally, the techniques available for the typing of both indel polymorphisms were until now not suitable for large-scale studies. In this report, we describe a simple, PCR-based typing method to distinguish rabbit CSN3 alleles. We analyse both ancient wild rabbit populations from the Iberian Peninsula and France, and the more recently derived English wild rabbits and domestic stocks. A new allele (C) showing another major indel (250 bp) in intron 1 was found, but exclusively detected in Iberian wild rabbits. In addition, our survey revealed the occurrence of new haplotypes in wild populations, suggesting that intragenic recombination is important in creating genetic diversity at this locus. This easy and low cost single-step PCR-based method results in an improvement over previous described techniques, can be easily set up in a routine molecular laboratory and would probably be a valuable tool in the management of rabbit domestic breeds.     

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.     

Haplotype diversity and the block structure of linkage disequilibrium

Several recent studies indicate that patterns of linkage disequilibrium in the human genome cannot be reconciled with a uniform distribution of recombination events, but crossovers appear to be localized in short hot-spots that separate longer stretches of DNA. Markers within these low-recombination blocks show increased levels of linkage disequilibrium and very low haplotype diversity. This could simplify study of the genetic basis of complex diseases if causal variants are common.     

The pH-dependent subunit dissociation and catalytic activity of bovine dopamine beta-hydroxylase

The soluble form of dopamine beta-hydroxylase from bovine adrenal medulla has previously been shown to exist as a tetrameric species of Mr = 290,000 composed of two disulfide-linked dimers. Here we report that this enzyme can also undergo a reversible tetramerdimer dissociation which is dependent on pH. Gel permeation chromatography of dopamine beta-hydroxylase at pH 5.0 demonstrates a Stokes radius of 5.8 nm. When the pH is shifted to 5.7, the Stokes radius changes to 6.9 nm. Sedimentation equilibrium analysis of the purified enzyme demonstrates that this change in molecular size is due to a change in molecular weight. At low protein concentration, the estimated Mr of the enzyme is 145,000 at pH 5.0 and at high protein concentration approaches 290,000 at pH 5.7. This change in Mr is consistent with the existence of a tetramer-dimer dissociation and a change in the equilibrium constant from 1.8 X 10(-6) M to 1.16 X 10(-9) M when the pH is increased from 5.0 to 5.7. This pH-dependent subunit dissociation is correlated with pH-dependent changes in enzyme activity. Purified bovine-soluble dopamine beta-hydroxylase activity is a hyperbolic function of tyramine concentration at pH 5.0. However, the hydroxylase activity displays non-hyperbolic kinetics at pH 6.0. The kinetic data obtained at pH 6.0 can be accounted for by fitting to a model containing two nonidentical catalytic forms of enzyme generated by the pH-dependent partial dissociation of tetrameric enzyme to dimeric subunits. The two catalytic forms have apparently identical maximal velocities; however, they differ in their Michaelis constants for the substrate; the dimeric form having a low Km and the tetrameric form having a high Km. Since the pH inside bovine adrenal medullary chromaffin granules is approximately 5.5, we conclude that the subunits of dopamine beta-hydroxylase are in dynamic dissociation in a physiologically important pH range.     

Linkage disequilibrium between an allele at the dopamine D4 receptor locus and Tourette syndrome, by the transmission-disequilibrium test.

Dopaminergic abnormalities are implicated in the pathogenesis of Tourette syndrome (TS) and chronic multiple tics. We used the transmission-disequilibrium test (TDT) method to test for linkage disequilibrium between a specific allele (the seven-repeat allele (DRD4*7R) of the exon 3 VNTR polymorphic site) at the D4 dopamine receptor locus (DRD4) and expression of chronic multiple tics and TS. This particular allele had been shown in functional studies to have different binding properties compared with the other common alleles in this DRD4 polymorphic system. We studied 64 family trios (consisting of an affected person and two parents, at least one heterozygous for DRD4*7R), including 12 nuclear family trios and 52 trios from four large TS kindreds. The DRD4*7R allele was transmitted significantly more frequently than expected (chi 2 TDT ranging from 8.47 [P < .004] to 10.80 [P = .001], depending on breadth of disease definition and inclusion or exclusion of inferred genotypes). Confirmation of this finding will depend on either replication in other samples or the identification of a transmitted functional mutation within this sample.     

Analysis of the MTHFR gene linkage disequilibrium structure and association of polymorphic gene variants with coronary atherosclerosis

Analysis of the genome-specific linkage disequilibrium patterns in certain populations is a highly promising approach to the identification of functional variants that underlie susceptibility to complex diseases. In the present study, the linkage disequilibrium patterns of the methylenetetrahydrofolate reductase gene (MTHFR) were examined in a group of patients with coronary atherosclerosis (coronary artery disease, CAD) and in a control sample from the Russian population. It was demonstrated that in the samples from one population, which were differentiated by the presence or absence of CAD, the MTHFR linkage disequilibrium patterns had similar features. Association of the MTHFR rs7533315 and rs2066462 polymorphisms with CAD was demonstrated. In addition, the evolution of the haplotypes and their role in the formation of CAD in the Russian population was reconstructed. The data on the association between genetic variability in the MTHFR locus and pathogenetically important indices of lipid metabolism were obtained. The high informativeness of the haplotype approach in case-control tests for associations with CAD was demonstrated.     

Extended haplotypes and linkage disequilibrium between 11 markers at the APOA1-C3-A4 gene cluster on chromosome 11

One hundred fifty-four unrelated French Caucasian subjects were typed for 11 RFLPs at or near the APOA1-C3-A4 gene cluster on the long arm of chromosome 11. All subjects belonged to families having lived in the Toulouse area (in the southwest of France) for over three generations. Allele frequencies for each RFLP were in agreement with previous studies in Caucasian populations for the APOA1/SstI marker. Pairwise linkage disequilibrium was determined. Among the 55 pairs studied, 30 are newly reported. Only three significant nonrandom associations were found: APOA1/MspI-3'APOC3/SstI, APOA1/MspI-3'APOA4/XbaI, and APOA4/DraI-APOA4/XbaI. Extended 11-marker haplotypes were constructed. Haplotype frequencies were estimated by the maximum-likelihood procedure and compared to expected frequencies calculated under the assumption of equilibrium. Among the 37 estimated haplotypes, seven containing at least four nonrandomly associated alleles showed markedly increased frequencies. These results, obtained in a geographically homogeneous population, confirm the existence of disequilibrium in the apolipoprotein cluster, but to a lower extent than previously reported in Caucasian populations, which were geographically more heterogeneous.     

Haplotypes and linkage disequilibrium at the phenylalanine hydroxylase locus, PAH, in a global representation of populations

Because defects in the phenylalanine hydroxylase gene (PAH) cause phenylketonuria (PKU), PAH was studied for normal polymorphisms and linkage disequilibrium soon after the gene was cloned. Studies in the 1980s concentrated on European populations in which PKU was common and showed that haplotype-frequency variation exists between some regions of the world. In European populations, linkage disequilibrium generally was found not to exist between RFLPs at opposite ends of the gene but was found to exist among the RFLPs clustered at each end. We have now undertaken the first global survey of normal variation and disequilibrium across the PAH gene. Four well-mapped single-nucleotide polymorphisms (SNPs) spanning approximately 75 kb, two near each end of the gene, were selected to allow linkage disequilibrium across most of the gene to be examined. These SNPs were studied as PCR-RFLP markers in samples of, on average, 50 individuals for each of 29 populations, including, for the first time, multiple populations from Africa and from the Americas. All four sites are polymorphic in all 29 populations. Although all but 5 of the 16 possible haplotypes reach frequencies >5% somewhere in the world, no haplotype was seen in all populations. Overall linkage disequilibrium is highly significant in all populations, but disequilibrium between the opposite ends is significant only in Native American populations and in one African population. This study demonstrates that the physical extent of linkage disequilibrium can differ substantially among populations from different regions of the world, because of both ancient genetic drift in the ancestor common to a large regional group of modern populations and recent genetic drift affecting individual populations.