Multiple trait multiple interval mapping of quantitative trait loci from inbred line crosses

ABSTRACT: BACKGROUND: Although many experiments have measurements on multiple traits, most studies performed the analysis of mapping of quantitative trait loci (QTL) for each trait separately using single trait analysis. Single trait analysis does not take advantage of possible genetic and environmental correlations between traits. In this paper, we propose a novel statistical method for multiple trait multiple interval mapping (MTMIM) of QTL for inbred line crosses. We also develop a novel score-based method for estimating genome-wide significance level of putative QTL effects suitable for the MTMIM model. The MTMIM method is implemented in the freely available and widely used Windows QTL Cartographer software. RESULTS: Throughout the paper, we provide compelling empirical evidences that: (1) the score-based threshold maintains proper type I error rate and tends to keep false discovery rate within an acceptable level; (2) the MTMIM method can deliver better parameter estimates and power than single trait multiple interval mapping method; (3) an analysis of Drosophila dataset illustrates how the MTMIM method can better extract information from datasets with measurements in multiple traits. CONCLUSIONS: The MTMIM method represents a convenient statistical framework to test hypotheses of pleiotropic QTL versus closely linked nonpleiotropic QTL, QTL by environment interaction, and to estimate the total genotypic variance-covariance matrix between traits and to decompose it in terms of QTL-specific variance-covariance matrices, therefore, providing more details on the genetic architecture of complex traits.     

Mapping quantitative trait loci for longitudinal traits in line crosses

Quantitative traits whose phenotypic values change over time are called longitudinal traits. Genetic analyses of longitudinal traits can be conducted using any of the following approaches: (1) treating the phenotypic values at different time points as repeated measurements of the same trait and analyzing the trait under the repeated measurements framework, (2) treating the phenotypes measured from different time points as different traits and analyzing the traits jointly on the basis of the theory of multivariate analysis, and (3) fitting a growth curve to the phenotypic values across time points and analyzing the fitted parameters of the growth trajectory under the theory of multivariate analysis. The third approach has been used in QTL mapping for longitudinal traits by fitting the data to a logistic growth trajectory. This approach applies only to the particular S-shaped growth process. In practice, a longitudinal trait may show a trajectory of any shape. We demonstrate that one can describe a longitudinal trait with orthogonal polynomials, which are sufficiently general for fitting any shaped curve. We develop a mixed-model methodology for QTL mapping of longitudinal traits and a maximum-likelihood method for parameter estimation and statistical tests. The expectation-maximization (EM) algorithm is applied to search for the maximum-likelihood estimates of parameters. The method is verified with simulated data and demonstrated with experimental data from a pseudobackcross family of Populus (poplar) trees.     

Mapping quantitative trait loci using multiple families of line crosses.

To avoid a loss in statistical power as a result of homozygous individuals being selected as parents of a mapping population, one can use multiple families of line crosses for quantitative trait genetic linkage analysis. Two strategies of combining data are investigated: the fixed-model and the random-model strategies. The fixed-model approach estimates and tests the average effect of gene substitution for each parent, while the random-model approach treats each effect of gene substitution as a random variable and directly estimates and tests the variance of gene substitution. Extensive Monte Carlo simulations verify that the two strategies perform equally well, although the random model is preferable in combining data from a large number of families. Simulations also show that there may be an optimal sampling strategy (number of families vs. number of individuals per family) in which QTL mapping reaches its maximum power and minimum estimation error. Deviation from the optimal strategy reduces the efficiency of the method.     

Using the mixed model for interval mapping of quantitative trait loci in outbred line crosses

Interval mapping by simple regression is a powerful method for the detection of quantitative trait loci (QTLs) in line crosses such as F2 populations. Due to the ease of computation of the regression approach, relatively complex models with multiple fixed effects, interactions between QTLs or between QTLs and fixed effects can easily be accommodated. However, polygenic effects, which are not targeted in QTL analysis, cannot be treated as random effects in a least squares analysis. In a cross between true inbred lines this is of no consequence, as the polygenic effect contributes just to the residual variance. In a cross between outbred lines, however, if a trait has high polygenic heritability, the additive polygenic effect has a large influence on variation in the population. Here we extend the fixed model for the regression interval mapping method to a mixed model using an animal model. This makes it possible to use not only the observations from progeny (e.g. F2), but also those from the parents (F1) to evaluate QTLs and polygenic effects. We show how the animal model using parental observations can be applied to an outbred cross and so increase the power and accuracy of QTL analysis. Three estimation methods, i.e. regression and an animal model either with or without parental observations, are applied to simulated data. The animal model using parental observations is shown to have advantages in estimating QTL position and additive genotypic value, especially when the polygenic heritability is large and the number of progeny per parent is small.     

