An empirical Bayes method for estimating epistatic effects of quantitative trait loci

Summary .   The genetic variance of a quantitative trait is often controlled by the segregation of multiple interacting loci. Linear model regression analysis is usually applied to estimating and testing effects of these quantitative trait loci (QTL). Including all the main effects and the effects of interaction (epistatic effects), the dimension of the linear model can be extremely high. Variable selection via stepwise regression or stochastic search variable selection (SSVS) is the common procedure for epistatic effect QTL analysis. These methods are computationally intensive, yet they may not be optimal. The LASSO (least absolute shrinkage and selection operator) method is computationally more efficient than the above methods. As a result, it has been widely used in regression analysis for large models. However, LASSO has never been applied to genetic mapping for epistatic QTL, where the number of model effects is typically many times larger than the sample size. In this study, we developed an empirical Bayes method (E-BAYES) to map epistatic QTL under the mixed model framework. We also tested the feasibility of using LASSO to estimate epistatic effects, examined the fully Bayesian SSVS, and reevaluated the penalized likelihood (PENAL) methods in mapping epistatic QTL. Simulation studies showed that all the above methods performed satisfactorily well. However, E-BAYES appears to outperform all other methods in terms of minimizing the mean-squared error (MSE) with relatively short computing time. Application of the new method to real data was demonstrated using a barley dataset.     

Shrinkage Estimation Method for Mapping Multiple Quantitative Trait Loci

In this article, shrinkage estimation method for multiple-marker analysis and for mapping multiple quantitative trait loci (QTL) was reviewed. For multiple-marker analysis, Xu (Genetics, 2003, 163:789-801) developed a Bayesian shrinkage estimation (BSE) method. The key to the success of this method is to allow each marker effect have its own variance parameter, which in turn has its own prior distribution so that the variance can be estimated from the data. Under this hierarchical model, a large number of markers can be handled although most of them may have negligible effects. Under epistatic genetic model, however, the running time is very long. To overcome this problem, a novel method of incorporating the idea described above into maximum likelihood, known as penalized likelihood method, was proposed. A simulated study showed that this method can handle a model with multiple effects, which are ten times larger than the sample size. For multiple QTL analysis, two modified versions for the BSE method were introduced: one is the fixed-interval method and another is the variable-interval method. The former deals with markers with intermediate density, and the latter can handle markers with extremely high density as well as model with epistatic effects. For the detection of epistatic effects, penalized likelihood method and the variable-interval approach of the BSE method are available.     

多QTL定位的压缩估计方法

本文综述了多标记分析和多QTL定位的压缩估计方法。对于前者,Xu（Genetics,2003,163：789—801）首先提出了Bayesian压缩估计方法。其关键在于让每个效应有一个特定的方差参数,而该方差又服从一定的先验分布,以致能从资料中估计之。由此,能够同时估计大量分子标记基因座的遗传效应,即使大多数标记的效应是可忽略的。然而,对于上位性遗传模型,其运算时间还是过长。为此,笔者将上述思想嵌入极大似然法,提出了惩罚最大似然方法。模拟研究显示：该方法能处理变量个数大于样本容量10倍左右的线性遗传模型。对于后者,本文详细介绍了基于固定区间和可变区间的Bayesian压缩估计方法。固定区间方法可处理中等密度的分子标记资料;可变区间方法则可分析高密度分子标记资料,甚至是上位性遗传模型。对于上位性检测,已介绍的惩罚最大似然方法和可变区间Bayesian压缩估计方法可供利用。应当指出,压缩估计方法在今后的eQTL和QTN定位以及基因互作网络分析等研究中也是有应用价值的。     

A penalized likelihood method for mapping epistatic quantitative trait Loci with one-dimensional genome searches

