QTL mapping and the genetic basis of adaptation: recent developments

Quantitative trait loci (QTL) mapping has been used in a number of evolutionary studies to study the genetic basis of adaptation by mapping individual QTL that explain the differences between differentiated populations and also estimating their effects and interaction in the mapping population. This analysis can provide clues about the evolutionary history of populations and causes of the population differentiation. QTL mapping analysis methods and associated computer programs provide us tools for such an inference on the genetic basis and architecture of quantitative trait variation in a mapping population. Current methods have the capability to separate and localize multiple QTL and estimate their effects and interaction on a quantitative trait. More recent methods have been targeted to provide a comprehensive inference on the overall genetic architecture of multiple traits in a number of environments. This development is important for evolutionary studies on the genetic basis of multiple trait variation, genotype by environment interaction, host–parasite interaction, and also microarray gene expression QTL analysis.     

Mapping quantitative trait loci in multiple populations of Arabidopsis thaliana identifies natural allelic variation for trichome density

The majority of biological traits are genetically complex. Mapping the quantitative trait loci (QTL) that determine these phenotypes is a powerful means for estimating many parameters of the genetic architecture for a trait and potentially identifying the genes responsible for natural variation. Typically, such experiments are conducted in a single mapping population and, therefore, have only the potential to reveal genomic regions that are polymorphic between the progenitors of the population. What remains unclear is how well the QTL identified in any one mapping experiment characterize the genetics that underlie natural variation in traits. Here we provide QTL mapping data for trichome density from four recombinant inbred mapping populations of Arabidopsis thaliana. By aligning the linkage maps for these four populations onto a common physical map, the results from each experiment were directly compared. Seven of the nine QTL identified are population specific while two were mapped in all four populations. Our results show that many lineage-specific alleles that either increase or decrease trichome density persist in natural populations and that most of this genetic variation is additive. More generally, these findings suggest that the use of multiple populations holds great promise for better understanding the genetic architecture of natural variation.     

Genomic dissection of variation in clutch size and egg mass in a wild great tit (Parus major) population

Clutch size and egg mass are life history traits that have been extensively studied in wild bird populations, as life history theory predicts a negative trade‐off between them, either at the phenotypic or at the genetic level. Here, we analyse the genomic architecture of these heritable traits in a wild great tit (Parus major) population, using three marker‐based approaches – chromosome partitioning, quantitative trait locus (QTL) mapping and a genome‐wide association study (GWAS). The variance explained by each great tit chromosome scales with predicted chromosome size, no location in the genome contains genome‐wide significant QTL, and no individual SNPs are associated with a large proportion of phenotypic variation, all of which may suggest that variation in both traits is due to many loci of small effect, located across the genome. There is no evidence that any regions of the genome contribute significantly to both traits, which combined with a small, nonsignificant negative genetic covariance between the traits, suggests the absence of genetic constraints on the independent evolution of these traits. Our findings support the hypothesis that variation in life history traits in natural populations is likely to be determined by many loci of small effect spread throughout the genome, which are subject to continued input of variation by mutation and migration, although we cannot exclude the possibility of an additional input of major effect genes influencing either trait.     

Ontogenetics of QTL: the genetic architecture of trichome density over time in Arabidopsis thaliana

Although much is known about the molecular genetic basis of trichome development in Arabidopsis thaliana, less is known about the underlying genetic basis of continuous variation in a trait known to be of adaptive importance: trichome density. The density of leaf trichomes is known to be a major determinant of herbivore damage in natural populations of A. thaliana and herbivores are a significant selective force on genetic variation for trichome density. A number of developmental changes occur during ontogeny in A. thaliana, including changes in trichome density. I used multiple interval mapping (MIM) analysis to identify QTL responsible for trichome density on both juvenile leaves and adult leaves in replicate, independent trials and asked whether those QTL changed with ontogeny. In both juvenile and adult leaves, I detected a single major QTL on chromosome 2 that explained much of the genetic variance. Although additional QTL were detected, there were no consistent differences in the genetic architecture of trichome density measured on juvenile and adult leaves. The finding of a single QTL of major effect for a trait of known adaptive importance suggests that genes of major effect may play an important role in adaptation.     

Gene networks in hexaploid wheat: interacting quantitative trait loci for grain protein content

