Power of quantitative trait locus mapping for polygenic binary traits using generalized and regression interval mapping in multi-family half-sib designs

A generalized interval mapping (GIM) method to map quantitative trait loci (QTL) for binary polygenic traits in a multi-family half-sib design is developed based on threshold theory and implemented using a Newton-Raphson algorithm. Statistical power and bias of QTL mapping for binary traits by GIM is compared with linear regression interval mapping (RIM) using simulation. Data on 20 paternal half-sib families were simulated with two genetic markers that bracketed an additive QTL. Data simulated and analysed were: (1) data on the underlying normally distributed liability (NDL) scale, (2) binary data created by truncating NDL data based on three thresholds yielding data sets with three different incidences, and (3) NDL data with polygenic and QTL effects reduced by a proportion equal to the ratio of the heritabilities on the binary versus NDL scale (reduced-NDL). Binary data were simulated with and without systematic environmental (herd) effects in an unbalanced design. GIM and RIM gave similar power to detect the QTL and similar estimates of QTL location, effects and variances. Presence of fixed effects caused differences in bias between RIM and GIM, where GIM showed smaller bias which was affected less by incidence. The original NDL data had higher power and lower bias in QTL parameter estimates than binary and reduced-NDL data. RIM for reduced-NDL and binary data gave similar power and estimates of QTL parameters, indicating that the impact of the binary nature of data on QTL analysis is equivalent to its impact on heritability.     

Application of the false discovery rate to quantitative trait loci interval mapping with multiple traits

Controlling the false discovery rate (FDR) has been proposed as an alternative to controlling the genome-wise error rate (GWER) for detecting quantitative trait loci (QTL) in genome scans. The objective here was to implement FDR in the context of regression interval mapping for multiple traits. Data on five traits from an F2 swine breed cross were used. FDR was implemented using tests at every 1 cM (FDR1) and using tests with the highest test statistic for each marker interval (FDRm). For the latter, a method was developed to predict comparison-wise error rates. At low error rates, FDR1 behaved erratically; FDRm was more stable but gave similar significance thresholds and number of QTL detected. At the same error rate, methods to control FDR gave less stringent significance thresholds and more QTL detected than methods to control GWER. Although testing across traits had limited impact on FDR, single-trait testing was recommended because there is no theoretical reason to pool tests across traits for FDR. FDR based on FDRm was recommended for QTL detection in interval mapping because it provides significance tests that are meaningful, yet not overly stringent, such that a more complete picture of QTL is revealed.     

Mapping QTLs and QTL × environment interaction for CIMMYT maize drought stress program using factorial regression and partial least squares methods

The study of QTL × environment interaction (QEI) is important for understanding genotype × environment interaction (GEI) in many quantitative traits. For modeling GEI and QEI, factorial regression (FR) models form a powerful class of models. In FR models, covariables (contrasts) defined on the levels of the genotypic and/or environmental factor(s) are used to describe main effects and interactions. In FR models for QTL expression, considerable numbers of genotypic covariables can occur as for each putative QTL an additional covariable needs to be introduced. For large numbers of genotypic and/or environmental covariables, least square estimation breaks down and partial least squares (PLS) estimation procedures become an attractive alternative. In this paper we develop methodology for analyzing QEI by FR for estimating effects and locations of QTLs and QEI and interpreting QEI in terms of environmental variables. A randomization test for the main effects of QTLs and QEI is presented. A population of F₂ derived F₃ families was evaluated in eight environments differing in drought stress and soil nitrogen content and the traits yield and anthesis silking interval (ASI) were measured. For grain yield, chromosomes 1 and 10 showed significant QEI, whereas in chromosomes 3 and 8 only main effect QTLs were observed. For ASI, QTL main effects were observed on chromosomes 1, 2, 6, 8, and 10, whereas QEI was observed only on chromosome 8. The assessment of the QEI at chromosome 1 for grain yield showed that the QTL main effect explained 35.8% of the QTL + QEI variability, while QEI explained 64.2%. Minimum temperature during flowering time explained 77.6% of the QEI. The QEI analysis at chromosome 10 showed that the QTL main effect explained 59.8% of the QTL + QEI variability, while QEI explained 40.2%. Maximum temperature during flowering time explained 23.8% of the QEI. Results of this study show the possibilities of using FR for mapping QTL and for dissecting QEI in terms of environmental variables. PLS regression is efficient in accounting for background noise produced by other QTLs.     

