Joint-multiple family linkage analysis predicts within-family variation better than single-family analysis of the maize nested association mapping population

Quantitative trait locus (QTL) mapping has been used to dissect the genetic architecture of complex traits and predict phenotypes for marker-assisted selection. Many QTL mapping studies in plants have been limited to one biparental family population. Joint analysis of multiple biparental families offers an alternative approach to QTL mapping with a wider scope of inference. Joint-multiple population analysis should have higher power to detect QTL shared among multiple families, but may have lower power to detect rare QTL. We compared prediction ability of single-family and joint-family QTL analysis methods with fivefold cross-validation for 6 diverse traits using the maize nested association mapping population, which comprises 25 biparental recombinant inbred families. Joint-family QTL analysis had higher mean prediction abilities than single-family QTL analysis for all traits at most significance thresholds, and was always better at more stringent significance thresholds. Most robust QTL (detected in >50% of data samples) were restricted to one family and were often not detected at high frequency by joint-family analysis, implying substantial genetic heterogeneity among families for complex traits in maize. The superior predictive ability of joint-family QTL models despite important genetic differences among families suggests that joint-family models capture sufficient smaller effect QTL that are shared across families to compensate for missing some rare large-effect QTL.     

Statistical power of QTL mapping methods applied to bacteria counts

Most QTL mapping methods assume that phenotypes follow a normal distribution, but many phenotypes of interest are not normally distributed, e.g. bacteria counts (or colony-forming units, CFU). Such data are extremely skewed to the right and can present a high amount of zero values, which are ties from a statistical point of view. Our objective is therefore to assess the efficiency of four QTL mapping methods applied to bacteria counts: (1) least-squares (LS) analysis, (2) maximum-likelihood (ML) analysis, (3) non-parametric (NP) mapping and (4) nested ANOVA (AN). A transformation based on quantiles is used to mimic observed distributions of bacteria counts. Single positions (1 marker, 1 QTL) as well as chromosome scans (11 markers, 1 QTL) are simulated. When compared with the analysis of a normally distributed phenotype, the analysis of raw bacteria counts leads to a strong decrease in power for parametric methods, but no decrease is observed for NP. However, when a mathematical transformation (MT) is applied to bacteria counts prior to analysis, parametric methods have the same power as NP. Furthermore, parametric methods, when coupled with MT, outperform NP when bacteria counts have a very high proportion of zeros (70.8%). Our results show that the loss of power is mainly explained by the asymmetry of the phenotypic distribution, for parametric methods, and by the existence of ties, for the non-parametric method. Therefore, mapping of QTL for bacterial diseases, as well as for other diseases assessed by a counting process, should focus on the occurrence of ties in phenotypes before choosing the appropriate QTL mapping method.     

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.     

Interval Mapping of Multiple Quantitative Trait Loci

The interval mapping method is widely used for the mapping of quantitative trait loci (QTLs) in segregating generations derived from crosses between inbred lines. The efficiency of detecting and the accuracy of mapping multiple QTLs by using genetic markers are much increased by employing multiple QTL models instead of the single QTL models (and no QTL models) used in interval mapping. However, the computational work involved with multiple QTL models is considerable when the number of QTLs is large. In this paper it is proposed to combine multiple linear regression methods with conventional interval mapping. This is achieved by fitting one QTL at a time in a given interval and simultaneously using (part of) the markers as cofactors to eliminate the effects of additional QTLs. It is shown that the proposed method combines the easy computation of the single QTL interval mapping method with much of the efficiency and accuracy of multiple QTL models.     

On the generalized poisson regression mixture model for mapping quantitative trait loci with count data

Statistical methods for mapping quantitative trait loci (QTL) have been extensively studied. While most existing methods assume normal distribution of the phenotype, the normality assumption could be easily violated when phenotypes are measured in counts. One natural choice to deal with count traits is to apply the classical Poisson regression model. However, conditional on covariates, the Poisson assumption of mean-variance equality may not be valid when data are potentially under- or overdispersed. In this article, we propose an interval-mapping approach for phenotypes measured in counts. We model the effects of QTL through a generalized Poisson regression model and develop efficient likelihood-based inference procedures. This approach, implemented with the EM algorithm, allows for a genomewide scan for the existence of QTL throughout the entire genome. The performance of the proposed method is evaluated through extensive simulation studies along with comparisons with existing approaches such as the Poisson regression and the generalized estimating equation approach. An application to a rice tiller number data set is given. Our approach provides a standard procedure for mapping QTL involved in the genetic control of complex traits measured in counts.     

