The Genetic Architecture of Quantitative Traits Cannot Be Inferred from Variance Component Analysis

Classical quantitative genetic analyses estimate additive and non-additive genetic and environmental components of variance from phenotypes of related individuals without knowing the identities of quantitative trait loci (QTLs). Many studies have found a large proportion of quantitative trait variation can be attributed to the additive genetic variance (V_A), providing the basis for claims that non-additive gene actions are unimportant. In this study, we show that arbitrarily defined parameterizations of genetic effects seemingly consistent with non-additive gene actions can also capture the majority of genetic variation. This reveals a logical flaw in using the relative magnitudes of variance components to indicate the relative importance of additive and non-additive gene actions. We discuss the implications and propose that variance component analyses should not be used to infer the genetic architecture of quantitative traits.     

Genetic effects at pleiotropic loci are context-dependent with consequences for the maintenance of genetic variation in populations

Context-dependent genetic effects, including genotype-by-environment and genotype-by-sex interactions, are a potential mechanism by which genetic variation of complex traits is maintained in populations. Pleiotropic genetic effects are also thought to play an important role in evolution, reflecting functional and developmental relationships among traits. We examine context-dependent genetic effects at pleiotropic loci associated with normal variation in multiple metabolic syndrome (MetS) components (obesity, dyslipidemia, and diabetes-related traits). MetS prevalence is increasing in Western societies and, while environmental in origin, presents substantial variation in individual response. We identify 23 pleiotropic MetS quantitative trait loci (QTL) in an F₁₆ advanced intercross between the LG/J and SM/J inbred mouse strains (Wustl:LG,SM-G16; n = 1002). Half of each family was fed a high-fat diet and half fed a low-fat diet; and additive, dominance, and parent-of-origin imprinting genotypic effects were examined in animals partitioned into sex, diet, and sex-by-diet cohorts. We examine the context-dependency of the underlying additive, dominance, and imprinting genetic effects of the traits associated with these pleiotropic QTL. Further, we examine sequence polymorphisms (SNPs) between LG/J and SM/J as well as differential expression of positional candidate genes in these regions. We show that genetic associations are different in different sex, diet, and sex-by-diet settings. We also show that over- or underdominance and ecological cross-over interactions for single phenotypes may not be common, however multidimensional synthetic phenotypes at loci with pleiotropic effects can produce situations that favor the maintenance of genetic variation in populations. Our findings have important implications for evolution and the notion of personalized medicine.     

Dissecting the loci underlying maturation timing in Atlantic salmon using haplotype and multi-SNP based association methods

Characterizing the role of different mutational effect sizes in the evolution of fitness-related traits has been a major goal in evolutionary biology for a century. Such characterization in a diversity of systems, both model and non-model, will help to understand the genetic processes underlying fitness variation. However, well-characterized genetic architectures of such traits in wild populations remain uncommon. In this study, we used haplotype-based and multi-SNP Bayesian association methods with sequencing data for 313 individuals from wild populations to test the mutational composition of known candidate regions for sea age at maturation in Atlantic salmon (Salmo salar). We detected an association at five loci out of 116 candidates previously identified in an aquaculture strain with maturation timing in wild Atlantic salmon. We found that at four of these five loci, variation explained by the locus was predominantly driven by a single SNP suggesting the genetic architecture of this trait includes multiple loci with simple, non-clustered alleles and a locus with potentially more complex alleles. This highlights the diversity of genetic architectures that can exist for fitness-related traits. Furthermore, this study provides a useful multi-SNP framework for future work using sequencing data to characterize genetic variation underlying phenotypes in wild populations.Subject terms: Evolutionary genetics, Genetic association study     

Distinct genetic control of parasite elimination, dissemination, and disease after Leishmania major infection

