The effect of epistasis between linked genes on quantitative trait locus analysis

The effect of epistasis between linked genes on quantitative trait locus (QTL) analysis was studied as a function of their contribution to the phenotypic variance and their genetic distance by simulation of F2 (at least 200 individuals) and recombinant inbred line (RIL) populations. Data sets were replicated 100 times. For F2 populations, the presence of epistasis improves the detection of QTLs having effects in opposite directions. Epistasis between linked QTLs (26.5 cM) was poorly detected even when its contribution was relatively high compared to the main effects, and was null for heritabilities lower than 0.10. The detection of false-positive main effects is strongly affected by the distance between epistatic QTLs. The closer they are (≤11.5 cM), the higher the probability of detecting false-positive main-effect QTLs and the lower the probability of detecting epistatic effects. In this case, the presence of main-effect QTLs is due to the deviation of the heterozygote from the homozygotes at each linked interacting QTL and is algebraically explained by the joint effect of the linkage and the additive-by-additive interaction, resulting in a heterosis at a single genomic region in the absence of simulated dominant genetic effects. The number of false-positive main effects only reached nominal levels at about 100 cM. For RIL populations, the number of false positives or the detection of existing epistasis does not depend on the distance, and the power to detect epistatic QTLs is much higher even with small sample sizes and low contributions to the trait. RIL populations are highly recommended to detect epistatic QTLs and to better infer the genetic architecture of a quantitative trait.     

High-resolution quantitative trait locus mapping reveals sign epistasis controlling ovariole number between two Drosophila species

Identifying the genes underlying genetically complex traits is of fundamental importance for medicine, agriculture, and evolutionary biology. However, the level of resolution offered by traditional quantitative trait locus (QTL) mapping is usually coarse. We analyze here a trait closely related to fitness, ovariole number. Our initial interspecific mapping between Drosophila sechellia (8 ovarioles/ovary) and D. simulans (15 ovarioles/ovary) identified a major QTL on chromosome 3 and a minor QTL on chromosome 2. To refine the position of the major QTL, we selected 1038 additional recombinants in the region of interest using flanking morphological markers (selective phenotyping). This effort generated approximately one recombination event per gene and increased the mapping resolution by approximately seven times. Our study thus shows that using visible markers to select for recombinants can efficiently increase the resolution of QTL mapping. We resolved the major QTL into two epistatic QTL, QTL3a and QTL3b. QTL3a shows sign epistasis: it has opposite effects in two different genetic backgrounds, the presence vs. the absence of the QTL3b D. sechellia allele. This property of QTL3a allows us to reconstruct the probable order of fixation of the QTL alleles during evolution.     

Selective genotyping with epistasis can be utilized for a major quantitative trait locus mapping in hypertension in rats

Epistasis used to be considered an obstacle in mapping quantitative trait loci (QTL) despite its significance. Numerous epistases have proved to be involved in quantitative genetics. We established a backcross model that demonstrates a major QTL for hypertension (Ht). Seventy-eight backcrossed rats (BC), derived from spontaneously hypertensive rats (SHR) and normotensive Fischer 344 rats, showed bimodal distribution of systolic blood pressure (BP) values and a phenotypic segregation ratio consistent with 1:1. In this backcross analysis, sarco(endo)plasmic reticulum Ca(2+)-dependent ATPase (Serca) II heterozygotes showed widespread bimodality in frequency distribution of BP values and obviously demonstrated Ht. First, in genome-wide screening, Mapmaker/QTL analysis mapped Ht at a locus between D1Mgh8 and D1Mit4 near Sa in all 78 BC. The peak logarithm of the odds (LOD) score reached 5.3. Second, Serca II heterozygous and homozygous BC were analyzed separately using Mapmaker/QTL. In the 35 Serca II heterozygous BC, the peak LOD score was 3.8 at the same locus whereas it did not reach statistical significance in the 43 Serca II homozygotes. Third, to map Ht efficiently, we selected 18 Serca II heterozygous BC with 9 highest and 9 lowest BP values. In these 18 BC, the peak LOD score reached 8.1. In 17 of the 18, D1Mgh8 genotypes (homo or hetero) qualitatively cosegregated with BP phenotypes (high or low) (P < 0.0001, by chi-square analysis). In conclusion, selective genotyping with epistasis can be utilized for a major QTL mapping near Sa on chromosome 1 in SHR.     

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.     

