Haplotyping a Quantitative Trait with a High-Density Map in Experimental Crosses

Background The ultimate goal of genetic mapping of quantitative trait loci (QTL) is the positional cloning of genes involved in any agriculturally or medically important phenotype. However, only a small portion (≤ 1%) of the QTL detected have been characterized at the molecular level, despite the report of hundreds of thousands of QTL for different traits and populations. Methods/Results We develop a statistical model for detecting and characterizing the nucleotide structure and organization of haplotypes that underlie QTL responsible for a quantitative trait in an F2 pedigree. The discovery of such haplotypes by the new model will facilitate the molecular cloning of a QTL. Our model is founded on population genetic properties of genes that are segregating in a pedigree, constructed with the mixture-based maximum likelihood context and implemented with the EM algorithm. The closed forms have been derived to estimate the linkage and linkage disequilibria among different molecular markers, such as single nucleotide polymorphisms, and quantitative genetic effects of haplotypes constructed by non-alleles of these markers. Results from the analysis of a real example in mouse have validated the usefulness and utilization of the model proposed. Conclusion The model is flexible to be extended to model a complex network of genetic regulation that includes the interactions between different haplotypes and between haplotypes and environments.


INTRODUCTION
The basic principle for quantitative trait locus (QTL) mapping is the cosegregation of the alleles at a QTL with those at one or a set of known polymorphic markers genotyped on a genome in an experimental cross [1,2]. If a QTL is cosegregating with molecular markers, the genetic effects of the QTL on a quantitative trait and its genomic position can be estimated from the marker genotypes and phenotypic values of the trait. This estimation process particularly assumes the QTL to be located within an interval constructed by a pair of flanking markers in which a test statistics calculated under the reduced (there is no QTL) and full model (there is a QTL) is used to test the existence of the QTL and estimate its position. This so-called interval mapping approach and its extensions [3][4][5] is robust and powerful for the detection of major QTL and presents the most efficient way to utilize marker information when marker maps are sparse [6]. However, interval mapping is limited by its incapacity to infer any information about the sequence structure and organization of the QTL. Partly for this reason, only a few QTL mapped from markers have been successfully cloned [7][8][9], despite a considerable number of QTL reported in the literature.
Interval QTL mapping also has an unsolved statistical difficulty when it is used with a high-density linkage map. With more markers genotyped, a genetic map for QTL identification has tended to be infinitely dense. For such an infinitely dense map in which markers are located everywhere over the genome, test statistics at nearby intervals are not independent any more. Thus, the critical threshold used to acclaim the existence of a QTL by interval mapping will be difficult to analytically determine. Although an empirical alternative based on permutation tests has been proposed for threshold determination [10], extensive computing may affect the use efficiency of interval mapping.
Despite its unsuitability for interval mapping of QTL, an infinitely dense map provides an important tool for characterizing genetic variants that contribute to quantitative variation via the analysis of haplotypes composed of non-alleles at a set of highly linked markers. Recent genetic studies suggest that a gene may determine a complex trait, such as body weight or drug response, through its haplotype rather than genotype [11,12]. The completion of the genome projects for several important organisms, Arabdopsis, chicken, human, mouse and poplar, has made massive amounts of DNA sequence data available. In particular, single nucleotide polymorphisms (SNPs), being the most common type of variant in the DNA sequence, provide a powerful means for genotyping the whole genome or any part of it. This facilitates the identification of specific SNP-constructed haplotypes which are responsible for quantitative traits. A set of SNPs that cause quantitative differences among individuals are called quantitative trait nucleotides (QTNs). Liu et al. [13] proposed a statistical model for estimating and testing haplotype effects at a QTN in a random sample drawn from a natural population. This model is based on the population genetic properties of gene segregation. Through the implementation of the EM algorithm, population genetic parameters of SNPs, such as haplotype frequencies, allele frequencies and linkage disequilibria, and quantitative genetic parameters, such as haplotype effects of a QTN, are estimated with closed forms.
The motivation of this work is to derive a statistical model for haplotype discovery responsible for quantitative variation in a mapping population derived an experimental cross. Unlike a natural population in which gene co-segregation analysis is based on linkage disequilibria [14], experimental crosses, such as the backcross or F 2 , have usually been analyzed in terms of the linkage between different markers and QTL. In this article, we will frame a general statistical model for estimating the linkage between different SNPs and testing haplotype effects within the context of linkage disequilibrium analysis in an F 2 pedigree. We show that the new model can test for the dependence of SNPs when a multipoint analysis is performed. We have derived closed forms for the EM algorithm to estimate a variety of genetic parameters. A worked example is used to validate the usefulness and utilization of the model.