In hexaploid wheat, single-locus and two-locus quantitative trait loci (QTL) analyses for grain protein content (GPC) were conducted using two different mapping populations (PI and PII). Main effect QTLs (M-QTLs), epistatic QTLs (E-QTLs) and QTL x environment interactions (QE, QQE) were detected using two-locus analyses in both the populations. Only a few QTLs were common in both the analyses, and the QTLs and the interactions detected in the two populations differed, suggesting the superiority of two-locus analysis and the need for using several mapping populations for QTL analysis. A sizable proportion of genetic variation for GPC was due to interactions (28.59% and 54.03%), rather than to M-QTL effects (7.24% and 7.22%), which are the only genetic effects often detected in the majority of QTL studies. Even E-QTLs made a marginal contribution to genetic variation (2.68% and 6.04%), thus suggesting that the major part of genetic variation is due to changes in gene networks rather than the presence or absence of specific genes. This is in sharp contrast to the genetic dissection of pre-harvest sprouting tolerance conducted by us earlier, where interaction effects were not substantial, suggesting that the nature of genetic variation also depends on the nature of the trait.     

Do large effect QTL fractionate? A case study at the maize domestication QTL teosinte branched1

Quantitative trait loci (QTL) mapping is a valuable tool for studying the genetic architecture of trait variation. Despite the large number of QTL studies reported in the literature, the identified QTL are rarely mapped to the underlying genes and it is usually unclear whether a QTL corresponds to one or multiple linked genes. Similarly, when QTL for several traits colocalize, it is usually unclear whether this is due to the pleiotropic action of a single gene or multiple linked genes, each affecting one trait. The domestication gene teosinte branched1 (tb1) was previously identified as a major domestication QTL with large effects on the differences in plant and ear architecture between maize and teosinte. Here we present the results of two experiments that were performed to determine whether the single gene tb1 explains all trait variation for its genomic region or whether the domestication QTL at tb1 fractionates into multiple linked QTL. For traits measuring plant architecture, we detected only one QTL per trait and these QTL all mapped to tb1. These results indicate that tb1 is the sole gene for plant architecture traits that segregates in our QTL mapping populations. For most traits related to ear morphology, we detected multiple QTL per trait in the tb1 genomic region, including a large effect QTL at tb1 itself plus one or two additional linked QTL. tb1 is epistatic to two of these additional QTL for ear traits. Overall, these results provide examples for both a major QTL that maps to a single gene, as well as a case in which a QTL fractionates into multiple linked QTL.     

Detection of QTLs for flowering date in three mapping populations of the model legume species Medicago truncatula

Adaptation to the environment and reproduction are dependent on the date of flowering in the season. The objectives of this paper were to evaluate the effect of photoperiod on flowering date of the model species for legume crops, Medicago truncatula and to describe genetic architecture of this trait in multiple mapping populations. The effect of photoperiod (12 and 18 h) was analysed on eight lines. Quantitative variation in three recombinant inbred lines (RILs) populations involving four parental lines was evaluated, and QTL detection was carried out. Flowering occurred earlier in long than in short photoperiods. Modelling the rate of progression to flowering with temperature and photoperiod gave high R(2), with line-specific parameters that indicated differential responses of the lines to both photoperiod and temperature. QTL detection showed a QTL on chromosome 7 that was common to all populations and seasons. Taking advantage of the multiple mapping populations, it was condensed into a single QTL with a support interval of only 0.9 cM. In a bioanalysis, six candidate genes were identified in this interval. This design also indicated other genomic regions that were involved in flowering date variation more specifically in one population or one season. The analysis on three different mapping populations detected more QTLs than on a single population, revealed more alleles and gave a more precise position of the QTLs that were common to several populations and/or seasons. Identification of candidate genes was a result of integration of QTL analysis and genomics in M. truncatula.     

Six new recombinant inbred populations for the study of quantitative traits in Arabidopsis thaliana

Quantitative approaches are now widely used to study the genetic architecture of complex traits. However, most studies have been conducted in single mapping populations, which sample only a fraction of the natural allelic variation available within a gene pool and can identify only a subset of the loci controlling the traits. To enable the progress towards an understanding of the global genetic architecture of a broad range of complex traits, we have developed and characterised six new Arabidopsis thaliana recombinant inbred populations. To evaluate the utility of these populations for integrating analyses from multiple populations, we identified quantitative trait loci (QTL) controlling flowering time in vernalized plants growing in 16 h days. We used the physical positions of markers to align the linkage maps of our populations with those of six existing populations. We identified seven QTL in genomic locations coinciding with those identified in previous studies and in addition a further eight QTL were identified. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. An erratum to this article can be found at     

Genetic architecture of complex traits in plants

Genetic architecture refers to the numbers and genome locations of genes that affect a trait, the magnitude of their effects, and the relative contributions of additive, dominant, and epistatic gene effects. Quantitative trait locus (QTL) mapping techniques are commonly used to investigate genetic architectures, but the scope of inferences drawn from QTL studies are often restricted by the limitations of the experimental designs. Recent advances in experimental and statistical procedures, including the simultaneous analysis of QTL that segregate in diverse germplasm, should improve genetic architecture studies. High-resolution QTL mapping methods are being developed that may define the specific DNA sequence variants underlying QTL. Studies of genetic architecture, combined with improved knowledge of the structure of plant populations, will impact our understanding of plant evolution and the design of crop improvement strategies.     

