Varying coefficient models for mapping quantitative trait loci using recombinant inbred intercrosses

There has been a great deal of interest in the development of methodologies to map quantitative trait loci (QTL) using experimental crosses in the last 2 decades. Experimental crosses in animal and plant sciences provide important data sources for mapping QTL through linkage analysis. The Collaborative Cross (CC) is a renewable mouse resource that is generated from eight genetically diverse founder strains to mimic the genetic diversity in humans. The recombinant inbred intercrosses (RIX) generated from CC recombinant inbred (RI) lines share similar genetic structures of F(2) individuals but with up to eight alleles segregating at any one locus. In contrast to F(2) mice, genotypes of RIX can be inferred from the genotypes of their RI parents and can be produced repeatedly. Also, RIX mice typically do not share the same degree of relatedness. This unbalanced genetic relatedness requires careful statistical modeling to avoid false-positive findings. Many quantitative traits are inherently complex with genetic effects varying with other covariates, such as age. For such complex traits, if phenotype data can be collected over a wide range of ages across study subjects, their dynamic genetic patterns can be investigated. Parametric functions, such as sigmoidal or logistic functions, have been used for such purpose. In this article, we propose a flexible nonparametric time-varying coefficient QTL mapping method for RIX data. Our method allows the QTL effects to evolve with time and naturally extends classical parametric QTL mapping methods. We model the varying genetic effects nonparametrically with the B-spline bases. Our model investigates gene-by-time interactions for RIX data in a very flexible nonparametric fashion. Simulation results indicate that the varying coefficient QTL mapping has higher power and mapping precision compared to parametric models when the assumption of constant genetic effects fails. We also apply a modified permutation procedure to control overall significance level.     

Bayesian quantitative trait loci mapping for multiple traits

