Using molecular markers to map multiple quantitative trait loci: models for backcross,recombinant inbred,and doubled haploid progeny

Summary To maximize parameter estimation efficiency and statistical power and to estimate epistasis, the parameters of multiple quantitative trait loci (QTLs) must be simultaneously estimated. If multiple QTL affect a trait, then estimates of means of QTL genotypes from individual locus models are statistically biased. In this paper, I describe methods for estimating means of QTL genotypes and recombination frequencies between marker and quantitative trait loci using multilocus backcross, doubled haploid, recombinant inbred, and testcross progeny models. Expected values of marker genotype means were defined using no double or multiple crossover frequencies and flanking markers for linked and unlinked quantitative trait loci. The expected values for a particular model comprise a system of nonlinear equations that can be solved using an interative algorithm, e.g., the Gauss-Newton algorithm. The solutions are maximum likelihood estimates when the errors are normally distributed. A linear model for estimating the parameters of unlinked quantitative trait loci was found by transforming the nonlinear model. Recombination frequency estimators were defined using this linear model. Certain means of linked QTLs are less efficiently estimated than means of unlinked QTLs.     

Mapping-Linked Quantitative Trait Loci Using Bayesian Analysis and Markov Chain Monte Carlo Algorithms

A Bayesian method for mapping linked quantitative trait loci (QTL) using multiple linked genetic markers is presented. Parameter estimation and hypothesis testing was implemented via Markov chain Monte Carlo (MCMC) algorithms. Parameters included were allele frequencies and substitution effects for two biallelic QTL, map positions of the QTL and markers, allele frequencies of the markers, and polygenic and residual variances. Missing data were polygenic effects and multi-locus marker-QTL genotypes. Three different MCMC schemes for testing the presence of a single or two linked QTL on the chromosome were compared. The first approach includes a model indicator variable representing two unlinked QTL affecting the trait, one linked and one unlinked QTL, or both QTL linked with the markers. The second approach incorporates an indicator variable for each QTL into the model for phenotype, allowing or not allowing for a substitution effect of a QTL on phenotype, and the third approach is based on model determination by reversible jump MCMC. Methods were evaluated empirically by analyzing simulated granddaughter designs. All methods identified correctly a second, linked QTL and did not reject the one-QTL model when there was only a single QTL and no additional or an unlinked QTL.     

A Random Model Approach to Interval Mapping of Quantitative Trait Loci

Mapping quantitative trait loci in outbred populations is important because many populations of organisms are noninbred. Unfortunately, information about the genetic architecture of the trait may not be available in outbred populations. Thus, the allelic effects of genes can not be estimated with ease. In addition, under linkage equilibrium, marker genotypes provide no information about the genotype of a QTL (our terminology for a single quantitative trait locus is QTL while multiple loci are referred to as QTLs). To circumvent this problem, an interval mapping procedure based on a random model approach is described. Under a random model, instead of estimating the effects, segregating variances of QTLs are estimated by a maximum likelihood method. Estimation of the variance component of a QTL depends on the proportion of genes identical-by-descent (IBD) shared by relatives at the locus, which is predicted by the IBD of two markers flanking the QTL. The marker IBD shared by two relatives are inferred from the observed marker genotypes. The procedure offers an advantage over the regression interval mapping in terms of high power and small estimation errors and provides flexibility for large sibships, irregular pedigree relationships and incorporation of common environmental and fixed effects.     

Genetic response from marker assisted selection in an outbred population for differing marker bracket sizes and with two identified quantitative trait loci.

Effect of flanking quantitative trait loci (QTL)-marker bracket size on genetic response to marker assisted selection in an outbred population was studied by simulation of a nucleus breeding scheme. In addition, genetic response with marker assisted selection (MAS) from two quantitative trait loci on the same and different chromosome(s) was investigated. QTL that explained either 5% or 10% of phenotypic variance were simulated. A polygenic component was simulated in addition to the quantitative trait loci. In total, 35% of the phenotypic variance was due to genetic factors. The trait was measured on females only. Having smaller marker brackets flanking the QTL increased the genetic response from MAS selection. This was due to the greater ability to trace the QTL transmission from one generation to the next with the smaller flanking QTL-marker bracket, which increased the accuracy of estimation of the QTL allelic effects. Greater negative covariance between effects at both QTL was observed when two QTL were located on the same chromosome compared to different chromosomes. Genetic response with MAS was greater when the QTL were on the same chromosome in the early generations and greater when they were on different chromosomes in the later generations of MAS.     

