Methodology and Accuracy of Estimation of Quantitative Trait Loci Parameters in a Half-Sib Design Using Maximum Likelihood

Maximum likelihood methods were developed for estimation of the six parameters relating to a marker-linked quantitative trait locus (QTL) segregating in a half-sib design, namely the QTL additive effect, the QTL dominance effect, the population mean, recombination between the marker and the QTL, the population frequency of the QTL alleles, and the within-family residual variance. The method was tested on simulated stochastic data with various family structures under two genetic models. A method for predicting the expected value of the likelihood was also derived and used to predict the lower bound sampling errors of the parameter estimates and the correlations between them. It was found that standard errors and confidence intervals were smallest for the population mean and variance, intermediate for QTL effects and allele frequency, and highest for recombination rate. Correlations among standard errors of the parameter estimates were generally low except for a strong negative correlation (r = -0.9) between the QTL's dominance effect and the population mean, and medium positive and negative correlations between the QTL's additive effect and, respectively, recombination rate (r = 0.5) and residual variance (r = -0.6). The implications for experimental design and method of analysis on power and accuracy of marker-QTL linkage experiments were discussed.     

Efficiency of triple test cross for detecting epistasis with marker information

The triple test cross (TTC) is an experimental design for detecting epistasis and estimating the components of genetic variance for quantitative traits. In this paper, we extend the analysis to include molecular information. The statistical power of the mating design was assessed under a model assuming that a finite number of loci affect the trait in question. Formulae are developed for the analysis with or without marker information relating to the recombination fraction between loci, the genetical properties of quantitative trait controlled by the quantitative trait loci (QTL), the linkage phases of the parents and population size. Application of these formulae showed that the recombination fraction between genes and the magnitude and the types of epistasis have important interactions in their effects on power. The results demonstrate that the TTC may have increased power to detect epistasis when marker information is present. However, the simulation experiments show that the standard deviation of the estimated expected mean square was higher with one marker than that with two, whereas the corresponding value without marker information was the lowest. In addition, we demonstrate that the relative position of QTL and markers and the number of markers can both affect the power of epistatic detection.     

Using molecular markers to map multiple quantitative trait loci: models for backcross,recombinant inbred,and doubled haploid progeny

Summary To maximize parameter estimation efficiency and statistical power and to estimate epistasis, the parameters of multiple quantitative trait loci (QTLs) must be simultaneously estimated. If multiple QTL affect a trait, then estimates of means of QTL genotypes from individual locus models are statistically biased. In this paper, I describe methods for estimating means of QTL genotypes and recombination frequencies between marker and quantitative trait loci using multilocus backcross, doubled haploid, recombinant inbred, and testcross progeny models. Expected values of marker genotype means were defined using no double or multiple crossover frequencies and flanking markers for linked and unlinked quantitative trait loci. The expected values for a particular model comprise a system of nonlinear equations that can be solved using an interative algorithm, e.g., the Gauss-Newton algorithm. The solutions are maximum likelihood estimates when the errors are normally distributed. A linear model for estimating the parameters of unlinked quantitative trait loci was found by transforming the nonlinear model. Recombination frequency estimators were defined using this linear model. Certain means of linked QTLs are less efficiently estimated than means of unlinked QTLs.     

The Effects of Selection on Linkage Analysis for Quantitative Traits

The effects of within-sample selection on the outcome of analyses detecting linkage between genetic markers and quantitative traits were studied. It was found that selection by truncation for the trait of interest significantly reduces the differences between marker genotype means thus reducing the power to detect linked quantitative trait loci (QTL). The size of this reduction is a function of proportion selected, the magnitude of the QTL effect, recombination rate between the marker locus and the QTL, and the allele frequency of the QTL. Proportion selected was the most influential of these factors on bias, e.g., for an allele substitution effect of one standard deviation unit, selecting the top 80%, 50% or 20% of the population required 2, 6 or 24 times the number of progeny, respectively, to offset the loss of power caused by this selection. The effect on power was approximately linear with respect to the size of gene effect, almost invariant to recombination rate, and a complex function of QTL allele frequency. It was concluded that experimental samples from animal populations which have been subjected to even minor amounts of selection will be inefficient in yielding information on linkage between markers and loci influencing the quantitative trait under selection.     

Detecting Marker-Qtl Linkage and Estimating Qtl Gene Effect and Map Location Using a Saturated Genetic Map