Epistasis is a common and important phenomenon, as indicated by results from a number of recent experiments. Unfortunately, the discovery of epistatic quantitative trait loci (QTL) is difficult since one must search for multiple QTL simultaneously in two or more dimensions. Such a multidimensional search necessitates many statistical tests, and a high statistical threshold must be adopted to avoid false positives. Furthermore, the large number of (interaction) parameters in comparison with the number of observations results in a serious danger of overfitting and overinterpretation of the data. In this article we present a new statistical framework for mapping epistasis in inbred line crosses. It is based on reducing the high dimensionality of the problem in two ways. First, epistatic QTL are mapped in a one-dimensional genome scan for high interactions between QTL and the genetic background. Second, the dimension of the search is bounded by penalized likelihood methods. We use simulated backcross data to illustrate the new approach.     

Bayesian model choice and search strategies for mapping interacting quantitative trait Loci

Most complex traits of animals, plants, and humans are influenced by multiple genetic and environmental factors. Interactions among multiple genes play fundamental roles in the genetic control and evolution of complex traits. Statistical modeling of interaction effects in quantitative trait loci (QTL) analysis must accommodate a very large number of potential genetic effects, which presents a major challenge to determining the genetic model with respect to the number of QTL, their positions, and their genetic effects. In this study, we use the methodology of Bayesian model and variable selection to develop strategies for identifying multiple QTL with complex epistatic patterns in experimental designs with two segregating genotypes. Specifically, we develop a reversible jump Markov chain Monte Carlo algorithm to determine the number of QTL and to select main and epistatic effects. With the proposed method, we can jointly infer the genetic model of a complex trait and the associated genetic parameters, including the number, positions, and main and epistatic effects of the identified QTL. Our method can map a large number of QTL with any combination of main and epistatic effects. Utility and flexibility of the method are demonstrated using both simulated data and a real data set. Sensitivity of posterior inference to prior specifications of the number and genetic effects of QTL is investigated.     

Mapping of epistatic quantitative trait loci in four-way crosses

Four-way crosses (4WC) involving four different inbred lines often appear in plant and animal commercial breeding programs. Direct mapping of quantitative trait loci (QTL) in these commercial populations is both economical and practical. However, the existing statistical methods for mapping QTL in a 4WC population are built on the single-QTL genetic model. This simple genetic model fails to take into account QTL interactions, which play an important role in the genetic architecture of complex traits. In this paper, therefore, we attempted to develop a statistical method to detect epistatic QTL in 4WC population. Conditional probabilities of QTL genotypes, computed by the multi-point single locus method, were used to sample the genotypes of all putative QTL in the entire genome. The sampled genotypes were used to construct the design matrix for QTL effects. All QTL effects, including main and epistatic effects, were simultaneously estimated by the penalized maximum likelihood method. The proposed method was confirmed by a series of Monte Carlo simulation studies and real data analysis of cotton. The new method will provide novel tools for the genetic dissection of complex traits, construction of QTL networks, and analysis of heterosis.     

Bias correction for estimated QTL effects using the penalized maximum likelihood method

A penalized maximum likelihood method has been proposed as an important approach to the detection of epistatic quantitative trait loci (QTL). However, this approach is not optimal in two special situations: (1) closely linked QTL with effects in opposite directions and (2) small-effect QTL, because the method produces downwardly biased estimates of QTL effects. The present study aims to correct the bias by using correction coefficients and shifting from the use of a uniform prior on the variance parameter of a QTL effect to that of a scaled inverse chi-square prior. The results of Monte Carlo simulation experiments show that the improved method increases the power from 25 to 88% in the detection of two closely linked QTL of equal size in opposite directions and from 60 to 80% in the identification of QTL with small effects (0.5% of the total phenotypic variance). We used the improved method to detect QTL responsible for the barley kernel weight trait using 145 doubled haploid lines developed in the North American Barley Genome Mapping Project. Application of the proposed method to other shrinkage estimation of QTL effects is discussed.     

Mapping epistatic quantitative trait loci with one-dimensional genome searches

