Allele frequency misspecification: effect on power and Type I error of model-dependent linkage analysis of quantitative traits under random ascertainment
© Mandal et al; licensee BioMed Central Ltd. 2006
Received: 14 December 2005
Accepted: 20 April 2006
Published: 20 April 2006
Studies of model-based linkage analysis show that trait or marker model misspecification leads to decreasing power or increasing Type I error rate. An increase in Type I error rate is seen when marker related parameters (e.g., allele frequencies) are misspecified and ascertainment is through the trait, but lod-score methods are expected to be robust when ascertainment is random (as is often the case in linkage studies of quantitative traits). In previous studies, the power of lod-score linkage analysis using the "correct" generating model for the trait was found to increase when the marker allele frequencies were misspecified and parental data were missing. An investigation of Type I error rates, conducted in the absence of parental genotype data and with misspecification of marker allele frequencies, showed that an inflation in Type I error rate was the cause of at least part of this apparent increased power. To investigate whether the observed inflation in Type I error rate in model-based LOD score linkage was due to sampling variation, the trait model was estimated from each sample using REGCHUNT, an automated segregation analysis program used to fit models by maximum likelihood using many different sets of initial parameter estimates.
The Type I error rates observed using the trait models generated by REGCHUNT were usually closer to the nominal levels than those obtained when assuming the generating trait model.
This suggests that the observed inflation of Type I error upon misspecification of marker allele frequencies is at least partially due to sampling variation. Thus, with missing parental genotype data, lod-score linkage is not as robust to misspecification of marker allele frequencies as has been commonly thought.
In model-based linkage analysis, one assumes that the causative mechanisms of both the trait and marker phenotypes are known without error, including the number of loci involved, the mode of inheritance and allele frequencies. With the advancement of genetic research in complex traits, we are dealing more with uncertain mode of inheritance and unknown parameters (allele frequency and degree of dominance) which ought to be known to apply model-dependent linkage analysis appropriately. Several studies have been reported in the literature about the implications of misspecifying the underlying genetic model of the trait in linkage analysis. Clerget-Darpoux et al.  quantified the effects of using wrong allele frequencies, penetrance and degree of dominance on the linkage test and on the recombination fraction estimate under the assu mption of a single locus model in lod-score linkage analysis. Goldin and Weeks  reported that a biased estimate of the recombination fraction is introduced in lod-score analysis by reducing the power to detect an existing linkage when underlying disease model parameters are incorrectly specified in the case of a common dichotomous trait. It is important to mention that the accurate estimation of marker allele frequencies is essential in the linkage analysis when reconstruction of the genotypes of key family members is needed due to the unavailability of biological samples (e.g., in the case of diseases with late age of onset, such as lung cancer, when parental data often are not available). So far, a few studies have been presented in the literature about the effect of using incorrect allele frequencies in linkage analysis. Ott  and Freimer et al.  reported an increase of the false positive rate in linkage when incorrect marker allele frequencies are assumed in the analyses and selection takes place through trait phenotypes. Misspecification of the assumed marker model also leads to an elevated Type I error rate when linkage analysis is performed on pedigrees ascertained with respect to the trait phenotypes . Asymptotically, bias is NOT introduced in the linkage analysis if there is a misspecification of the trait related parameters and the ascertainment is through the trait phenotypes or if there is a misspecification of the marker related parameters and the ascertainment is through the marker phenotypes .
In previous studies, [7, 8] the power of lod-score linkage analysis using the correct model for the trait appeared to increase when either some or all of the parental marker information was missing and the marker allele frequencies were misspecified in the analysis. An investigation of Type I error rates was then conducted in the absence of parental data and with misspecification of marker allele frequencies. This showed that an inflation in Type I error rate was the cause of at least part of this apparent increased power [8, 9]. The objective of the current study is to investigate, using a different analytic strategy, whether the observed inflation in Type I error rate in model-based lod-score linkage is due to sampling variation.
Simulation of data
Computer simulations generating a quantitative trait and marker data in nuclear families were performed using G.A.S.P. V3.3 . The trait simulated was due to an additive major locus with two equifrequent alleles. The heritability of the trait due to the single trait locus was 90%, with the remaining phenotypic variation (10%) due to a residual environmental component. Data were simulated for highly polymorphic five-allele marker loci linked to the trait locus with recombination fractions of 0.0 and 0.01, as well as an unlinked marker, and with varying combinations of allele frequencies. The simulated marker allele frequencies were 0.5, 0.3, 0.2, 0.1, 0.05, 0.01 and 0.001 for the "first" allele, with the other four alleles being equally frequent. For each marker model, 10,000 samples of 300 nuclear families with sibship size of two (300 independent sibpairs) were simulated. Marker data for all parents were suppressed (missing), in order to provide the situation in which marker allele frequencies are used in lod-score linkage analysis. These samples were then analyzed in two different ways, as described below.
LOD-score linkage analysis assuming the "true" generating model for the trait
Allele frequencies (first, rest) in a five-allele marker used in the study. The frequency of the most common marker allele is underlined. All combinations of true and assumed allele frequencies were tested.
True allele frequencies
Assumed allele frequencies
LOD-score linkage analysis assuming a trait model estimated from the data
The same simulated data were reanalyzed using LODLINK as described above except that instead of assuming the "true" simulation model, parameters estimated by segregation analysis of the trait data were assumed for the trait model. The segregation analysis was performed on each replicate sample using REGCHUNT , an automated quantitative trait segregation analysis program that fits models using the method of maximum likelihood using many different sets of initial estimates. REGCHUNT incorporates the functionality of REGC  running multiple times. This program generates the initial values randomly from a uniform distribution. It discards invalid initial estimates and aborts maximizations in which an allele frequency goes to a bound. For this study, 100 sets of initial estimates in a Mendelian three-distribution model were used for model fitting; and the best fitting model obtained from each sample was used in the subsequent linkage analysis of that sample.
Power and Type I error rates
For both methods of analysis, power was determined as the proportion of 10,000 samples in which analysis of a marker linked to the trait produced maximum LOD scores equal to or exceeding thresholds corresponding to nominal significance levels of 0.01, 0.001 and 0.0001. Type I error rate was determined in a similar manner based on analysis of the unlinked marker, using nominal significance levels of 0.01 and 0.001, since 10,000 replications are not adequate to accurately estimate Type I error at the 0.0001 level. Note that commonly used lod-score thresholds of 1.0, 2.0 and 3.0 correspond to p-values of approximately 0.0159, 0.0012, and 0.0001, respectively [13, 14].
Observed power to detect linkage using the generating trait model
Observed power to detect linkage using the estimated trait model from REGCHUNT
Observed false-positive error rates using the generating trait model
When the most common allele(s) were misspecified as even more common in the analysis, the observed Type I error rate tended to decrease slightly, unless the misspecification was relatively large.
Observed false-positive error rates using the estimated trait model from REGCHUNT
When performing lod-score linkage analysis of a randomly sampled quantitative trait with missing parental marker information, the power of lod-score linkage analysis often appears to increase as the marker allele frequencies are misspecified even when the "correct" generating model for the trait is used in the analysis. The effect of misspecification on power is more pronounced when the nominal significance level is more significant. The apparent increase in power appears to be due at least in part to an increase in Type I error rates. A severe increase in the Type I error rate was observed for lod-score linkage analysis in a sample size of 300 independent sibpairs without parental data when marker allele frequencies were misspecified in the analysis by underestimating the most common marker allele. The inflation may be caused by sampling variation resulting in some extreme samples for which the "true" generating model is not a good fit. It is also possible that this inflation in Type I error rate is observed because these sample sizes were not large enough to display asymptotic properties. Different proportions of available parental information have a significant effect on the power of lod-score linkage analysis depending upon the sample sizes . Omission of all parental marker data has a substantial deleterious effect on the power observed in lod-score linkage analysis for small or moderate sample sizes (100 or 300 independent sib pairs). In the case of larger sample sizes (500 independent sib pairs) this effect is observed only at nominal significance levels more extreme than 0.0001. Use of marker information for at least one parent in the analysis reduces the observed Type I error rate to the nominal level or close to the nominal level in many situations. In addition, use of multiple sibpairs from the same sibship has a similar effect by allowing parental genotypes to be imputed and used in the analysis .
In an attempt to investigate the effect of misspecification of marker allele frequencies on trait heritability, similar analysis was done on a simulated trait with 50% heritability due to an additive major locus with two equifrequent alleles and the other one due completely to a random environmental effect. LODLINK was used to test for linkage in each of 10,000 replicates of 300 families with sibship size two. When the trait was completely environmentally determined, the Type I error rate of the lod-score linkage analysis was found to be robust to misspecification of marker allele frequencies, even if the allele frequencies were severely misspecified (± 0.8 of true frequencies). However, for the trait with heritability of 50%, the Type I error rate was sometimes as much as 21 times the nominal significance level of 0.001 when the allele frequency of a common allele in a five-allele polymorphic marker was misspecified in the analysis as rare (data not shown). In general, the Type I error rates for the trait with a heritability of 90% were the highest among the three traits studied. Thus, when the parental genotype data are missing, the model-based lod-score test of linkage is robust with respect to Type I error for environmentally determined traits but not for the traits which are largely genetically determined. The results for the trait with heritability of 50% were similar to that of the trait heritability of 90% but less dramatic .
In an attempt to resolve the problem of inflated Type I error rates, the model was estimated from each sample using REGCHUNT followed by linkage analysis with LODLINK. The data obtained show that the observed power and Type I error rates of model-based LOD score linkage are often inflated when marker allele frequencies are misspecified, even when using this analytic strategy. However, the Type I error rates observed in this case are usually closer to the nominal levels than those obtained when assuming the generating trait model. This suggests that the observed inflation of Type I error upon misspecification of marker allele frequencies is at least partially due to sampling variation. However, sampling variation is inadequate to completely explain the pattern of inflation of Type I error rates due to marker allele misspecification. The observed Type I error rates are possibly due to not having asymptotic sample sizes. The effect of sibship size is expected to be small relative to the effect of locus heterozygosity on the ability to infer missing parental genotypes from sibship data. If there are two alleles in a marker, no sibship size can determine with certainty the parental mating type. On the other hand, parental mating type can often be determined with certainty (the higher the locus heterozygosity, the more frequently this is possible) in a sibship of size two if there are four or more alleles. From our findings, we have observed that, for two-allele markers, sib-pair sample sizes of 100 and 300 with sibships of size five (i.e., number of families is 10 and 30, respectively) always produce a better power than sib-pair sizes of 100 and 300 with sibships of size two (i.e., number of families is 100 and 300, respectively) when "correct" allele frequencies are used in the analysis. However, as we increased the number of marker alleles to five, the power using 300 sib-pairs with sibships of size five (i.e., the number of families is 30) is not always better than the power achieved using 300 sib-pairs with a sibship size of two (i.e., number of families is 300). However, when only 100 sibpairs are available for study, sibships of size five (i.e., number of families is 10) yield better power than sibships of size two (i.e., number of families is 100) even with a five-allele marker, since power is not optimal for a sample size of only 100 sib-pairs, so the extra power that is obtained by being able to use the multiple sibs in a large sibship to impute the missing parental genotype is important. In general, an increase in the marker heterozygosity leads to more complete inference of parental genotypes . The power with 300 sibpairs is found to be quite high even when we do not have much information to impute the parental genotypes. However, in the case of a disease for which genetic heterogeneity exists, the sample sizes of 300 pairs with 30 families might be more powerful than 300 pairs in 300 families, since power will also not be optimal in this situation. This issue should be explored in future studies.
It is a common practice in simulation studies to use the generating model as the "best" model. The data from the present study suggest that the estimated model generally performs better, with respect to Type I error, than the generating model. Therefore, we need to use caution in interpreting simulation study results obtained using the generating model. These results also suggest that in order to avoid spurious inferences of linkage when performing linkage analysis of randomly ascertained quantitative traits with missing parental data, a very large sample size is needed. Moreover, both the trait model parameters and marker allele frequencies should be estimated from the sample data for randomly ascertained quantitative traits to attempt to reduce the false positive rates.
A thorough investigation was previously conducted on power and Type I error using the model-free Haseman-Elston (H-E) sibpair linkage method [15, 16]. Power results were reported on three types of traits with varying portion of random environmental effect (10%, 30% and 50%). The H-E sib-pair linkage method was found to be robust in situations considered with misspecifications of allele frequencies, except for a slight decrease in power when sample size was small and when the marker was not very polymorphic.
The estimates of marker allele frequencies are irrelevant to linkage analysis when marker genotypes are known for all family members. However, if parental genotypes are not available, accurate estimates of allele frequencies for highly polymorphic DNA marker loci with high heterozygosities are often not available in the underlying study population. If the genotypes at a marker locus are not known and cannot be imputed for a family member, then inaccurate estimates of marker allele frequencies may lead to false positive evidence for linkage if model-based linkage analysis is used. However, if a quantitative trait is well characterized and marker allele frequencies are accurately estimated the type I error rate for the model-based linkage method is not substantially greater than the nominal rate.
In conclusion, the model-based lod-score method of linkage analysis is not robust to misspecification of marker allele frequencies when the families are randomly ascertained and parental marker data are missing, especially when the true trait model is assumed, rather than one estimated from the data. There is some indication, however, that a slight overestimate of the most common marker allele would be safer to use than an underestimate.
The results of this paper were obtained by using the program package S.A.G.E., which is supported by a U.S. Public Health Service Resource Grant (RR03655) from the National Center for Research Resources. The research was supported in part by NIH Contract N01-HG-65404 and by the Intramural Research Program of the National Human Genome Research Institute, National Institutes of Health.
- Clerget-Darpoux F, Bonaiti-Pellie C, Hochez J: Effects of misspecifying genetic parameters in lod score analysis. Biometrics. 1986, 42: 393-399.PubMedView ArticleGoogle Scholar
- Goldin LR, Weeks DE: Two-locus models of disease: comparison of likelihood and nonparametric linkage methods. Am J Hum Genet. 1993, 53: 908-915.PubMedPubMed CentralGoogle Scholar
- Ott J: Strategies for characterizing highly polymorphic markers in human gene mapping. Am J Hum Genet. 1992, 51: 283-290.PubMedPubMed CentralGoogle Scholar
- Freimer NB, Sandkuijl LA, Blower SM: Incorrect specification of marker allele frequencies: effects on linkage analysis. Am J Hum Genet. 1993, 52: 1102-1110.PubMedPubMed CentralGoogle Scholar
- Knapp M, Seuchter SA, Baur MP: The effect of misspecifying allele frequencies in incompletely typed families. Genet Epidemiol. 1993, 10: 413-418. 10.1002/gepi.1370100614.PubMedView ArticleGoogle Scholar
- Amos CI, Williamson JA: Robustness of the maximum-likelihood (LOD) method for detecting linkage. Am J Hum Genet. 1993, 52: 213-214.PubMedPubMed CentralGoogle Scholar
- Mandal DM: A comparison of the power of model dependent and model independent linkage analysis using computer simulation. PhD Dissertation. 1996, Louisiana State University Health Sciences CenterGoogle Scholar
- Mandal DM, Wilson AF, Keats BJB, Bailey-Wilson JE: Factors affecting inflation of Type I error of model-based linkage under random ascertainment. Am J Hum Genet. 1998, 63 (A298):
- Mandal DM, Wilson AF, Sorant AJM, Keats BJB, Bailey-Wilson JE: Comparison of the Type I error rate of model-dependent linkage analysis for quantitative traits under random ascertainment. Am J Hum Genet. 1999, 65 (A389):
- Wilson AF, Bailey-Wilson JE, Pugh EW, Sorant AJM: The Genometric Analysis Simulation Program (G.A.S.P.) A software tool for testing and investigating methods in statistical genetics. Am J Hum Genet. 1996, 59 (A193):
- S.A.G.E. (Release 2.2). S.A.G.E. Statistical Analysis for Genetic Epidemiology, Release 2.2. LSU Medical Center, New Orleans. 1994
- Atwood LD, Kammerer CM, Mitchell BD: Exploring the HDL likelihood surface. Genet Epidemiol. 1993, 10: 641-645. 10.1002/gepi.1370100652.PubMedView ArticleGoogle Scholar
- Morton NE: Sequential tests for the detection of linkage. Am J Hum Genet. 1955, 7: 277-318.PubMedPubMed CentralGoogle Scholar
- Ott J: Methods of analysis and resources available for genetic trait mapping. J Hered. 1999, 90: 68-70. 10.1093/jhered/90.1.68.PubMedView ArticleGoogle Scholar
- Mandal DM, Wilson AF, Elston RC, Weissbecker K, Keats BJ, Bailey-Wilson JE: Effects of misspecification of allele frequencies on the type I error rate of model-free linkage analysis. Hum Hered. 2000, 50: 126-132. 10.1159/000022900.PubMedView ArticleGoogle Scholar
- Mandal DM, Wilson AF, Bailey-Wilson JE: Effects of misspecification of allele frequencies on the power of Haseman-Elston sib-pair linkage method for quantitative traits. Am J Med Genet. 2001, 103: 308-313. 10.1002/ajmg.1566.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.