Detection of quantitative trait loci in outbred populations with incomplete marker data.

Augmentation of marker genotypes for ungenotyped individuals is implemented in a Bayesian approach via the use of Markov chain Monte Carlo techniques. Marker data on relatives and phenotypes are combined to compute conditional posterior probabilities for marker genotypes of ungenotyped individuals. The presented procedure allows the analysis of complex pedigrees with ungenotyped individuals to detect segregating quantitative trait loci (QTL). Allelic effects at the QTL were assumed to follow a normal distribution with a covariance matrix based on known QTL position and identity by descent probabilities derived from flanking markers. The Bayesian approach estimates variance due to the single QTL, together with polygenic and residual variance. The method was empirically tested through analyzing simulated data from a complex granddaughter design. Ungenotyped dams were related to one or more sons or grandsires in the design. Heterozygosity of the marker loci and size of QTL were varied. Simulation results indicated a significant increase in power when ungenotyped dams were included in the analysis.     

Interval mapping of quantitative trait loci for reproductive,morphological, and seed traits of soybean (Glycine max L.)

Quantitative trait loci (QTL) were mapped in segregating progeny from a cross between two soybean (Glycine max (L.) Merr.) cultivars: Minsoy (PI 27.890) and Noir 1 (PI 290.136). The 15 traits analyzed included reproductive, morphological, and seed traits, seed yield and carbon isotope discrimination ratios (¹³C/¹²C). Genetic variation was detected for all of the traits, and transgressive segregation was a common phenomenon. One hundred and thirty-two linked genetic markers and 24 additional unlinked markers were used to locate QTL by interval mapping and one-way analysis of variance, respectively. Quantitative trait loci controlling 11 of the 15 traits studied were localized to intervals in 6 linkage groups. Quantitative trait loci for developmental and morphological traits (R1, R5, R8, plant height, canopy height, leaf area, etc.) tended to be clustered in three intervals, two of which were also associated with seed yield. Quantitative trait loci for seed oil were separated from all the other QTL. Major QTL for maturity and plant height were linked to RFLP markers R79 (31% variation) and G173 (53% variation). Quantitative trait loci associated with unlinked markers included possible loci for seed protein and weight. Linkage between QTL is discussed in relation to the heritability and genetic correlation of the traits.     

Selection mapping of loci for quantitative disease resistance in a diverse maize population

The selection response of a complex maize population improved primarily for quantitative disease resistance to northern leaf blight (NLB) and secondarily for common rust resistance and agronomic phenotypes was investigated at the molecular genetic level. A tiered marker analysis with 151 simple sequence repeat (SSR) markers in 90 individuals of the population indicated that on average six alleles per locus were available for selection. An improved test statistic for selection mapping was developed, in which quantitative trait loci (QTL) are identified through the analysis of allele-frequency shifts at mapped multiallelic loci over generations of selection. After correcting for the multiple tests performed, 25 SSR loci showed evidence of selection. Many of the putatively selected loci were unlinked and dispersed across the genome, which was consistent with the diffuse distribution of previously published QTL for NLB resistance. Compelling evidence for selection was found on maize chromosome 8, where several putatively selected loci colocalized with published NLB QTL and a race-specific resistance gene. Analysis of F(2) populations derived from the selection mapping population suggested that multiple linked loci in this chromosomal segment were, in part, responsible for the selection response for quantitative resistance to NLB.     

Further insights of the variance component method for detecting QTL in livestock and aquacultural species: relaxing the assumption of additive effects

Complex traits may show some degree of dominance at the gene level that may influence the statistical power of simple models, i.e. assuming only additive effects to detect quantitative trait loci (QTL) using the variance component method. Little has been published on this topic even in species where relatively large family sizes can be obtained, such as poultry, pigs, and aquacultural species. This is important, when the idea is to select regions likely to be harbouring dominant QTL or in marker assisted selection. In this work, we investigated the empirical power and accuracy to both detect and localise dominant QTL with or without incorporating dominance effects explicitly in the model of analysis. For this purpose, populations with variable family sizes and constant population size and different values for dominance variance were simulated. The results show that when using only additive effects there was little loss in power to detect QTL and estimates of position, using or not using dominance, were empirically unbiased. Further, there was little gain in accuracy of positioning the QTL with most scenarios except when simulating an overdominant QTL.     