The discovery of epistatically interacting QTL is hampered by the intractability and low power to detect QTL in multidimensional genome searches. We describe a new method that maps epistatic QTL by identifying loci of high QTL by genetic background interaction. This approach allows detection of QTL involved not only in pairwise but also higher-order interaction, and does so with one-dimensional genome searches. The approach requires large populations derived from multiple related inbred-line crosses as is more typically available for plants. Using maximum likelihood, the method contrasts models in which QTL allelic values are either nested within, or fixed over, populations. We apply the method to simulated doubled-haploid populations derived from a diallel among three inbred parents and illustrate the power of the method to detect QTL of different effect size and different levels of QTL by genetic background interaction. Further, we show how the method can be used in conjunction with standard two-locus QTL detection models that use two-dimensional genome searches and find that the method may double the power to detect first-order epistasis.     

Bayesian analysis for genetic architecture of dynamic traits

The dissection of the genetic architecture of quantitative traits, including the number and locations of quantitative trait loci (QTL) and their main and epistatic effects, has been an important topic in current QTL mapping. We extend the Bayesian model selection framework for mapping multiple epistatic QTL affecting continuous traits to dynamic traits in experimental crosses. The extension inherits the efficiency of Bayesian model selection and the flexibility of the Legendre polynomial model fitting to the change in genetic and environmental effects with time. We illustrate the proposed method by simultaneously detecting the main and epistatic QTLs for the growth of leaf age in a doubled-haploid population of rice. The behavior and performance of the method are also shown by computer simulation experiments. The results show that our method can more quickly identify interacting QTLs for dynamic traits in the models with many numbers of genetic effects, enhancing our understanding of genetic architecture for dynamic traits. Our proposed method can be treated as a general form of mapping QTL for continuous quantitative traits, being easier to extend to multiple traits and to a single trait with repeat records.     

Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse

The extent to which epistasis contributes to adaptation, population differentiation, and speciation is a long-standing and important problem in evolutionary genetics. Using recombinant inbred (RI) lines of Arabidopsis thaliana grown under natural field conditions, we have examined the genetic architecture of fitness-correlated traits with respect to epistasis; we identified both single-locus additive and two-locus epistatic QTL for natural variation in fruit number, germination, and seed length and width. For fruit number, we found seven significant epistatic interactions, but only two additive QTL. For seed germination, length, and width, there were from two to four additive QTL and from five to eight epistatic interactions. The epistatic interactions were both positive and negative. In each case, the magnitude of the epistatic effects was roughly double that of the effects of the additive QTL, varying from -41% to +29% for fruit number and from -5% to +4% for seed germination, length, and width. A number of the QTL that we describe participate in more than one epistatic interaction, and some loci identified as additive also may participate in an epistatic interaction; the genetic architecture for fitness traits may be a network of additive and epistatic effects. We compared the map positions of the additive and epistatic QTL for germination, seed width, and seed length from plants grown in both the field and the greenhouse. While the total number of significant additive and epistatic QTL was similar under the two growth conditions, the map locations were largely different. We found a small number of significant epistatic QTL x environment effects when we tested directly for them. Our results support the idea that epistatic interactions are an important part of natural genetic variation and reinforce the need for caution in comparing results from greenhouse-grown and field-grown plants.     

Bayesian model selection for genome-wide epistatic quantitative trait loci analysis

The problem of identifying complex epistatic quantitative trait loci (QTL) across the entire genome continues to be a formidable challenge for geneticists. The complexity of genome-wide epistatic analysis results mainly from the number of QTL being unknown and the number of possible epistatic effects being huge. In this article, we use a composite model space approach to develop a Bayesian model selection framework for identifying epistatic QTL for complex traits in experimental crosses from two inbred lines. By placing a liberal constraint on the upper bound of the number of detectable QTL we restrict attention to models of fixed dimension, greatly simplifying calculations. Indicators specify which main and epistatic effects of putative QTL are included. We detail how to use prior knowledge to bound the number of detectable QTL and to specify prior distributions for indicators of genetic effects. We develop a computationally efficient Markov chain Monte Carlo (MCMC) algorithm using the Gibbs sampler and Metropolis-Hastings algorithm to explore the posterior distribution. We illustrate the proposed method by detecting new epistatic QTL for obesity in a backcross of CAST/Ei mice onto M16i.     

