Generalized linear model for interval mapping of quantitative trait loci

We developed a generalized linear model of QTL mapping for discrete traits in line crossing experiments. Parameter estimation was achieved using two different algorithms, a mixture model-based EM (expectation–maximization) algorithm and a GEE (generalized estimating equation) algorithm under a heterogeneous residual variance model. The methods were developed using ordinal data, binary data, binomial data and Poisson data as examples. Applications of the methods to simulated as well as real data are presented. The two different algorithms were compared in the data analyses. In most situations, the two algorithms were indistinguishable, but when large QTL are located in large marker intervals, the mixture model-based EM algorithm can fail to converge to the correct solutions. Both algorithms were coded in C++ and interfaced with SAS as a user-defined SAS procedure called PROC QTL.     

Mapping quantitative trait loci for binary traits using a heterogeneous residual variance model: an application to Marek's disease susceptibility in chickens

A typical problem in mapping quantitative trait loci (QTLs) comes from missing QTL genotype. A routine method for parameter estimation involving missing data is the mixture model maximum likelihood method. We developed an alternative QTL mapping method that describes a mixture of several distributions by a single model with a heterogeneous residual variance. The two methods produce similar results, but the heterogeneous residual variance method is computationally much faster than the mixture model approach. In addition, the new method can automatically generate sampling variances of the estimated parameters. We derive the new method in the context of QTL mapping for binary traits in a F2 population. Using the heterogeneous residual variance model, we identified a QTL on chromosome IV that controls Marek's disease susceptibility in chickens. The QTL alone explains 7.2% of the total disease variation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mapping quantitative trait loci in tetraploid populations

Knowledge of quantitative trait locus (QTL) mapping in polyploids is almost void, albeit many exquisite strategies of QTL mapping have been proposed and extensive investigations have been carried out in diploid animals and plants. In this paper we develop a simple algorithm which uses an iteratively reweighted least square method to map QTLs in tetraploid populations. The method uses information from all markers in a linkage group to infer the probability distribution of QTL genotype under the assumption of random chromosome segregation. Unlike QTL mapping in diploid species, here we estimate and test the compound 'gametic effect', which consists of the composite 'genic effect' of alleles and higher-order gene interactions. The validity and efficiency of the proposed method are investigated through simulation studies. Results show that the method can successfully locate QTLs and separates different sources (e.g. additive and dominance) of variance components contributed by the QTLs.     

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     

QTL mapping in outbred half-sib families using Bayesian model selection

In this article, we propose a model selection method, the Bayesian composite model space approach, to map quantitative trait loci (QTL) in a half-sib population for continuous and binary traits. In our method, the identity-by-descent-based variance component model is used. To demonstrate the performance of this model, the method was applied to map QTL underlying production traits on BTA6 in a Chinese half-sib dairy cattle population. A total of four QTLs were detected, whereas only one QTL was identified using the traditional least square (LS) method. We also conducted two simulation experiments to validate the efficiency of our method. The results suggest that the proposed method based on a multiple-QTL model is efficient in mapping multiple QTL for an outbred half-sib population and is more powerful than the LS method based on a single-QTL model.     

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.     

A mixed model QTL analysis for sugarcane multiple-harvest-location trial data

Sugarcane-breeding programs take at least 12 years to develop new commercial cultivars. Molecular markers offer a possibility to study the genetic architecture of quantitative traits in sugarcane, and they may be used in marker-assisted selection to speed up artificial selection. Although the performance of sugarcane progenies in breeding programs are commonly evaluated across a range of locations and harvest years, many of the QTL detection methods ignore two- and three-way interactions between QTL, harvest, and location. In this work, a strategy for QTL detection in multi-harvest-location trial data, based on interval mapping and mixed models, is proposed and applied to map QTL effects on a segregating progeny from a biparental cross of pre-commercial Brazilian cultivars, evaluated at two locations and three consecutive harvest years for cane yield (tonnes per hectare), sugar yield (tonnes per hectare), fiber percent, and sucrose content. In the mixed model, we have included appropriate (co)variance structures for modeling heterogeneity and correlation of genetic effects and non-genetic residual effects. Forty-six QTLs were found: 13 QTLs for cane yield, 14 for sugar yield, 11 for fiber percent, and 8 for sucrose content. In addition, QTL by harvest, QTL by location, and QTL by harvest by location interaction effects were significant for all evaluated traits (30 QTLs showed some interaction, and 16 none). Our results contribute to a better understanding of the genetic architecture of complex traits related to biomass production and sucrose content in sugarcane.     

Joint mapping of quantitative trait Loci for multiple binary characters