Derivation of the shrinkage estimates of quantitative trait locus effects

The shrinkage estimate of a quantitative trait locus (QTL) effect is the posterior mean of the QTL effect when a normal prior distribution is assigned to the QTL. This note gives the derivation of the shrinkage estimate under the multivariate linear model. An important lemma regarding the posterior mean of a normal likelihood combined with a normal prior is introduced. The lemma is then used to derive the Bayesian shrinkage estimates of the QTL effects.     

Estimation of quantitative trait locus allele frequency via a modified granddaughter design

A method is described on the basis of a modification of the granddaughter design to obtain estimates of quantitative trait loci (QTL) allele frequencies in dairy cattle populations and to determine QTL genotypes for both homozygous and heterozygous grandsires. The method is based on determining the QTL allele passed from grandsires to their maternal granddaughters using haplotypes consisting of several closely linked genetic markers. This method was applied to simulated data of 10 grandsire families, each with 500 granddaughters, and a QTL with a substitution effect of 0.4 phenotypic standard deviations and to actual data for a previously analyzed QTL in the center of chromosome 6, with substitution effect of 1 phenotypic standard deviation on protein percentage. In the simulated data the standard error for the estimated QTL substitution effect with four closely linked multiallelic markers was only 7% greater than the expected standard error with completely correct identification of QTL allele origin. The method estimated the population QTL allelic frequency as 0.64 +/- 0.07, compared to the simulated value of 0.7. In the actual data, the frequency of the allele that increases protein percentage was estimated as 0.63 +/- 0.06. In both data sets the hypothesis of equal allelic frequencies was rejected at P < 0.05.     

More about quantitative trait locus mapping with diallel designs

We present a general regression-based method for mapping quantitative trait loci (QTL) by combining different populations derived from diallel designs. The model expresses, at any map position, the phenotypic value of each individual as a function of the specific-mean of the population to which the individual belongs, the additive and dominance effects of the alleles carried by the parents of that population and the probabilities of QTL genotypes conditional on those of neighbouring markers. Standard linear model procedures (ordinary or iteratively reweighted least-squares) are used for estimation and test of the parameters.     

A haplotype-based algorithm for multilocus linkage disequilibrium mapping of quantitative trait loci with epistasis

For tightly linked loci, cosegregation may lead to nonrandom associations between alleles in a population. Because of its evolutionary relationship with linkage, this phenomenon is called linkage disequilibrium. Today, linkage disequilibrium-based mapping has become a major focus of recent genome research into mapping complex traits. In this article, we present a new statistical method for mapping quantitative trait loci (QTL) of additive, dominant, and epistatic effects in equilibrium natural populations. Our method is based on haplotype analysis of multilocus linkage disequilibrium and exhibits two significant advantages over current disequilibrium mapping methods. First, we have derived closed-form solutions for estimating the marker-QTL haplotype frequencies within the maximum-likelihood framework implemented by the EM algorithm. The allele frequencies of putative QTL and their linkage disequilibria with the markers are estimated by solving a system of regular equations. This procedure has significantly improved the computational efficiency and the precision of parameter estimation. Second, our method can detect marker-QTL disequilibria of different orders and QTL epistatic interactions of various kinds on the basis of a multilocus analysis. This can not only enhance the precision of parameter estimation, but also make it possible to perform whole-genome association studies. We carried out extensive simulation studies to examine the robustness and statistical performance of our method. The application of the new method was validated using a case study from humans, in which we successfully detected significant QTL affecting human body heights. Finally, we discuss the implications of our method for genome projects and its extension to a broader circumstance. The computer program for the method proposed in this article is available at the webpage http://www.ifasstat.ufl.edu/genome/~LD.     

棉花数量性状基因定位研究进展

棉花的许多重要性状,包括产量、纤维品质、株型、抗病抗逆性、生理生化等都是数量性状,受遗传和环境因子的共同作用。近年来,随着分子生物学技术的进步,棉花基因组研究得到迅速发展,为棉花数量性状基因（quantitative trait locus,QTL）定位奠定了坚实的基础。概述了近十几年来棉花QTL定位研究及分子标记辅助选择的进展,结合研究实践指出了棉花QTL定位及标记辅助选择存在的问题,并对其发展方向做出了初步探讨。     