A simulation study was carried out on a backcross population in order to determine the effect of marker spacing, gene effect and population size on the power of marker-quantitative trait loci (QTL) linkage experiments and on the standard error of maximum likelihood estimates (MLE) of QTL gene effect and map location. Power of detecting a QTL was virtually the same for a marker spacing of 10 cM as for an infinite number of markers and was only slightly decreased for marker spacing of 20 or even 50 cM. The advantage of using interval mapping as compared to single-marker analysis was slight. ``Resolving power' of a marker-QTL linkage experiment was defined as the 95% confidence interval for the QTL map location that would be obtained when scoring an infinite number of markers. It was found that reducing marker spacing below the resolving power did not add appreciably to narrowing the confidence interval. Thus, the 95% confidence interval with infinite markers sets the useful marker spacing for estimating QTL map location for a given population size and estimated gene effect.     

A comparison of regression interval mapping and multiple interval mapping for linked QTL

Regression interval mapping and multiple interval mapping are compared with regard to mapping linked quantitative trait loci (QTL) in inbred-line cross experiments. For that purpose, a simulation study was performed using genetic models with two linked QTL. Data were simulated for F(2) populations of different sizes and with all QTL and marker alleles fixed for alternative alleles in the parental lines. The criteria for comparison are power of QTL identification and the accuracy of the QTL position and effect estimates. Further, the estimates of the relative QTL variance are assessed. There are distinct differences in the QTL position estimates between the two methods. Multiple interval mapping tends to be more powerful as compared to regression interval mapping. Multiple interval mapping further leads to more accurate QTL position and QTL effect estimates. The superiority increased with wider marker intervals and larger population sizes. If QTL are in repulsion, the differences between the two methods are very pronounced. For both methods, the reduction of the marker interval size from 10 to 5 cM increases power and greatly improves QTL parameter estimates. This contrasts with findings in the literature for single QTL scenarios, where a marker density of 10 cM is generally considered as sufficient. The use of standard (asymptotic) statistical theory for the computation of the standard errors of the QTL position and effect estimates proves to give much too optimistic standard errors for regression interval mapping as well as for multiple interval mapping.     

QTL analysis: unreliability and bias in estimation procedures

Several statistical methods which employ multiple marker data are currently available for the analysis of quantitative trait loci (QTL) in experimental populations. Although comparable estimates of QTL location and effects have been obtained by these methods, using simulated and real data sets, their accuracy and reliability have not been extensively investigated. The present study specifically examines the merit of using F₂ and doubled haploid populations for locating QTL and estimating their effects. Factors which may affect accuracy and reliability of QTL mapping, such as the number and position of the markers available, the accuracy of the marker locations and the size of the experimental population used, are considered. These aspects are evaluated for QTL of differing heritabilities and locations along the chromosome.A population of 300 F₂ individuals and 150 doubled haploid lines gave estimates of QTL position and effect which were comparable, albeit extremely unreliable. Even for a QTL of high heritability (10%), the confidence interval was 35 cM. There was little increase in reliability to be obtained from using 300, rather than 200, F₂ individuals and 100 doubled haploid lines gave similar results to 150. QTL estimates were not significantly improved either by using the expected, rather than the observed, marker positions or by using a dense map of markers rather than a sparse map. A QTL which was asymmetrically located in the linkage group resulted in inaccurate estimates of QTL position which were seriously biassed at low heritability of the QTL. In a population of 300 F₂ individuals the bias increased from 4 cM to 20 cM, for a QTL with 10% and 2% heritability respectively.     

Conclusion of LOD-score analysis for family data generated under two-locus models.

The power to detect linkage by the LOD-score method is investigated here for diseases that depend on the effects of two genes. The classical strategy is, first, to detect a major-gene (MG) effect by segregation analysis and, second, to seek for linkage with genetic markers by the LOD-score method using the MG parameters. We already showed that segregation analysis can lead to evidence for a MG effect for many two-locus models, with the estimates of the MG parameters being very different from those of the two genes involved in the disease. We show here that use of these MG parameter estimates in the LOD-score analysis may lead to a failure to detect linkage for some two-locus models. For these models, use of the sib-pair method gives a non-negligible increase of power to detect linkage. The linkage-homogeneity test among subsamples differing for the familial disease distribution provides evidence of parameter misspecification, when the MG parameters are used. Moreover, for most of the models, use of the MG parameters in LOD-score analysis leads to a large bias in estimation of the recombination fraction and sometimes also to a rejection of linkage for the true recombination fraction. A final important point is that a strong evidence of an MG effect, obtained by segregation analysis, does not necessarily imply that linkage will be detected for at least one of the two genes, even with the true parameters and with a close informative marker.     