Mapping and analysis of quantitative trait loci in experimental populations

Simple statistical methods for the study of quantitative trait loci (QTL), such as analysis of variance, have given way to methods that involve several markers and high-resolution genetic maps. As a result, the mapping community has been provided with statistical and computational tools that have much greater power than ever before for studying and locating multiple and interacting QTL. Apart from their immediate practical applications, the lessons learnt from this evolution of QTL methodology might also be generally relevant to other types of functional genomics approach that are aimed at the dissection of complex phenotypes, such as microarray assessment of gene expression.     

Functional mapping of quantitative trait loci underlying the character process: a theoretical framework

Unlike a character measured at a finite set of landmark points, function-valued traits are those that change as a function of some independent and continuous variable. These traits, also called infinite-dimensional characters, can be described as the character process and include a number of biologically, economically, or biomedically important features, such as growth trajectories, allometric scalings, and norms of reaction. Here we present a new statistical infrastructure for mapping quantitative trait loci (QTL) underlying the character process. This strategy, termed functional mapping, integrates mathematical relationships of different traits or variables within the genetic mapping framework. Logistic mapping proposed in this article can be viewed as an example of functional mapping. Logistic mapping is based on a universal biological law that for each and every living organism growth over time follows an exponential growth curve (e.g., logistic or S-shaped). A maximum-likelihood approach based on a logistic-mixture model, implemented with the EM algorithm, is developed to provide the estimates of QTL positions, QTL effects, and other model parameters responsible for growth trajectories. Logistic mapping displays a tremendous potential to increase the power of QTL detection, the precision of parameter estimation, and the resolution of QTL localization due to the small number of parameters to be estimated, the pleiotropic effect of a QTL on growth, and/or residual correlations of growth at different ages. More importantly, logistic mapping allows for testing numerous biologically important hypotheses concerning the genetic basis of quantitative variation, thus gaining an insight into the critical role of development in shaping plant and animal evolution and domestication. The power of logistic mapping is demonstrated by an example of a forest tree, in which one QTL affecting stem growth processes is detected on a linkage group using our method, whereas it cannot be detected using current methods. The advantages of functional mapping are also discussed.     

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

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

Controlling the Type I and Type II Errors in Mapping Quantitative Trait Loci

Although the interval mapping method is widely used for mapping quantitative trait loci (QTLs), it is not very well suited for mapping multiple QTLs. Here, we present the results of a computer simulation to study the application of exact and approximate models for multiple QTLs. In particular, we focus on an automatic two-stage procedure in which in the first stage ``important' markers are selected in multiple regression on markers. In the second stage a QTL is moved along the chromosomes by using the pre-selected markers as cofactors, except for the markers flanking the interval under study. A refined procedure for cases with large numbers of marker cofactors is described. Our approach will be called MQM mapping, where MQM is an acronym for ``multiple-QTL models' as well as for ``marker-QTL-marker.' Our simulation work demonstrates the great advantage of MQM mapping compared to interval mapping in reducing the chance of a type I error (i.e., a QTL is indicated at a location where actually no QTL is present) and in reducing the chance of a type II error (i.e., a QTL is not detected).     

A simple regression method for mapping quantitative trait loci in line crosses using flanking markers