Defining the assumptions underlying modeling of epistatic QTL using variance component methods

Variance component models are commonly used to detect quantitative trait loci (QTL) in general pedigrees. The variance-covariance structure of the random QTL effect is given by the identity by descent (IBD) between genotypes. Epistatic effects have previously been modeled, both for unlinked and linked loci, as a random effect with a variance-covariance structure given by the Hadamard product between the IBD matrices of the direct QTL effects. In the original papers, the model was given but not derived. Here, we identify the underlying assumptions of this previously proposed model. It assumes that either an unlinked QTL or a fully informative marker (i.e., all marker alleles are unique in the base generation) is located between the loci. We discuss the need of developing a general algorithm to estimate the variance-covariance structure of the random epistatic effect for linked loci.     

Dynamics of molecular markers linked to the resistance loci in a mosquito-Plasmodium system

Models on the evolution of resistance to parasitism generally assume fitness tradeoffs between the costs of being parasitized and the costs associated with resistance. This study tested this assumption using the yellow fever mosquito Aedes aegypti and malaria parasite Plasmodium gallinaceum system. Experimental mosquito populations were created by mixing susceptible and resistant strains in equal proportions, and then the dynamics of markers linked to loci for Plasmodium resistance and other unlinked neutral markers were determined over 12 generations. We found that when the mixed population was maintained under parasite-free conditions, the frequencies of alleles specific to the susceptible strain at markers closely linked to the loci for resistance (QTL markers) as well as other unlinked markers increased significantly in the first generation and then fluctuated around equilibrium frequencies for all six markers. However, when the mixed population was exposed to an infected blood meal every generation, allele frequencies at the QTL markers for resistance were not significantly changed. Small population size caused significant random fluctuations of allele frequencies at all marker loci. Consistent allele frequency changes in the QTL markers and other unlinked markers suggest that the reduced fitness in the resistant population has a genome-wide effect on the genetic makeup of the mixed population. Continuous exposure to parasites promoted the maintenance of alleles from the resistant Moyo-R strain in the mixed population. The results are discussed in relation to the proposed malaria control strategy through genetic disruption of vector competence.     

Precision Mapping of Quantitative Trait Loci

Adequate separation of effects of possible multiple linked quantitative trait loci (QTLs) on mapping QTLs is the key to increasing the precision of QTL mapping. A new method of QTL mapping is proposed and analyzed in this paper by combining interval mapping with multiple regression. The basis of the proposed method is an interval test in which the test statistic on a marker interval is made to be unaffected by QTLs located outside a defined interval. This is achieved by fitting other genetic markers in the statistical model as a control when performing interval mapping. Compared with the current QTL mapping method (i.e., the interval mapping method which uses a pair or two pairs of markers for mapping QTLs), this method has several advantages. (1) By confining the test to one region at a time, it reduces a multiple dimensional search problem (for multiple QTLs) to a one dimensional search problem. (2) By conditioning linked markers in the test, the sensitivity of the test statistic to the position of individual QTLs is increased, and the precision of QTL mapping can be improved. (3) By selectively and simultaneously using other markers in the analysis, the efficiency of QTL mapping can be also improved. The behavior of the test statistic under the null hypothesis and appropriate critical value of the test statistic for an overall test in a genome are discussed and analyzed. A simulation study of QTL mapping is also presented which illustrates the utility, properties, advantages and disadvantages of the method.     

Identification of selective sweeps using a dynamically adjusted number of linked microsatellites

There is currently large interest in distinguishing the signatures of genetic variation produced by demographic events from those produced by natural selection. We propose a simple multilocus statistical test to identify candidate sites of selective sweeps with high power. The test is based on the variability profile measured in an array of linked microsatellites. We also show that the analysis of flanking markers drastically reduces the number of false positives among the candidates that are identified in a genomewide survey of unlinked loci and find that this property is maintained in many population-bottleneck scenarios. However, for a certain range of intermediately severe population bottlenecks we find genomic signatures that are very similar to those produced by a selective sweep. While in these worst-case scenarios the power of the proposed test remains high, the false-positive rate reaches values close to 50%. Hence, selective sweeps may be hard to identify even if multiple linked loci are analyzed. Nevertheless, the integration of information from multiple linked loci always leads to a considerable reduction of the false-positive rate compared to a genome scan of unlinked loci. We discuss the application of this test to experimental data from Drosophila melanogaster.     

