Individual genotyping remains prohibitively expensive for a genome scan. One method to reduce cost is to pool DNA from individuals with extreme phenotypic values and to measure the allele frequency difference between pools (Barcellos et al., 1997 ; Daniels et al., 1998; Fisher et al., 1999; Hill et al., 1999; Shaw et al., 1998; Stockton et al., 1998; Suzuki et al., 1998). Initial attention focused on pooled designs for dichotomous traits and case- control studies (Risch and Teng 1998). More recently, pooled tests have been discussed for quantitative traits, a more appropriate model for diseases such as obesity and hypertension.

In the absence of experimental error, the optimal design for an unrelated population is to compare frequencies between pools of the most extreme 27% of individuals ranked by phenotypic value, retaining 80% of the information of individual genotyping (Bader et al., 2001). Experimental sources of error, primarily allele frequency measurement error, degrade the test power (Jawaid et al., 2002).

Population stratification poses a second challenge to practical use of pooled tests for human populations. Genomic control methods, developed to reduce stratification effects in genotype-based association tests (Devlin and Roeder 1999; Pritchard and Rosenberg 1999; Pritchard et al. 2001 ; Zhang and Zhou, 2001), are not directly applicable to pooled tests.

Here we present optimized pooled DNA test designs, including family-based tests robust to stratification. Estimates of test power explicitly include allele frequency measurement error. This distinguishes our treatment from prior theoretical work, permits the optimization of test design as a function of known parameters, and provides a bridge to

experimentalists seeking practical guidance for whether to attempt and how to perform pooled association tests.

SUMMARY OF THE INVENTION The invention is drawn to a method for detecting an association in a population of unrelated individuals between a genetic locus and a quantitative phenotype, wherein two or more alleles occur at the locus, and wherein the phenotype is expressed using a numerical phenotypic value whose range falls within a first numerical limit and a second numerical limit. This method comprises the steps of : a) obtaining the phenotypic value for each individual in the population; b) determining the minimum number of individuals from the population required for detecting the association using a non-centrality parameter; c) selecting a first subpopulation of individuals having phenotypic values that are higher than a predetermined lower limit and pooling DNA from the individuals in the first subpopulation to provide an upper pool; d) selecting a second subpopulation of individuals having phenotypic values that are lower than a predetermined upper limit and pooling DNA from the individuals in the second subpopulation to provide a lower pool; e) for one or more genetic loci, measuring the frequency of occurrence of each allele at said locus in the upper pool and the lower pool; f) for a particular genetic locus, measuring the difference in frequency of occurrence of a specified allele between the upper pool and the lower pool; and g) determining that an association exists if the allele frequency difference between the pools is larger than a predetermined value.

In one embodiment of the invention, the difference in frequency of occurrence of the specified allele has associated with it an error of measurement. In one aspect of the invention the error of measurement is 0.04. In another, the error of measurement is 0.01.

In another embodiment of the invention, the predetermined lower limit is set so that the upper pool ranges from including the highest 37% of the population to including the highest 19% of the population and the predetermined upper limit is set so that the lower pool ranges from including the lowest 37% of the population to including the lowest 19% of the population. In another aspect of the invention, the predetermined lower limit is set so that the upper pool includes the highest 27% of the population and the predetermined upper limit is set so that the lower pool includes the lowest 27% of the population.

In another embodiment of the invention, the genetic locus has two alleles.

In another embodiment of the invention, the population includes individuals who may be classified into classes. In one aspect of the invention, the classes are based on an age group, gender, race or ethnic origin. In another aspect of the invention, the members of a class are included in the pools.

In another embodiment of the invention the method is used for determining the genetic basis of disease predisposition.

In another embodiment of the invention, the genetic locus which is analyzed for determining the genetic basis of disease predisposition contains a single nucleotide polymorphism.