The use of flanking marker methods has proved to be a powerful tool for the mapping of quantitative trait loci (QTL) in the segregating generations derived from crosses between inbred lines. Methods to analyse these data, based on maximum-likelihood, have been developed and provide good estimates of QTL effects in some situations. Maximum-likelihood methods are, however, relatively complex and can be computationally slow. In this paper we develop methods for mapping QTL based on multiple regression which can be applied using any general statistical package. We use the example of mapping in an F(2) population and show that these regression methods produce very similar results to those obtained using maximum likelihood. The relative simplicity of the regression methods means that models with more than a single QTL can be explored and we give examples of two lined loci and of two interacting loci. Other models, for example with more than two QTL, with environmental fixed effects, with between family variance or for threshold traits, could be fitted in a similar way. The ease, speed of application and generality of regression methods for flanking marker analysis, and the good estimates they obtain, suggest that they should provide the method of choice for the analysis of QTL mapping data from inbred line crosses.     

Multiple Trait Analysis of Genetic Mapping for Quantitative Trait Loci

We present in this paper models and statistical methods for performing multiple trait analysis on mapping quantitative trait loci (QTL) based on the composite interval mapping method. By taking into account the correlated structure of multiple traits, this joint analysis has several advantages, compared with separate analyses, for mapping QTL, including the expected improvement on the statistical power of the test for QTL and on the precision of parameter estimation. Also this joint analysis provides formal procedures to test a number of biologically interesting hypotheses concerning the nature of genetic correlations between different traits. Among the testing procedures considered are those for joint mapping, pleiotropy, QTL by environment interaction, and pleiotropy vs. close linkage. The test of pleiotropy (one pleiotropic QTL at a genome position) vs. close linkage (multiple nearby nonpleiotropic QTL) can have important implications for our understanding of the nature of genetic correlations between different traits in certain regions of a genome and also for practical applications in animal and plant breeding because one of the major goals in breeding is to break unfavorable linkage. Results of extensive simulation studies are presented to illustrate various properties of the analyses.     

MCQTL: multi-allelic QTL mapping in multi-cross design

The aim of the MCQTL software package is to perform QTL mapping in multi-cross designs. It allows the analysis of the usual populations derived from inbred lines and can link the families by assuming that the QTL locations are the same in all of them. Moreover, a diallel modelling of the QTL genotypic effects is allowed in multiple related families. The implemented model is a linear regression model. A composite interval mapping and an iterative QTL mapping are implemented to deal with multiple QTL models. Marker cofactor selections by forward or backward stepwise methods are implemented as well as computation of threshold test value by permutation. AVAILABILITY: The program is available on request after signing a licence agreement; free of charge for academic and non-profit organizations at http://www.genoplante.org (Bioinformatics products).     

An algorithm for molecular dissection of tumor progression

The volumetric growth of tumor cells as a function of time is most often likely to be a complex trait, controlled by the combined influences of multiple genes and environmental influences. Genetic mapping has proven to be a powerful tool for detecting and identifying specific genes affecting complex traits, i.e., quantitative trait loci (QTL), based on polymorphic markers. In this article, we present a novel statistical model for genetic mapping of QTL governing tumor growth trajectories in humans. In principle, this model is a combination of functional mapping proposed to map function-valued traits and linkage disequilibrium mapping designed to provide high resolution mapping of QTL by making use of recombination events created at a historic time. We implement an EM-simplex hybrid algorithm for parameter estimation, in which a closed-form solution for the EM algorithm is derived to estimate the population genetic parameters of QTL including the allele frequencies and the coefficient of linkage disequilibrium, and the simplex algorithm incorporated to estimate the curve parameters describing the dynamic changes of cancer cells for different QTL genotypes. Extensive simulations are performed to investigate the statistical properties of our model. Through a number of hypothesis tests, our model allows for cutting-edge studies aimed to decipher the genetic mechanisms underlying cancer growth, development and differentiation. The implications of our model in gene therapy for cancer research are discussed.     

Zero-inflated generalized Poisson regression mixture model for mapping quantitative trait loci underlying count trait with many zeros