A correlation method for detecting and estimating linkage between a marker locus and a quantitative trait locus using inbred lines

The advent of molecular genetic markers has stimulated interest in detecting linkage between a marker locus and a quantitative trait locus (QTL) because the marker locus, even without direct effect on the quantitative trait, could be useful in increasing the response to selection. A correlation method for detecting and estimating linkage between a marker locus and a QTL is described using selfing and sib-mating populations. Computer simulations were performed to estimate the power of the method, the sample size (N) needed to detect linkage, and the recombination value (r). The power of this method was a function of the expected recombination value E(r), the standardized difference (d) between the QTL genotypic means, and N. The power was highest at complete linkage, decreased with an increase in E(r), and then increased at E(r)=0.5. A larger d and N led to a higher power. The sample size needed to detect linkage was dependent upon E(r) and d. The sample size had a minimum value at E(r)=0, increased with an increase in E(r) and a decrease in d. In general, the r was overestimated. With an increase in d, the r was closer to its expectation. Detection of linkage by the proposed method under incomplete linkage was more efficient than estimation of recombination values. The correlation method and the method of comparison of marker-genotype means have a similar power when there is linkage, but the former has a slightly higher power than the latter when there is no linkage.     

Linkage disequilibrium fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

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.     

Bayesian analysis of linkage between genetic markers and quantitative trait loci. II. Combining prior knowledge with experimental evidence

Summary A Bayesian method was developed for identifying genetic markers linked to quantitative trait loci (QTL) by analyzing data from daughter or granddaughter designs and single markers or marker pairs. Traditional methods may yield unrealistic results because linkage tests depend on number of markers and QTL gene effects associated with selected markers are overestimated. The Bayesian or posterior probability of linkage combines information from a daughter or granddaughter design with the prior probability of linkage between a marker locus and a QTL. If the posterior probability exceeds a certain quantity, linkage is declared. Upon linkage acceptance, Bayesian estimates of marker-QTL recombination rate and QTL gene effects and frequencies are obtained. The Bayesian estimates of QTL gene effects account for different amounts of information by shrinking information from data toward the mean or mode of a prior exponential distribution of gene effects. Computation of the Bayesian analysis is feasible. Exact results are given for biallelic QTL, and extensions to multiallelic QTL are suggested.     

Optimal use of regression models in genome-wide association studies

The performance of linear regression models in genome-wide association studies is influenced by how marker information is parameterized in the model. Considering the impact of parameterization is especially important when using information from multiple markers to test for association. Properties of the population, such as linkage disequilibrium (LD) and allele frequencies, will also affect the ability of a model to provide statistical support for an underlying quantitative trait locus (QTL). Thus, for a given location in the genome, the relationship between population properties and model parameterization is expected to influence the performance of the model in providing evidence for the position of a QTL. As LD and allele frequencies vary throughout the genome and between populations, understanding the relationship between these properties and model parameterization is of considerable importance in order to make optimal use of available genomic data. Here, we evaluate the performance of regression-based association models using genotype and haplotype information across the full spectrum of allele frequency and LD scenarios. Genetic marker data from 200 broiler chickens were used to simulate genomic conditions by selecting individual markers to act as surrogate QTL (sQTL) and then investigating the ability of surrounding markers to estimate sQTL genotypes and provide statistical support for their location. The LD and allele frequencies of markers and sQTL are shown to have a strong effect on the performance of models relative to one another. Our results provide an indication of the best choice of model parameterization given certain scenarios of marker and QTL LD and allele frequencies. We demonstrate a clear advantage of haplotype-based models, which account for phase uncertainty over other models tested, particularly for QTL with low minor allele frequencies. We show that the greatest advantage of haplotype models over single-marker models occurs when LD between markers and the causal locus is low. Under these situations, haplotype models have a greater accuracy of predicting the location of the QTL than other models tested.     

A simulation study on the accuracy of position and effect estimates of linked QTL and their asymptotic standard deviations using multiple interval mapping in an F2 scheme