Multi-environment QTL mixed models for drought stress adaptation in wheat

Many quantitative trait loci (QTL) detection methods ignore QTL-by-environment interaction (QEI) and are limited in accommodation of error and environment-specific variance. This paper outlines a mixed model approach using a recombinant inbred spring wheat population grown in six drought stress trials. Genotype estimates for yield, anthesis date and height were calculated using the best design and spatial effects model for each trial. Parsimonious factor analytic models best captured the variance–covariance structure, including genetic correlations, among environments. The 1RS.1BL rye chromosome translocation (from one parent) which decreased progeny yield by 13.8 g m⁻² was explicitly included in the QTL model. Simple interval mapping (SIM) was used in a genome-wide scan for significant QTL, where QTL effects were fitted as fixed environment-specific effects. All significant environment-specific QTL were subsequently included in a multi-QTL model and evaluated for main and QEI effects with non-significant QEI effects being dropped. QTL effects (either consistent or environment-specific) included eight yield, four anthesis, and six height QTL. One yield QTL co-located (or was linked) to an anthesis QTL, while another co-located with a height QTL. In the final multi-QTL model, only one QTL for yield (6 g m⁻²) was consistent across environments (no QEI), while the remaining QTL had significant QEI effects (average size per environment of 5.1 g m⁻²). Compared to single trial analyses, the described framework allowed explicit modelling and detection of QEI effects and incorporation of additional classification information about genotypes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A mixed-model quantitative trait loci (QTL) analysis for multiple-environment trial data using environmental covariables for QTL-by-environment interactions, with an example in maize

Complex quantitative traits of plants as measured on collections of genotypes across multiple environments are the outcome of processes that depend in intricate ways on genotype and environment simultaneously. For a better understanding of the genetic architecture of such traits as observed across environments, genotype-by-environment interaction should be modeled with statistical models that use explicit information on genotypes and environments. The modeling approach we propose explains genotype-by-environment interaction by differential quantitative trait locus (QTL) expression in relation to environmental variables. We analyzed grain yield and grain moisture for an experimental data set composed of 976 F(5) maize testcross progenies evaluated across 12 environments in the U.S. corn belt during 1994 and 1995. The strategy we used was based on mixed models and started with a phenotypic analysis of multi-environment data, modeling genotype-by-environment interactions and associated genetic correlations between environments, while taking into account intraenvironmental error structures. The phenotypic mixed models were then extended to QTL models via the incorporation of marker information as genotypic covariables. A majority of the detected QTL showed significant QTL-by-environment interactions (QEI). The QEI were further analyzed by including environmental covariates into the mixed model. Most QEI could be understood as differential QTL expression conditional on longitude or year, both consequences of temperature differences during critical stages of the growth.     

Mapping Quantitative Trait Loci for Complex Binary Diseases Using Line Crosses

A composite interval gene mapping procedure for complex binary disease traits is proposed in this paper. The binary trait of interest is assumed to be controlled by an underlying liability that is normally distributed. The liability is treated as a typical quantitative character and thus described by the usual quantitative genetics model. Translation from the liability into a binary (disease) phenotype is through the physiological threshold model. Logistic regression analysis is employed to estimate the effects and locations of putative quantitative trait loci (our terminology for a single quantitative trait locus is QTL while multiple loci are referred to as QTLs). Simulation studies show that properties of this mapping procedure mimic those of the composite interval mapping for normally distributed data. Potential utilization of the QTL mapping procedure for resolving alternative genetic models (e.g., single- or two-trait-locus model) is discussed.     

A comparison of interval mapping procedures including the QTL-cluster mapping