Bayesian joint mapping of quantitative trait loci for Gaussian and categorical characters in line crosses

Methodology for joint mapping of quantitative trait loci (QTL) affecting continuous and binary characters in experimental crosses is presented. The procedure consists of a Bayesian Gaussian-threshold model implemented via Markov chain Monte Carlo, which bypasses bottlenecks due to high-dimensional integrals required in maximum likelihood approaches. The method handles multiple binary traits and multiple QTL. Modeling of ordered categorical traits is discussed as well. Features of the method are illustrated using simulated datasets representing a backcross design, and the data are analyzed using mixed-trait and single-trait models. The mixed-trait analysis provides greater detection power of a QTL than a single-trait analysis when the QTL affects two or more traits. The number of QTL inferred in the mixed-trait analysis does not pertain to a specific trait, but the roles of each QTL on specific traits can be assessed from estimates of its effects. The impacts of varying incidence level and sample size on the mixed-trait QTL mapping analysis are investigated as well.     

A general mixture model for mapping quantitative trait loci by using molecular markers

Summary In a segregating population a quantitative trait may be considered to follow a mixture of (normal) distributions, the mixing proportions being based on Mendelian segregation rules. A general and flexible mixture model is proposed for mapping quantitative trait loci (QTLs) by using molecular markers. A method is discribed to fit the model to data. The model makes it possible to (1) analyse non-normally distributed traits such as lifetimes, counts or percentages in addition to normally distributed traits, (2) reduce environmental variation by taking into account the effects of experimental design factors and interaction between genotype and environment, (3) reduce genotypic variation by taking into account the effects of two or more QTLs simultaneously, (4) carry out a (combined) analysis of different population types, (5) estimate recombination frequencies between markers or use known marker distances, (6) cope with missing marker observations, (7) use markers as covariables in detection and mapping of QTLs, and finally to (8) implement the mapping in standard statistical packages.     

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     

A simple linear regression method for quantitative trait loci linkage analysis with censored observations

Standard quantitative trait loci (QTL) mapping techniques commonly assume that the trait is both fully observed and normally distributed. When considering survival or age-at-onset traits these assumptions are often incorrect. Methods have been developed to map QTL for survival traits; however, they are both computationally intensive and not available in standard genome analysis software packages. We propose a grouped linear regression method for the analysis of continuous survival data. Using simulation we compare this method to both the Cox and Weibull proportional hazards models and a standard linear regression method that ignores censoring. The grouped linear regression method is of equivalent power to both the Cox and Weibull proportional hazards methods and is significantly better than the standard linear regression method when censored observations are present. The method is also robust to the proportion of censored individuals and the underlying distribution of the trait. On the basis of linear regression methodology, the grouped linear regression model is computationally simple and fast and can be implemented readily in freely available statistical software.     

Linkage analysis of quantitative trait loci in multiple line crosses

Simple line crosses, for example, backcross and F₂, are commonly used in mapping quantitative trait loci (QTL). However, these simple crosses are rarely used alone in commercial plant breeding; rather, crosses involving multiple inbred lines or several simple crosses but connected by shared inbred lines may be common in plant breeding. Mapping QTL using crosses of multiple lines is more relevant to plant breeding. Unfortunately, current statistical methods and computer programs of QTL mapping are all designed for simple line crosses or multiple line crosses but under a regular mating system. It is not straightforward to extend the existing methods to handle multiple line crosses under irregular and complicated mating designs. The major hurdle comes from irregular inbreeding, multiple generations, and multiple alleles. In this study, we develop a Bayesian method implemented via the Markov chain Monte Carlo (MCMC) algorithm for mapping QTL using complicated multiple line crosses. With the MCMC algorithm, we are able to draw a complete path of the gene flow from founder alleles to their descendents via a recursive process. This has greatly simplified the problem caused by irregular mating and inbreeding in the mapping population. Adopting the reversible jump MCMC algorithm, we are able to simultaneously search for multiple QTL along the genome. We can even infer the posterior distribution of the number of QTL, one of the most important parameters in QTL study. Application of the new MCMC based QTL mapping procedure is demonstrated using two different mating designs. Design I involves two inbred lines and their derived F₁, F₂, and BC populations. Design II is a half-diallel cross involving three inbred lines. The two designs appear different, but can be handled with the same robust computer program.     

Advances in Bayesian multiple quantitative trait loci mapping in experimental crosses