Elimination of pathogens is the basis of host resistance to infections; however, relationship between persisting pathogens and disease has not been clarified. Leishmania major infection in mice is an important model of host–pathogen relationship. Infected BALB/c mice exhibit high parasite numbers in lymph nodes and spleens, and a chronic disease with skin lesions, splenomegaly, and hepatomegaly, increased serum IgE levels and cytokine imbalance. Although numerous gene loci affecting these disease symptoms have been reported, genes controlling parasites’ elimination or dissemination have never been mapped. We therefore compared genetics of the clinical and immunologic symptomatology with parasite load in (BALB/c?×?CcS-11) F₂ hybrids and mapped five loci, two of which control parasite elimination or dissemination. Lmr5 influences parasite loads in spleens (and skin lesions, splenomegaly, and serum IgE, IL-4, and IFNγ levels), and Lmr20 determines parasite numbers in draining lymph nodes (and serum levels of IgE and IFNγ), but no skin or visceral pathology. Three additional loci do not affect parasite numbers but influence significantly the disease phenotype—Lmr21: skin lesions and IFNγ levels, Lmr22: IL-4 levels, Lmr23: IFNγ levels, indicating that development of L. major-caused disease includes critical regulations additional to control of parasite spread.     

Genomic regions of pufferfishes responsible for host specificity of a monogenean parasite, Heterobothrium okamotoi

The genetic mechanisms underlying host specificity of parasitic infections are largely unknown. After hatching, the larvae of the monogenean parasite, Heterobothrium okamotoi, attach to the gill filaments of hosts and the post-larval worms develop there by consuming nutrients from the host. The susceptibility to H. okamotoi infection differs markedly among fish species. While this parasite can grow on tiger pufferfish (also called fugu), Takifugu rubripes, it appears to be rejected by a close congener, grass pufferfish, Takifugu niphobles, after initial attachment to the gills. To determine the genetic architecture of the pufferfish responsible for this host specificity, we performed genome-wide quantitative trait loci analysis. We raised second generation (F₂) hybrids of the two pufferfish species and experimentally infected them with the monogenean in vivo. To assess possible differences in host mechanisms between early and later periods of infection, we sampled fish three h and 21 days after exposure. Genome scanning of fish from the 3 h infection trial revealed suggestive quantitative trait loci on linkage groups 2 and 14, which affected the number of parasites on the gill. However, analysis of fish 21 days p.i. detected a significant quantitative trait locus on linkage group 9 and three other suggestive quantitative trait loci on linkage groups 7, 18 and 22. These results indicated the polygenic nature of the host mechanisms involved in the infection/rejection of H. okamotoi. Moreover the analyses suggested that host factors may play a more important role during the growth period of the parasite than during initial host recognition at the time of attachment. Within the 95% confidence interval of the linkage group 9 quantitative trait locus in the fugu genome, there were 214 annotated protein-coding genes, including immunity-related genes such as Irak4, Muc2 and Muc5ac.     

Genetic association with boldness and maternal performance in a free-ranging population of grey seals (Halichoerus grypus)

Individual variation in quantitative traits clearly influence many ecological and evolutionary processes. Moderate to high heritability estimates of personality and life-history traits suggest some level of genetic control over these traits. Yet, we know very little of the underlying genetic architecture of phenotypic variation in the wild. In this study, we used a candidate gene approach to investigate the association of genetic variants with repeated measures of boldness and maternal performance traits (weaning mass and lactation duration) collected over an 11- and 28-year period, respectively, in a free-ranging population of grey seals on Sable Island National Park Reserve, Canada. We isolated and re-sequenced five genes: dopamine receptor D4 (DRD4), serotonin transporter (SERT), oxytocin receptor (OXTR), and melanocortin receptors 1 (MC1R) and 5 (MC5R). We discovered single nucleotide polymorphisms (SNPs) in each gene; and, after accounting for loci in linkage disequilibrium and filtering due to missing data, we were able to test for genotype-phenotype relationships at seven loci in three genes (DRD4, SERT, and MC1R). We tested for association between these loci and traits of 180 females having extreme shy-bold phenotypes using mixed-effects models. One locus within SERT was significantly associated with boldness (effect size = 0.189) and a second locus within DRD4 with weaning mass (effect size = 0.232). Altogether, genotypes explained 6.52–13.66% of total trait variation. Our study substantiates SERT and DRD4 as important determinants of behaviour, and provides unique insight into the molecular mechanisms underlying maternal performance variation in a marine predator.Subject terms: Behavioural ecology, Evolutionary genetics, Behavioural genetics, Genetic association study, Quantitative trait     