Epistatic interaction is an important genetic basis of grain yield and its components in maize

A population of 294 recombinant inbred lines (RIL) derived from Yuyu22, an elite maize hybrid extending broadly in China, has been constructed to investigate the genetic basis of grain yield, and associated yield components in maize. The main-effect quantitative trait loci (QTL), digenic epistatic interactions, and their interactions with the environment for grain yield and its three components were identified by using the mixed linear model approach. Thirty-two main-effect QTL and forty-four pairs of digenic epistatic interactions were detected for the four measured traits in four environments. Our results suggest that both additive effects and epistasis (additive × additive) effects are important genetic bases of grain yield and its components in the RIL population. Only 30.4% of main-effect QTL for ear length were involved in epistatic interactions. This implies that many loci in epistatic interactions may not have significant effects for traits alone but may affect trait expression by epistatic interaction with the other loci.     

Identifying quantitative trait locus by genetic background interactions in association studies

Association studies are designed to identify main effects of alleles across a potentially wide range of genetic backgrounds. To control for spurious associations, effects of the genetic background itself are often incorporated into the linear model, either in the form of subpopulation effects in the case of structure or in the form of genetic relationship matrices in the case of complex pedigrees. In this context epistatic interactions between loci can be captured as an interaction effect between the associated locus and the genetic background. In this study I developed genetic and statistical models to tie the locus by genetic background interaction idea back to more standard concepts of epistasis when genetic background is modeled using an additive relationship matrix. I also simulated epistatic interactions in four-generation randomly mating pedigrees and evaluated the ability of the statistical models to identify when a biallelic associated locus was epistatic to other loci. Under additive-by-additive epistasis, when interaction effects of the associated locus were quite large (explaining 20% of the phenotypic variance), epistasis was detected in 79% of pedigrees containing 320 individuals. The epistatic model also predicted the genotypic value of progeny better than a standard additive model in 78% of simulations. When interaction effects were smaller (although still fairly large, explaining 5% of the phenotypic variance), epistasis was detected in only 9% of pedigrees containing 320 individuals and the epistatic and additive models were equally effective at predicting the genotypic values of progeny. Epistasis was detected with the same power whether the overall epistatic effect was the result of a single pairwise interaction or the sum of nine pairwise interactions, each generating one ninth of the epistatic variance. The power to detect epistasis was highest (94%) at low QTL minor allele frequency, fell to a minimum (60%) at minor allele frequency of about 0.2, and then plateaued at about 80% as alleles reached intermediate frequencies. The power to detect epistasis declined when the linkage disequilibrium between the DNA marker and the functional polymorphism was not complete.     

Modeling quantitative trait Loci and interpretation of models

A quantitative genetic model relates the genotypic value of an individual to the alleles at the loci that contribute to the variation in a population in terms of additive, dominance, and epistatic effects. This partition of genetic effects is related to the partition of genetic variance. A number of models have been proposed to describe this relationship: some are based on the orthogonal partition of genetic variance in an equilibrium population. We compare a few representative models and discuss their utility and potential problems for analyzing quantitative trait loci (QTL) in a segregating population. An orthogonal model implies that estimates of the genetic effects are consistent in a full or reduced model in an equilibrium population and are directly related to the partition of the genetic variance in the population. Linkage disequilibrium does not affect the estimation of genetic effects in a full model, but would in a reduced model. Certainly linkage disequilibrium would complicate the detection of QTL and epistasis. Using different models does not influence the detection of QTL and epistasis. However, it does influence the estimation and interpretation of genetic effects.     

Simultaneous mapping of epistatic QTL in DU6i × DBA/2 mice