Many complex human diseases and traits of biological and/or economic importance are determined by interacting networks of multiple quantitative trait loci (QTL) and environmental factors. Mapping QTL is critical for understanding the genetic basis of complex traits, and for ultimate identification of genes responsible. A variety of sophisticated statistical methods for QTL mapping have been developed. Among these developments, the evolution of Bayesian approaches for multiple QTL mapping over the past decade has been remarkable. Bayesian methods can jointly infer the number of QTL, their genomic positions and their genetic effects. Here, we review recently developed and still developing Bayesian methods and associated computer software for mapping multiple QTL in experimental crosses. We compare and contrast these methods to clearly describe the relationships among different Bayesian methods. We conclude this review by highlighting some areas of future research.     

Investigating the probability of sign inconsistency in the regression coefficients of markers flanking quantitative trait loci

Estimates of the locations and effects of quantitative trait loci (QTL) can be obtained by regressing phenotype on marker genotype. Under certain basic conditions, the signs of regression coefficients flanking QTL must be the same. There is no guarantee, however, that the signs of the regression coefficient estimates will be the same. We use sign inconsistency to describe the situation in which there is disagreement between the signs of the estimated regression coefficients flanking QTL. The presence of sign inconsistency can undermine the effectiveness of QTL mapping strategies that presume intervals whose markers have regression coefficient estimates of differing sign to be devoid of QTL. This article investigates the likelihood of sign inconsistency under various conditions. We derive an analytic expression for the approximate probability of sign inconsistency in the single-QTL case. We also examine sign inconsistency probabilities when multiple QTL are present through simulation. We have discovered that the probability of sign inconsistency can be unacceptably high, even when the conditions for QTL detection are otherwise quite favorable.     

Detection and parameter estimation for quantitative trait loci using regression models and multiple markers

A strategy of multi-step minimal conditional regression analysis has been developed to determine the existence of statistical testing and parameter estimation for a quantitative trait locus (QTL) that are unaffected by linked QTLs. The estimation of marker-QTL recombination frequency needs to consider only three cases: 1) the chromosome has only one QTL, 2) one side of the target QTL has one or more QTLs, and 3) either side of the target QTL has one or more QTLs. Analytical formula was derived to estimate marker-QTL recombination frequency for each of the three cases. The formula involves two flanking markers for case 1), two flanking markers plus a conditional marker for case 2), and two flanking markers plus two conditional markers for case 3). Each QTL variance and effect, and the total QTL variance were also estimated using analytical formulae. Simulation data show that the formulae for estimating marker-QTL recombination frequency could be a useful statistical tool for fine QTL mapping. With 1 000 observations, a QTL could be mapped to a narrow chromosome region of 1.5 cM if no linked QTL is present, and to a 2.8 cM chromosome region if either side of the target QTL has at least one linked QTL.     

Joint tests for quantitative trait loci in experimental crosses

Selective genotyping is common because it can increase the expected correlation between QTL genotype and phenotype and thus increase the statistical power of linkage tests (i.e., regression-based tests). Linkage can also be tested by assessing whether the marginal genotypic distribution conforms to its expectation, a marginal-based test. We developed a class of joint tests that, by constraining intercepts in regression-based analyses, capitalize on the information available in both regression-based and marginal-based tests. We simulated data corresponding to the null hypothesis of no QTL effect and the alternative of some QTL effect at the locus for a backcross and an F2 intercross between inbred strains. Regression-based and marginal-based tests were compared to corresponding joint tests. We studied the effects of random sampling, selective sampling from a single tail of the phenotypic distribution, and selective sampling from both tails of the phenotypic distribution. Joint tests were nearly as powerful as all competing alternatives for random sampling and two-tailed selection under both backcross and F2 intercross situations. Joint tests were generally more powerful for one-tailed selection under both backcross and F2 intercross situations. However, joint tests cannot be recommended for one-tailed selective genotyping if segregation distortion is suspected.     

Genomewide markers as cofactors for precision mapping of quantitative trait loci

In composite interval mapping of quantitative trait loci (QTL), subsets of background markers are used to account for the effects of QTL outside the marker interval being tested. Here, I propose a QTL mapping approach (called G model) that utilizes genomewide markers as cofactors. The G model involves backward elimination on a given chromosome after correcting for genomewide marker effects, calculated under a random effects model, at all the other chromosomes. I simulated a trait controlled by 15 or 30 QTL, mapping populations of N = 96, 192, and 384 recombinant inbreds, and N _M = 192 and 384 evenly spaced markers. In the C model, which utilized subsets of background markers, the number of QTL detected and the number of false positives depended on the number of cofactors used, with five cofactors being too few with N = 384 and 20–40 cofactors being too many with N = 96. A window size of 0 cM for excluding cofactors maintained the number of true QTL detected while decreasing the number of false positives. The number of true QTL detected was generally higher with the G model than with the C model, and the G model led to good control of the type I error rate in simulations where the null hypothesis of no marker–QTL linkage was true. Overall, the results indicated that the G model is useful in QTL mapping because it is less subjective and has equal, if not better, performance when compared with the traditional approach of using subsets of markers to account for background QTL.     