Approaches like multiple interval mapping using a multiple-QTL model for simultaneously mapping QTL can aid the identification of multiple QTL, improve the precision of estimating QTL positions and effects, and are able to identify patterns and individual elements of QTL epistasis. Because of the statistical problems in analytically deriving the standard errors and the distributional form of the estimates and because the use of resampling techniques is not feasible for several linked QTL, there is the need to perform large-scale simulation studies in order to evaluate the accuracy of multiple interval mapping for linked QTL and to assess confidence intervals based on the standard statistical theory. From our simulation study it can be concluded that in comparison with a monogenetic background a reliable and accurate estimation of QTL positions and QTL effects of multiple QTL in a linkage group requires much more information from the data. The reduction of the marker interval size from 10 cM to 5 cM led to a higher power in QTL detection and to a remarkable improvement of the QTL position as well as the QTL effect estimates. This is different from the findings for (single) interval mapping. The empirical standard deviations of the genetic effect estimates were generally large and they were the largest for the epistatic effects. These of the dominance effects were larger than those of the additive effects. The asymptotic standard deviation of the position estimates was not a good criterion for the accuracy of the position estimates and confidence intervals based on the standard statistical theory had a clearly smaller empirical coverage probability as compared to the nominal probability. Furthermore the asymptotic standard deviation of the additive, dominance and epistatic effects did not reflect the empirical standard deviations of the estimates very well, when the relative QTL variance was smaller/equal to 0.5. The implications of the above findings are discussed.     

The effects of selective genotyping on estimates of proportion of recombination between linked quantitative trait loci

Selective genotyping is the marker assay of only the more extreme phenotypes for a quantitative trait and is intended to increase the efficiency of quantitative trait loci (QTL) mapping. We show that selective genotyping can bias estimates of the recombination frequency between linked QTLs — upwardly when QTLs are in repulsion phase, and downwardly when QTLs are in coupling phase. We examined these biases under simple models involving two QTLs segregating in a backcross or F₂ population, using both analytical models and computer simulations. We found that bias is a function of the proportion selected, the magnitude of QTL effects, distance between QTLs and the dominance of QTLs. Selective genotyping thus may decrease the power of mapping multiple linked QTLs and bias the construction of a marker map. We suggest a large proportion than previously suggested (50%) or the entire population be genotyped if linked QTLs of large effects (explain > 10% phenotypic variance) are evident. New models need to be developed to explicitly incorporate selection into QTL map construction.     

Interval mapping of quantitative trait loci with selective DNA pooling data

Selective DNA pooling is an efficient method to identify chromosomal regions that harbor quantitative trait loci (QTL) by comparing marker allele frequencies in pooled DNA from phenotypically extreme individuals. Currently used single marker analysis methods can detect linkage of markers to a QTL but do not provide separate estimates of QTL position and effect, nor do they utilize the joint information from multiple markers. In this study, two interval mapping methods for analysis of selective DNA pooling data were developed and evaluated. One was based on least squares regression (LS-pool) and the other on approximate maximum likelihood (ML-pool). Both methods simultaneously utilize information from multiple markers and multiple families and can be applied to different family structures (half-sib, F2 cross and backcross). The results from these two interval mapping methods were compared with results from single marker analysis by simulation. The results indicate that both LS-pool and ML-pool provided greater power to detect the QTL than single marker analysis. They also provide separate estimates of QTL location and effect. With large family sizes, both LS-pool and ML-pool provided similar power and estimates of QTL location and effect as selective genotyping. With small family sizes, however, the LS-pool method resulted in severely biased estimates of QTL location for distal QTL but this bias was reduced with the ML-pool.     

Combined linkage and association mapping of quantitative trait loci by multiple markers

Using multiple diallelic markers, variance component models are proposed for high-resolution combined linkage and association mapping of quantitative trait loci (QTL) based on nuclear families. The objective is to build a model that may fully use marker information for fine association mapping of QTL in the presence of prior linkage. The measures of linkage disequilibrium and the genetic effects are incorporated in the mean coefficients and are decomposed into orthogonal additive and dominance effects. The linkage information is modeled in variance-covariance matrices. Hence, the proposed methods model both association and linkage in a unified model. On the basis of marker information, a multipoint interval mapping method is provided to estimate the proportion of allele sharing identical by descent (IBD) and the probability of sharing two alleles IBD at a putative QTL for a sib-pair. To test the association between the trait locus and the markers, both likelihood-ratio tests and F-tests can be constructed on the basis of the proposed models. In addition, analytical formulas of noncentrality parameter approximations of the F-test statistics are provided. Type I error rates of the proposed test statistics are calculated to show their robustness. After comparing with the association between-family and association within-family (AbAw) approach by Abecasis and Fulker et al., it is found that the method proposed in this article is more powerful and advantageous based on simulation study and power calculation. By power and sample size comparison, it is shown that models that use more markers may have higher power than models that use fewer markers. The multiple-marker analysis can be more advantageous and has higher power in fine mapping QTL. As an application, the Genetic Analysis Workshop 12 German asthma data are analyzed using the proposed methods.     