Methods for multiple-marker mapping of quantitative trait loci in half-sib populations

In this paper we consider the detection of individual loci controlling quantitative traits of interest (quantitative trait loci or QTLs) in the large half-sib family structure found in some species. Two simple approaches using multiple markers are proposed, one using least squares and the other maximum likelihood. These methods are intended to provide a relatively fast screening of the entire genome to pinpoint regions of interest for further investigation. They are compared with a more traditional single-marker least-squares approach. The use of multiple markers is shown to increase power and has the advantage of providing an estimate for the location of the QTL. The maximum-likelihood and the least-squares approaches using multiple markers give similar power and estimates for the QTL location, although the likelihood approach also provides estimates of the QTL effect and sire heterozygote frequency. A number of assumptions have been made in order to make the likelihood calculations feasible, however, and computationally it is still more demanding than the least-squares approach. The least-squares approach using multiple markers provides a fast method that can easily be extended to include additional effects.     

Permutation Tests for Multiple Loci Affecting a Quantitative Character

The problem of detecting minor quantitative trait loci (QTL) responsible for genetic variation not explained by major QTL is of importance in the complete dissection of quantitative characters. Two extensions of the permutation-based method for estimating empirical threshold values are presented. These methods, the conditional empirical threshold (CET) and the residual empirical threshold (RET), yield critical values that can be used to construct tests for the presence of minor QTL effects while accounting for effects of known major QTL. The CET provides a completely nonparametric test through conditioning on markers linked to major QTL. It allows for general nonadditive interactions among QTL, but its practical application is restricted to regions of the genome that are unlinked to the major QTL. The RET assumes a structural model for the effect of major QTL, and a threshold is constructed using residuals from this structural model. The search space for minor QTL is unrestricted, and RET-based tests may be more powerful than the CET-based test when the structural model is approximately true.     

Detection of putative quantitative trait loci in line crosses under infinitesimal genetic models

Quantitative trait locus (QTL) mapping studies often employ segregating generations derived from a cross between genetically divergent inbred lines. In the analysis of such data it is customary to fit a single QTL and use a null hypothesis which assumes that the genomic region under study contributes no genetic variance. To explore the situation in which multiple linked genes contribute to the genetic variance, we simulated an F₂-mapping experiment in which the genetic difference between the two original inbred strains was caused by a large number of loci, each having equal effect on the quantitative trait. QTLs were either in coupling, dispersion or repulsion phase in the base population of inbred lines, with the expected F₂ genetic variance explained by the QTLs being equivalent in the three models. Where QTLs were in coupling phase, one inbred line was fixed for all plus alleles, and the other line was fixed for minus alleles. Where QTLs were in dispersion phase, they were assumed to be randomly fixed for one or other allele (as if the inbred lines had evolved from a common ancestor by random drift). Where QTLs were in repulsion phase alleles within an inbred line were alternating plus and minus at adjacent loci, and alternative alleles were fixed in the two inbred lines. In all these genetic models a standard interval mapping test statistic used to determine whether there is a QTL of large effect segregating in the population was inflated on average. Furthermore, the use of a threshold for QTL detection derived under the assumption that no QTLs were segregating would often lead to spurious conclusions regards the presence of genes of large effects (i.e. type I errors). The employment of an alternative model for the analysis, including linked markers as cofactors in the analysis of a single interval, reduced the problem of type I error rate, although test statistics were still inflated relative to the case of no QTLs. It is argued that in practice one should take into account the difference between the strains or the genetic variance in the F₂ population when setting significance thresholds. In addition, tests designed to probe the adequacy of a single-QTL model or of an alternative infinitesimal coupling model are described. Such tests should be applied in QTL mapping studies to help dissect the true nature of genetic variation.     

A simple and rapid method for calculating identity-by-descent matrices using multiple markers