Using the deterministic sampling, patterns of the log-likelihood surfaces expected in some interval mapping procedures for estimating the position of, and the effect for, QTL(s) were investigated for the situations where a single QTL or closely linked QTLs are contained in a chromosome segment bracketed with two markers. The mapping procedures compared were the conventional, likelihood-based interval mapping (IM), the regression interval mapping (RIM), and the QTL-cluster mapping (CM) in which the conditional probabilities of transmission of the whole segment marked by the flanking markers were taken into consideration. The half-sib design was used here, and several cases of the true genetic model were considered, differing in the number of QTLs contained in the marker interval, the linkage phase for the sire, and the magnitude of the QTL(s) effect. For the true genetic models where a single QTL or closely linked QTLs being in coupling phase are contained in the interval, with (R)IM, clear global maxima of the log-likelihood were observed within the range of the marker interval. It was shown that the estimates of the QTL(s) effect at the marked segment level are expected to be unbiased. On the other hand, in a setting where the linkage phase for the linked QTLs located in the interval was different from coupling and repulsion, there was found a ridge along the interval for the log-likelihood surface with (R)IM, indicating the dependency between the estimates of the position of, and the effect for, the putative QTL. For this case, it was found that the position of the putative QTL could be estimated as that of one of the flanking markers, and the estimate of the QTL effect be biased. In contrast, it was revealed that CM is expected to provide the unbiased estimate of the QTL(s)-effect at the segment level for any case of the true genetic models considered here. If the aim is for marker-assisted selection rather than mapping closely linked QTLs individually, then the CM approach is considered to be useful.     

Mapping QTLs with epistatic effects and QTL×environment interactions by mixed linear model approaches

A new methodology based on mixed linear models was developed for mapping QTLs with digenic epistasis and QTL×environment (QE) interactions. Reliable estimates of QTL main effects (additive and epistasis effects) can be obtained by the maximum-likelihood estimation method, while QE interaction effects (additive×environment interaction and epistasis×environment interaction) can be predicted by the-best-linear-unbiased-prediction (BLUP) method. Likelihood ratio and t statistics were combined for testing hypotheses about QTL effects and QE interactions. Monte Carlo simulations were conducted for evaluating the unbiasedness, accuracy, and power for parameter estimation in QTL mapping. The results indicated that the mixed-model approaches could provide unbiased estimates for both positions and effects of QTLs, as well as unbiased predicted values for QE interactions. Additionally, the mixed-model approaches also showed high accuracy and power in mapping QTLs with epistatic effects and QE interactions. Based on the models and the methodology, a computer software program (QTLMapper version 1.0) was developed, which is suitable for interval mapping of QTLs with additive, additive×additive epistasis, and their environment interactions. Received: 23 October 1998 / Accepted: 11 May 1999     

On the use of mathematically-derived traits in QTL mapping

Mathematically-derived traits from two or more component traits, either by addition, subtraction, multiplication, or division, have been frequently used in genetics and breeding. When used in quantitative trait locus (QTL) mapping, derived traits sometimes show discrepancy with QTL identified for the component traits. We used three QTL distributions and three genetic effects models, and an actual maize mapping population, to investigate the efficiency of using derived traits in QTL mapping, and to understand the genetic and biological basis of derived-only QTL, i.e., QTL identified for a derived trait but not for any component trait. Results indicated that the detection power of the four putative QTL was consistently greater than 90% for component traits in simulated populations, each consisting of 200 recombinant inbred lines. Lower detection power and higher false discovery rate (FDR) were observed when derived traits were used. In an actual maize population, simulations were designed based on the observed QTL distributions and effects. When derived traits were used, QTL detected for both component and derived traits had comparable power, but those detected for component traits but not for derived traits had low detection power. The FDR from subtraction and division in the maize population were higher than the FDR from addition and multiplication. The use of derived traits increased the gene number, caused higher-order gene interactions than observed in component traits, and possibly complicated the linkage relationship between QTL as well. The increased complexity of the genetic architecture with derived traits may be responsible for the reduced detection power and the increased FDR. Derived-only QTL identified in practical genetic populations can be explained either as minor QTL that are not significant in QTL mapping of component traits, or as false positives.     

Quantitative trait loci for grain moisture at harvest and field grain drying rate in maize (Zea mays, L.)