Combining data from multiple inbred line crosses improves the power and resolution of quantitative trait loci mapping

Rodent inbred line crosses are widely used to map genetic loci associated with complex traits. This approach has proven to be powerful for detecting quantitative trait loci (QTL); however, the resolution of QTL locations, typically approximately 20 cM, means that hundreds of genes are implicated as potential candidates. We describe analytical methods based on linear models to combine information available in two or more inbred line crosses. Our strategy is motivated by the hypothesis that common inbred strains of the laboratory mouse are derived from a limited ancestral gene pool and thus QTL detected in multiple crosses are likely to represent shared ancestral polymorphisms. We demonstrate that the combined-cross analysis can improve the power to detect weak QTL, can narrow support intervals for QTL regions, and can be used to separate multiple QTL that colocalize by chance. Moreover, combined-cross analysis can establish the allelic states of a QTL among a set of parental lines, thus providing critical information for narrowing QTL regions by haplotype analysis.     

Mapping quantitative trait loci in noninbred mosquito crosses

The identification of genes that affect quantitative traits has been of great interest to geneticists for many decades, and many statistical methods have been developed to map quantitative trait loci (QTL). Most QTL mapping studies in experimental organisms use purely inbred lines, where the two homologous chromosomes in each individual are identical. As a result, many existing QTL mapping methods developed for experimental organisms are applicable only to genetic crosses between inbred lines. However, it may be difficult to obtain inbred lines for certain organisms, e.g., mosquitoes. Although statistical methods for QTL mapping in outbred populations, e.g., humans, can be applied for such crosses, these methods may not fully take advantage of the uniqueness of these crosses. For example, we can generally assume that the two grandparental lines are homozygous at the QTL of interest, but such information is not be utilized through methods developed for outbred populations. In addition, mating types and phases can be relatively easy to establish through the analysis of adjacent markers due to the large number of offspring that can be collected, substantially simplifying the computational need. In this article, motivated by a mosquito intercross experiment involving two selected lines that are not genetically homozygous across the genome, we develop statistical methods for QTL mapping for genetic crosses involving noninbred lines. In our procedure, we first infer parental mating types and use likelihood-based methods to infer phases in each parent on the basis of genotypes of offspring and one parent. A hidden Markov model is then employed to estimate the number of high-risk alleles at marker positions and putative QTL positions between markers in each offspring, and QTL mapping is finally conducted through the inferred QTL configuration across all offspring in all crosses. The performance of the proposed methods is assessed through simulation studies, and the usefulness of this method is demonstrated through its application to a mosquito data set.     

Bayesian multiple quantitative trait loci mapping for complex traits using markers of the entire genome

A Bayesian methodology has been developed for multiple quantitative trait loci (QTL) mapping of complex binary traits that follow liability threshold models. Unlike most QTL mapping methods where only one or a few markers are used at a time, the proposed method utilizes all markers across the genome simultaneously. The outperformance of our Bayesian method over the traditional single-marker analysis and interval mapping has been illustrated via simulations and real data analysis to identify candidate loci associated with colorectal cancer.     

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.     

A powerful and robust method for mapping quantitative trait loci in general pedigrees

The variance-components model is the method of choice for mapping quantitative trait loci in general human pedigrees. This model assumes normally distributed trait values and includes a major gene effect, random polygenic and environmental effects, and covariate effects. Violation of the normality assumption has detrimental effects on the type I error and power. One possible way of achieving normality is to transform trait values. The true transformation is unknown in practice, and different transformations may yield conflicting results. In addition, the commonly used transformations are ineffective in dealing with outlying trait values. We propose a novel extension of the variance-components model that allows the true transformation function to be completely unspecified. We present efficient likelihood-based procedures to estimate variance components and to test for genetic linkage. Simulation studies demonstrated that the new method is as powerful as the existing variance-components methods when the normality assumption holds; when the normality assumption fails, the new method still provides accurate control of type I error and is substantially more powerful than the existing methods. We performed a genomewide scan of monoamine oxidase B for the Collaborative Study on the Genetics of Alcoholism. In that study, the results that are based on the existing variance-components method changed dramatically when three outlying trait values were excluded from the analysis, whereas our method yielded essentially the same answers with or without those three outliers. The computer program that implements the new method is freely available.