Haplotype and diplotype
A haplotype represents a linear arrangement of nucleotides (alleles) at different SNPs on a single chromosome, or part of a chromosome. The pair of haplotypes is called a diplotype. The observed phenotype of a diplotype is called a genotype. A diplotype is always constructed by two haplotypes, one from the maternal parent and the other from the paternal parent. Suppose there are two different SNPs on the same genomic region, one with two alleles A and a and the other with two alleles B and b, respectively. Allele A from SNP 1 and allele B from SNP 2 are located on the first homologous chromosome, whereas allele a from SNP 1 and allele b from SNP 2 located on the second homologous chromosome. , its growth will be poorer. The statistical model being developed will be used to determine which diplotype is associated with better growth in experimental crosses.

Linkage disequilibrium in the F 2 intercross
A general model: Haplotype analysis in the backcross is straightforward because the diplotype is determined for all the backcross genotype. Simple analysis of variance can be used to detect haplotype effects on a quantitative trait. In the F 2 , this is not a case in which the double heterozygote is a mixture of two possible diplotypes.
Suppose many SNPs are genotyped each of which is segregating in a 1:2:1 Mendelian ratio in the F 2 population. As seen in the human genome [15], these SNPs are divided into different haplotype blocks. For a given block, there are a particular number of representative SNPs or htSNPs that uniquely identify the common haplotypes in this block or QTN. Several algorithms have been developed to identify a minimal subset of htSNPs that can characterize the most common haplotypes [16][17][18]. Consider a QTN that contains L htSNPs among which there exist linkage disequilibria of different orders. The two alleles, 1 and 0, at each of these SNPs are symbolized by r 1 ,…,r L , respectively. For a cross initiated with two inbred parents, the allele frequencies for each of these htSNPs should be 1/2. A haplotype frequency, denoted as p r1r2ÁÁÁrL , is decomposed into the following components: where D's are the linkage disequilibria of different orders among particular SNPs.  Totally, L SNPs form 2 L haplotypes expressed as [r 1 …r L ], 2 L21 (2 L +1) diplotypes, i.e., a pair of maternally-(m) and paternally-derived haplotypes (p), expressed as [r 1 m …r L m ] [r 1 p …r L p ] (r 1 m , r 1 p ,…;r L m ,r L p = 1,0) and 3 L genotypes expressed as r 1 r9 1 /…/r L r9 L (r 1 $r9 1, …,r L $r9 L = 1,0). Only genotypes can be observed. The number of diplotypes is smaller than the number of genotypes because the genotypes that are heterozygous at two or more SNPs contain multiple different diplotypes. Diplotype (and therefore genotype) frequencies can be expressed in terms of haplotype frequencies. We use P ½r m 1 :::r m L ½r p 1 :::r p L and P r1 r 0 1 =...=rL r 0 L to denote the diplotype and genotype frequencies, respectively, and n r1 r 0 1 =...=rL r 0 L to denote genotype observation.
A special case: Two-point linkage disequilibrium: For two given SNPs (S 1 and S 2 ), there are four different haplotypes in a cross population. According to the definition given above, these four haplotypes are denoted as [11], [10], [01] and [00], whose frequencies in a cross population are, respectively, expressed as Assume that the two SNPs are linked with a recombination fraction r. The haplotype frequencies can be expressed in terms of r, i.e., p 11~1 2 (1{r), p 10~1 2 r, p 01~1 2 r and p 00~1 2 (1{r). Combining equation (2), this establishes the relation between the linkage disequilibrium and recombination fraction as The first three first-order linkage disequilibria can be used to describe the linkage between different SNPs and crossover interference, whereas the last second-order linkage disequilibrium is thought to be associated with chromatid interference.

Haplotyping a trait with two SNPs
Our interest is to search for the haplotype diversity that can explain phenotypic variation in a complex trait. The association between haplotype diversity and phenotypic variation has been detected in several studies of drug responses [11,12]. This allows us to assume that a particular haplotype is different from other haplotypes for a given trait. Here, our focus will be on modelling haplotype effects in experimental crosses. Although haplotypes (comprising diplotypes) can be directly observed in the backcross, this is not possible for the F 2 because their heterozygous genotypes are not concordant with diplotypes or haplotypes. For the F 2 population, the effects of different haplotypes on the phenotype need be postulated from observed zygotic genotypes. The inference of diplotypes for a particular genotype is statistically a missing data problem that can be formulated by a finite mixture model.
Mixture model: The statistical method for the genomewide scan of QTN is formulated on the basis of a finite mixture model. The mixture model assumes that each observation comes from one of an assumed set of distributions. The mixture model derived to detect haplotype effects on a quantitative trait based on SNP genotype data contains three major parts: (1) the mixture proportions of each distribution, denoted as the relative frequencies of different diplotypes for the same SNP genotype, (2) the mean for each diplotype in the density function, and (3) the residual variance common to all diplotypes.
For simplicity, we consider a QTN that is composed of only two SNPs each with two alleles designated as 1 and 0. These two SNPs segregating in the F 2 population form four haplotypes whose frequencies are arrayed in vector Hp = (p 11 , p 10 , p 01 , p 00 ). All the genotypes are consistent with diplotypes, except for the double heterozygote, 10/10, that contains two different diplotypes [11][00] with a frequency of 2 p 11 p 00 and [10][01] with a frequency of 2 p 10 p 01 ( Table 1). The relative frequencies of different diplotypes for the double heterozygote are a function of haplotype frequencies.
A total of n individuals in the F 2 are classified into 9 genotypes for the two SNPs, each genotype with observation generally expressed as n r1 r 0 1 =r2 r 0 2 (r 1 $r9 1 ,r 2 $r9 2 ,r 3 $r9 3 = 1,0). The frequency of each genotype can be expressed in terms of haplotype frequencies (Table 1). Considering a quantitative trait controlled by diplotype (rather than genotype) diversity, the phenotypic value of the trait (y i ) for individual i is expressed by a linear model, i.e., where j i is the indicator variable defined as 1 if a diplotype considered is compatible with subject i and as 0 otherwise, u ½r m 1 r m 2 ½r p 1 r p 2 ~u½r p 1 r p 2 ½r m 1 r m 2 is the genotypic value for diplotype ½r m 1 r m 2 ½r p 1 r p 2 , and e i is the residual error distributed as N(0,s 2 ). Assume that this QTN triggers an effect on the trait because at least one haplotype is different from the remaining seven. Without loss of generality, let [11] be such a distinct haplotype, called risk haplotype, designated as A. All the other non-risk haplotypes, [10], [01] and [00], are collectively expressed as Ā . The risk and nonrisk haplotypes form three composite diplotypes AA (2), AĀ (1) and Ā Ā (0). Let m 2 , m 1 and m 0 be the genotypic value of the three composite diplotypes, respectively ( Table 1). The means for different composite diplotypes and residual variance are arrayed by a quantitative genetic parameter vector H q = (m 2, m 1, m 0, s 2 ).
Likelihoods: With the above notation, we construct two likelihoods, one for haplotype frequencies (H p ) based on SNP data (S) and the other for quantitative genetic parameters (H q ) based on haplotype frequencies (H p ), phenotypic (y) and SNP data (S). They are, respectively, expressed as where f j (y i ) is a normal distribution density function of composite diplotype j (j = 2,1,0), i.e., It can be seen from the above likelihood functions that, although most zygote genotypes contain a single component (diplotype), the double heterozygote is the mixture of two possible diplotypes weighted by w and 1-w, expressed as w~p 11 p 00 p 11 p 00 zp 10 p 01 , ð9Þ which represents the relative frequency of diplotype [11][00] for the double heterozygote.
It should be noted that L(H p , H q | y, S) relies on the haplotype frequencies defined in L(H p |S) and, thus, the latter is thought to be nested within the former. The estimates of parameters that maximize L(H p |S) can also maximize the L(H p , H q | y, S).
The EM algorithm: A closed-form solution for the EM algorithm has been derived to estimate the unknown parameters that maximize the two likelihoods of (26) [13]. The estimates of haplotype frequencies are based on the log-likelihood function L(H p |M), whereas the estimates of diplotype genotypic means and residual variance are based on the log-likelihood function L(H p , H q | y, M). These two different types of parameters can be estimated using a two-stage hierarchical EM algorithm.
At a higher hierarchy of the EM algorithm, the E step is aimed to calculate the relative frequency (w) of diplotype [11][00] in the double heterozygote is calculated by equation (9). The M step is aimed to estimate the haplotype frequencies based on the probabilities calculated in the previous iteration usinĝ  At a lower hierarchy of the EM algorithm, the E step is derived to calculate the posterior probability (V [ 11 ] [00]i ) of individual i with the double heterozygous genotype to be diplotype [11][00] by Note that for all the other genotypes, such posterior probabilities do not exist. By assuming that [11] is a risk haplotype, the M step is derived to estimate the genotypic values (m j ) for each composite diplotype and the residual variance based on the calculated posterior probabilities bŷ where _ n n~n 11=10 zn 10=11 , € n n~n 11=00 zn 10=00 zn 01=01 zn 01=00 zn 00=00 : Iterations including the E and M steps are repeated at the higher hierarchy between equations (9) and (10) and at the lower hierarchy among equations (12) and (13) With the genotypic means of composite diplotypes, we can estimate the overall mean (m) and additive (a) and dominant genetic effects (d) due to the QTN detected, respectively, bŷ Model selection: The likelihood L(H p , H q | y, S) is formulated by assuming that haplotype [11][11] is a risk haplotype. However, a real risk haplotype is unknown from raw data (y, S). An additional step for the choice of the most likely risk haplotype should be implemented. The simplest way to do so is to calculate the likelihood values by assuming that any one of the four haplotypes can be a risk haplotype (Table 1). Thus, we obtain four possible likelihood values under different risk haplotypes; that is, (1) L 1 (Ĥ H p ,Ĥ H 1q jy,S) for [11], (2) L 2 (Ĥ H p ,Ĥ H 2q jy,S) for [10] 1,…,4). The largest likelihood value calculated is thought to correspond to the most likely risk haplotype.
In practice, it is also possible that there exist more than one risk haplotypes for a QTN. Relative to the bi-''allelic'' QTN with one risk haplotype, such a QTN is called a multi-''allelic'' QTN. If there are two risk haplotypes, we will have six composite diplotypes. Assuming that [11] (denoted by A 1 ) and [10] (denoted by A 2 ) are risk haplotypes and the remaining haplotypes [10] and [01] are non-risk haplotypes (denoted by A 3 ), then six composite diplotypes, expressed as can be specified according to the diplotype distribution as shown in Table 1. Totally, there are six such haplotype combinations for a two-SNP QTL, each combination corresponding to a likelihood value. Based on the calculated likelihoods, we can determine a most likely risk and non-risk haplotype combination. If there are three risk haplotypes, we will have 10 different composite diplotypes. The optimal risk and non-risk haplotype combination will be selected from three combinations based on the likelihoods.
The likelihood can be used as a criterion to select the optimal risk and non-risk haplotype combination when the number of risk haplotype is the same. However, when the number of risk haplotype is different, an AIC-or BIC-based model selection strategy [20] should be used because of different numbers of parameters being estimated in this case.
Hypothesis tests: We can test two major hypotheses in the following sequence: (1) the association between two SNPs by testing their linkage disequilibrium, and (2) the difference of a given haplotype from the remaining haplotypes by testing the significance of haplotype additive and dominant effects on the trait. The linkage disequilibrium between two given SNPs can be tested using two alternative hypotheses: The log-likelihood ratio test statistic for the significance of LD is calculated by comparing the likelihood values under the H 1 (full model) and H 0 (reduced model) using The LR 1 is considered to asymptotically follow a x 2 distribution with one degree of freedom. Diplotype or haplotype effects on the trait, i.e., the existence of a QTN, can be tested using the following hypotheses expressed as The log-likelihood ratio test statistic (LR 2 ) under these two hypotheses can be similarly calculated, where the tildes and hats denote the MLEs of parameters under the null and alternative hypotheses of (17), respectively. Although the critical threshold for determining the existence of a QTN can be based on empirical permutation tests, the LR 2 may asymptotically follow a x 2 distribution with two degrees of freedom, so that the threshold can be obtained from the x 2 distribution table.

Haplotyping a trait with multiple SNPs
Haplotype structure: The statistical method for QTN mapping is exemplified by a set of three SNPs, S 1 -S 3 , for a QTN. Two alleles 1 and 0 at each SNP are symbolized by r 1 , r 2 and r 3 , respectively. Eight haplotypes, [ In the F 2 population, there are 27 genotypes for the three SNPs. Let n r1 r 0 1 =r2 r 0 2 =r3 r 0 3 (r 1 $r9 1 ,r 2 $r9 2 ,r 3 $r9 3 = 1,0) be the number of offspring for a genotype. As seen in Table 2, the frequency of each genotype is expressed in terms of haplotype frequencies. Similar to equation (25), the phenotypic value of the trait for individual i is expressed, at the diplotype level, as where j i is the indicator variable defined as 1 if a diplotype considered is compatible with subject i and as 0 otherwise, u , and e i is the residual error distributed as N(0,s 2 ). Note that m and p stand for the maternally and paternally derived alleles, respectively.
By assuming [111] as a risk haplotype (labelled by A) and all the others as non-risk haplotypes (labelled by Ā ), Table 2 provides the formulation of genotypic values for three composite diplotypes, m 2 for AA, m 1 for AĀ and m 0 for Ā Ā . The haplotype effect parameters and residual covariance matrix are arrayed by a quantitative genetic parameter vector H q = (m 2 , m 1 , m 0 ,s 2 ).
Likelihoods and algorithms:With the above notation, we construct two likelihoods, one for haplotype frequencies (H p ) based on SNP data (S) and the other for quantitative genetic parameters (H q ) based on haplotype frequencies (H p ), phenotypic (y) and SNP data (S). They are, respectively, expressed as logL(HpjS)~constant logL(HqjHp,y,S)z where w . 's (w . = 12w) are defined below, and f j ( y j ) (j = 2, 1, 0) is a normal distribution density function of composite diplotype j. A two-stage hierarchical EM algorithm is derived to estimate haplotype frequencies and quantitative genetic parameters. At the higher hierarchy of the EM framework, we calculate the proportions of a particular diplotype within double or triple heterozygous genotypes (E step) by The calculated relative proportions by equation (21) were used to estimate the haplotype frequencies with  At the lower hierarchy of the EM framework, we calculate the posterior probabilities of a double or triple heterozygous individual i to be a particular diplotype (AĀ ) (E step), for which where [111] is assumed as the risk haplotype, expressed as With the calculated posterior probabilities by the above equation (23), we then estimate the quantitative genetic parameters, H q , based on the log-likelihood equations. These equations have similar, but more complicated, forms like equations (12) and (13). Hypothesis tests can be made for linkage disequilibria among three SNPs and haplotype effects. Four different linkage disequilibria, D 12 , D 13 , D 23 and D 123 , that describe the linkage among three SNPs can each be tested using the null hypotheses described by equation (21). The log-likelihood ratios for each hypothesis are thought to follow a x 2 distribution. R-SNP model: The idea for haplotyping a quantitative trait is described for two-and three-SNP models. It is possible that these models are too simple to characterize genetic variants for quantitative variation. With the analytical line for the two-and three-SNP sequencing model, a model can be developed to include an arbitrary number of SNPs whose sequences are associated with the phenotypic variation. A key issue for the multi-SNP sequencing model is how to distinguish among 2 r21 different diplotypes for the same genotype heterozygous at r loci. The relative frequencies of these diplotypes can be expressed in terms of haplotype frequencies. The integrative EM algorithm can be employed to estimate the MLEs of haplotype frequencies. A general formula for estimating haplotype frequencies can be derived.

RESULTS
The statistical model described above can be used to map and identify QTNs for a quantitative trait in an F 2 population. Because the marker data we have for mouse are microsatellites rather than SNPs, we use these microsatellite markers as a surrogate of SNPs for the purpose to demonstrate the utility of the model. Our marker data were from Vaughn et al.'s [21] study in which a linkage map composed of 19 chromosomes was constructed with 96 microsatellite markers for 502 F 2 mice (259 males and 243 females) derived from two strains, the Large (LG/J) and Small (SM/J). This map has a total map distance of ,1780 cM (in Haldane's units) and an average interval length of ,23 cM. The F 2 progeny was measured for their body mass at 10 weekly intervals starting at age 7 days. The raw weights were corrected for the effects of each covariate due to dam, litter size at birth, parity and sex [21]. Here, only adult body weights at week 10 are used for ''QTN'' analysis.
For each F 2 mouse, the parental origin of alleles at each marker can be discerned in molecular studies. Let L and S be the alleles inherited from the Large (LG/J) and Small (SM/J) strains, respectively. For any pair of markers, there are four different haplotypes, LL, LS, SL and SS, whose frequencies are accordingly denoted as By assuming all the four haplotypes as a risk haplotype, respectively, the above model allows for the estimates of haplotype frequencies by the EM iteration at the higher hierarchy and of composite genotypic values by the EM iteration at the lower hierarchy. The estimated haplotype frequencies are used to estimate linkage disequilibrium based on equation (14) and the recombination fraction (r) based on equation (4). This estimation process is moved from the first (M 1 -M 2 ) to last pair of markers (M 6 -M 7 ) on chromosome 1 and then from chromosome 1 to 19. Table 3 tabulates the results of the MLEs of haplotype frequencies and log-likelihoods under the assumptions of different risk haplotypes. A total of 96 markers are sparsely located on 19 mouse chromosomes, with the estimated recombination fractions from the linkage disequilibrium model [8] consistent with those obtained from the linkage model [21]. Significant likelihood ratios for testing haplotype effects were determined by critical values obtained from the x 2 -square distribution with two degrees of freedom with a Bonferroni adjustment to the type I error. The adjusted critical values for the two-and three-marker QTN models are 18.20 and 18.76, respectively, at the 5% significance level. Significant haplotype effects are detected for a total of eight marker pairs (Table 3), which include one pair on chromosome 4, two consecutive pairs on chromosome 6, four consecutive pairs on chromosome 7 and one pair on chromosome 14. For some pairs, multiple significant risk haplotypes were detected. Risk haplotypes purely composed of alleles inherited from the LG/J or SM/J parent exert a positive or negative additive effect on body weight, respectively. Based on the relative values of estimated additive and dominant effects, the significant marker pairs detected display partial dominant effects ( Table 3).
The results from the three-marker model are basically consistent with those from the two-marker model ( Table 4). The advantage of the three-marker model is that it incorporates the interferences between adjacent marker intervals into the estimation process and, thus, can potentially increase the estimation precision of haplotype effects.

DISCUSSION
Quantitative trait locus (QTL) mapping aims to identify narrow chromosomal segments for a quantitative trait by using a statistical method, and has proven its value to study the genetic architecture of the trait in a variety of species [6][7][8]. The limitations of this technique lie in its inability to characterize the structure and organization of DNA sequences and statistical difficulty in deriving the distribution of test statistics under the null hypothesis of no QTL [22]. At least partly for these reasons, despite thousands of QTL reported for different traits and populations, a very small portion of them have been cloned [9]. With the completion of the genome projects for several important organisms, a new line of thought in the post genomic era has begun to emerge for the identification of specific combinations of nucleotides or haplotypes that contribute to a complex quantitative trait [13,23].
Theory and methods for haplotype discovery have well been established for natural populations [13] in which the non-random association among different single nucleotide polymorphims (SNP), specified by the coefficients of linkage disequilibria, lays a foundation for the mixture model of haplotyping a quantitative trait. In this article, we derived a statistical model for detecting haplotypes and estimating their effects on quantitative variation of a trait in experimental crosses. We used the principle of linkage disequilibrium analysis to characterize the linkage among different markers that is usually described by the recombination fractions in a commonly used F 2 population, initiated with two inbred lines. We established an interchangeable relationship between the linkage and linkage disequilibrium. The merit of this relationship in trait haplotyping includes the incorporation of interferences between adjacent marker intervals into the estimation and test of haplotype effects when multiple markers are modelled simultaneously.
The haplotyping model developed in this article was used to analyze a published F 2 population of mouse [21], but we used microsatellite markers as the surrogate of SNPs so that we can detect the effects of haplotypes constructed by microsatellite alleles. The whole-genome of mouse was scanned for haplotype effects on body weight by a two-and multi-marker model, respectively. Consistent results were observed from the two models, which suggests that four regions in mouse chromosomes 4, 6, 7, and 14 contribute to variation in body weight. These findings are in a good agreement with those from traditional interval QTL mapping [21]. But our haplotype discovery is more informative in terms of the characterization of specific haplotype structure and organization responsible for trait variation.
We have proposed a new model for haplotyping a quantitative trait in the F 2 progeny population. The tenet of the model can be extended to haplotype a complicated trans-generational pedigree, founded with multiple original parents and involving individuals with different relatedness. The model can also be modified to dissect the epistatic effects of different genes [23] and the interaction of genes with environment. For these extensions, haplotype selection aimed to detect the risk haplotypes that are expressed differently from the others present many challenges, but is crucial for the facilitation of the process of detecting the association between haplotype diversity and phenotypic variation.
Our haplotyping model offers a powerful tool for positional cloning of QTL that are important for a complex trait. Flint et al. [9] reviewed the potential of currently available cloning strategies, such as probabilistic ancestral haplotype reconstruction, Yin-Yang crosses and in silico analysis of sequence variants, to identify genes that underlie QTL in rodents. Our model, in conjunction with these strategies, may open a new gateway for the illustration of a detailed picture of the genetic architecture for a complex trait.