Joint mapping for multiple quantitative traits has shed new light on genetic mapping by pinpointing pleiotropic effects and close linkage. Joint mapping also can improve statistical power of QTL detection. However, such a joint mapping procedure has not been available for discrete traits. Most disease resistance traits are measured as one or more discrete characters. These discrete characters are often correlated. Joint mapping for multiple binary disease traits may provide an opportunity to explore pleiotropic effects and increase the statistical power of detecting disease loci. We develop a maximum-likelihood method for mapping multiple binary traits. We postulate a set of multivariate normal disease liabilities, each contributing to the phenotypic variance of one disease trait. The underlying liabilities are linked to the binary phenotypes through some underlying thresholds. The new method actually maps loci for the variation of multivariate normal liabilities. As a result, we are able to take advantage of existing methods of joint mapping for quantitative traits. We treat the multivariate liabilities as missing values so that an expectation-maximization (EM) algorithm can be applied here. We also extend the method to joint mapping for both discrete and continuous traits. Efficiency of the method is demonstrated using simulated data. We also apply the new method to a set of real data and detect several loci responsible for blast resistance in rice.     

Mapping of QTLs governing agronomic and yield traits in chickpea

Chickpea is one of the most important leguminous cool season food crops, cultivated prevalently in South Asia and Middle East. The main objective of this study was to identify quantitative trait loci (QTLs) associated with seven agronomic and yield traits in two recombinant inbred line populations of chickpea derived from the crosses JG62 × Vijay (JV population) and Vijay × ICC4958 (VI population) from at least three environments. Single locus QTL analysis involved composite interval mapping (CIM) for individual traits and multiple-trait composite interval mapping (MCIM) for correlated traits to detect pleiotropic QTLs. Two-locus analysis was conducted to identify the main effect QTLs (M-QTLs), epistatic QTLs (E-QTLs) and QTL × environment interactions. Through CIM analysis, a total of 106 significant QTLs (41 in JV and 65 in VI populations) were identified for the seven traits, of which one QTL each for plant height and days to maturity was common in both the populations. Six pleiotropic QTLs that were consistent over the environments were also identified. LG2 in JV and LG1a in VI contained at least one QTL for each trait. Hence, concentrating on these LGs in molecular breeding programs is most likely to bring simultaneous improvement in these traits.     

Mitochondrial genetic effects on latent class variables associated with susceptibility to alcoholism

We report the results of statistical genetic analyses of data from the Collaborative Study on the Genetics of Alcoholism prepared for the Genetic Analysis Workshop 14 to detect and characterize maternally inherited mitochondrial genetic effects on variation in latent class psychiatric/behavioral variables employed in the diagnosis of alcoholism. Using published extensions to variance decomposition methods for statistical genetic analysis of continuous and discrete traits we: 1) estimated the proportion of the variance in each trait due to the effects of mitochondrial DNA (mtDNA), 2) tested for pleiotropy, both mitochondrial genetic and residual additive genetic, between trait pairs, and 3) evaluated whether the simultaneous estimation of mitochondrial genetic effects on these traits improves our ability to detect and localize quantitative trait loci (QTL) in the nuclear genome. After correction for multiple testing, we find significant (p < 0.009) mitochondrial genetic contributions to the variance for two latent class variables. Although we do detect significant residual additive genetic correlations between the two traits, there is no evidence of a residual mitochondrial genetic correlation between them. Evidence for autosomal QTL for these traits is improved when linkage screens are conditioned on significant mitochondrial genetic effects. We conclude that mitochondrial genes may contribute to variation in some latent class psychiatric/behavioral variables associated with alcoholism.     

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.     

Study on mapping Quantitative Trait Loci for animal complex binary traits using Bayesian-Markov chain Monte Carlo approach

It is a challenging issue to map Quantitative Trait Loci (QTL) underlying complex discrete traits,which usually show discontinuous distribution and less information,using conventional statisti-cal methods. Bayesian-Markov chain Monte Carlo (Bayesian-MCMC) approach is the key procedure in mapping QTL for complex binary traits,which provides a complete posterior distribution for QTL parameters using all prior information. As a consequence,Bayesian estimates of all interested vari-ables can be obtained straightforwardly basing on their posterior samples simulated by the MCMC algorithm. In our study,utilities of Bayesian-MCMC are demonstrated using simulated several ani-mal outbred full-sib families with different family structures for a complex binary trait underlied by both a QTL and polygene. Under the Identity-by-Descent-Based variance component random model,three samplers basing on MCMC,including Gibbs sampling,Metropolis algorithm and reversible jump MCMC,were implemented to generate the joint posterior distribution of all unknowns so that the QTL parameters were obtained by Bayesian statistical inferring. The results showed that Bayesian-MCMC approach could work well and robust under different family structures and QTL effects. As family size increases and the number of family decreases,the accuracy of the parameter estimates will be im-proved. When the true QTL has a small effect,using outbred population experiment design with large family size is the optimal mapping strategy.     

Study on mapping Quantitative Trait Loci for animal complex binary traits using Bayesian-Markov chain Monte Carlo approach