Phenotypes measured in counts are commonly observed in nature. Statistical methods for mapping quantitative trait loci (QTL) underlying count traits are documented in the literature. The majority of them assume that the count phenotype follows a Poisson distribution with appropriate techniques being applied to handle data dispersion. When a count trait has a genetic basis, “naturally occurring” zero status also reflects the underlying gene effects. Simply ignoring or miss-handling the zero data may lead to wrong QTL inference. In this article, we propose an interval mapping approach for mapping QTL underlying count phenotypes containing many zeros. The effects of QTLs on the zero-inflated count trait are modelled through the zero-inflated generalized Poisson regression mixture model, which can handle the zero inflation and Poisson dispersion in the same distribution. We implement the approach using the EM algorithm with the Newton-Raphson algorithm embedded in the M-step, and provide a genome-wide scan for testing and estimating the QTL effects. The performance of the proposed method is evaluated through extensive simulation studies. Extensions to composite and multiple interval mapping are discussed. The utility of the developed approach is illustrated through a mouse F₂ intercross data set. Significant QTLs are detected to control mouse cholesterol gallstone formation.     

Bayesian mapping of genomewide interacting quantitative trait loci for ordinal traits

Development of statistical methods and software for mapping interacting QTL has been the focus of much recent research. We previously developed a Bayesian model selection framework, based on the composite model space approach, for mapping multiple epistatic QTL affecting continuous traits. In this study we extend the composite model space approach to complex ordinal traits in experimental crosses. We jointly model main and epistatic effects of QTL and environmental factors on the basis of the ordinal probit model (also called threshold model) that assumes a latent continuous trait underlies the generation of the ordinal phenotypes through a set of unknown thresholds. A data augmentation approach is developed to jointly generate the latent data and the thresholds. The proposed ordinal probit model, combined with the composite model space framework for continuous traits, offers a convenient way for genomewide interacting QTL analysis of ordinal traits. We illustrate the proposed method by detecting new QTL and epistatic effects for an ordinal trait, dead fetuses, in a F(2) intercross of mice. Utility and flexibility of the method are also demonstrated using a simulated data set. Our method has been implemented in the freely available package R/qtlbim, which greatly facilitates the general usage of the Bayesian methodology for genomewide interacting QTL analysis for continuous, binary, and ordinal traits in experimental crosses.     

Complex genetic effects in quantitative trait locus identification: a computationally tractable random model for use in F(2) populations

Methodology for mapping quantitative trait loci (QTL) has focused primarily on treating the QTL as a fixed effect. These methods differ from the usual models of genetic variation that treat genetic effects as random. Computationally expensive methods that allow QTL to be treated as random have been explicitly developed for additive genetic and dominance effects. By extending these methods with a variance component method (VCM), multiple QTL can be mapped. We focused on an F(2) crossbred population derived from inbred lines and estimated effects for each individual and their corresponding marker-derived genetic covariances. We present extensions to pairwise epistatic effects, which are computationally intensive because a great many individual effects must be estimated. But by replacing individual genetic effects with average genetic effects for each marker class, genetic covariances are approximated. This substantially reduces the computational burden by reducing the dimensions of covariance matrices of genetic effects, resulting in a remarkable gain in the speed of estimating the variance components and evaluating the residual log-likelihood. Preliminary results from simulations indicate competitiveness of the reduced model with multiple-interval mapping, regression interval mapping, and VCM with individual genetic effects in its estimated QTL positions and experimental power.     

Robust Bayesian mapping of quantitative trait loci using Student-t distribution for residual

In most quantitative trait loci (QTL) mapping studies, phenotypes are assumed to follow normal distributions. Deviations from this assumption may affect the accuracy of QTL detection, leading to detection of false positive QTL. To improve the robustness of QTL mapping methods, we replace the normal distribution assumption for residuals in a multiple QTL model with a Student-t distribution that is able to accommodate residual outliers. A Robust Bayesian mapping strategy is proposed on the basis of the Bayesian shrinkage analysis for QTL effects. The simulations show that Robust Bayesian mapping approach can substantially increase the power of QTL detection when the normality assumption does not hold and applying it to data already normally distributed does not influence the result. The proposed QTL mapping method is applied to mapping QTL for the traits associated with physics–chemical characters and quality in rice. Similarly to the simulation study in the real data case the robust approach was able to detect additional QTLs when compared to the traditional approach. The program to implement the method is available on request from the first or the corresponding author. Xin Wang and Zhongze Piao contributed equally to this study.