Quantitative trait locus mapping in natural populations: progress, caveats and future directions

Over the last 15 years quantitative trait locus (QTL) mapping has become a popular method for understanding the genetic basis of continuous variation in a variety of systems. For example, the technique is now an integral tool in medical genetics, livestock production, plant breeding and population genetics of model organisms. Ten years ago, it was suggested that the method could be used to understand continuous variation in natural populations. In this review I: (i) clarify what is meant by natural population in the QTL context, (ii) discuss whether evolutionary biologists have successfully mapped QTL in natural populations, (iii) highlight some of the questions that have been addressed by QTL mapping in natural populations, (iv) describe how QTL mapping can be conducted in unmanipulated natural populations, (v) highlight some of the limitations of QTL mapping and (vi) try to predict some future directions for QTL mapping in natural populations.     

Analysis of quantitative trait loci that influence animal behavior

Behavioral differences between inbred strains of mice and rats have a genetic basis that can now be dissected using quantitative trait locus (QTL) analysis. Over the last 10 years, a large number of genetic loci that influence behavior have been mapped. In this article I review what that information has revealed about the genetic architecture of behavior. I show that most behaviors are influenced by QTL of small effect, each contributing to less than 10% of the variance of a behavioral trait. The small effect of each QTL on behavioral variation suggests that the mutational spectrum is different from that which results in Mendelian disorders. Regions of DNA should be appropriately prioritized to find the molecular variants, for instance by looking at sequences that control the level of gene expression rather than variants in coding regions. While the number of allelic loci that can contribute to a trait is large, this is not necessarily the case: the analysis of selected strains shows that a remarkably small number of QTL can explain the bulk of the genetic variation in behavior. I conclude by arguing that genetic mapping has more to offer than a starting point for positional cloning projects. With advances in multivariate analyses, mapping can also test hypotheses about the psychological processes that give rise to behavioral variation.     

Quantitative Trait Loci Differentiating the Outbreeding Mimulus Guttatus from the Inbreeding M. Platycalyx

Theoretical predictions about the evolution of selfing depend on the genetic architecture of loci controlling selfing (monogenic vs. polygenic determination, large vs. small effect of alleles, dominance vs. recessiveness), and studies of such architecture are lacking. We inferred the genetic basis of mating system differences between the outbreeding Mimulus guttatus and the inbreeding M. platycalyx by quantitative trait locus (QTL) mapping using random amplified polymorphic DNA and isozyme markers. One to three QTL were detected for each of five mating system characters, and each QTL explained 7.6-28.6% of the phenotypic variance. Taken together, QTL accounted for up to 38% of the variation in mating system characters, and a large proportion of variation was unaccounted for. Inferred QTL often affected more than one trait, contributing to the genetic correlation between those traits. These results are consistent with the hypothesis that quantitative variation in plant mating system characters is primarily controlled by loci with small effect.     

Detection of QTL for flowering time in multiple families of elite maize

Flowering time is a fundamental quantitative trait in maize that has played a key role in the postdomestication process and the adaptation to a wide range of climatic conditions. Flowering time has been intensively studied and recent QTL mapping results based on diverse founders suggest that the genetic architecture underlying this trait is mainly based on numerous small-effect QTL. Here, we used a population of 684 progenies from five connected families to investigate the genetic architecture of flowering time in elite maize. We used a joint analysis and identified nine main effect QTL explaining approximately 50?% of the genotypic variation of the trait. The QTL effects were small compared with the observed phenotypic variation and showed strong differences between families. We detected no epistasis with the genetic background but four digenic epistatic interactions in a full 2-dimensional genome scan. Our results suggest that flowering time in elite maize is mainly controlled by main effect QTL with rather small effects but that epistasis may also contribute to the genetic architecture of the trait.     

Genetic architecture of a selection response in Arabidopsis thaliana

Quantitative trait locus (QTL) mapping has become an established and effective method for studying the genetic architecture of complex traits. In this report, we use a QTL mapping approach in combination with data from a large selection experiment in Arabidopsis thaliana to explore a response to selection of experimental populations with differentiated genetic backgrounds. Experimental populations with genetic backgrounds derived from ecotypes Landsberg and Niederzenz were exposed to multiple generations of fertility and viability selection. This selection resulted in phenotypic shifts in a number of life-history and fitness-related characters including early development time, flowering time, dry biomass, longevity, and fruit production. Quantitative trait loci were mapped for these traits and their positions were compared to previously characterized allele frequency changes in the experimental populations (Ungerer et al. 2003). Quantitative trait locus positions largely colocalized with genomic regions under strong and consistent selection in populations with differentiated genetic backgrounds, suggesting that alleles for these traits were selected similarly in differentiated genetic backgrounds. However, one QTL region exhibited a more variable response; being positively selected on one genetic background but apparently neutral in another. This study demonstrates how QTL mapping approaches can be combined with map-based population genetic data to study how selection acts on standing genetic variation in populations.     

