Bayesian robust analysis for genetic architecture of quantitative traits

MOTIVATION: In most quantitative trait locus (QTL) mapping studies, phenotypes are assumed to follow normal distributions. Deviations from this assumption may affect the accuracy of QTL detection and lead to detection of spurious QTLs. To improve the robustness of QTL mapping methods, we replaced the normal distribution for residuals in multiple interacting QTL models with the normal/independent distributions that are a class of symmetric and long-tailed distributions and are able to accommodate residual outliers. Subsequently, we developed a Bayesian robust analysis strategy for dissecting genetic architecture of quantitative traits and for mapping genome-wide interacting QTLs in line crosses. RESULTS: Through computer simulations, we showed that our strategy had a similar power for QTL detection compared with traditional methods assuming normal-distributed traits, but had a substantially increased power for non-normal phenotypes. When this strategy was applied to a group of traits associated with physical/chemical characteristics and quality in rice, more main and epistatic QTLs were detected than traditional Bayesian model analyses under the normal assumption.     

Mapping multiple QTL of different effects: comparison of a simple sequential testing strategy and multiple QTL mapping

The aim of this study was to explore, by computer simulation, the mapping of QTLs in a realistic but complex situation of many (linked) QTLs with different effects, and to compare two QTL mapping methods. A novel method to dissect genetic variation on multiple chromosomes using molecular markers in backcross and F₂ populations derived from inbred lines was suggested, and its properties tested using simulations. The rationale for this sequential testing method was to explicitly test for alternative genetic models. The method consists of a series of four basic statistical tests to decide whether variance was due to a single QTL, two QTLs, multiple QTLs, or polygenes, starting with a test to detect genetic variance associated with a particular chromosome. The method was able to distinguish between different QTL configurations, in that the probability to `detect' the correct model was high, varying from 0.75 to 1. For example, for a backcross population of 200 and an overall heritability of 50%, in 78% of replicates a polygenic model was detected when that was the underlying true model. To test the method for multiple chromosomes, QTLs were simulated on 10 chromosomes, following a geometric series of allele effects, assuming positive alleles were in coupling in the founder lines For these simulations, the sequential testing method was compared to the established Multiple QTL Mapping (MQM) method. For a backcross population of 400 individuals, power to detect genetic variance was low with both methods when the heritability was 0.40. For example, the power to detect genetic variation on a chromosome on which 6 QTLs explained 12.6% of the genetic variance, was less than 60% for both methods. For a large heritability (0.90), the power of MQM to detect genetic variance and to dissect QTL configurations was generally better, due to the simultaneous fitting of markers on all chromosomes. It is concluded that when testing different QTL configurations on a single chromosome using the sequential testing procedure, regions of other chromosomes which explain a significant amount of variation should be fitted in the model of analysis. This study reinforces the need for large experiments in plants and other species if the aim of a genome scan is to dissect quantitative genetic variation.     

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.     

Limitations of a two-step moment method for mapping linked multiple QTL

As an alternative to multiple-interval mapping a two-step moment method was recently proposed to map linked multiple quantitative trait loci (QTLs). The advantage of this moment method was supposed to be its simplicity and computational efficiency, especially in detecting closely linked QTLs within a marker bracket, but also in mapping QTLs in different marker intervals. Using simulations it is shown that the two-step moment method may give poor results compared with multiple-interval mapping, irrespective of whether the QTLs are in the same or in different marker intervals, especially if linked QTLs are in repulsion. The criteria of comparison are number of identified QTLs, likelihood ratio test statistics, means and empirical standard errors of the QTL position and QTL effects estimates, and the accuracy of the residual variance estimates. Further, the joint conditional probabilities of QTL genotypes for two putative QTLs within a marker interval were derived and compared with the unmodified approach ignoring the non-independence of the conditional probabilities.     

Mapping quantitative trait loci for binary traits using a heterogeneous residual variance model: an application to Marek's disease susceptibility in chickens

A typical problem in mapping quantitative trait loci (QTLs) comes from missing QTL genotype. A routine method for parameter estimation involving missing data is the mixture model maximum likelihood method. We developed an alternative QTL mapping method that describes a mixture of several distributions by a single model with a heterogeneous residual variance. The two methods produce similar results, but the heterogeneous residual variance method is computationally much faster than the mixture model approach. In addition, the new method can automatically generate sampling variances of the estimated parameters. We derive the new method in the context of QTL mapping for binary traits in a F2 population. Using the heterogeneous residual variance model, we identified a QTL on chromosome IV that controls Marek's disease susceptibility in chickens. The QTL alone explains 7.2% of the total disease variation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Use of randomization testing to detect multiple epistatic QTLs

Here, we describe a randomization testing strategy for mapping interacting quantitative trait loci (QTLs). In a forward selection strategy, non-interacting QTLs and simultaneously mapped interacting QTL pairs are added to a total genetic model. Simultaneous mapping of epistatic QTLs increases the power of the mapping strategy by allowing detection of interacting QTL pairs where none of the QTL can be detected by their marginal additive and dominance effects. Randomization testing is used to derive empirical significance thresholds for every model selection step in the procedure. A simulation study was used to evaluate the statistical properties of the proposed randomization tests and for which types of epistasis simultaneous mapping of epistatic QTLs adds power. Least squares regression was used for QTL parameter estimation but any other QTL mapping method can be used. A genetic algorithm was used to search for interacting QTL pairs, which makes the proposed strategy feasible for single processor computers. We believe that this method will facilitate the evaluation of the importance at epistatic interaction among QTLs controlling multifactorial traits and disorders.     

Power analysis for quantitative trait locus mapping in populations derived by multiple backcrosses

 Populations derived by multiple backcrosses are potentially useful for quantitative trait locus (QTL) mapping studies. Comparisons of relative power to detect QTL using populations derived by multiple back-crosses are needed to make decisions when mapping projects are initiated. The objective of this study was to theoretically compare the power to detect QTL in populations derived by multiple backcrosses relative to mapping in a recombinant inbred population of equal size. Backcrossing results in a reduction in genetic variance with each generation and also results in an increasing frequency of the recurrent parent marker genotype. The relevant outcome for QTL mapping is a reduction in genetic variance to partition between marker genotype classes and increasing unbalance of the number of individuals contributing to the mean of the marker genotypes. Both of these factors lead to a decrease in the power to detect a QTL as the number of backcross generations increases. Experimental error was held constant with the populations compared. From a theoretical standpoint, backcross-derived populations offer few advantages for QTL detection. If, however, a backcrossing approach is the most efficient method to achieve a desired breeding objective and if QTL detection is an objective of equal or less importance, backcross-derived populations are a reasonable approach to QTL detection. Received: 4 August 1996 / Accepted: 4 April 1997     

A simple method for calculating the statistical power for detecting a QTL located in a marker interval

We developed a simple method for calculating the statistical power for detecting a QTL located in an interval flanked by two markers. The statistical method for QTL detection is assumed to be the Haley and Knott's simple regression method of interval mapping. This method allows us to answer one of the fundamental questions in designing a QTL mapping experiment: What is the minimum marker density required to detect a QTL explaining a certain heritable proportion of the phenotypic variance (denoted by h(2)) with a power gamma under a Type I error alpha in an F(2) or other mating designs with a sample size n? Computing the statistical power only requires the ability to evaluate a non-central F-distribution function and the inverse function of this distribution.     

A modified algorithm for the improvement of composite interval mapping

Composite interval mapping (CIM) is the most commonly used method for mapping quantitative trait loci (QTL) with populations derived from biparental crosses. However, the algorithm implemented in the popular QTL Cartographer software may not completely ensure all its advantageous properties. In addition, different background marker selection methods may give very different mapping results, and the nature of the preferred method is not clear. A modified algorithm called inclusive composite interval mapping (ICIM) is proposed in this article. In ICIM, marker selection is conducted only once through stepwise regression by considering all marker information simultaneously, and the phenotypic values are then adjusted by all markers retained in the regression equation except the two markers flanking the current mapping interval. The adjusted phenotypic values are finally used in interval mapping (IM). The modified algorithm has a simpler form than that used in CIM, but a faster convergence speed. ICIM retains all advantages of CIM over IM and avoids the possible increase of sampling variance and the complicated background marker selection process in CIM. Extensive simulations using two genomes and various genetic models indicated that ICIM has increased detection power, a reduced false detection rate, and less biased estimates of QTL effects.     

On the differences between maximum likelihood and regression interval mapping in the analysis of quantitative trait loci

The differences between maximum-likelihood (ML) and regression (REG) interval mapping in the analysis of quantitative trait loci (QTL) are investigated analytically and numerically by simulation. The analytical investigation is based on the comparison of the solution sets of the ML and REG methods in the estimation of QTL parameters. Their differences are found to relate to the similarity between the conditional posterior and conditional probabilities of QTL genotypes and depend on several factors, such as the proportion of variance explained by QTL, relative QTL position in an interval, interval size, difference between the sizes of QTL, epistasis, and linkage between QTL. The differences in mean squared error (MSE) of the estimates, likelihood-ratio test (LRT) statistics in testing parameters, and power of QTL detection between the two methods become larger as (1) the proportion of variance explained by QTL becomes higher, (2) the QTL locations are positioned toward the middle of intervals, (3) the QTL are located in wider marker intervals, (4) epistasis between QTL is stronger, (5) the difference between QTL effects becomes larger, and (6) the positions of QTL get closer in QTL mapping. The REG method is biased in the estimation of the proportion of variance explained by QTL, and it may have a serious problem in detecting closely linked QTL when compared to the ML method. In general, the differences between the two methods may be minor, but can be significant when QTL interact or are closely linked. The ML method tends to be more powerful and to give estimates with smaller MSEs and larger LRT statistics. This implies that ML interval mapping can be more accurate, precise, and powerful than REG interval mapping. The REG method is faster in computation, especially when the number of QTL considered in the model is large. Recognizing the factors affecting the differences between REG and ML interval mapping can help an efficient strategy, using both methods in QTL mapping to be outlined.     

Bayesian mapping of quantitative trait loci for multiple complex traits with the use of variance components

Complex traits important for humans are often correlated phenotypically and genetically. Joint mapping of quantitative-trait loci (QTLs) for multiple correlated traits plays an important role in unraveling the genetic architecture of complex traits. Compared with single-trait analysis, joint mapping addresses more questions and has advantages for power of QTL detection and precision of parameter estimation. Some statistical methods have been developed to map QTLs underlying multiple traits, most of which are based on maximum-likelihood methods. We develop here a multivariate version of the Bayes methodology for joint mapping of QTLs, using the Markov chain-Monte Carlo (MCMC) algorithm. We adopt a variance-components method to model complex traits in outbred populations (e.g., humans). The method is robust, can deal with an arbitrary number of alleles with arbitrary patterns of gene actions (such as additive and dominant), and allows for multiple phenotype data of various types in the joint analysis (e.g., multiple continuous traits and mixtures of continuous traits and discrete traits). Under a Bayesian framework, parameters--including the number of QTLs--are estimated on the basis of their marginal posterior samples, which are generated through two samplers, the Gibbs sampler and the reversible-jump MCMC. In addition, we calculate the Bayes factor related to each identified QTL, to test coincident linkage versus pleiotropy. The performance of our method is evaluated in simulations with full-sib families. The results show that our proposed Bayesian joint-mapping method performs well for mapping multiple QTLs in situations of either bivariate continuous traits or mixed data types. Compared with the analysis for each trait separately, Bayesian joint mapping improves statistical power, provides stronger evidence of QTL detection, and increases precision in estimation of parameter and QTL position. We also applied the proposed method to a set of real data and detected a coincident linkage responsible for determining bone mineral density and areal bone size of wrist in humans.     

Regression-based multi-trait QTL mapping using a structural equation model

Quantitative trait loci (QTL) mapping often results in data on a number of traits that have well-established causal relationships. Many multi-trait QTL mapping methods that account for the correlation among multiple traits have been developed to improve the statistical power and the precision of QTL parameter estimation. However, none of these methods are capable of incorporating the causal structure among the traits. Consequently, genetic functions of the QTL may not be fully understood. Structural equation modeling (SEM) allows researchers to explicitly characterize the causal structure among the variables and to decompose effects into direct, indirect, and total effects. In this paper, we developed a multi-trait SEM method of QTL mapping that takes into account the causal relationships among traits related to grain yield. Performance of the proposed method is evaluated by simulation study and applied to data from a wheat experiment. Compared with single trait analysis and the multi-trait least-squares analysis, our multi-trait SEM improves statistical power of QTL detection and provides important insight into how QTLs regulate traits by investigating the direct, indirect, and total QTL effects. The approach also helps build biological models that more realistically reflect the complex relationships among QTL and traits and is more precise and efficient in QTL mapping than single trait analysis.     

Empirical Bayesian elastic net for multiple quantitative trait locus mapping

In multiple quantitative trait locus (QTL) mapping, a high-dimensional sparse regression model is usually employed to account for possible multiple linked QTLs. The QTL model may include closely linked and thus highly correlated genetic markers, especially when high-density marker maps are used in QTL mapping because of the advancement in sequencing technology. Although existing algorithms, such as Lasso, empirical Bayesian Lasso (EBlasso) and elastic net (EN) are available to infer such QTL models, more powerful methods are highly desirable to detect more QTLs in the presence of correlated QTLs. We developed a novel empirical Bayesian EN (EBEN) algorithm for multiple QTL mapping that inherits the efficiency of our previously developed EBlasso algorithm. Simulation results demonstrated that EBEN provided higher power of detection and almost the same false discovery rate compared with EN and EBlasso. Particularly, EBEN can identify correlated QTLs that the other two algorithms may fail to identify. When analyzing a real dataset, EBEN detected more effects than EN and EBlasso. EBEN provides a useful tool for inferring high-dimensional sparse model in multiple QTL mapping and other applications. An R software package ‘EBEN'' implementing the EBEN algorithm is available on the Comprehensive R Archive Network (CRAN).     

Detection of putative quantitative trait loci in line crosses under infinitesimal genetic models

Quantitative trait locus (QTL) mapping studies often employ segregating generations derived from a cross between genetically divergent inbred lines. In the analysis of such data it is customary to fit a single QTL and use a null hypothesis which assumes that the genomic region under study contributes no genetic variance. To explore the situation in which multiple linked genes contribute to the genetic variance, we simulated an F₂-mapping experiment in which the genetic difference between the two original inbred strains was caused by a large number of loci, each having equal effect on the quantitative trait. QTLs were either in coupling, dispersion or repulsion phase in the base population of inbred lines, with the expected F₂ genetic variance explained by the QTLs being equivalent in the three models. Where QTLs were in coupling phase, one inbred line was fixed for all plus alleles, and the other line was fixed for minus alleles. Where QTLs were in dispersion phase, they were assumed to be randomly fixed for one or other allele (as if the inbred lines had evolved from a common ancestor by random drift). Where QTLs were in repulsion phase alleles within an inbred line were alternating plus and minus at adjacent loci, and alternative alleles were fixed in the two inbred lines. In all these genetic models a standard interval mapping test statistic used to determine whether there is a QTL of large effect segregating in the population was inflated on average. Furthermore, the use of a threshold for QTL detection derived under the assumption that no QTLs were segregating would often lead to spurious conclusions regards the presence of genes of large effects (i.e. type I errors). The employment of an alternative model for the analysis, including linked markers as cofactors in the analysis of a single interval, reduced the problem of type I error rate, although test statistics were still inflated relative to the case of no QTLs. It is argued that in practice one should take into account the difference between the strains or the genetic variance in the F₂ population when setting significance thresholds. In addition, tests designed to probe the adequacy of a single-QTL model or of an alternative infinitesimal coupling model are described. Such tests should be applied in QTL mapping studies to help dissect the true nature of genetic variation.     

Simultaneous mapping of epistatic QTL in chickens reveals clusters of QTL pairs with similar genetic effects on growth

We used simultaneous mapping of interacting quantitative trait locus (QTL) pairs to study various growth traits in a chicken F2 intercross. The method was shown to increase the number of detected QTLs by 30 % compared with a traditional method detecting QTLs by their marginal genetic effects. Epistasis was shown to be an important contributor to the genetic variance of growth, with the largest impact on early growth (before 6 weeks of age). There is also evidence for a discrete set of interacting loci involved in early growth, supporting the previous findings of different genetic regulation of early and late growth in chicken. The genotype-phenotype relationship was evaluated for all interacting QTL pairs and 17 of the 21 evaluated QTL pairs could be assigned to one of four clusters in which the pairs in a cluster have very similar genetic effects on growth. The genetic effects of the pairs indicate commonly occurring dominance-by-dominance, heterosis and multiplicative interactions. The results from this study clearly illustrate the increase in power obtained by using this novel method for simultaneous detection of epistatic QTL, and also how visualization of genotype-phenotype relationships for epistatic QTL pairs provides new insights to biological mechanisms underlying complex traits.     

Bayesian functional mapping of dynamic quantitative traits

Without consideration of other linked QTLs responsible for dynamic trait, original functional mapping based on a single QTL model is not optimal for analyzing multiple dynamic trait loci. Despite that composite functional mapping incorporates the effects of genetic background outside the tested QTL in mapping model, the arbitrary choice of background markers also impact on the power of QTL detection. In this study, we proposed Bayesian functional mapping strategy that can simultaneously identify multiple QTL controlling developmental patterns of dynamic traits over the genome. Our proposed method fits the change of each QTL effect with the time by Legendre polynomial and takes the residual covariance structure into account using the first autoregressive equation. Also, Bayesian shrinkage estimation was employed to estimate the model parameters. Especially, we specify the gamma distribution as the prior for the first-order auto-regressive coefficient, which will guarantee the convergence of Bayesian sampling. Simulations showed that the proposed method could accurately estimate the QTL parameters and had a greater statistical power of QTL detection than the composite functional mapping. A real data analysis of leaf age growth in rice is used for the demonstration of our method. It shows that our Bayesian functional mapping can detect more QTLs as compared to composite functional mapping.     

Mapping quantitative trait loci with dominant markers in four-way crosses

 A common problem in mapping quantitative trait loci (QTLs) is that marker data are often incomplete. This includes missing data, dominant markers, and partially informative markers, arising in outbred populations. Here we briefly present an iteratively re-weighted least square method (IRWLS) to incorporate dominant and missing markers for mapping QTLs in four-way crosses under a heterogeneous variance model. The algorithm uses information from all markers in a linkage group to infer the QTL genotype. Monte Carlo simulations indicate that with half dominant markers, QTL detection is almost as efficient as with all co-dominant markers. However, the precision of the estimated QTL parameters generally decreases as more markers become missing or dominant. Notable differences are observed on the standard deviation of the estimated QTL position for varying levels of marker information content. The method is relatively simple so that more complex models including multiple QTLs or fixed effects can be fitted. Finally, the method can be readily extended to QTL mapping in full-sib families. Received: 16 June 1998 / Accepted: 29 September 1998     

Generalized Linear Model for Mapping Discrete Trait Loci Implemented with LASSO Algorithm

Generalized estimating equation (GEE) algorithm under a heterogeneous residual variance model is an extension of the iteratively reweighted least squares (IRLS) method for continuous traits to discrete traits. In contrast to mixture model-based expectation–maximization (EM) algorithm, the GEE algorithm can well detect quantitative trait locus (QTL), especially large effect QTLs located in large marker intervals in the manner of high computing speed. Based on a single QTL model, however, the GEE algorithm has very limited statistical power to detect multiple QTLs because of ignoring other linked QTLs. In this study, the fast least absolute shrinkage and selection operator (LASSO) is derived for generalized linear model (GLM) with all possible link functions. Under a heterogeneous residual variance model, the LASSO for GLM is used to iteratively estimate the non-zero genetic effects of those loci over entire genome. The iteratively reweighted LASSO is therefore extended to mapping QTL for discrete traits, such as ordinal, binary, and Poisson traits. The simulated and real data analyses are conducted to demonstrate the efficiency of the proposed method to simultaneously identify multiple QTLs for binary and Poisson traits as examples.     

The Statistical Power of Inclusive Composite Interval Mapping in Detecting Digenic Epistasis Showing Common F2 Segregation Ratios

Epistasis is a commonly observed genetic phenomenon and an important source of variation of complex traits,which could maintain additive variance and therefore assure the long-term genetic gain in breeding.Inclusive composite interval mapping（ICIM） is able to identify epistatic quantitative trait loci（QTLs） no matter whether the two interacting QTLs have any additive effects.In this article,we conducted a simulation study to evaluate detection power and false discovery rate（FDR） of ICIM epistatic mapping,by considering F₂ and doubled haploid（DH） populations,different F₂ segregation ratios and population sizes.Results indicated that estimations of QTL locations and effects were unbiased,and the detection power of epistatic mapping was largely affected by population size,heritability of epistasis,and the amount and distribution of genetic effects.When the same likelihood of odd（LOD） threshold was used,detection power of QTL was higher in F₂ population than power in DH population;meanwhile FDR in F₂ was also higher than that in DH.The increase of marker density from 10 cM to 5 cM led to similar detection power but higher FDR.In simulated populations,ICIM achieved better mapping results than multiple interval mapping（MIM） in estimation of QTL positions and effect.At the end,we gave epistatic mapping results of ICIM in one actual population in rice（Oryza sativa L.）.     

A Sequential Quantitative Trait Locus Fine-Mapping Strategy Using Recombinant-Derived Progeny

A thorough understanding of the quantitative trait loci(QTLs)that underlie agronomically important traits in crops would greatly increase agricultural productivity.Although advances have been made in QTL cloning,the majority of QTLs remain unknown because of their low heritability and minor contributions to phenotypic performance.Here we summarize the key advantages and disadvantages of current QTL fine-mapping methodologies,and then introduce a sequential QTL fine-mapping strategy based on both genotypes and phenotypes of progeny derived from recombinants.With this mapping strategy,experimental errors could be dramatically diminished so as to reveal the authentic genetic effect of target QTLs.The number of progeny required to detect QTLs atvarious R~2 values was calculated,and the backcross generation suitable to start QTL fine-mapping was also estimated.This mapping strategy has proved to be very powerful in narrowing down QTL regions,particularly minor-effect QTLs,as revealed by fine-mapping of various resistance QTLs in maize.Application of this sequential QTL mapping strategy should accelerate cloning of agronomically important QTLs,which is currently a substantial challenge in crops.