Mapping quantitative trait loci with epistatic effects

Epistatic variance can be an important source of variation for complex traits. However, detecting epistatic effects is difficult primarily due to insufficient sample sizes and lack of robust statistical methods. In this paper, we develop a Bayesian method to map multiple quantitative trait loci (QTLs) with epistatic effects. The method can map QTLs in complicated mating designs derived from the cross of two inbred lines. In addition to mapping QTLs for quantitative traits, the proposed method can even map genes underlying binary traits such as disease susceptibility using the threshold model. The parameters of interest are various QTL effects, including additive, dominance and epistatic effects of QTLs, the locations of identified QTLs and even the number of QTLs. When the number of QTLs is treated as an unknown parameter, the dimension of the model becomes a variable. This requires the reversible jump Markov chain Monte Carlo algorithm. The utility of the proposed method is demonstrated through analysis of simulation data.     

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.     

Quantitative trait locus analysis of longitudinal quantitative trait data in complex pedigrees

There is currently considerable interest in genetic analysis of quantitative traits such as blood pressure and body mass index. Despite the fact that these traits change throughout life they are commonly analyzed only at a single time point. The genetic basis of such traits can be better understood by collecting and effectively analyzing longitudinal data. Analyses of these data are complicated by the need to incorporate information from complex pedigree structures and genetic markers. We propose conducting longitudinal quantitative trait locus (QTL) analyses on such data sets by using a flexible random regression estimation technique. The relationship between genetic effects at different ages is efficiently modeled using covariance functions (CFs). Using simulated data we show that the change in genetic effects over time can be well characterized using CFs and that including parameters to model the change in effect with age can provide substantial increases in power to detect QTL compared with repeated measure or univariate techniques. The asymptotic distributions of the methods used are investigated and methods for overcoming the practical difficulties in fitting CFs are discussed. The CF-based techniques should allow efficient multivariate analyses of many data sets in human and natural population genetics.     

Complex genetic interactions in a quantitative trait locus

Whether in natural populations or between two unrelated members of a species, most phenotypic variation is quantitative. To analyze such quantitative traits, one must first map the underlying quantitative trait loci. Next, and far more difficult, one must identify the quantitative trait genes (QTGs), characterize QTG interactions, and identify the phenotypically relevant polymorphisms to determine how QTGs contribute to phenotype. In this work, we analyzed three Saccharomyces cerevisiae high-temperature growth (Htg) QTGs (MKT1, END3, and RHO2). We observed a high level of genetic interactions among QTGs and strain background. Interestingly, while the MKT1 and END3 coding polymorphisms contribute to phenotype, it is the RHO2 3′UTR polymorphisms that are phenotypically relevant. Reciprocal hemizygosity analysis of the Htg QTGs in hybrids between S288c and ten unrelated S. cerevisiae strains reveals that the contributions of the Htg QTGs are not conserved in nine other hybrids, which has implications for QTG identification by marker-trait association. Our findings demonstrate the variety and complexity of QTG contributions to phenotype, the impact of genetic background, and the value of quantitative genetic studies in S. cerevisiae.     

The X chromosome in quantitative trait locus mapping

The X chromosome requires special treatment in the mapping of quantitative trait loci (QTL). However, most QTL mapping methods, and most computer programs for QTL mapping, have focused exclusively on autosomal loci. We describe a method for appropriate treatment of the X chromosome for QTL mapping in experimental crosses. We address the important issue of formulating the null hypothesis of no linkage appropriately. If the X chromosome is treated like an autosome, a sex difference in the phenotype can lead to spurious linkage on the X chromosome. Further, the number of degrees of freedom for the linkage test may be different for the X chromosome than for autosomes, and so an X chromosome-specific significance threshold is required. To address this issue, we propose a general procedure to obtain chromosome-specific significance thresholds that controls the genomewide false positive rate at the desired level. We apply our methods to data on gut length in a large intercross of mice carrying the Sox10Dom mutation, a model of Hirschsprung disease. We identified QTL contributing to variation in gut length on chromosomes 5 and 18. We found suggestive evidence of linkage to the X chromosome, which would be viewed as strong evidence of linkage if the X chromosome was treated as an autosome. Our methods have been implemented in the package R/qtl.     