We have mapped epistatic quantitative trait loci (QTL) in an F₂ cross between DU6i × DBA/2 mice. By including these epistatic QTL and their interaction parameters in the genetic model, we were able to increase the genetic variance explained substantially (8.8%–128.3%) for several growth and body composition traits. We used an analysis method based on a simultaneous search for epistatic QTL pairs without assuming that the QTL had any effect individually. We were able to detect several QTL that could not be detected in a search for marginal QTL effects because the epistasis cancelled out the individual effects of the QTL. In total, 23 genomic regions were found to contain QTL affecting one or several of the traits and eight of these QTL did not have significant individual effects. We identified 44 QTL pairs with significant effects on the traits, and, for 28 of the pairs, an epistatic QTL model fit the data significantly better than a model without interactions. The epistatic pairs were classified by the significance of the epistatic parameters in the genetic model, and visual inspection of the two-locus genotype means identified six types of related genotype–phenotype patterns among the pairs. Five of these patterns resembled previously published patterns of QTL interactions.     

Modeling epistasis of quantitative trait loci using Cockerham's model

We use the orthogonal contrast scales proposed by Cockerham to construct a genetic model, called Cockerham's model, for studying epistasis between genes. The properties of Cockerham's model in modeling and mapping epistatic genes under linkage equilibrium and disequilibrium are investigated and discussed. Because of its orthogonal property, Cockerham's model has several advantages in partitioning genetic variance into components, interpreting and estimating gene effects, and application to quantitative trait loci (QTL) mapping when compared to other models, and thus it can facilitate the study of epistasis between genes and be readily used in QTL mapping. The issues of QTL mapping with epistasis are also addressed. Real and simulated examples are used to illustrate Cockerham's model, compare different models, and map for epistatic QTL. Finally, we extend Cockerham's model to multiple loci and discuss its applications to QTL mapping.     

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.     

A Model Selection Approach for the Identification of Quantitative Trait Loci in Experimental Crosses, Allowing Epistasis

The identification of quantitative trait loci (QTL) and their interactions is a crucial step toward the discovery of genes responsible for variation in experimental crosses. The problem is best viewed as one of model selection, and the most important aspect of the problem is the comparison of models of different sizes. We present a penalized likelihood approach, with penalties on QTL and pairwise interactions chosen to control false positive rates. This extends the work of Broman and Speed to allow for pairwise interactions among QTL. A conservative version of our penalized LOD score provides strict control over the rate of extraneous QTL and interactions; a more liberal criterion is more lenient on interactions but seeks to maintain control over the rate of inclusion of false loci. The key advance is that one needs only to specify a target false positive rate rather than a prior on the number of QTL and interactions. We illustrate the use of our model selection criteria as exploratory tools; simulation studies demonstrate reasonable power to detect QTL. Our liberal criterion is comparable in power to two Bayesian approaches.     

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.     

Evidence for widespread epistatic interactions influencing Marek's disease virus viremia levels in chicken

Marek's disease (MD), a T cell lymphoma induced by the Marek's disease virus (MDV), is the main chronic infectious disease concern threatening the poultry industry. Enhancing genetic resistance to MD in commercial poultry is an attractive method to augment MD vaccines, which is currently the control method of choice. In order to implement this control strategy through marker-assisted selection (MAS), it is necessary to identify quantitative trait loci (QTL) or genes that influence MD incidence. Previous studies have demonstrated that it is possible to identify QTL that confer MD resistance in both experimental and commercial chickens. With the advent of the chicken genome sequence and new genomic tools, and evidence that interactions are important in understanding complex traits, the line 6 x 7 F(2) experimental resource population was re-evaluated with finer resolution for epistatic interactions. The F(2) population, consisting of 272 individuals and previously genotyped with 133 genetic markers, was combined along with 576 additional single nucleotide polymorphisms (SNPs) genotyped on 80 individuals in each of the distribution tails for MD and other associated traits, and tested for the presence of main effects and two-way epistatic interactions accounting for MD incidence, viremia titers, and length of survival. Main effects were generally not significant but a large number of highly significant interactions, involving loci located throughout the genome, were identified that account for MDV viremia titers in infected birds. These results suggest that resistance to MD is highly complex and will require the incorporation of epistatic interaction analyses and functional genomic approaches to reveal the underlying genetic basis.