BRIEF DESCRIPTION OF THE DRAWINGS Figure 1. The information retained by the between-family pooled test design, expressed as a fraction of the information from individual genotyping followed by a between- family test, is depicted sibships of size 4,2, and 1, each population having 1000 total individuals. The optimal pooling fraction, indicated by an arrow, shifts to lower values as the number of sibs per family decreases. The optimal fraction and corresponding information retained also shift to lower values as the minor allele frequency decreases, with results shown for frequencies 0.1 and 0.01. The raw measurement error is 0.01.

Figure 2. The optimal number of sibs to select from each family (top panel) and the information retained relative to individual genotyping (bottom panel) are shown for sibship sizes 2-5,6,8,16, and 32 as a function of the scaled measurement error K. For sibships through 5, it is always optimal to select just the highest and lowest sib.

Figure 3. The optimal fraction of families to select (top panel) and information retained (lower panel) are displayed for sibships of size 2 through 6 as a function of the scaled measurement error K.

Figure 4. The optimal pooling fraction (top panel) and information retained (bottom panel) for between-family and within-family tests of a population of 500 sib-pairs are shown as a function of raw measurement error for marker frequencies 0.5 and 0.01. The within-family tests include pre-selection of discordant-like families.

Figure 5. The optimal pooling fraction (top panel) and the information retained (bottom panel) from exact numerical calculations (solid line) and an analytical fit (dashed line) are displayed as a function of the normalized measurement error K. The fit coincides with the exact results for the information retained.

DETAILED DESCRIPTION We present optimized designs for pooled DNA tests conducted on a population of Nls families, each a sibship of size s (N total individuals). The genotypic correlation within a sibship is denoted r, with typical values of 1/4,1/2, and 1 for half-sibs, full-sibs, and monozygotic twins. Sibships may also represent inbred lines; in this case, r is the genetic correlation within each line. Sibs in different families are assumed to have uncorrelated genotypes.

To conduct a pooled DNA test for association of a particular allele Al with a quantitative trait, individuals are selected for an upper pool, comprising higher phenotypic values, and a lower pool, comprising lower phenotypic value, using designs reminiscent of selection strategies for optimizing breeding value and for QTL mapping (Hill 1971 ; Kimura and Crow 1978; Ollivier et al. 1997). We restrict attention to balanced designs in which each pool has fN individuals, with f < 0.5 defined as the pooling fraction. Balanced designs are favored when high and low phenotypes are treated symmetrically ; asymmetry can favor unbalanced designs (Jawaid et al., 2002).

We consider four designs: (i) unrelated individuals (s = 1), in which thefiv individuals having highest and lowest phenotypic values are selected for the upper and lower pools respectively; (ii) between-family, in which all s sibs from the fN/s families having highest and lowest mean phenotypic values are selected for the upper and lower pools; (iii) within-family, in which the s'sibs having highest and lowest phenotypic values within each family are selected for the upper and lower pools, yielding a pooling fraction f = s'ls ; (iv) within-family with pre-selection of discordant families, in which a fractionJ'of families with greatest within-family phenotypic variance are selected, where Xs is the phenotype of sib s and X is the family mean, then the extreme high and low sib within each selected family are selected for the upper and lower pool for a final pooling fraction f = f lN.

A suitable statistic for a two-sided test for each design is <BR> <BR> <BR> 2 (_ p) 2<BR> <BR> <BR> Var (-)' where the estimated frequencies of allele A1 in the upper and lower pools are denoted Pu and PL. The variance is the sum of three terms, Var (pu-PL) = VS + VC + VM- The sampling variance Vs represents the unavoidable error in estimating the population

frequency from a finite sample. The concentration variance Vc arises from sample-to-sample DNA concentration variance within a pool. The measurement variance is VM = vs 2, where s is the experimental allele frequency measurement error for each pool. We assume that the three sources of variation are independent, which should be justified when individual and pooled DNA samples are treated uniformly. In an ideal experiment, Vc and VM vanish, and the total variance is Vs.

Under the null hypothesis, Z2 has a % 2 diskibution with one degree of freedom.