Hybrids with low grain moisture (GM) at harvest are specially required in mid- to short-season environments. One of the most important factors determining this trait is field grain drying rate (FDR). To produce hybrids with low GM at harvest, inbred lines can be obtained through selection for either GM or FDR. Thus, a single-cross population (181 F _2:3-generation plants) of two divergent inbred lines was evaluated to locate QTL affecting GM at harvest and FDR as a starting point for marker assisted selection (MAS). Moisture measurements were made with a hand-held moisture meter. Detection of QTL was facilitated with interval mapping in one and two dimensions including an interaction term, and a genetic linkage map of 122 SSR loci covering 1,557.8 cM. The markers were arranged in ten linkage groups. QTL mapping was made for the mean trait performance of the F _2:3 population across years. Ten QTL and an interaction were associated with GM. These QTL accounted for 54.8 and 65.2% of the phenotypic and genotypic variation, respectively. Eight QTL and two interactions were associated with FDR accounting for 35.7 and 45.2% of the phenotypic and genotypic variation, respectively. Two regions were in common between traits. The interaction between QTL for GM at harvest had practical implications for MAS. We conclude that MAS per se will not be an efficient method for reducing GM at harvest and/or increasing FDR. A selection index including both molecular marker information and phenotypic values, each appropriately weighted, would be the best selection strategy.     

Novel and major QTL for branch angle detected by using DH population from an exotic introgression in rapeseed (<Emphasis Type="Italic">Brassica napus</Emphasis> L.)

Key message

A high-density SNP map was constructed and several novel QTL for branch angle across six environments in Brassica napus were identified.

Abstract

Branch angle is a major determinant for the ideotype of a plant, while the mechanisms underlying this trait in Brassica napus remain elusive. Herein, we developed one doubled haploid population from a cross involving one Capsella bursa-pastoris derived B. napus intertribal introgression line with the compressed branches and wooden stems, and constructed a high-density SNP map covering the genetic distance of 2242.14 cM, with an average marker interval of 0.73 cM. After phenotypic measurements across six environments, the inclusive composite interval mapping algorithm was conducted to analyze the QTL associated with branch angle. In single-environment analysis, a total of 17 QTL were detected and mainly distributed on chromosomes A01, A03, A09 and C03. Of these, three major QTL, qBA.A03-2, qBA.C03-3 and qBA.C03-4 were steadily expressed, each explaining more than 10% of the phenotypic variation in at least two environments. Compared with other results on rapeseed branch angle, these major QTL were newly detected. In QTL by environment interactions (QEI) mapping, 10 QTL were identified, and the QTL average effect and QEI effect were estimated. Of these, 7 QTL were detected in both single-environment analysis and QEI mapping. Based on the physical positions of SNPs and the functional annotation of the Arabidopsis thaliana genome, 27 genes within the QTL regions were selected as candidate genes, including early auxin-responsive genes, small auxin-up RNA, auxin/indoleacetic acid and gretchenhagen-3. These results may pave the way for deciphering the genetic control of branch angle in B. napus.

     

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.     

Mapping QTL main and interaction influences on milling quality in elite US rice germplasm

Rice (Oryza sativa L.) head-rice yield (HR) is a key export and domestic quality trait whose genetic control is poorly understood. With the goal of identifying genomic regions influencing HR, quantitative-trait-locus (QTL) mapping was carried out for quality-related traits in recombinant inbred lines (RILs) derived from crosses of common parent Cypress, a high-HR US japonica cultivar, with RT0034, a low-HR indica line (129 RILs) and LaGrue, a low-HR japonica cultivar (298 RILs), grown in two US locations in 2005–2007. Early heading increased HR in the Louisiana (LA) but not the Arkansas (AR) location. Fitting QTL-mapping models to separate QTL main and QTL × environment interaction (QEI) effects and identify epistatic interactions revealed six main-effect HR QTLs in the two crosses, at four of which Cypress contributed the increasing allele. Multi-QTL models accounted for 0.36 of genetic and 0.21 of genetic × environment interaction of HR in MY1, and corresponding proportions of 0.25 and 0.37 in MY2. The greater HR advantage of Cypress in LA than in AR corresponded to a genomewide pattern of opposition of HR-increasing QTL effects by AR-specific effects, suggesting a selection strategy for improving this cultivar for AR. Treating year–location combinations as independent environments resulted in underestimation of QEI effects, evidently owing to lower variation among years within location than between location. Identification of robust HR QTLs in elite long-grain germplasm is suggested to require more detailed attention to the interaction of plant and grain development parameters with environmental conditions than has been given to date.     

Parallel computing in interval mapping of quantitative trait loci.

