A comparison between methods for linkage disequilibrium fine mapping of quantitative trait loci

We present a maximum likelihood method for mapping quantitative trait loci that uses linkage disequilibrium information from single and multiple markers. We made paired comparisons between analyses using a single marker, two markers and six markers. We also compared the method to single marker regression analysis under several scenarios using simulated data. In general, our method outperformed regression (smaller mean square error and confidence intervals of location estimate) for quantitative trait loci with dominance effects. In addition, the method provides estimates of the frequency and additive and dominance effects of the quantitative trait locus.     

Robust variance-components approach for assessing genetic linkage in pedigrees.

To assess evidence for genetic linkage from pedigrees, I developed a limited variance-components approach. In this method, variability among trait observations from individuals within pedigrees is expressed in terms of fixed effects from covariates and effects due to an unobservable trait-affecting major locus, random polygenic effects, and residual nongenetic variance. The effect attributable to a locus linked to a marker is a function of the additive and dominance components of variance of the locus, the recombination fraction, and the proportion of genes identical by descent at the marker locus for each pair of sibs. For unlinked loci, the polygenic variance component depends only on the relationship between the relative pair. Parameters can be estimated by either maximum-likelihood methods or quasi-likelihood methods. The forms of quasi-likelihood estimators are provided. Hypothesis tests derived from the maximum-likelihood approach are constructed by appeal to asymptotic theory. A simulation study showed that the size of likelihood-ratio tests was appropriate but that the monogenic component of variance was generally underestimated by the likelihood approach.     

DIFFERENTIAL DOMINANCE OF PLEIOTROPIC LOCI FOR MOUSE SKELETAL TRAITS

The term "differential dominance" describes the situation in which the dominance effects at a pleiotropic locus vary between traits. Directional selection on the phenotype can lead to balancing selection on differentially dominant pleiotropic loci. Even without any individual overdominant traits, some linear combination of traits will display overdominance at a locus displaying differential dominance. Multivariate overdominance may be responsible, in part, for high levels of heterozygosity found in natural populations. We examine differential dominance of 70 mouse skeletal traits at 92 quantitative trait loci (QTL). Our results indicate moderate to strong additive and dominance effects at pleiotropic loci, low levels of individual-trait overdominance, and universal multivariate overdominance. Multivariate overdominance affects a range of 6% to 81% of morphospace, with a mean of 32%. Multivariate overdominance tends to affect a larger percentage of morphospace at pleiotropic loci with antagonistic effects on multiple traits (42%). We conclude that multivariate overdominance is common and should be considered in models and in empirical studies of the role of genetic variation in evolvability.     

Use of Multiple Genetic Markers in Prediction of Breeding Values

Genotypes at a marker locus give information on transmission of genes from parents to offspring and that information can be used in predicting the individuals' additive genetic value at a linked quantitative trait locus (MQTL). In this paper a recursive method is presented to build the gametic relationship matrix for an autosomal MQTL which requires knowledge on recombination rate between the marker locus and the MQTL linked to it. A method is also presented to obtain the inverse of the gametic relationship matrix. This information can be used in a mixed linear model for simultaneous evaluation of fixed effects, gametic effects at the MQTL and additive genetic effects due to quantitative trait loci unlinked to the marker locus (polygenes). An equivalent model can be written at the animal level using the numerator relationship matrix for the MQTL and a method for obtaining the inverse of this matrix is presented. Information on several unlinked marker loci, each of them linked to a different locus affecting the trait of interest, can be used by including an effect for each MQTL. The number of equations per animal in this case is 2m + 1 where m is the number of MQTL. A method is presented to reduce the number of equations per animal to one by combining information on all MQTL and polygenes into one numerator relationship matrix. It is illustrated how the method can accommodate individuals with partial or no marker information. Numerical examples are given to illustrate the methods presented. Opportunities to use the presented model in constructing genetic maps are discussed.     

Derivation of single-locus relationship coefficients conditional on marker information