Genome-Wide Association Analysis of Adaptation Using Environmentally Predicted Traits

Current methods for studying the genetic basis of adaptation evaluate genetic associations with ecologically relevant traits or single environmental variables, under the implicit assumption that natural selection imposes correlations between phenotypes, environments and genotypes. In practice, observed trait and environmental data are manifestations of unknown selective forces and are only indirectly associated with adaptive genetic variation. In theory, improved estimation of these forces could enable more powerful detection of loci under selection. Here we present an approach in which we approximate adaptive variation by modeling phenotypes as a function of the environment and using the predicted trait in multivariate and univariate genome-wide association analysis (GWAS). Based on computer simulations and published flowering time data from the model plant Arabidopsis thaliana, we find that environmentally predicted traits lead to higher recovery of functional loci in multivariate GWAS and are more strongly correlated to allele frequencies at adaptive loci than individual environmental variables. Our results provide an example of the use of environmental data to obtain independent and meaningful information on adaptive genetic variation.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

Dissecting adaptive clinal variation: markers, inversions and size/stress associations in Drosophila melanogaster from a central field population

Many organisms show latitudinal variation for quantitative traits that is assumed to be due to climatic adaptation. These clines provide an opportunity to study the genetics of the adaptive process both at the phenotypic and the underlying molecular levels. Yet researchers rarely try to link variation in quantitative traits to their underlying molecular genetic basis. We describe a novel approach for exploring the genetic basis for clinal variation in size and stress traits in Drosophila melanogaster. We look for associations between genetic markers and traits that exhibit clinal patterns on the east coast of Australia using a single, geographically central population. There are strong associations between markers found within In(3R)Payne and variation in size, suggesting that this inversion explains much of the clinal variation in this trait. We also find that development time is associated with the Adh allozyme locus, cold resistance is negatively associated with the In(3L)Payne inversion and a genetic marker for Hsp70, a heat‐shock protein, is associated with heat resistance. Finally we discuss the importance of inversions in clinal variation for quantitative traits and for identifying quantitative trait loci.     

Molecular-Marker-Facilitated Investigations of Quantitative-Trait Loci in Maize. I. Numbers, Genomic Distribution and Types of Gene Action

Individual genetic factors which underlie variation in quantitative traits of maize were investigated in each of two F2 populations by examining the mean trait expressions of genotypic classes at each of 17-20 segregating marker loci. It was demonstrated that the trait expression of marker locus classes could be interpreted in terms of genetic behavior at linked quantitative trait loci (QTLs). For each of 82 traits evaluated, QTLs were detected and located to genomic sites. The numbers of detected factors varied according to trait, with the average trait significantly influenced by almost two-thirds of the marked genomic sites. Most of the detected associations between marker loci and quantitative traits were highly significant, and could have been detected with fewer than the 1800-1900 plants evaluated in each population. The cumulative, simple effects of marker-linked regions of the genome explained between 8 and 40% of the phenotypic variation for a subset of 25 traits evaluated. Single marker loci accounted for between 0.3% and 16% of the phenotypic variation of traits. Individual plant heterozygosity, as measured by marker loci, was significantly associated with variation in many traits. The apparent types of gene action at the QTLs varied both among traits and between loci for given traits, although overdominance appeared frequently, especially for yield-related traits. The prevalence of apparent overdominance may reflect the effects of multiple QTLs within individual marker-linked regions, a situation which would tend to result in overestimation of dominance. Digenic epistasis did not appear to be important in determining the expression of the quantitative traits evaluated. Examination of the effects of marked regions on the expression of pairs of traits suggests that genomic regions vary in the direction and magnitudes of their effects on trait correlations, perhaps providing a means of selecting to dissociate some correlated traits. Marker-facilitated investigations appear to provide a powerful means of examining aspects of the genetic control of quantitative traits. Modifications of the methods employed herein will allow examination of the stability of individual gene effects in varying genetic backgrounds and environments.     

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host–parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host–parasite interactions. In the QTL panel used here, Daphnia''s resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host–parasite systems. Only the Pasteuria–Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium–Daphnia system remains unclear.     