It is a challenging issue to map Quantitative Trait Loci (QTL) underlying complex discrete traits, which usually show discontinuous distribution and less information, using conventional statistical methods. Bayesian-Markov chain Monte Carlo (Bayesian-MCMC) approach is the key procedure in mapping QTL for complex binary traits, which provides a complete posterior distribution for QTL parameters using all prior information. As a consequence, Bayesian estimates of all interested variables can be obtained straightforwardly basing on their posterior samples simulated by the MCMC algorithm. In our study, utilities of Bayesian-MCMC are demonstrated using simulated several animal outbred full-sib families with different family structures for a complex binary trait underlied by both a QTL and polygene. Under the Identity-by-Descent-Based variance component random model, three samplers basing on MCMC, including Gibbs sampling, Metropolis algorithm and reversible jump MCMC, were implemented to generate the joint posterior distribution of all unknowns so that the QTL parameters were obtained by Bayesian statistical inferring. The results showed that Bayesian-MCMC approach could work well and robust under different family structures and QTL effects. As family size increases and the number of family decreases, the accuracy of the parameter estimates will be improved. When the true QTL has a small effect, using outbred population experiment design with large family size is the optimal mapping strategy.     

Multi-trait and multi-environment QTL analyses of yield and a set of physiological traits in pepper

Key message

A mixed model framework was defined for QTL analysis of multiple traits across multiple environments for a RIL population in pepper. Detection power for QTLs increased considerably and detailed study of QTL by environment interactions and pleiotropy was facilitated.

Abstract

For many agronomic crops, yield is measured simultaneously with other traits across multiple environments. The study of yield can benefit from joint analysis with other traits and relations between yield and other traits can be exploited to develop indirect selection strategies. We compare the performance of three multi-response QTL approaches based on mixed models: a multi-trait approach (MT), a multi-environment approach (ME), and a multi-trait multi-environment approach (MTME). The data come from a multi-environment experiment in pepper, for which 15 traits were measured in four environments. The approaches were compared in terms of number of QTLs detected for each trait, the explained variance, and the accuracy of prediction for the final QTL model. For the four environments together, the superior MTME approach delivered a total of 47 regions containing putative QTLs. Many of these QTLs were pleiotropic and showed quantitative QTL by environment interaction. MTME was superior to ME and MT in the number of QTLs, the explained variance and accuracy of predictions. The large number of model parameters in the MTME approach was challenging and we propose several guidelines to help obtain a stable final QTL model. The results confirmed the feasibility and strengths of novel mixed model QTL methodology to study the architecture of complex traits.     

Natural variation for carbohydrate content in Arabidopsis. Interaction with complex traits dissected by quantitative genetics

Besides being a metabolic fuel, carbohydrates play important roles in plant growth and development, in stress responses, and as signal molecules. We exploited natural variation in Arabidopsis (Arabidopsis thaliana) to decipher the genetic architecture determining carbohydrate content. A quantitative trait locus (QTL) approach in the Bay-0 x Shahdara progeny grown in two contrasting nitrogen environments led to the identification of 39 QTLs for starch, glucose, fructose, and sucrose contents representing at least 14 distinct polymorphic loci. A major QTL for fructose content (FR3.4) and a QTL for starch content (ST3.4) were confirmed in heterogeneous inbred families. Several genes associated with carbon (C) metabolism colocalize with the identified QTL. QTLs for senescence-related traits, and for flowering time, water status, and nitrogen-related traits, previously detected with the same genetic material, colocalize with C-related QTLs. These colocalizations reflect the complex interactions of C metabolism with other physiological processes. QTL fine-mapping and cloning could thus lead soon to the identification of genes potentially involved in the control of different connected physiological processes.     

Log-normal variation belts for growth curves

Prediction (confidence) or tolerance belts compound the uncertainty of sample estimates with the estimated extent of individual variation. The latter is therefore better described by variation belts, in which sample estimates are simply substituted for population parameters. Variation belts can provide valuable graphical indications concerning the goodness of fit of postulated error models. While multiplicative least-squares (MLS) methods appear appropriate in principle for biological growth, they are unsatisfactory in practice when logarithmically transformed data are heteroscedastic. Heteroscedastic multiplicative error models can be fitted by iteratively reweighted multiplicative least squares (IRMLS), but unacceptable negative or infinite residual variance estimates and unreasonably wide variation belts are occasionally obtained. These difficulties can be prevented by constrained iteratively reweighted multiplicative least squares (CIRMLS). Examples are presented concerning the metabolic allometry of white rats, the somatic growth of male elephant seals, and the growth of an experimental population of Paramecium caudatum.     

QTL mapping for yield and yield contributing traits in two mapping populations of bread wheat