A fast, partly recursive deterministic method for calculating Identity-by-Descent (IBD) probabilities was developed with the objective of using IBD in Quantitative Trait Locus (QTL) mapping. The method combined a recursive method for a single marker locus with a method to estimate IBD between sibs using multiple markers. Simulated data was used to compare the deterministic method developed in the present paper with a stochastic method (LOKI) for precision in estimating IBD probabilities and performance in the task of QTL detection with the variance component approach. This comparison was made in a variety of situations by varying family size and degree of polymorphism among marker loci. The following were observed for the deterministic method relative to MCMC: (i) it was an order of magnitude faster; (ii) its estimates of IBD probabilities were found to agree closely, even though it does not extract information when haplotypes are not known with certainty; (iii) the shape of the profile for the QTL test statistic as a function of location was similar, although the magnitude of the test statistic was slightly smaller; and (iv) the estimates of QTL variance was similar. It was concluded that the method proposed provided a rapid means of calculating the IBD matrix with only a small loss in precision, making it an attractive alternative to the use of stochastic MCMC methods. Furthermore, developments in marker technology providing denser maps would enhance the relative advantage of this method.     

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

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

Genetic models to estimate additive and non-additive effects of marker-associated QTL using multiple regression techniques

Summary The development of molecular markers has recently raised expectations for their application in selection programs. However, some questions related to quantitative trait loci (QTL) identification are still unanswered. The objectives of this paper are (1) to develop statistical genetic models for detecting and locating on the genome multi-QTL with additive, dominance and epistatic effects using multiple linear regression analysis in the backcross and F_n generations from the cross of two inbred lines; and (2) to discuss the bias caused by linked and unlinked QTL on the genetic estimates. Non-linear models were developed for different backcross and F_n generations when both epistasis and no epistasis were assumed. Generation analysis of marked progenies is suggested as a way of increasing the number of observations for the estimates without additional cost for molecular scoring. Some groups of progenies can be created in different generations from the same scored individuals. The non-linear models were transformed into approximate multivariate linear models to which combined stepwise and standard regression analysis could be applied. Expressions for the biases of the marker classes from linked QTL were obtained when no epistasis was assumed. When epistasis was assumed, these expressions increased in complexity, and the biases were caused by both linked and unlinked QTL.     

Molecular markers for seed colour in Brassica juncea.

A detailed RFLP map was used to map QTLs associated with seed colour in Brassica juncea using a doubled-haploid population derived from a cross between a black/brown-seeded cultivar and a yellow-seeded breeding line. Segregation analysis suggested that seed colour was under control of 2 unlinked loci with duplicate gene action. However, QTL analysis revealed 3 QTLs, SC-B4, SC-A10 and SC-A6, affecting seed colour. The QTLs were consistent across environments, and individually explained 43%, 31%, and 16%, respectively, and collectively 62% of the phenotypic variation in the population. Digenic interaction analysis showed that closest flanking locus of QTL SC-B4, wg7b6cNM, had strong epistasis with the locus wg5a1a, which is tightly linked to QTL SC-A6. The interaction of these 2 loci explained 27% of the phenotypic variation in the population, while the whole model explained 84%. In a multiple regression model, the effects of QTL SC-A10, as well as its interaction with other loci, were non-significant, whereas the effects of loci wg7b6cNM and wg5a1a and their interaction were significant. Ninety-eight percent of the DH lines carried the expected alleles of loci wg7b6cNM and wg5a1a for seed colour, confirming that only these 2 loci were linked to seed colour in B. juncea. Four additional digenic interactions significantly affected seed colour, and all 5 digenic interactions were consistent across environments.     

Methodologies for segregation analysis and QTL mapping in plants

Most characters of biological interest and economic importance are quantitative traits. To uncover the genetic architecture of quantitative traits, two approaches have become popular in China. One is the establishment of an analytical model for mixed major-gene plus polygenes inheritance and the other the discovery of quantitative trait locus (QTL). Here we review our progress employing these two approaches. First, we proposed joint segregation analysis of multiple generations for mixed major-gene plus polygenes inheritance. Second, we extended the multilocus method of Lander and Green (1987), Jiang and Zeng (1997) to a more generalized approach. Our methodology handles distorted, dominant and missing markers, including the effect of linked segregation distortion loci on the estimation of map distance. Finally, we developed several QTL mapping methods. In the Bayesian shrinkage estimation (BSE) method, we suggested a method to test the significance of QTL effects and studied the effect of the prior distribution of the variance of QTL effect on QTL mapping. To reduce running time, a penalized maximum likelihood method was adopted. To mine novel genes in crop inbred lines generated in the course of normal crop breeding work, three methods were introduced. If a well-documented genealogical history of the lines is available, two-stage variance component analysis and multi-QTL Haseman-Elston regression were suggested; if unavailable, multiple loci in silico mapping was proposed.