Mapping quantitative trait loci by genotyping haploid tissues.

Mapping strategies based on a half- or full-sib family design have been developed to map quantitative trait loci (QTL) for outcrossing species. However, these strategies are dependent on controlled crosses where marker-allelic frequency and linkage disequilibrium between the marker and QTL may limit their application. In this article, a maximum-likelihood method is developed to map QTL segregating in an open-pollinated progeny population using dominant markers derived from haploid tissues from single meiotic events. Results from the haploid-based mapping strategy are not influenced by the allelic frequencies of markers and their linkage disequilibria with QTL, because the probabilities of QTL genotypes conditional on marker genotypes of haploid tissues are independent of these population parameters. Parameter estimation and hypothesis testing are implemented via expectation/conditional maximization algorithm. Parameters estimated include the additive effect, the dominant effect, the population mean, the chromosomal location of the QTL in the interval, and the residual variance within the QTL genotypes, plus two population parameters, outcrossing rate and QTL-allelic frequency. Simulation experiments show that the accuracy and power of parameter estimates are affected by the magnitude of QTL effects, heritability levels of a trait, and sample sizes used. The application and limitation of the method are discussed.     

Multipoint identity-by-descent prediction using dense markers to map quantitative trait loci and estimate effective population size

A novel multipoint method, based on an approximate coalescence approach, to analyze multiple linked markers is presented. Unlike other approximate coalescence methods, it considers all markers simultaneously but only two haplotypes at a time. We demonstrate the use of this method for linkage disequilibrium (LD) mapping of QTL and estimation of effective population size. The method estimates identity-by-descent (IBD) probabilities between pairs of marker haplotypes. Both LD and combined linkage and LD mapping rely on such IBD probabilities. The method is approximate in that it considers only the information on a pair of haplotypes, whereas a full modeling of the coalescence process would simultaneously consider all haplotypes. However, full coalescence modeling is computationally feasible only for few linked markers. Using simulations of the coalescence process, the method is shown to give almost unbiased estimates of the effective population size. Compared to direct marker and haplotype association analyses, IBD-based QTL mapping showed clearly a higher power to detect a QTL and a more realistic confidence interval for its position. The modeling of LD could be extended to estimate other LD-related parameters such as recombination rates.     

Detection of quantitative trait Loci influencing recombination using recombinant inbred lines

The genetic basis of variation in recombination in higher plants is polygenic and poorly understood, despite its theoretical and practical importance. Here a method of detecting quantitative trait loci (QTL) influencing recombination in recombinant inbred lines (RILs) is proposed that relies upon the fact that genotype data within RILs carry the signature of past recombination. Behavior of the segregational genetic variance in numbers of chromosomal crossovers (recombination) over generations is described for self-, full-sib-, and half-sib-generated RILs with no dominance in true crossovers. This genetic variance, which as a fraction of the total phenotypic variance contributes to the statistical power of the method, was asymptotically greatest with half sibbing, less with sibbing, and least with selfing. The statistical power to detect a recombination QTL declined with diminishing QTL effect, genome target size, and marker density. For reasonably tight marker linkage power was greater with less intense inbreeding for later generations and vice versa for early generations. Generational optima for segregation variance and statistical power were found, whose onset and narrowness varied with marker density and mating design, being more pronounced for looser marker linkage. Application of this method to a maize RIL population derived from inbred lines Mo17 and B73 and developed by selfing suggested two putative QTL (LOD > 2.4) affecting certain chromosomes, and using a canonical transformation another putative QTL was detected. However, permutation tests failed to support their presence (experimentwise alpha = 0.05). Other populations with more statistical power and chosen specifically for recombination QTL segregation would be more effective.