The coefficient of relationship is defined as the correlation between the additive genetic values of two individuals. This coefficient can be defined specifically for a single quantitative trait locus (QTL) and may deviate considerably from the overall expectation if it is taken conditional on information from linked marker loci. Conditional halfsib correlations are derived under a simple genetic model with a biallelic QTL linked to a biallelic marker locus. The conditional relationship coefficients are shown to depend on the recombination rate between the marker and the QTL and the population frequency of the marker alleles, but not on parameters of the QTL, i.e. number and frequency of QTL alleles, degree of dominance etc., nor on the (usually unknown) QTL genotype of the sire. Extensions to less simplified cases (multiple alleles at the marker locus and the QTL, two marker loci flanking the QTL) are given. For arbitrary pedigrees, conditional relationship coefficients can also be derived from the conditional gametic covariance matrix suggested by Fernando and Grossman (1989). The connection of these two approaches is discussed. The conditional relationship coefficient can be used for marker-assisted genetic evaluation as well as for the detection of QTL and the estimation of their effects.     

Bayesian analysis of linkage between genetic markers and quantitative trait loci. I. Prior knowledge

Summary Prior information on gene effects at individual quantitative trait loci (QTL) and on recombination rates between marker loci and QTL is derived. The prior distribution of QTL gene effects is assumed to be exponential with major effects less likely than minor ones. The prior probability of linkage between a marker and another single locus is a function of the number and length of chromosomes, and of the map function relating recombination rate to genetic distance among loci. The prior probability of linkage between a marker locus and a quantitative trait depends additionally on the number of detectable QTL, which may be determined from total additive genetic variance and minimum detectable QTL effect. The use of this prior information should improve linkage tests and estimates of QTL effects.     

Genetic evaluation methods for populations with dominance and inbreeding

Summary The effect of inbreeding on mean and genetic covariance matrix for a quantitative trait in a population with additive and dominance effects is shown. This genetic covariance matrix is a function of five relationship matrices and five genetic parameters describing the population. Elements of the relationship matrices are functions of Gillois (1964) identity coefficients for the four genes at a locus in two individuals. The equivalence of the path coefficient method (Jacquard 1966) and the tabular method (Smith and Mäki-Tanila 1990) to compute the covariance matrix of additive and dominance effects in a population with inbreeding is shown. The tabular method is modified to compute relationship matrices rather than the covariance matrix, which is trait dependent. Finally, approximate and exact Best Linear Unbiased Predictions (BLUP) of additive and dominance effects are compared using simulated data with inbreeding but no directional selection. The trait simulated was affected by 64 unlinked biallelic loci with equal effect and complete dominance. Simulated average inbreeding levels ranged from zero in generation one to 0.35 in generation five. The approximate method only accounted for the effect of inbreeding on mean and additive genetic covariance matrix, whereas the exact accounted for all of the changes in mean and genetic covariance matrix due to inbreeding. Approximate BLUP, which is computable for large populations where exact BLUP is not feasible, yielded unbiased predictions of additive and dominance effects in each generation with only slightly reduced accuracies relative to exact BLUP.     

A Bayesian Approach to Detect Quantitative Trait Loci Using Markov Chain Monte Carlo

Markov chain Monte Carlo (MCMC) techniques are applied to simultaneously identify multiple quantitative trait loci (QTL) and the magnitude of their effects. Using a Bayesian approach a multi-locus model is fit to quantitative trait and molecular marker data, instead of fitting one locus at a time. The phenotypic trait is modeled as a linear function of the additive and dominance effects of the unknown QTL genotypes. Inference summaries for the locations of the QTL and their effects are derived from the corresponding marginal posterior densities obtained by integrating the likelihood, rather than by optimizing the joint likelihood surface. This is done using MCMC by treating the unknown QTL genotypes, and any missing marker genotypes, as augmented data and then by including these unknowns in the Markov chain cycle along with the unknown parameters. Parameter estimates are obtained as means of the corresponding marginal posterior densities. High posterior density regions of the marginal densities are obtained as confidence regions. We examine flowering time data from double haploid progeny of Brassica napus to illustrate the proposed method.     

Effects of genetic drift on variance components under a general model of epistasis