Unifying Genetic Canalization,Genetic Constraint,and Genotype-by-Environment Interaction: QTL by Genomic Background by Environment Interaction of Flowering Time in Boechera stricta

Natural populations exhibit substantial variation in quantitative traits. A quantitative trait is typically defined by its mean and variance, and to date most genetic mapping studies focus on loci altering trait means but not (co)variances. For single traits, the control of trait variance across genetic backgrounds is referred to as genetic canalization. With multiple traits, the genetic covariance among different traits in the same environment indicates the magnitude of potential genetic constraint, while genotype-by-environment interaction (GxE) concerns the same trait across different environments. While some have suggested that these three attributes of quantitative traits are different views of similar concepts, it is not yet clear, however, whether they have the same underlying genetic mechanism. Here, we detect quantitative trait loci (QTL) influencing the (co)variance of phenological traits in six distinct environments in Boechera stricta, a close relative of Arabidopsis. We identified nFT as the QTL altering the magnitude of phenological trait canalization, genetic constraint, and GxE. Both the magnitude and direction of nFT''s canalization effects depend on the environment, and to our knowledge, this reversibility of canalization across environments has not been reported previously. nFT''s effects on trait covariance structure (genetic constraint and GxE) likely result from the variable and reversible canalization effects across different traits and environments, which can be explained by the interaction among nFT, genomic backgrounds, and environmental stimuli. This view is supported by experiments demonstrating significant nFT by genomic background epistatic interactions affecting phenological traits and expression of the candidate gene, FT. In contrast to the well-known canalization gene Hsp90, the case of nFT may exemplify an alternative mechanism: Our results suggest that (at least in traits with major signal integrators such as flowering time) genetic canalization, genetic constraint, and GxE may have related genetic mechanisms resulting from interactions among major QTL, genomic backgrounds, and environments.     

QTL analysis of floral traits in Louisiana iris hybrids

The formation of hybrid zones between nascent species is a widespread phenomenon. The evolutionary consequences of hybridization are influenced by numerous factors, including the action of natural selection on quantitative trait variation. Here we examine how the genetic basis of floral traits of two species of Louisiana Irises affects the extent of quantitative trait variation in their hybrids. Quantitative trait locus (QTL) mapping was used to assess the size (magnitude) of phenotypic effects of individual QTL, the degree to which QTL for different floral traits are colocalized, and the occurrence of mixed QTL effects. These aspects of quantitative genetic variation would be expected to influence (1) the number of genetic steps (in terms of QTL substitutions) separating the parental species phenotypes; (2) trait correlations; and (3) the potential for transgressive segregation in hybrid populations. Results indicate that some Louisiana Iris floral trait QTL have large effects and QTL for different traits tend to colocalize. Transgressive variation was observed for six of nine traits, despite the fact that mixed QTL effects influence few traits. Overall, our QTL results imply that the genetic basis of floral morphology and color traits might facilitate the maintenance of phenotypic divergence between Iris fulva and Iris brevicaulis, although a great deal of phenotypic variation was observed among hybrids.     

Accounting for a quantitative trait locus for plasma triglyceride levels: utilization of variants in multiple genes

Background

For decades, research efforts have tried to uncover the underlying genetic basis of human susceptibility to a variety of diseases. Linkage studies have resulted in highly replicated findings and helped identify quantitative trait loci (QTL) for many complex traits; however identification of specific alleles accounting for linkage remains elusive. The purpose of this study was to determine whether with a sufficient number of variants a linkage signal can be fully explained.