Quantitative trait loci (QTL) detection in multicross inbred designs: recovering QTL identical-by-descent status information from marker data

Mapping quantitative trait loci in plants is usually conducted using a population derived from a cross between two inbred lines. The power of such QTL detection and the parameter estimates depend largely on the choice of the two parental lines. Thus, the QTL detected in such populations represent only a small part of the genetic architecture of the trait. In addition, the effects of only two alleles are characterized, which is of limited interest to the breeder, while common pedigree breeding material remains unexploited for QTL mapping. In this study, we extend QTL mapping methodology to a generalized framework, based on a two-step IBD variance component approach, applicable to any type of breeding population obtained from inbred parents. We then investigate with simulated data mimicking conventional breeding programs the influence of different estimates of the IBD values on the power of QTL detection. The proposed method would provide an alternative to the development of specifically designed recombinant populations, by utilizing the genetic variation actually managed by plant breeders. The use of these detected QTL in assisting breeding would thus be facilitated.     

CONNECTING QTLS TO THE G-MATRIX OF EVOLUTIONARY QUANTITATIVE GENETICS

Evolutionary quantitative genetics has recently advanced in two distinct streams. Many biologists address evolutionary questions by estimating phenotypic selection and genetic (co)variances ( G matrices). Simultaneously, an increasing number of studies have applied quantitative trait locus (QTL) mapping methods to dissect variation. Both conceptual and practical difficulties have isolated these two foci of quantitative genetics. A conceptual integration follows from the recognition that QTL allele frequencies are the essential variables relating the G -matrix to marker-based mapping experiments. Breeding designs initiated from randomly selected parental genotypes can be used to estimate QTL-specific genetic (co)variances. These statistics appropriately distill allelic variation and provide an explicit population context for QTL mapping estimates. Within this framework, one can parse the G -matrix into a set of mutually exclusive genomic components and ask whether these parts are similar or dissimilar in their respective features, for example the magnitude of phenotypic effects and the extent and nature of pleiotropy. As these features are critical determinants of sustained response to selection, the integration of QTL mapping methods into G -matrix estimation can provide a concrete, genetically based experimental program to investigate the evolutionary potential of natural populations.     

Quantitative trait loci mapping in five new large recombinant inbred line populations of Arabidopsis thaliana genotyped with consensus single-nucleotide polymorphism markers

Quantitative approaches conducted in a single mapping population are limited by the extent of genetic variation distinguishing the parental genotypes. To overcome this limitation and allow a more complete dissection of the genetic architecture of complex traits, we built an integrated set of 15 new large Arabidopsis thaliana recombinant inbred line (RIL) populations optimized for quantitative trait loci (QTL) mapping, having Columbia as a common parent crossed to distant accessions. Here we present 5 of these populations that were validated by investigating three traits: flowering time, rosette size, and seed production as an estimate of fitness. The large number of RILs in each population (between 319 and 377 lines) and the high density of evenly spaced genetic markers scored ensure high power and precision in QTL mapping even under a minimal phenotyping framework. Moreover, the use of common markers across the different maps allows a direct comparison of the QTL detected within the different RIL sets. In addition, we show that following a selective phenotyping strategy by performing QTL analyses on genotypically chosen subsets of 164 RILs (core populations) does not impair the power of detection of QTL with phenotypic contributions >7%.     

Alleles versus mutations: Understanding the evolution of genetic architecture requires a molecular perspective on allelic origins

Perspectives on the role of large‐effect quantitative trait loci (QTL) in the evolution of complex traits have shifted back and forth over the past few decades. Different sets of studies have produced contradictory insights on the evolution of genetic architecture. I argue that much of the confusion results from a failure to distinguish mutational and allelic effects, a limitation of using the Fisherian model of adaptive evolution as the lens through which the evolution of adaptive variation is examined. A molecular‐based perspective reveals that allelic differences can involve the cumulative effects of many mutations plus intragenic recombination, a model that is supported by extensive empirical evidence. I discuss how different selection regimes could produce very different architectures of allelic effects under a molecular‐based model, which may explain conflicting insights on genetic architecture from studies of variation within populations versus between divergently selected populations. I address shortcomings of genome‐wide association study (GWAS) practices in light of more suitable models of allelic evolution, and suggest alternate GWAS strategies to generate more valid inferences about genetic architecture. Finally, I discuss how adopting more suitable models of allelic evolution could help redirect research on complex trait evolution toward addressing more meaningful questions in evolutionary biology.     

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.