Linear regression analysis is considered the least computationally demanding method for mapping quantitative trait loci (QTL). However, simultaneous search for multiple QTL, the use of permutations to obtain empirical significance thresholds, and larger experimental studies significantly increase the computational demand. This report describes an easily implemented parallel algorithm, which significantly reduces the computing time in both QTL mapping and permutation testing. In the example provided, the analysis time was decreased to less than 15% of a single processor system by the use of 18 processors. We indicate how the efficiency of the analysis could be improved by distributing the computations more evenly to the processors and how other ways of distributing the data facilitate the use of more processors. The use of parallel computing in QTL mapping makes it possible to routinely use permutations to obtain empirical significance thresholds for multiple traits and multiple QTL models. It could also be of use to improve the computational efficiency of the more computationally demanding QTL analysis methods.     

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.     

Interactions between markers can be caused by the dominance effect of quantitative trait loci

F₂ populations are commonly used in genetic studies of animals and plants. For simplicity, most quantitative trait locus or loci (QTL) mapping methods have been developed on the basis of populations having two distinct genotypes at each polymorphic marker or gene locus. In this study, we demonstrate that dominance can cause the interactions between markers and propose an inclusive linear model that includes marker variables and marker interactions so as to completely control both additive and dominance effects of QTL. The proposed linear model is the theoretical basis for inclusive composite-interval QTL mapping (ICIM) for F₂ populations, which consists of two steps: first, the best regression model is selected by stepwise regression, which approximately identifies markers and marker interactions explaining both additive and dominance variations; second, the interval mapping approach is applied to the phenotypic values adjusted by the regression model selected in the first step. Due to the limited mapping population size, the large number of variables, and multicollinearity between variables, coefficients in the inclusive linear model cannot be accurately determined in the first step. Interval mapping is necessary in the second step to fine tune the QTL to their true positions. The efficiency of including marker interactions in mapping additive and dominance QTL was demonstrated by extensive simulations using three QTL distribution models with two population sizes and an actual rice F₂ population.     

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.     

Quantitative trait locus mapping with background control in genetic populations of clonal F1 and double cross

In this study, we considered five categories of molecular markers in clonal F₁ and double cross populations, based on the number of distinguishable alleles and the number of distinguishable genotypes at the marker locus. Using the completed linkage maps, incomplete and missing markers were imputed as fully informative markers in order to simplify the linkage mapping approaches of quantitative trait genes. Under the condition of fully informative markers, we demonstrated that dominance effect between the female and male parents in clonal F₁ and double cross populations can cause the interactions between markers. We then developed an inclusive linear model that includes marker variables and marker interactions so as to completely control additive effects of the female and male parents, as well as the dominance effect between the female and male parents. The linear model was finally used for background control in inclusive composite interval mapping (ICIM) of quantitative trait locus (QTL). The efficiency of ICIM was demonstrated by extensive simulations and by comparisons with simple interval mapping, multiple‐QTL models and composite interval mapping. Finally, ICIM was applied in one actual double cross population to identify QTL on days to silking in maize.     

Strategies for genetic mapping of categorical traits

The search for efficient and powerful statistical methods and optimal mapping strategies for categorical traits under various experimental designs continues to be one of the main tasks in genetic mapping studies. Methodologies for genetic mapping of categorical traits can generally be classified into two groups, linear and non-linear models. We develop a method based on a threshold model, termed mixture threshold model to handle ordinal (or binary) data from multiple families. Monte Carlo simulations are done to compare its statistical efficiencies and properties of the proposed non-linear model with a linear model for genetic mapping of categorical traits using multiple families. The mixture threshold model has notably higher statistical power than linear models. There may be an optimal sampling strategy (family size vs number of families) in which genetic mapping reaches its maximal power and minimal estimation errors. A single large-sibship family does not necessarily produce the maximal power for detection of quantitative trait loci (QTL) due to genetic sampling of QTL alleles. The QTL allelic model has a marked impact on efficiency of genetic mapping of categorical traits in terms of statistical power and QTL parameter estimation. Compared with a fixed number of QTL alleles (two or four), the model with an infinite number of QTL alleles and normally distributed allelic effects results in loss of statistical power. The results imply that inbred designs (e.g. F₂ or four-way crosses) with a few QTL alleles segregating or reducing number of QTL alleles (e.g. by selection) in outbred populations are desirable in genetic mapping of categorical traits using data from multiple families. This revised version was published online in July 2006 with corrections to the Cover Date.