Method

We used comprehensive fine-mapping using a dense set of single nucleotide polymorphisms (SNPs) across the entire quantitative trait locus (QTL) on human chromosome 7q36 linked to plasma triglyceride levels. Analyses included measured genotype and combined linkage association analyses.

Results

Screening this linkage region, we found an over representation of nominally significant associations in five genes (MLL3, DPP6, PAXIP1, HTR5A, INSIG1). However, no single genetic variant was sufficient to account for the linkage. On the other hand, multiple variants capturing the variation in these five genes did account for the linkage at this locus. Permutation analyses suggested that this reduction in LOD score was unlikely to have occurred by chance (p = 0.008).

Discussion

With recent findings, it has become clear that most complex traits are influenced by a large number of genetic variants each contributing only a small percentage to the overall phenotype. We found that with a sufficient number of variants, the linkage can be fully explained. The results from this analysis suggest that perhaps the failure to identify causal variants for linkage peaks may be due to multiple variants under the linkage peak with small individual effect, rather than a single variant of large effect.     

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

The importance of context to the genetic architecture of diabetes-related traits is revealed in a genome-wide scan of a LG/J?��?SM/J murine model

Variations in diabetic phenotypes are caused by complex interactions of genetic effects, environmental factors, and the interplay between the two. We tease apart these complex interactions by examining genome-wide genetic and epigenetic effects on diabetes-related traits among different sex, diet, and sex-by-diet cohorts in a Mus musculus model. We conducted a genome-wide scan for quantitative trait loci that affect serum glucose and insulin levels and response to glucose stress in an F₁₆ Advanced Intercross Line of the LG/J and SM/J intercross (Wustl:LG,SM-G16). Half of each sibship was fed a high-fat diet and half was fed a relatively low-fat diet. Context-dependent genetic (additive and dominance) and epigenetic (parent-of-origin imprinting) effects were characterized by partitioning animals into sex, diet, and sex-by-diet cohorts. We found that different cohorts often have unique genetic effects at the same loci, and that genetic signals can be masked or erroneously assigned to specific cohorts if they are not considered individually. Our data demonstrate that the effects of genes on complex trait variation are highly context-dependent and that the same genomic sequence can affect traits differently depending on an individual??s sex and/or dietary environment. Our results have important implications for studies of complex traits in humans.     

Accounting for Genetic Architecture Improves Sequence Based Genomic Prediction for a Drosophila Fitness Trait

The ability to predict quantitative trait phenotypes from molecular polymorphism data will revolutionize evolutionary biology, medicine and human biology, and animal and plant breeding. Efforts to map quantitative trait loci have yielded novel insights into the biology of quantitative traits, but the combination of individually significant quantitative trait loci typically has low predictive ability. Utilizing all segregating variants can give good predictive ability in plant and animal breeding populations, but gives little insight into trait biology. Here, we used the Drosophila Genetic Reference Panel to perform both a genome wide association analysis and genomic prediction for the fitness-related trait chill coma recovery time. We found substantial total genetic variation for chill coma recovery time, with a genetic architecture that differs between males and females, a small number of molecular variants with large main effects, and evidence for epistasis. Although the top additive variants explained 36% (17%) of the genetic variance among lines in females (males), the predictive ability using genomic best linear unbiased prediction and a relationship matrix using all common segregating variants was very low for females and zero for males. We hypothesized that the low predictive ability was due to the mismatch between the infinitesimal genetic architecture assumed by the genomic best linear unbiased prediction model and the true genetic architecture of chill coma recovery time. Indeed, we found that the predictive ability of the genomic best linear unbiased prediction model is markedly improved when we combine quantitative trait locus mapping with genomic prediction by only including the top variants associated with main and epistatic effects in the relationship matrix. This trait-associated prediction approach has the advantage that it yields biologically interpretable prediction models.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.