We analyze the changes in the mean and variance components of a quantitative trait caused by changes in allele frequencies, concentrating on the effects of genetic drift. We use a general representation of epistasis and dominance that allows an arbitrary relation between genotype and phenotype for any number of diallelic loci. We assume initial and final Hardy-Weinberg and linkage equilibrium in our analyses of drift-induced changes. Random drift generates transient linkage disequilibria that cause correlations between allele frequency fluctuations at different loci. However, we show that these have negligible effects, at least for interactions among small numbers of loci. Our analyses are based on diffusion approximations that summarize the effects of drift in terms of F, the inbreeding coefficient, interpreted as the expected proportional decrease in heterozygosity at each locus. For haploids, the variance of the trait mean after a population bottleneck is var(delta(z)) = sigma(n)k=1 FkV(A(k)), where n is the number of loci contributing to the trait variance, V(A(1)) = V(A) is the additive genetic variance, and V(A(k)) is the kth-order additive epistatic variance. The expected additive genetic variance after the bottleneck, denoted (V*(A)), is closely related to var(delta(z)); (V*(A)) = (1 - F) sigma(n)k=1 kFk-1V(A(k)). Thus, epistasis inflates the expected additive variance above V(A)(1 - F), the expectation under additivity. For haploids (and diploids without dominance), the expected value of every variance component is inflated by the existence of higher order interactions (e.g., third-order epistasis inflates (V*(AA. This is not true in general with diploidy, because dominance alone can reduce (V*(A)) below V(A)(1 - F) (e.g., when dominant alleles are rare). Without dominance, diploidy produces simple expressions: var(delta(z)) = sigma(n)k=1 (2F)kV(A(k)) and (V(A)) = (1 - F) sigma(n)k=1 k(2F)k-1V(A(k)). With dominance (and even without epistasis), var(delta(z)) and (V*(A)) no longer depend solely on the variance components in the base population. For small F, the expected additive variance simplifies to (V*(A)) approximately equal to (1 - F)V(A) + 4FV(AA) + 2FV(D) + 2FC(AD), where C(AD) is a sum of two terms describing covariances between additive effects and dominance and additive X dominance interactions. Whether population bottlenecks lead to expected increases in additive variance depends primarily on the ratio of nonadditive to additive genetic variance in the base population, but dominance precludes simple predictions based solely on variance components. We illustrate these results using a model in which genotypic values are drawn at random, allowing extreme and erratic epistatic interactions. Although our analyses clarify the conditions under which drift is expected to increase V(A), we question the evolutionary importance of such increases.     

Interactions between markers can be caused by the dominance effect of quantitative trait loci

F₂ populations are commonly used in genetic studies of animals and plants. For simplicity, most quantitative trait locus or loci (QTL) mapping methods have been developed on the basis of populations having two distinct genotypes at each polymorphic marker or gene locus. In this study, we demonstrate that dominance can cause the interactions between markers and propose an inclusive linear model that includes marker variables and marker interactions so as to completely control both additive and dominance effects of QTL. The proposed linear model is the theoretical basis for inclusive composite-interval QTL mapping (ICIM) for F₂ populations, which consists of two steps: first, the best regression model is selected by stepwise regression, which approximately identifies markers and marker interactions explaining both additive and dominance variations; second, the interval mapping approach is applied to the phenotypic values adjusted by the regression model selected in the first step. Due to the limited mapping population size, the large number of variables, and multicollinearity between variables, coefficients in the inclusive linear model cannot be accurately determined in the first step. Interval mapping is necessary in the second step to fine tune the QTL to their true positions. The efficiency of including marker interactions in mapping additive and dominance QTL was demonstrated by extensive simulations using three QTL distribution models with two population sizes and an actual rice F₂ population.     

QTL analysis: a simple ‘marker-regression’ approach

A method to locate quantitative trait loci (QTL) on a chromosome and to estimate their additive and dominance effects is described. It applies to generations derived from an F₁ by selfing or backcrossing and to doubled haploid lines, given that marker genotype information (RFLP, RAPD, etc.) and quantitative trait data are available. The method involves regressing the additive difference between marker genotype means at a locus against a function of the recombination frequency between that locus and a putative QTL. A QTL is located, as by other regression methods, at that point where the residual mean square is minimised. The estimates of location and gene effects are consistent and as reliable as conventional flanking-marker methods. Further applications include the ability to test for the presence of two, or more, linked QTL and to compare different crosses for the presence of common QTL. Furthermore, the technique is straightforward and may be programmed using standard pc-based statistical software.     

Heterosis for biomass-related traits in Arabidopsis investigated by quantitative trait loci analysis of the triple testcross design with recombinant inbred lines

Arabidopsis thaliana has emerged as a leading model species in plant genetics and functional genomics including research on the genetic causes of heterosis. We applied a triple testcross (TTC) design and a novel biometrical approach to identify and characterize quantitative trait loci (QTL) for heterosis of five biomass-related traits by (i) estimating the number, genomic positions, and genetic effects of heterotic QTL, (ii) characterizing their mode of gene action, and (iii) testing for presence of epistatic effects by a genomewide scan and marker x marker interactions. In total, 234 recombinant inbred lines (RILs) of Arabidopsis hybrid C24 x Col-0 were crossed to both parental lines and their F1 and analyzed with 110 single-nucleotide polymorphism (SNP) markers. QTL analyses were conducted using linear transformations Z1, Z2, and Z3 calculated from the adjusted entry means of TTC progenies. With Z1, we detected 12 QTL displaying augmented additive effects. With Z2, we mapped six QTL for augmented dominance effects. A one-dimensional genome scan with Z3 revealed two genomic regions with significantly negative dominance x additive epistatic effects. Two-way analyses of variance between marker pairs revealed nine digenic epistatic interactions: six reflecting dominance x dominance effects with variable sign and three reflecting additive x additive effects with positive sign. We conclude that heterosis for biomass-related traits in Arabidopsis has a polygenic basis with overdominance and/or epistasis being presumably the main types of gene action.     

Quantitative trait locus mapping with background control in genetic populations of clonal F1 and double cross

In this study, we considered five categories of molecular markers in clonal F₁ and double cross populations, based on the number of distinguishable alleles and the number of distinguishable genotypes at the marker locus. Using the completed linkage maps, incomplete and missing markers were imputed as fully informative markers in order to simplify the linkage mapping approaches of quantitative trait genes. Under the condition of fully informative markers, we demonstrated that dominance effect between the female and male parents in clonal F₁ and double cross populations can cause the interactions between markers. We then developed an inclusive linear model that includes marker variables and marker interactions so as to completely control additive effects of the female and male parents, as well as the dominance effect between the female and male parents. The linear model was finally used for background control in inclusive composite interval mapping (ICIM) of quantitative trait locus (QTL). The efficiency of ICIM was demonstrated by extensive simulations and by comparisons with simple interval mapping, multiple‐QTL models and composite interval mapping. Finally, ICIM was applied in one actual double cross population to identify QTL on days to silking in maize.     

Pleiotropic effects on mandibular morphology I. Developmental morphological integration and differential dominance

Pleiotropy refers to a single genetic locus that affects more than one phenotypic trait. Pleiotropic effects of genetic loci are thought to play an important role in evolution, reflecting functional and developmental relationships among phenotypes. In a previous study, we examined pleiotropic effects displayed by quantitative trait loci (QTLs) on murine mandibular morphology in relation to mandibular structure and function. In replicating most of our previous QTLs and increasing our sample size, this study strengthens and extends our earlier results. As in our previous study, we find that QTL effects tend to be restricted to developmentally or functionally related traits. In addition, we examine patterns of differential dominance for pleiotropic QTL effects. Differential dominance occurs when dominance patterns for a single locus vary among traits. We find that multivariate overdominance is a common and substantial phenomenon, and may potentially provide an explanation for the persistence of heterozygosity in natural populations.     

Genetic effects at pleiotropic loci are context-dependent with consequences for the maintenance of genetic variation in populations

Context-dependent genetic effects, including genotype-by-environment and genotype-by-sex interactions, are a potential mechanism by which genetic variation of complex traits is maintained in populations. Pleiotropic genetic effects are also thought to play an important role in evolution, reflecting functional and developmental relationships among traits. We examine context-dependent genetic effects at pleiotropic loci associated with normal variation in multiple metabolic syndrome (MetS) components (obesity, dyslipidemia, and diabetes-related traits). MetS prevalence is increasing in Western societies and, while environmental in origin, presents substantial variation in individual response. We identify 23 pleiotropic MetS quantitative trait loci (QTL) in an F₁₆ advanced intercross between the LG/J and SM/J inbred mouse strains (Wustl:LG,SM-G16; n = 1002). Half of each family was fed a high-fat diet and half fed a low-fat diet; and additive, dominance, and parent-of-origin imprinting genotypic effects were examined in animals partitioned into sex, diet, and sex-by-diet cohorts. We examine the context-dependency of the underlying additive, dominance, and imprinting genetic effects of the traits associated with these pleiotropic QTL. Further, we examine sequence polymorphisms (SNPs) between LG/J and SM/J as well as differential expression of positional candidate genes in these regions. We show that genetic associations are different in different sex, diet, and sex-by-diet settings. We also show that over- or underdominance and ecological cross-over interactions for single phenotypes may not be common, however multidimensional synthetic phenotypes at loci with pleiotropic effects can produce situations that favor the maintenance of genetic variation in populations. Our findings have important implications for evolution and the notion of personalized medicine.     

Molecular-Marker-Facilitated Investigations of Quantitative-Trait Loci in Maize. I. Numbers, Genomic Distribution and Types of Gene Action

Individual genetic factors which underlie variation in quantitative traits of maize were investigated in each of two F2 populations by examining the mean trait expressions of genotypic classes at each of 17-20 segregating marker loci. It was demonstrated that the trait expression of marker locus classes could be interpreted in terms of genetic behavior at linked quantitative trait loci (QTLs). For each of 82 traits evaluated, QTLs were detected and located to genomic sites. The numbers of detected factors varied according to trait, with the average trait significantly influenced by almost two-thirds of the marked genomic sites. Most of the detected associations between marker loci and quantitative traits were highly significant, and could have been detected with fewer than the 1800-1900 plants evaluated in each population. The cumulative, simple effects of marker-linked regions of the genome explained between 8 and 40% of the phenotypic variation for a subset of 25 traits evaluated. Single marker loci accounted for between 0.3% and 16% of the phenotypic variation of traits. Individual plant heterozygosity, as measured by marker loci, was significantly associated with variation in many traits. The apparent types of gene action at the QTLs varied both among traits and between loci for given traits, although overdominance appeared frequently, especially for yield-related traits. The prevalence of apparent overdominance may reflect the effects of multiple QTLs within individual marker-linked regions, a situation which would tend to result in overestimation of dominance. Digenic epistasis did not appear to be important in determining the expression of the quantitative traits evaluated. Examination of the effects of marked regions on the expression of pairs of traits suggests that genomic regions vary in the direction and magnitudes of their effects on trait correlations, perhaps providing a means of selecting to dissociate some correlated traits. Marker-facilitated investigations appear to provide a powerful means of examining aspects of the genetic control of quantitative traits. Modifications of the methods employed herein will allow examination of the stability of individual gene effects in varying genetic backgrounds and environments.     

The Contribution of Quantitative Trait Loci and Neutral Marker Loci to the Genetic Variances and Covariances among Quantitative Traits in Random Mating Populations

Using Cockerham's approach of orthogonal scales, we develop genetic models for the effect of an arbitrary number of multiallelic quantitative trait loci (QTLs) or neutral marker loci (NMLs) upon any number of quantitative traits. These models allow the unbiased estimation of the contributions of a set of marker loci to the additive and dominance variances and covariances among traits in a random mating population. The method has been applied to an analysis of allozyme and quantitative data from the European oyster. The contribution of a set of marker loci may either be real, when the markers are actually QTLs, or apparent, when they are NMLs that are in linkage disequilibrium with hidden QTLs. Our results show that the additive and dominance variances contributed by a set of NMLs are always minimum estimates of the corresponding variances contributed by the associated QTLs. In contrast, the apparent contribution of the NMLs to the additive and dominance covariances between two traits may be larger than, equal to or lower than the actual contributions of the QTLs. We also derive an expression for the expected variance explained by the correlation between a quantitative trait and multilocus heterozygosity. This correlation explains only a part of the genetic variance contributed by the markers, i.e., in general, a combination of additive and dominance variances and, thus, provides only very limited information relative to the method supplied here.     

Dissecting total genetic variance into additive and dominance components of purebred and crossbred pig traits

The partition of the total genetic variance into its additive and non-additive components can differ from trait to trait, and between purebred and crossbred populations. A quantification of these genetic variance components will determine the extent to which it would be of interest to account for dominance in genomic evaluations or to establish mate allocation strategies along different populations and traits. This study aims at assessing the contribution of the additive and dominance genomic variances to the phenotype expression of several purebred Piétrain and crossbred (Piétrain × Large White) pig performances. A total of 636 purebred and 720 crossbred male piglets were phenotyped for 22 traits that can be classified into six groups of traits: growth rate and feed efficiency, carcass composition, meat quality, behaviour, boar taint and puberty. Additive and dominance variances estimated in univariate genotypic models, including additive and dominance genotypic effects, and a genomic inbreeding covariate allowed to retrieve the additive and dominance single nucleotide polymorphism variances for purebred and crossbred performances. These estimated variances were used, together with the allelic frequencies of the parental populations, to obtain additive and dominance variances in terms of genetic breeding values and dominance deviations. Estimates of the Piétrain and Large White allelic contributions to the crossbred variance were of about the same magnitude in all the traits. Estimates of additive genetic variances were similar regardless of the inclusion of dominance. Some traits showed relevant amount of dominance genetic variance with respect to phenotypic variance in both populations (i.e. growth rate 8%, feed conversion ratio 9% to 12%, backfat thickness 14% to 12%, purebreds-crossbreds). Other traits showed higher amount in crossbreds (i.e. ham cut 8% to 13%, loin 7% to 16%, pH semimembranosus 13% to 18%, pH longissimus dorsi 9% to 14%, androstenone 5% to 13% and estradiol 6% to 11%, purebreds-crossbreds). It was not encountered a clear common pattern of dominance expression between groups of analysed traits and between populations. These estimates give initial hints regarding which traits could benefit from accounting for dominance for example to improve genomic estimated breeding value accuracy in genetic evaluations or to boost the total genetic value of progeny by means of assortative mating.     

The quantitative genetics of fluctuating asymmetry

Fluctuating asymmetry (subtle departures from identical expression of a trait across an axis of symmetry) in many taxa is under stabilizing selection for reduced asymmetry. However, lack of reliable estimates of genetic parameters for asymmetry variation hampers our ability to predict the evolutionary outcome of this selection. Here we report on a study, based on analysis of variation within and between isofemale lines and of generation means (line-cross analysis), designed to dissect in detail the quantitative genetics of positional fluctuating asymmetry (PFA) in bristle number in natural populations of Drosophila falleni. PFA is defined as the difference between the two sides of the body in the placement or position of components of a meristic trait. Heritability (measured at 25 degrees C) of two related measures of PFA were 13% and 21%, both of which differed significantly from zero. In contrast, heritability estimates for fluctuating asymmetry in the total number of anterior (0.7%) and transverse (2.4%) sternopleural bristles were smaller, not significant, and in quantitative agreement with previously published estimates. Heritabilities for bristle number (trait size) were considerably greater than that for any asymmetry measure. The experimental design controlled for the potentially confounding effects of common familial environment, and repeated testing revealed that PFA differences between lines were genetically stable for up to 16 generations in the laboratory at 25 degrees C. We performed line cross analysis between strains at the extremes of the PFA distribution (highest and lowest values); parental strains, F1, F1r (reciprocal), F2, backcross, and backcross reciprocal generations were represented. The inheritance of PFA was described best by additive and dominance effects localized to the X-chromosomes, whereas autosomal dominance effects were also detected. Epistatic, maternal, and cytoplasmic effects were not detected. The inheritance of trait size was notably more complex and involved significant autosomal additive, dominance, and epistatic effects; maternal dominance effects; and additive and dominance effects localized to the X-chromosomes. The additive genetic correlation between PFA and its associated measure of trait size was negative (-0.049), but not statistically significant, indicating that the loci contributing additive genetic effects to these traits are probably different. It is suggested that PFA may be a sensitive measure of developmental instability because PFA taps the ability of an organism to integrate interconnected developmental pathways.     

Association studies under general disease models

There is great expectation that the levels of association found between genetic markers and disease status will play a role in the location of disease genes. This expectation follows from regarding association as being proportional to linkage disequilibrium and therefore inversely related to recombination value. For disease genes with more than two alleles, the association measure is instead a weighted average of linkage disequilibria, with the weights depending on allele frequencies and genotype susceptibilities at the disease loci. There is no longer a simple relationship, even in expectation, with recombination. We adopt a general framework to examine association mapping methods which helps to clarify the nature of case-control and transmission/disequilibrium-type tests and reveals the relationship between measures of association and coefficients of linkage disequilibrium. In particular, we can show the consequences of additive and nonadditive effects at the trait locus on the behavior of these tests. These concepts have a natural extension to marker haplotypes. The association of two-locus marker haplotypes with disease phenotype depends on a weighted average of three-locus disequilibria (two markers with each disease locus). It is likely that these two-marker analyses will provide additional information in association mapping studies.