Under the alternate hypothesis, the tested marker is assumed to be a bi-allelic quantitative trait locus (QTL) with alleles Al and A2 occurring at frequencies p and (1-p) _ q. For between-family tests, the alleles are assumed to be in Hardy-Weinberg equilibrium and the population is assumed to be random mating; these assumptions may be relaxed for within- family tests. The variance of the allele frequency per individual is a p = pq l 2. For each design, the allele frequency is estimated as p = (#U + PL)/2. The estimated variance of the allele frequency per individual is denoted 6 P and equals p (l-p)/2.

The mean phenotypic effects are mG = a, d, and-a for genotypes G = AlAl, AlA2, and A2A2, respectively. The dominance ratio d/a describes the inheritance mode with typical values-1,0, and 1 for pure recessive, additive, or dominant inheritance. The proportion of trait variance accounted for by the QTL is denoted, The mean QTL effect is m = (p-q) a + 2pqd. Phenotypic values are assumed to be normally distributed for each genotype with mean µG = mG - m and residual variance o-2 =1-s Q arising from all genetic and environmental factors other than the QTL. The distribution of phenotypic values in the population is a mixture of three normal distributions with overall mean 0 and variance 1. The total phenotypic correlation between sibs from genetic factors (including the QTL) and environmental factors is termed t.

1 The non-centrality parameter (NCP), NCP = [E(#U-#L)]2/Var(#U-#L),

measures the information provided by a pooled DNA test. The notation E (0) is the expectation of an observable O. The approach followed below is to evaluate the numerator of the NCP as a function of the model parameters, providing accurate analytical results when possible and simulation results otherwise. For the denominator of the NCP, analytical results are obtained for the null hypothesis. For the alternative hypothesis, the expected allele frequencies for each pool have offseting changes from p to p 6p (see Methods for derivation), and the value of the denominator decreases by a small value proportional to (6plu)'. We make a conservative approximation by ignoring the change and using the null hypothesis denominator for the alternative hypothesis as well. In this case, the NCP equals (zα/2-z1-ß)2, where a and (3 are the type I and II error rates for the two-sided test.

Maximizing the NCP optimizes the test.

The denominator of the NCP is shown in the Methods to have the form where r is the coefficient of variation for DNA concentration. The constant G depends only on the family structure and equals 1 for pools of unrelated individuals, sR for the between-family design, and (1-r) for both within-family designs; the standard notation R relates the sib genotypic correlation r to family-based variance components, R =-. #[1+(s-1)r]. s Typically-c is less than 10% ; 12 may usually be ignored relative to G. The term K2 is used as shorthand, Referred to as the scaled measurement error, K represents the raw measurement error, E, scaled by the remaining sources of error in the allele frequency difference. In practice, K can be calculated prior to pooling because it depends on known quantities.

The numerator of the NCP is shown in the Methods to have the form

where (z) is the normal density (27r)-1/2 exp (-Z2/2), <) (z) is the cumulative normal probability and (D-' (z) its functional inverse. The constant F equals 1 for pools of unrelated individuals, R21T for between-family pools, and (l-R) 2/(l-T) for within-family pools without pre-selection. For the within-family design using discordant-like pre-selection, F= (1- r) 2/2 (1-t) for sib-pairs (expressions for larger sibships are unwieldy). The term R has the same definition as before, and T is the standard factor relating the sib phenotypic correlation t to family-based variance components, T = 1#[1+(s-1)t]. s Combining terms, the analytical result for the NCP, valid for small QTL effect, is The first factor is identical to the NCP for an association test performed by individual genotyping on a population of N unrelated individuals; the second factor, withr = 0, is the correction for individual genotyping a population of Nls families each having s sibs and then performing either a between-family test, with F/G = RlsT, or a within-family test, with FIG = (s-1) Rls (1-T). The third factor represents the fraction of information retained when the association test is performed by pooling instead of individual genotyping, and maximizing this factor with respect to the pooling fractionfprovides the optimal pool size.

When the measurement error s = 0, tests are optimized withf= 0.27 and 80% of the information is retained (Bader et al. 2001). As s increases, the maximum information that can be retained is determined entirely by the single collective term K.

Expressions for F, G, and K2 are summarized in Table I, and we now provide examples of each family-based design. Information retained by the between-family design is depicted in Fig. 1, with results for 3 sibship sizes: sib-quads, sib-pairs, and unrelated individuals, each population having 1000 total individuals. The optimal pooling fraction, indicated by an arrow, shifts to lower values as the number of sibs per family decreases. The optimal fraction and corresponding information retained also shift to lower values as the minor allele frequency decreases, with results shown for frequencies 0.1 and 0.01. The raw measurement error s = 0.01 in this example, and the pooling fraction and information retained would decrease for larger s (see Fig. 4 for examples of changing s).

In Fig. 2, the optimal number of sibs to select from each family (top panel) and the information retained relative to individual genotyping (bottom panel) are shown as a function of the scaled measurement error K for sibship sizes of 2-5,6,8,16, and 32. For sibships through 5, it is always optimal to select just the highest and lowest sib. For larger families and small measurement error, the top and bottom quarters of the sibs are pooled and 80% of the information is retained. The pooling fraction and information decrease as the measurement error increase.

Within-family tests can be improved by pre-selection of discordant-like families, as shown in Fig. 3. The optimal fraction of families to select (top panel) and information retained (bottom panel) are displayed for sibships of size 2 through 6 as a function of the scaled measurement error K (results determined by computer simulation). The fraction of families and information retained both decrease as K increases. Discordant pre-selection has the greatest benefit for sib-pairs: for the smallest values of K, only 56% of families are selected, retaining 80% of the information; had all families been used, only 60% of the information would have been retained. Pre-selection is less important for trios and larger sibships.

In Fig. 4, the optimal pooling fraction (top panel) and information retained (bottom panel) using between-family pools and using within-family pools with discordant-like pre- selection are displayed for a population of 500 sib-pairs (1000 individuals) as a function of the raw measurement error s. Results are shown marker frequencies 0.5 and 0.01. With no measurement error, the optimal pooling fraction of 0.27 retains 80% of the information in each case. As measurement error increases, the optimal pooling fraction decreases, as does the information retained.

The information loss increases for rarer alleles and is worse for the within-family test than for the between-family test. This behavior can be deduced from the scaled error K 2, which is inversely proportional to the allele frequency sampling variance. Since the sampling variance is 3x smaller within-family vs. between-family, K z is 3x larger, 4Ns Zpl-p) vs.

4Ns 2/3p (1-p), and more information is lost. The inverse dependence of K 2 on the allele frequency explains the decrease in power for rare alleles.

Because the allele frequency difference between sibs is uncorrelated from their allele frequency mean, the between-family and within-family tests are independent estimators of CA even when individuals contribute their DNA under both designs. The NCP of a combined

test is the sum of the NCPs for each test and it too follows a % 2 distribution with 1 degree of freedom. In practice, estimates for CSA may obtained by inverting the expressions for E (PU-p/) provided in Table I, then weighting each estimator by the inverse of its variance.

Population stratification may be indicated by a difference between the estimates for CA from a between-family and within-family test. In the absence of stratification, the difference follows a normal distribution with variance where the"+"and"-"subscripts refers to the between-family and within-family designs respectively, Yi = #[#-1(1-f~)], and V represents the total variance, Vs + Yc + VM, for each design. When stratification is indicated, the between-family estimate of GA may be unreliable but the within-family estimate remains robust.

In Fig. 5, the optimal pooling fraction (top panel) and the information retained (bottom panel) are displayed as a function of the scaled measurement error K. The information retained is calculated assuming no concentration variance. In addition to the numerically calculated results, an accurate fit is shown using the functional form f = 1-d) [A- (3/A) ln A-0. 067], with A justification for this functional form is provided in the Methods. The greatest deviations for the pooling fraction are at K = 0.5, where the fit yields a pooling fraction that is 0.006 too high, and at K = 3.5, where the fit is 0.01 too low. The information retained using the analytical value for the pooling fraction coincides with the exact numerical results on the scale of the figure. The experimental measurement error s corresponding to the scaled error K depends on the population structure and marker frequency. For example, for a population of 500 cases, 500 matched unrelated controls, and 10% marker frequency, s = 0.0067K is the raw error corresponding to K.

Based on the pooled designs described above, we outline a prospective study using 100,000 markers to detect QTLs with a 1% effect. If 100 false-positives are permitted from pooled tests (the false-positives may be resolved using individual genotyping) and 80% power is required, the NCP is 17. We assume pooling of discordant sib-pairs to protect against stratification effects. At the scaled error K = 1 where the pooled tests are still close to

maximum power, the pooling fraction would be 21%, 65% of the information of a population would be retained, and a population of 2600 individuals would be required. The raw measurement error corresponding to K = 1 for this population size is 0.005 for an allele with 50% frequency and 0.002 for an allele with 5% frequency, 5x to 10 x more precise than achieved by current-day instrumentation.

We can account for current-day precision by setting K = 10, which from Fig. 5 is seen to retain 7.7% of the information and corresponds to a pooling fraction of 1.6% of a total population of 22,000. In this case, the precision required for a pooled test is 0.017 for an allele with 50% frequency and 0.007 for an allele with 5% frequency. These precisions are within the range of current performance, especially if repeated measures are used to decrease the effective measurement error. The cost to collect and score such a population for multiple disease-related phenotypes would be under $50 million. Selection schemes could then be applied to generate pools for each phenotype in turn.

As noted previously, pooled tests perform worse for within-family tests and rare alleles, and may therefore be difficult to apply to disease-risk variants under negative selection pressure. The loss of power may be less severe for pharmacogenetic studies of variants affecting drug response, where selection pressure is absent, and for test crosses of model organisms (Grupe et al. 2001) or agricultural species whose marker frequencies are under experimental control.

The analysis provided here for quantitative traits may be extended to threshold characters yielding dichotomous classifications of a population. For case-control classification, the disease prevalence corresponds to the pooling fraction. When the quantitative character is available for measurement, it is approximately 4x more efficient to compare unrelated individuals with extremely high vs. extremely low characters than to compare the derived cases vs. controls (Bader et al. 2001).

In summary, we have derived the optimal pooling fractions for within-family and between-family tests of association. With ideal instrumentation, 80% of the information is retained and the optimal pooling fraction is 27%. As allele frequency measurement error increases, the optimal pooling fraction and the information retained both decreases. The information loss is more severe for low-frequency alleles and for within-family tests. The optimal pooling fraction depends on a single parameter representing the measurement error, and optimized pooling designs are provided as a function of this parameter.

EXAMPLES Example 1: Sampling variance and concentration variance Let pi represent the frequency of allele A, for individual i, either 0,1/2, or 1, and cl represent the concentration of DNA contributed by this individual to a pool of n individuals.

Neglecting measurement error, the allele frequency p for the pool is

which defines the relative concentration error #ci'. The terms 6pi and oc'are uncorrelated, and each has expectation zero. Furthermore, the sum of the 8terms is constrained to be zero. The variance of p* is

We have used with the concentration coefficient of variation defined as # # [Var(ci)]1/2/c0 and the genotypic correlation between a pair of individuals defined as rjj.

For the between-family design, a pool of n individuals contains n/s sibships of size s and genotypic correlation r. The result for Var(p*) is with R = (l/s) [l + (s-l) r]. Since the individuals in the upper and lower pools are unrelated, Vs + Vc = 2 Var (p).

For a within-family design, the allele frequency difference between pools is where i and j label individuals in the upper and lower pools respectively. The variance is Example 2: Expected allele frequency difference and non-centrality parameter The genotype-dependent phenotype distribution is defined using a variance components model, Xki=Yk+Yki+µki.

Family and individual effects are normally distributed with mean zero and variance The family index is k, the sib index is i, and the individual phenotypes XkZ'are the sum of Y., the family effect excluding the QTL, Yli, the individual effect excluding the QTL, and Slkí the QTL effect u (Gk.) for sib i. The total phenotypic correlation between sibs is t. Both r and u relate to the genetic background shared between sibs, r being the genotypic correlation (1 for monozygotic twins, 1/2 for full sibs, 1/4 for half sibs) and u being the shared genotype expectation (1 for monozygotic twins, 1/4 for full sibs, 0 for half sibs) (Falconer and Mackay 1996).

The observed phenotypes Xki are re-expressed as family means and individual deviations from family means, #Xki = Xki - Xk# Similar quantities are defined for the QTL effects, #µki = µki - µk#, and the variances of the observed quantities excluding QTL effects are Var(#Xki-#µki)=(1-T)#R2 When the QTL effects are small, T (l/) [l + (s-l) t] is an accurate approximation.

The probability that sibling 1 from family k with genotypes G = (G1, G2,..., Gs) is selected for the upper pool is 1-#[(X'-µG)/#], where # (z) is the cumulative normal probability. The variable under selection, denoted, is either Xk. (between-family pools) or AXkl (within-family pools); µG is either µk# (between-family pools) or #µkl (within-family pools); the variance of X-µG is #2, either To R (between-family pools) or (1-T) R2 (within-family pools); and is the selection threshold applied to X. Because the labeling of sibs is arbitrary, the fractionf of individuals selected for pooling is equal to the probability that sib 1 is selected, i. e. the probability that X is greater than the selection threshold, where Pr (G) is the probability of observing the sibship genotypes G.

To calculate the allele frequency of the selected individuals, the threshold is required as a function off Numerical inversion may be applied to the above equation.

Alternatively, when the QTL effect is small (ju. G < #), the linear approximation #[(X'-µG)/#]##(X'/#)-(µG/#)#(X'/#) is accurate, where (z) = d (z)/dz is the normal probability density. The terms linear in LG cancel in the sum over G, yielding f =1-#(X'/#).

The expected allele frequency of the resulting pool is

where PG represents the allele frequency of sib 1. Using the linear expansion for #[(X'-µG)/#] yields An analogous expression for the lower pools gives a symmetric result, yielding

where X'/# has been replaced by #-1(1-f).

The expectation of the correlation between p and µ for an individual is E(pµ) = p2[a-(p-q)a-2pqd]-2pq#1/2#[d-(p-q)a-2pqd] = pq[a-(p-q)d] =#p#A Similarly, the correlation between sibs i and j is E(piµj)=rij#p#A, where rij is their genotypic correlation. Summing over sibs yields either R#p#A (between-family pools) or (1-R) #p#A (within-family pools) for E (PG PG), with R = (l/s) [l + (s-l) r] as before.

Selecting discordant-like sib-pairs is equivalent to selection based on ##Xki#, and the within-family analytical results are directly applicable. For larger families, discordant-like families are pre-selected in decreasing rank order of the within-family phenotypic variance summed over siblings s.

We have ascertained that the analytical results for the NCP are virtually indistinguishable from exact numerical results when the QTL effect is 5% or less of the trait variance. For larger effects, roughly when the effect size aA2 approaches the minor allele

frequency, the genotype-dependent phenotype distributions become resolved, transforming a complex trait into Mendelian trait amenable to traditional linkage analysis.

Example 3: Analytical fit for the optimal pooling fraction Optimizing the pooling fraction is equivalent to maximizing the objective function I = 2y2/(f+f2#2), where y is shorthand for #[#-1(1-f)]. Writing f as 1-#(z) and optimizing using dildo = 0 yields y#(1+2f#2)-2zf#(1+f#2)=0.

We have used y = + (z), dy/dz =-yz, and df/dz =-y.

When K2 is large, z is also large, andfmay be replaced by its asymptotic expansion for large z, f = y- (Z-1-Z-3). With this substitution, the optimum satisfies z3 =1.

2y#2 Taking the natural logarithm of both sides and equating exponents, When K and z are both large, the term 3 In z is asymptotically small, giving An improved fit is obtained by perturbation theory by writing z=B(#)[1+b(#)], where Substituting this expression for z into J (z) and simplifying, B2b+3ln[B(1+b)]=0, which gives the asymptotic form b = (3/B2) 1nB, or z-B- (3/B) InB.

For clarity, the functional dependence of B and b on K has been suppressed.

The asymptotic form provides a good fit when K is much larger than 1 but not for smaller values. Since the asymptotic behavior for large K is not affected by introducing terms of lower order in K, the fit can be improved for small K without degrading the fit at large K by writing z = a-(3/A)lnA+a1, where

The constants al, a2, and a3 are then selected to fit the exact numerical results at particular values of K. Fitting the results z = 0.612 at # = 0 and z = 0.8047 at K = 1 provides the particular parameters al =-0. 067, a2 = 2, a3 = 3.

Table I. The non-centrality parameter for family-based pooled DNA designs a Design F G Unrelated individuals 1 1 Between-family R2lT sR Within-family (1-R)2/(1-T) 1-r <BR> <BR> Within-family, discordant (1-r)2/2(1-t)<BR> <BR> <BR> <BR> pre-selection a The non-centrality parameter (NCP) is [E(#U-#L)]2/Var(#U-#L). The A p R), where F is provided for each design, f is the pooling fraction, #A2 and #R2 are the additive and residual variance for a QTL with allele frequency p, #p2 is p(1-p)/2, #(z) is the normal probability density and (D (z) is the cumulative normal probability. The denominator of the NCP is [2 (G p INf] + 2, where G is provided for each design, # is the coefficient of variation for DNA sample concentrations in the pool, N is the total number of individuals before selection, and s is the raw measurement error. The combined expression for the NCP is where #2 is N#2/[(G+#2)#p2] and K is termed the scaled error. Each sibship has s sibs with genotypic correlation r and phenotypic correlation t ; R and T are(1/s)[1+(s-1)r] and (1/s)[1 + (s-1) t], respectively. b Analytical results are for sib-pairs only. For larger families see numerical results (Fig. 3).

References Abecasis GR, Noguchi E, Heinzmann A, Traherne JA, Bhattacharyya A, Leaves NI, Anderson GG, Zhang Y, Lench NJ, Carey A, Cardon LR, Moffatt MF, Cookson OC (2001) Extent and distribution of linkage disequilibrium in three genomic regions. Am J Hum Gen 68: 191-197 Ardlie KG, Kruglyak L, Seielstad M (2002) Patterns of linkage disequilibrium in the human genome. Nat Rev Genet 3: 299-309 Bader JS, Bansal A, and Sham P (2001) Efficient SNP-based tests of association for quantitative phonotypes using pooled DNA. Genescreen (in press) Barcellos LF, Klitz W, Field LL, Tobias R, Bowcock AM, Wilson R, Nelson MP, Nagatomi J, Thomson G (1997) Association mapping of disease loci, by use of a pooled DNA genomic screen. Am J Hum Gen 61: 734-747 Collins FS, Guyer MS, Chakarvarti A (1997) Variations on a theme: cataloging human DNA sequence variation. Science 274: 1580-1581 Daniels J, Holmans P, Williams N, Turic D, McGuffin P, Plomin R, Owen MJ (1998) A simple method for analysing microsatellite allele image patterns generated from DNA pools and its applications to allelic association studies. American Journal of Human Genetics 62: 1189-97 Devlin B, Roeder K (1999) Genomic control for association studies. Biometrics 55: 788-808 Falconer (1965) The inheritance of liability to certain diseases estimated from the incidence among relatives. Ann Hum Gen 51: 227-33 Falconer DS, MacKay TFC (1996) Introduction to quantitative genetics. Boston: Addison-Wesley Fisher PJ, Turic D, Williams NM, McGuffin P, Asherson P, Ball D, Craig I, Eley T, Hill L, Chorney K, Chorney MJ, Benbow CP, Lubinski D, Plomin R, Owen MJ (1999) DNA pooling identifies QTLs on chromosome 4 for general cognitive ability in children. Hum Mol Gen 8: 915-22 Grupe A, Germer S, Usuka J, Aud D, Belknap JK, Klein RF, Ahluwalia MK, Higuchi R, Peltz G (2001) In silico mapping of complex disease-related traits in mice. Science 292: 1915-1918 Hill, WG. (1971) Design and efficiency of selection experiments for estimating genetic parameters. Biometrics 27 : 293-311 Hill L, Craig IW, Asherson P, Ball D, Eley T, Ninomiya T, Fisher PJ, Turic D, McGuffin P, Owen MJ, Chorney K, Chorney MJ, Benbow CP, Lubinski D, Thompson LA,

Plomin R (1999) DNA pooling and dense marker maps: a systematic search for genes for cognitive ability. Neuroreport 10: 843-848 Jawaid A, Bader JS, Purcell S, Cherry SS, Sham P (2002) Optimal selection strategies for QTL mapping using pooled DNA samples. European Journal of Human Genetics (in press) Kimura M, Crow JF (1978) Effect of overall phenotypic selection on genetic change at individual loci. Proc Natl Acad Sci USA 75: 6168-6171 Ollivier L, Messer LA, Rothschild MF, Legault C (1997) The use of selection experiments for detecting quantitative trait loci. Genet Res, Camb 69: 227-232 Ott J (1999) Analysis of Human Genetic Linkage. Third edition. Johns Hopkins University Press, Baltimore Pritchard JK, Stephens M, Rosenberg NA, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945-959 Pritchard JK, Rosenberg NA (1999) Use of unlinked genetic markers to detect population stratification in association studies. Am J Hum Gen 65: 220-228 Reich DE, Cargill M, Bolk S, Ireland J, Sabeti PC, Richter DJ, Lavery T, Kouyoumjian R, Farhadian SF, Ward R, Lander ES (2001) Linkage disequilibrium in the human genome. Nature 411: 199-204 Risch N and Teng J (1998) The relative power of family-based and case-control designs for linkage diequilibrium studies of complex human diseases 1. DNA pooling.

Genome Res 8: 1273 Risch N, Merikangas K (1996) The future of genetic studies of complex human diseases. Science 273: 1516-1517 Satten GA, Flanders DW, and Yang Q (2001) Accounting for unmeasured population substructure inb case-control studies of genetic association using a novel latent-class model.

Am J Hum Gen 68: 466-477 Schork NJ, Nath SK, Fallin D, Chakarvati A (2000) Linkage disequilibrium analysis of biallelic DNA markers, human quantitative trait loci, and threshold-defined case and control Subjects. Am J Hum Gen 67: 1208-1218 Sham PC, SS Cherny, S Purcell, and JK Hewitt (2000) Power of linkage versus association analyses of quantitative traits, by use of variance-components models, for sibship data. Am J Hum Gen 66: 1616-1630 Shaw SH, Carrasquillo MM, Kashuk C, Puffenberger EG, Chakravarti A (1998) Allele frequency distributions in pooled DNA samples: applications to mapping complex disease genes. Genome Res 8: 111-123

Stockton DW, Lewis RA, Abboud EB, Al Rajhi A, Jabak M, Anderson KL, Lupski JR (1998) A novel locus for Leber congenital amaurosis on chromosome 14q24. Human Genetics 103: 328-333 Suzuki K, Bustos T, Spritz RA (1998) Linkage disequilibrium mapping of the gene for Margarita Island ectodermal dysplasia (ED4) to 1 lq23. American Journal of Human Genetics 63: 1102-1107 Zhang S, Zhao H (2001) Quantitative similarity-based association tests using population samples. American Journal of Human Genetics 69 : 601-614