Most quantitative trait loci (QTL) mapping experiments typically collect phenotypic data on multiple correlated complex traits. However, there is a lack of a comprehensive genomewide mapping strategy for correlated traits in the literature. We develop Bayesian multiple-QTL mapping methods for correlated continuous traits using two multivariate models: one that assumes the same genetic model for all traits, the traditional multivariate model, and the other known as the seemingly unrelated regression (SUR) model that allows different genetic models for different traits. We develop computationally efficient Markov chain Monte Carlo (MCMC) algorithms for performing joint analysis. We conduct extensive simulation studies to assess the performance of the proposed methods and to compare with the conventional single-trait model. Our methods have been implemented in the freely available package R/qtlbim (http://www.qtlbim.org), which greatly facilitates the general usage of the Bayesian methodology for unraveling the genetic architecture of complex traits.     

Bayesian mapping of multiple quantitative trait loci from incomplete inbred line cross data.

A novel fine structure mapping method for quantitative traits is presented. It is based on Bayesian modeling and inference, treating the number of quantitative trait loci (QTLs) as an unobserved random variable and using ideas similar to composite interval mapping to account for the effects of QTLs in other chromosomes. The method is introduced for inbred lines and it can be applied also in situations involving frequent missing genotypes. We propose that two new probabilistic measures be used to summarize the results from the statistical analysis: (1) the (posterior) QTL intensity, for estimating the number of QTLs in a chromosome and for localizing them into some particular chromosomal regions, and (2) the locationwise (posterior) distributions of the phenotypic effects of the QTLs. Both these measures will be viewed as functions of the putative QTL locus, over the marker range in the linkage group. The method is tested and compared with standard interval and composite interval mapping techniques by using simulated backcross progeny data. It is implemented as a software package. Its initial version is freely available for research purposes under the name Multimapper at URL http://www.rni.helsinki.fi/mjs.     

Bayesian LASSO for quantitative trait loci mapping

The mapping of quantitative trait loci (QTL) is to identify molecular markers or genomic loci that influence the variation of complex traits. The problem is complicated by the facts that QTL data usually contain a large number of markers across the entire genome and most of them have little or no effect on the phenotype. In this article, we propose several Bayesian hierarchical models for mapping multiple QTL that simultaneously fit and estimate all possible genetic effects associated with all markers. The proposed models use prior distributions for the genetic effects that are scale mixtures of normal distributions with mean zero and variances distributed to give each effect a high probability of being near zero. We consider two types of priors for the variances, exponential and scaled inverse-chi(2) distributions, which result in a Bayesian version of the popular least absolute shrinkage and selection operator (LASSO) model and the well-known Student's t model, respectively. Unlike most applications where fixed values are preset for hyperparameters in the priors, we treat all hyperparameters as unknowns and estimate them along with other parameters. Markov chain Monte Carlo (MCMC) algorithms are developed to simulate the parameters from the posteriors. The methods are illustrated using well-known barley data.     

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.     

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.     

Bayesian mapping of multiple quantitative trait loci from incomplete outbred offspring data

A general fine-scale Bayesian quantitative trait locus (QTL) mapping method for outcrossing species is presented. It is suitable for an analysis of complete and incomplete data from experimental designs of F2 families or backcrosses. The amount of genotyping of parents and grandparents is optional, as well as the assumption that the QTL alleles in the crossed lines are fixed. Grandparental origin indicators are used, but without forgetting the original genotype or allelic origin information. The method treats the number of QTL in the analyzed chromosome as a random variable and allows some QTL effects from other chromosomes to be taken into account in a composite interval mapping manner. A block-update of ordered genotypes (haplotypes) of the whole family is sampled once in each marker locus during every round of the Markov Chain Monte Carlo algorithm used in the numerical estimation. As a byproduct, the method gives the posterior distributions for linkage phases in the family and therefore it can also be used as a haplotyping algorithm. The Bayesian method is tested and compared with two frequentist methods using simulated data sets, considering two different parental crosses and three different levels of available parental information. The method is implemented as a software package and is freely available under the name Multimapper/outbred at URL http://www.rni.helsinki.fi/mjs/.     

Bayesian mapping of quantitative trait loci for complex binary traits

A complex binary trait is a character that has a dichotomous expression but with a polygenic genetic background. Mapping quantitative trait loci (QTL) for such traits is difficult because of the discrete nature and the reduced variation in the phenotypic distribution. Bayesian statistics are proved to be a powerful tool for solving complicated genetic problems, such as multiple QTL with nonadditive effects, and have been successfully applied to QTL mapping for continuous traits. In this study, we show that Bayesian statistics are particularly useful for mapping QTL for complex binary traits. We model the binary trait under the classical threshold model of quantitative genetics. The Bayesian mapping statistics are developed on the basis of the idea of data augmentation. This treatment allows an easy way to generate the value of a hypothetical underlying variable (called the liability) and a threshold, which in turn allow the use of existing Bayesian statistics. The reversible jump Markov chain Monte Carlo algorithm is used to simulate the posterior samples of all unknowns, including the number of QTL, the locations and effects of identified QTL, genotypes of each individual at both the QTL and markers, and eventually the liability of each individual. The Bayesian mapping ends with an estimation of the joint posterior distribution of the number of QTL and the locations and effects of the identified QTL. Utilities of the method are demonstrated using a simulated outbred full-sib family. A computer program written in FORTRAN language is freely available on request.     

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.     

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.     

Generalized linear mixed models for mapping multiple quantitative trait loci

Many biological traits are discretely distributed in phenotype but continuously distributed in genetics because they are controlled by multiple genes and environmental variants. Due to the quantitative nature of the genetic background, these multiple genes are called quantitative trait loci (QTL). When the QTL effects are treated as random, they can be estimated in a single generalized linear mixed model (GLMM), even if the number of QTL may be larger than the sample size. The GLMM in its original form cannot be applied to QTL mapping for discrete traits if there are missing genotypes. We examined two alternative missing genotype-handling methods: the expectation method and the overdispersion method. Simulation studies show that the two methods are efficient for multiple QTL mapping (MQM) under the GLMM framework. The overdispersion method showed slight advantages over the expectation method in terms of smaller mean-squared errors of the estimated QTL effects. The two methods of GLMM were applied to MQM for the female fertility trait of wheat. Multiple QTL were detected to control the variation of the number of seeded spikelets.     

Empirical Bayesian elastic net for multiple quantitative trait locus mapping

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

Bayesian mapping of quantitative trait loci under complicated mating designs

Quantitative trait loci (QTL) are easily studied in a biallelic system. Such a system requires the cross of two inbred lines presumably fixed for alternative alleles of the QTL. However, development of inbred lines can be time consuming and cost ineffective for species with long generation intervals and severe inbreeding depression. In addition, restriction of the investigation to a biallelic system can sometimes be misleading because many potentially important allelic interactions do not have a chance to express and thus fail to be detected. A complicated mating design involving multiple alleles mimics the actual breeding system. However, it is difficult to develop the statistical model and algorithm using the classical maximum-likelihood method. In this study, we investigate the application of a Bayesian method implemented via the Markov chain Monte Carlo (MCMC) algorithm to QTL mapping under arbitrarily complicated mating designs. We develop the method under a mixed-model framework where the genetic values of founder alleles are treated as random and the nongenetic effects are treated as fixed. With the MCMC algorithm, we first draw the gene flows from the founders to the descendants for each QTL and then draw samples of the genetic parameters. Finally, we are able to simultaneously infer the posterior distribution of the number, the additive and dominance variances, and the chromosomal locations of all identified QTL.     

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

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

New Arabidopsis recombinant inbred line populations genotyped using SNPWave and their use for mapping flowering-time quantitative trait loci

The SNPWave marker system, based on SNPs between the reference accessions Colombia-0 and Landsberg erecta (Ler), was used to distinguish a set of 92 Arabidopsis accessions from various parts of the world. In addition, we used these markers to genotype three new recombinant inbred line populations for Arabidopsis, having Ler as a common parent that was crossed with the accessions Antwerp-1, Kashmir-2, and Kondara. The benefit of using multiple populations that contain many similar markers and the fact that all markers are linked to the physical map of Arabidopsis facilitates the quantitative comparison of maps. Flowering-time variation was analyzed in the three recombinant inbred line populations. Per population, four to eight quantitative trait loci (QTL) were detected. The comparison of the QTL positions related to the physical map allowed the estimate of 12 different QTL segregating for flowering time for which Ler has an allele different from one, two, or three of the other accessions.     

Model for mapping imprinted quantitative trait loci in an inbred F2 design

The role of imprinting in shaping development has been ubiquitously observed in plants, animals, and humans. However, a statistical method that can detect and estimate the effects of imprinted quantitative trait loci (iQTL) over the genome has not been extensively developed. In this article, we propose a maximum likelihood approach for testing and estimating the imprinted effects of iQTL that contribute to variation in a quantitative trait. This approach, implemented with the EM algorithm, allows for a genome-wide scan for the existence of iQTL. This approach was used to reanalyze published data in an F(2) family derived from the LG/S and SM/S mouse strains. Several iQTL that regulate the growth of body weight by expressing paternally inherited alleles were identified. Our approach provides a standard procedure for testing the statistical significance of iQTL involved in the genetic control of complex traits.     

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.     

High-resolution mapping of quantitative trait loci by selective recombinant genotyping

Selective recombinant genotyping (SRG) is a three-stage procedure for high-resolution mapping of a QTL that has previously been mapped to a known confidence interval (target C.I.). In stage 1, a large mapping population is accessed and phenotyped, and a proportion, P, of the high and low tails is selected. In stage 2, the selected individuals are genotyped for a pair of markers flanking the target C.I., and a group of R individuals carrying recombinant chromosomes in the target interval are identified. In stage 3, the recombinant individuals are genotyped for a set of M markers spanning the target C.I. Extensive simulations showed that: (1) Standard error of QTL location (SEQTL) decreased when QTL effect (d) or population size (N) increased, but was constant for given "power factor" (PF = d(2)N); (2) increasing the proportion selected in the tails beyond 0.25 had only a negligible effect on SEQTL; and (3) marker spacing in the target interval had a remarkably powerful effect on SEQTL, yielding a reduction of up to 10-fold in going from highest (24 cM) to lowest (0.29 cM) spacing at given population size and QTL effect. At the densest marker spacing, SEQTL of 1.0-0.06 cM were obtained at PF = 500-16,000. Two new genotyping procedures, the half-section algorithm and the golden section/half-section algorithm, allow the equivalent of complete haplotyping of the target C.I. in the recombinant individuals to be achieved with many fewer data points than would be required by complete individual genotyping.     

A general mixture model approach for mapping quantitative trait loci from diverse cross designs involving multiple inbred lines

Most current statistical methods developed for mapping quantitative trait loci (QTL) based on inbred line designs apply to crosses from two inbred lines. Analysis of QTL in these crosses is restricted by the parental genetic differences between lines. Crosses from multiple inbred lines or multiple families are common in plant and animal breeding programmes, and can be used to increase the efficiency of a QTL mapping study. A general statistical method using mixture model procedures and the EM algorithm is developed for mapping QTL from various cross designs of multiple inbred lines. The general procedure features three cross design matrices, W, that define the contribution of parental lines to a particular cross and a genetic design matrix, D, that specifies the genetic model used in multiple line crosses. By appropriately specifying W matrices, the statistical method can be applied to various cross designs, such as diallel, factorial, cyclic, parallel or arbitrary-pattern cross designs with two or multiple parental lines. Also, with appropriate specification for the D matrix, the method can be used to analyse different kinds of cross populations, such as F2 backcross, four-way cross and mixed crosses (e.g. combining backcross and F2). Simulation studies were conducted to explore the properties of the method, and confirmed its applicability to diverse experimental designs.     

Mapping quantitative trait loci using the experimental designs of recombinant inbred populations

In the data collection of the QTL experiments using recombinant inbred (RI) populations, when individuals are genotyped for markers in a population, the trait values (phenotypes) can be obtained from the genotyped individuals (from the same population) or from some progeny of the genotyped individuals (from the different populations). Let Fu be the genotyped population and Fv (v>or=u) be the phenotyped population. The experimental designs that both marker genotypes and phenotypes are recorded on the same populations can be denoted as (Fu/Fv, u=v) designs and that genotypes and phenotypes are obtained from the different populations can be denoted as (Fu/Fv, v>u) designs. Although most of the QTL mapping experiments have been conducted on the backcross and F2(F2/F2) designs, the other (Fu/Fv, v>or=u) designs are also very popular. The great benefits of using the other (Fu/Fv, v>or=u) designs in QTL mapping include reducing cost and environmental variance by phenotyping several progeny for the genotyped individuals and taking advantages of the changes in population structures of other RI populations. Current QTL mapping methods including those for the (Fu/Fv, u=v) designs, mostly for the backcross or F2/F2 design, and for the F2/F3 design based on a one-QTL model are inadequate for the investigation of the mapping properties in the (Fu/Fv, uor=u) designs. In addition, the QTL mapping properties of the proposed and approximate methods in different designs are discussed. Simulations were performed to evaluate the performance of the proposed and approximate methods. The proposed method is proven to be able to correct the problems of the approximate and current methods for improving the resolution of genetic architecture of quantitative traits and can serve as an effective tool to explore the QTL mapping study in the system of RI populations.