Rank-based statistical methodologies for quantitative trait locus mapping

This article addresses the identification of genetic loci (QTL and elsewhere) that influence nonnormal quantitative traits with focus on experimental crosses. QTL mapping is typically based on the assumption that the traits follow normal distributions, which may not be true in practice. Model-free tests have been proposed. However, nonparametric estimation of genetic effects has not been studied. We propose an estimation procedure based on the linear rank test statistics. The properties of the new procedure are compared with those of traditional likelihood-based interval mapping and regression interval mapping via simulations and a real data example. The results indicate that the nonparametric method is a competitive alternative to the existing parametric methodologies.     

Estimation of recombination parameters between a quantitative trait locus (QTL) and two marker gene loci

Summary A new method is described to obtain maximum likelihood estimates of recombination frequencies between quantitative trait loci (QTL) and marker gene loci; it is based on Fisher's method of scoring and numerical differentiation. The method is applied to data from chromosome-doubled monoploid lines of barley originating from the F₁ generation of a cross between two well-adapted barley varieties. The lines segregated for marker gene loci ddt (DDT resistance) and s (short rachilla hairs) on chromosome 7. The quantitative trait of single-kernel weight was found statistically significantly associated with locus s, but not with locus ddt. The association is ascribed to a QTL designated Kw1. It could not be ascribed to pleiotropism at locus s since the recombination frequency between s and Kw1 (0.26±0.09) differed significantly from zero. The recombination frequencies between Kw1 and ddt and between ddt and s were 0.42±0.07 and 0.31±0.03, respectively, suggesting the locus order ddt, s, Kw1. The segregation ratio for alleles in locus Kw1 was estimated to be 4357, which is not significantly different from a 11 ratio. Means and standard deviations of single-kernel weight for lines with either of the two Kw1 alleles were estimated; the Kw1 locus accounted for 25% of the variance of the single kernel weight.     

Efficient algorithms for quantitative trait loci mapping problems.

Rapid advances in molecular genetics push the need for efficient data analysis. Advanced algorithms are necessary for extracting all possible information from large experimental data sets. We present a general linear algebra framework for quantitative trait loci (QTL) mapping, using both linear regression and maximum likelihood estimation. The formulation simplifies future comparisons between and theoretical analyses of the methods. We show how the common structure of QTL analysis models can be used to improve the kernel algorithms, drastically reducing the computational effort while retaining the original analysis results. We have evaluated our new algorithms on data sets originating from two large F(2) populations of domestic animals. Using an updating approach, we show that 1-3 orders of magnitude reduction in computational demand can be achieved for matrix factorizations. For interval-mapping/composite-interval-mapping settings using a maximum likelihood model, we also show how to use the original EM algorithm instead of the ECM approximation, significantly improving the convergence and further reducing the computational time. The algorithmic improvements makes it feasible to perform analyses which have previously been deemed impractical or even impossible. For example, using the new algorithms, it is reasonable to perform permutation testing using exhaustive search on populations of 200 individuals using an epistatic two-QTL model.     

The effects of selected sampling on the transmission disequilibrium test of a quantitative trait locus

We investigate how sampling of parents or children based on their extreme phenotypic values selected from clinical databases would affect the power of identification of quantitative trait loci (QTL) by a transmission disequilibrium test (TDT). We consider three selective sampling schemes based on the selection of phenotypic values of parents or children in nuclear families: (1) two children, one of extreme value, the other random; (2) two children extremely discordant; (3) one parent of extreme value. Other family members not specified will be recruited randomly with regard to phenotypic values. Our study shows that the second sampling scheme can always enhance the power for QTL identification, sometimes dramatically so. The increase in the statistical power of the TDT is particularly dramatic when h2 at the QTL under test is small or intermediate (e.g. 0.05 or 0.10). For the other two sampling schemes, under dominant effects at the QTL, the power is always increased relative to random sampling; however, under recessive or additive genetic effects, the power gain is generally minor or even decreased a little sometimes. Allele frequencies at the QTL and the selection stringency are important for determining the effect of selective sampling on the power of QTL identification. Our study is useful as a practical guideline on how to perform the TDT efficiently in practice by taking advantage of the extensive databases accumulated that are enriched with people of extreme phenotypic values.