In bread wheat, single-locus and two-locus QTL analyses were conducted for seven yield and yield contributing traits using two different mapping populations (P I and P II). Single-locus QTL analyses involved composite interval mapping (CIM) for individual traits and multiple-trait composite interval mapping (MCIM) for correlated yield traits to detect the pleiotropic QTLs. Two-locus analyses were conducted to detect main effect QTLs (M-QTLs), epistatic QTLs (E-QTLs) and QTL × environment interactions (QE and QQE). Only a solitary QTL for spikelets per spike was common between the above two populations. HomoeoQTLs were also detected, suggesting the presence of triplicate QTLs in bread wheat. Relatively fewer QTLs were detected in P I than in P II. This may be partly due to low density of marker loci on P I framework map (173) than in P II (521) and partly due to more divergent parents used for developing P II. Six QTLs were important which were pleiotropic/coincident involving more than one trait and were also consistent over environments. These QTLs could be utilized efficiently for marker assisted selection (MAS).     

Distribution of genomic regions differentiating oak species assessed by QTL detection

Pedunculate oak and sessile oak are two sympatric interfertile species that exhibit leaf morphological differences. We aimed to detect quantitative trait loci (QTLs) of these traits in order to locate genomic regions involved in species differentiation. A total of 15 leaf morphological traits were assessed in a mixed forest stand composed of Quercus petraea and Q. robur and in a full-sib pedigree of Q. robur. The progeny of the full-sib family were vegetatively propagated in two successive experiments comprising 174 and 216 sibs, and assessments were made on two leaves collected on each of the 1080 and 1530 cuttings corresponding to the two experiments. Traits that exhibited strong species differences in the mixed stand tended also to have higher repeatability values in the mapping population, thus indicating higher genetic control. A genetic map was constructed for QTL detection. Composite interval mapping with the one QTL model was used for QTL detection. From one to three QTLs were detected for 13 traits. In-depth analysis of the QTLs, controlling the five morphological traits that exhibited the highest interspecific differences in the mixed stand, indicated that they were distributed on six linkage groups, with two clusters comprising QTLs of at least two discriminant traits. These results were reinforced when error 1 for QTL detection was set at 5% at the chromosome level, as up to nine clusters could be identified. In conclusion, traits involved in interspecific differentiation of oaks are under polygenic control and widespread in clusters across the genome.     

QTLs for days to silking in a recombinant inbred line maize population subjected to high and low nitrogen regimes

Days to silking (DTS) is one of the most important traits in maize (Zea mays). To investigate its genetic basis, a recombinant inbred line population was subjected to high and low nitrogen (N) regimes to detect quantitative trait loci (QTLs) associated with DTS. Three QTLs were identified under the high N regime; these explained 25.4% of the phenotypic variance. Due to additive effects, the QTL on chromosome 6 increased DTS up to 0.66 days; while the other two QTLs mapped on chromosome 9 (one linked with Phi061 and the other linked with Nc134) decreased DTS 0.89 and 0.91 days, respectively. Under low N regime, two QTLs were mapped on chromosomes 6 and 9, which accounted for 25.9% of the phenotypic variance. Owing to additive effects, the QTL on chromosome 6 increased DTS 0.67 days, while the other QTL on chromosome 9 decreased it 1.48 days. The QTL on chromosome 6, flanked by microsatellite markers Bnlg1600 and Phi077, was detected under both N regimes. In conclusion, we identified four QTLs, one on chromosome 6 and three on chromosome 9. These results contribute to our understanding of the genetic basis of DTS and will be useful for developing marker-assisted selection in maize breeding programs.     

A New Method to Infer Causal Phenotype Networks Using QTL and Phenotypic Information

In the context of genetics and breeding research on multiple phenotypic traits, reconstructing the directional or causal structure between phenotypic traits is a prerequisite for quantifying the effects of genetic interventions on the traits. Current approaches mainly exploit the genetic effects at quantitative trait loci (QTLs) to learn about causal relationships among phenotypic traits. A requirement for using these approaches is that at least one unique QTL has been identified for each trait studied. However, in practice, especially for molecular phenotypes such as metabolites, this prerequisite is often not met due to limited sample sizes, high noise levels and small QTL effects. Here, we present a novel heuristic search algorithm called the QTL+phenotype supervised orientation (QPSO) algorithm to infer causal directions for edges in undirected phenotype networks. The two main advantages of this algorithm are: first, it does not require QTLs for each and every trait; second, it takes into account associated phenotypic interactions in addition to detected QTLs when orienting undirected edges between traits. We evaluate and compare the performance of QPSO with another state-of-the-art approach, the QTL-directed dependency graph (QDG) algorithm. Simulation results show that our method has broader applicability and leads to more accurate overall orientations. We also illustrate our method with a real-life example involving 24 metabolites and a few major QTLs measured on an association panel of 93 tomato cultivars. Matlab source code implementing the proposed algorithm is freely available upon request.