Quantitative genetic analysis of natural variation in body size in Drosophila melanogaster

Life-History Evolution and the Genetics of Fitness Components in Drosophila melanogaster

Life-history traits or "fitness components"-such as age and size at maturity, fecundity and fertility, age-specific rates of survival, and life span-are the major phenotypic determinants of Darwinian fitness. Analyzing the evolution and genetics of these phenotypic targets of selection is central to our understanding of adaptation. Due to its simple and rapid life cycle, cosmopolitan distribution, ease of maintenance in the laboratory, well-understood evolutionary genetics, and its versatile genetic toolbox, the "vinegar fly" Drosophila melanogaster is one of the most powerful, experimentally tractable model systems for studying "life-history evolution." Here, I review what has been learned about the evolution and genetics of life-history variation in D. melanogaster by drawing on numerous sources spanning population and quantitative genetics, genomics, experimental evolution, evolutionary ecology, and physiology. This body of work has contributed greatly to our knowledge of several fundamental problems in evolutionary biology, including the amount and maintenance of genetic variation, the evolution of body size, clines and climate adaptation, the evolution of senescence, phenotypic plasticity, the nature of life-history trade-offs, and so forth. While major progress has been made, important facets of these and other questions remain open, and the D. melanogaster system will undoubtedly continue to deliver key insights into central issues of life-history evolution and the genetics of adaptation.

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

Genetic basis of thermal plasticity variation in Drosophila melanogaster body size

Body size is a quantitative trait that is closely associated to fitness and under the control of both genetic and environmental factors. While developmental plasticity for this and other traits is heritable and under selection, little is known about the genetic basis for variation in plasticity that can provide the raw material for its evolution. We quantified genetic variation for body size plasticity in Drosophila melanogaster by measuring thorax and abdomen length of females reared at two temperatures from a panel representing naturally segregating alleles, the Drosophila Genetic Reference Panel (DGRP). We found variation between genotypes for the levels and direction of thermal plasticity in size of both body parts. We then used a Genome-Wide Association Study (GWAS) approach to unravel the genetic basis of inter-genotype variation in body size plasticity, and used different approaches to validate selected QTLs and to explore potential pleiotropic effects. We found mostly “private QTLs”, with little overlap between the candidate loci underlying variation in plasticity for thorax versus abdomen size, for different properties of the plastic response, and for size versus size plasticity. We also found that the putative functions of plasticity QTLs were diverse and that alleles for higher plasticity were found at lower frequencies in the target population. Importantly, a number of our plasticity QTLs have been targets of selection in other populations. Our data sheds light onto the genetic basis of inter-genotype variation in size plasticity that is necessary for its evolution.

Environmental conditions can influence development and lead to the production of phenotypes adjusted to the conditions adults will live in. This developmental plasticity, which can help organisms cope with environmental heterogeneity, is heritable and under selection. Its evolution will depend on available genetic variation. Using a panel of D. melanogaster flies representing naturally segregating alleles, we identified DNA sequence variants associated to variation in thermal plasticity for body size. We found that these variants correspond to a diverse set of functions and that their effects differ between body parts and properties of the thermal response. Our results shed new light onto the long discussed genes for plasticity.

Genomic regions and genes related to inter-population differences in body size in the ground beetle Carabus japonicus

Body size is a key trait in diversification among animal species, and revealing the gene regions responsible for body size diversification among populations or related species is important in evolutionary biology. We explored the genomic regions associated with body size differences in Carabus japonicus ground beetle populations by quantitative trait locus (QTL) mapping of F2 hybrids from differently sized parents from two populations using restriction site-associated DNA sequencing and de novo assembly of the beetle whole genome. The assembled genome had a total length of 191 Mb with a scaffold N50 of 0.73 Mb; 14,929 protein-coding genes were predicted. Three QTLs on different linkage groups had major effects on the overall size, which is composed chiefly of elytral length. In addition, we found QTLs on autosomal and X chromosomal linkage groups that affected head length and width, thoracic width, and elytral width. We determined the gene loci potentially related to control of body size in scaffolds of the genome sequence, which contained the QTL regions. The genetic basis of body size variation based on a small number of major loci would promote differentiation in body size in response to selection pressures related to variations in environmental conditions and inter-specific interactions.

Quantitative Trait Loci and Antagonistic Associations for Two Developmentally Related Traits in the Drosophila Head

In insects, some developmentally related traits are negatively correlated. Here, we mapped Quantitative Trait Loci (QTL) for traits of eye size and head capsule, in an intercontinental set of recombinant inbred lines (RILs) of Drosophila melanogaster Composite interval mapping identified QTL on all major chromosomes. Two negatively correlated traits (size of eyes and between-eyes distance) were influenced by one QTL that appeared to be antagonistic between the traits (QTL cytological range is 25F5-30A6), consistent with a negative genetic correlation between these traits of the head capsule. Comparisons of QTL across traits indicated a nonrandom distribution over the genome, with a considerable overlap between some QTL across traits. Developmentally-related traits were influenced by QTL in a pattern that is consistent both with 1) the sign of the genetic correlation between the traits and 2) a constraint in the micro-evolutionary differentiation in the traits.

Natural Genetic Variation and Candidate Genes for Morphological Traits in Drosophila melanogaster

Body size is a complex character associated to several fitness related traits that vary within and between species as a consequence of environmental and genetic factors. Latitudinal and altitudinal clines for different morphological traits have been described in several species of Drosophila and previous work identified genomic regions associated with such variation in D. melanogaster. However, the genetic factors that orchestrate morphological variation have been barely studied. Here, our main objective was to investigate genetic variation for different morphological traits associated to the second chromosome in natural populations of D. melanogaster along latitudinal and altitudinal gradients in Argentina. Our results revealed weak clinal signals and a strong population effect on morphological variation. Moreover, most pairwise comparisons between populations were significant. Our study also showed important within-population genetic variation, which must be associated to the second chromosome, as the lines are otherwise genetically identical. Next, we examined the contribution of different candidate genes to natural variation for these traits. We performed quantitative complementation tests using a battery of lines bearing mutated alleles at candidate genes located in the second chromosome and six second chromosome substitution lines derived from natural populations which exhibited divergent phenotypes. Results of complementation tests revealed that natural variation at all candidate genes studied, invected, Fasciclin 3, toucan, Reticulon-like1, jing and CG14478, affects the studied characters, suggesting that they are Quantitative Trait Genes for morphological traits. Finally, the phenotypic patterns observed suggest that different alleles of each gene might contribute to natural variation for morphological traits. However, non-additive effects cannot be ruled out, as wild-derived strains differ at myriads of second chromosome loci that may interact epistatically with mutant alleles.

Partially repeatable genetic basis of benthic adaptation in threespine sticklebacks

The extent to which convergent adaptation to similar ecological niches occurs by a predictable genetic basis remains a fundamental question in biology. Threespine stickleback fish have undergone an adaptive radiation in which ancestral oceanic populations repeatedly colonized and adapted to freshwater habitats. In multiple lakes in British Columbia, two different freshwater ecotypes have evolved: a deep-bodied benthic form adapted to forage near the lake substrate, and a narrow-bodied limnetic form adapted to forage in open water. Here, we use genome-wide linkage mapping in marine × benthic F2 genetic crosses to test the extent of shared genomic regions underlying benthic adaptation in three benthic populations. We identify at least 100 Quantitative Trait Loci (QTL) harboring genes influencing skeletal morphology. The majority of QTL (57%) are unique to one cross. However, four genomic regions affecting eight craniofacial and armor phenotypes are found in all three benthic populations. We find that QTL are clustered in the genome and overlapping QTL regions are enriched for genomic signatures of natural selection. These findings suggest that benthic adaptation has occurred via both parallel and nonparallel genetic changes.

Parallel effects of the inversion In(3R)Payne on body size across the North American and Australian clines in Drosophila melanogaster

Chromosomal inversions are thought to play a major role in climatic adaptation. In D. melanogaster, the cosmopolitan inversion In(3R)Payne exhibits latitudinal clines on multiple continents. As many fitness traits show similar clines, it is tempting to hypothesize that In(3R)P underlies observed clinal patterns for some of these traits. In support of this idea, previous work in Australian populations has demonstrated that In(3R)P affects body size but not development time or cold resistance. However, similar data from other clines of this inversion are largely lacking; finding parallel effects of In(3R)P across multiple clines would considerably strengthen the case for clinal selection. Here, we have analysed the phenotypic effects of In(3R)P in populations originating from the endpoints of the latitudinal cline along the North American east coast. We measured development time, egg-to-adult survival, several size-related traits (femur and tibia length, wing area and shape), chill coma recovery, oxidative stress resistance and triglyceride content in homokaryon lines carrying In(3R)P or the standard arrangement. Our central finding is that the effects of In(3R)P along the North American cline match those observed in Australia: standard arrangement lines were larger than inverted lines, but the inversion did not influence development time or cold resistance. Similarly, In(3R)P did not affect egg-to-adult survival, oxidative stress resistance and lipid content. In(3R)P thus seems to specifically affect size traits in populations from both continents. This parallelism strongly suggests an adaptive pattern, whereby the inversion has captured alleles associated with growth regulation and clinal selection acts on size across both continents.

Genomic Evidence for Adaptive Inversion Clines in Drosophila melanogaster

Clines in chromosomal inversion polymorphisms-presumably driven by climatic gradients-are common but there is surprisingly little evidence for selection acting on them. Here we address this long-standing issue in Drosophila melanogaster by using diagnostic single nucleotide polymorphism (SNP) markers to estimate inversion frequencies from 28 whole-genome Pool-seq samples collected from 10 populations along the North American east coast. Inversions In(3L)P, In(3R)Mo, and In(3R)Payne showed clear latitudinal clines, and for In(2L)t, In(2R)NS, and In(3R)Payne the steepness of the clinal slopes changed between summer and fall. Consistent with an effect of seasonality on inversion frequencies, we detected small but stable seasonal fluctuations of In(2R)NS and In(3R)Payne in a temperate Pennsylvanian population over 4 years. In support of spatially varying selection, we observed that the cline in In(3R)Payne has remained stable for >40 years and that the frequencies of In(2L)t and In(3R)Payne are strongly correlated with climatic factors that vary latitudinally, independent of population structure. To test whether these patterns are adaptive, we compared the amount of genetic differentiation of inversions versus neutral SNPs and found that the clines in In(2L)t and In(3R)Payne are maintained nonneutrally and independent of admixture. We also identified numerous clinal inversion-associated SNPs, many of which exhibit parallel differentiation along the Australian cline and reside in genes known to affect fitness-related traits. Together, our results provide strong evidence that inversion clines are maintained by spatially-and perhaps also temporally-varying selection. We interpret our data in light of current hypotheses about how inversions are established and maintained.

Genome-Wide Analysis Reveals Novel Regulators of Growth in Drosophila melanogaster

Organismal size depends on the interplay between genetic and environmental factors. Genome-wide association (GWA) analyses in humans have implied many genes in the control of height but suffer from the inability to control the environment. Genetic analyses in Drosophila have identified conserved signaling pathways controlling size; however, how these pathways control phenotypic diversity is unclear. We performed GWA of size traits using the Drosophila Genetic Reference Panel of inbred, sequenced lines. We find that the top associated variants differ between traits and sexes; do not map to canonical growth pathway genes, but can be linked to these by epistasis analysis; and are enriched for genes and putative enhancers. Performing GWA on well-studied developmental traits under controlled conditions expands our understanding of developmental processes underlying phenotypic diversity.

Genetic studies in Drosophila have elucidated conserved signaling pathways and environmental factors that together control organismal size. In humans, hundreds of genes are associated with height variation, but these associations have not been performed in a controlled environment. As a result we are still lacking an understanding of the mechanisms creating size variability within a species. Here, under carefully controlled environmental conditions, we identify naturally occurring genetic variants that are associated with size diversity in Drosophila. We identify a cluster of associations close to the kek1 locus, a well-characterized growth regulator, but otherwise find that most variants are located in or close to genes that do not belong to the conserved pathways but may interact with these in a biological network. We validate 33 novel growth regulatory genes that participate in diverse cellular processes, most notably cellular metabolism and cell polarity. This study is the first genome-wide association analysis of natural variants underlying size in Drosophila and our results complement the knowledge we have accumulated on this trait from mutational studies of single genes.

The FlyCatwalk: a high-throughput feature-based sorting system for artificial selection in Drosophila

Experimental evolution is a powerful tool for investigating complex traits. Artificial selection can be applied for a specific trait and the resulting phenotypically divergent populations pool-sequenced to identify alleles that occur at substantially different frequencies in the extreme populations. To maximize the proportion of loci that are causal to the phenotype among all enriched loci, population size and number of replicates need to be high. These requirements have, in fact, limited evolution studies in higher organisms, where the time investment required for phenotyping is often prohibitive for large-scale studies. Animal size is a highly multigenic trait that remains poorly understood, and an experimental evolution approach may thus aid in gaining new insights into the genetic basis of this trait. To this end, we developed the FlyCatwalk, a fully automated, high-throughput system to sort live fruit flies (Drosophila melanogaster) based on morphometric traits. With the FlyCatwalk, we can detect gender and quantify body and wing morphology parameters at a four-old higher throughput compared with manual processing. The phenotyping results acquired using the FlyCatwalk correlate well with those obtained using the standard manual procedure. We demonstrate that an automated, high-throughput, feature-based sorting system is able to avoid previous limitations in population size and replicate numbers. Our approach can likewise be applied for a variety of traits and experimental settings that require high-throughput phenotyping.

Complexity of the genetic basis of ageing in nature revealed by a clinal study of lifespan and methuselah, a gene for ageing, in Drosophila from eastern Australia

Clinal studies are a powerful tool for understanding the genetic basis of climatic adaptation. However, while clines in quantitative traits and genetic polymorphisms have been observed within and across continents, few studies have attempted to demonstrate direct links between them. The gene methuselah in Drosophila has been shown to have a major effect on stress response and longevity phenotypes based largely on laboratory studies of induced mutations in the mth gene. Clinal patterns in the most common mth haplotype and for lifespan (both increasing with latitude) have been observed in North American populations of D. melanogaster, implicating climatic selection. While these clinal patterns have led some to suggest that mth influences ageing in natural populations, limited evidence on the association between the two has so far been collected. Here, we describe a significant cline in the mth haplotype in eastern Australian D. melanogaster populations that parallel the cline in North America. We also describe a cline in mth gene expression. These findings further support the idea that mth is itself under selection. In contrast, we show that lifespan has a strong nonlinear clinal pattern, increasing southwards from the tropics, but then decreasing again from mid-latitudes. Furthermore, in association studies, we find no evidence for a direct link between mth haplotype and lifespan. Thus, while our data support a role for mth variation being under natural selection, we found no link to naturally occurring variation in lifespan and ageing in Australian populations of D. melanogaster. Our results indicate that the mth locus likely has genetic background and environment-specific effects.

Primary genome scan for complex body shape-related traits in the common carp Cyprinus carpio

To identify quantitative trait loci (QTL) that affect body shape in common carp Cyprinus carpio, a linkage map, 2159·23 cM long, was constructed with a total of 307 markers covering 51 linkage groups (LG). The map included 167 new single nucleotide polymorphism (SNP) markers derived from expressed sequence tags (EST) together with 140 microsatellite markers reported earlier. A primary genome scan was conducted for QTL for standard length (L(S)), head length (L(H)), body height (H(B)), body width (W(B)) and tail length (L(TAIL)) in an F1 line containing 92 offspring. A total of 15 suggestive QTL on six LGs were found to associate with L(S), L(H), H(B), W(B) and L(TAIL) which explained 10·7-17·4% of the variance. Five significant QTL were detected for body-shape related traits and located for LGs (lg1, 12 and 20). These QTL included: one associated with L(S) (21·1% variance explained), three for H(B) (almost 20% variance explained) and one for W(B) (20·7% variance explained).

The genetic architecture of growth rate in juvenile Takifugu species

Closely related species have often evolved dramatic differences in body size. Takifugu rubripes (fugu) is a large marine pufferfish whose genome has been sequenced, whereas T. niphobles is the smallest species among Takifugu. We show that, unsurprisingly, the juvenile growth rate of T. rubripes is higher than that of T. niphobles in a laboratory setting. We produced F(2) progenies of their F(1) hybrids and found one quantitative trait locus (QTL) significantly associated with variation in juvenile body size. This QTL region (3.5 Mb) contains no known genes directly related to growth phenotype (such as IGFs) except Fgf21, which inhibits growth hormone signaling in mouse. The QTL in Takifugu spp. is distinct from the region previously known to control body size variations in stickleback or tilapia. Our results suggest that in the fish tested herein, genomic regions underlying body size evolution might have different genetic origins. They also suggest that many diverse traits in Takifugu spp. are amenable to genetic mapping.

Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Genome-wide transcription analysis of clinal genetic variation in Drosophila

Clinal variation in quantitative traits is widespread, but its genetic basis awaits identification. Drosophila melanogaster shows adaptive, clinal variation in traits such as body size along latitudinal gradients on multiple continents. To investigate genome wide transcription differentiation between North and South that might contribute to the clinal phenotypic variation, we compared RNA expression patterns during development of D. melanogaster from tropical northern and temperate southern populations using whole genome tiling arrays. We found that genes that were differentially expressed between the cline ends were generally associated with metabolism and growth, and experimental alteration of expression of a sample of them generally resulted in altered body size in the predicted direction, sometimes significantly so. We further identified the serpent (srp) transcription factor binding sites to be enriched near genes up-regulated in expression in the south. Analysis of clinal populations revealed a significant cline in the expression level of srp. Experimental over-expression of srp increased body size, as predicted from its clinal expression pattern, suggesting that it may be involved in regulating adaptive clinal variation in Drosophila. This study identified a handful of genes that contributed to clinal phenotypic variation through altered gene expression level, yet misexpression of individual gene led to modest body size change.

Molecular basis of adaptive shift in body size in Drosophila melanogaster: functional and sequence analyses of the Dca gene

Latitudinal body size clines in animals conforming to Bergmann's rule occur on many continents but isolating their underlying genetic basis remains a challenge. In Drosophila melanogaster, the gene Dca accounts for approximately 5-10% of the natural wing size variation (McKechnie SW, Blacket MJ, Song SV, Rako L, Carroll X, Johnson TK, Jensen LT, Lee SF, Wee CW, Hoffmann AA. 2010. A clinally varying promoter polymorphism associated with adaptive variation in wing size in Drosophila. Mol Ecol. 19:775-784). We present here functional evidence that Dca is a negative regulator of wing size. A significant negative latitudinal cline of Dca gene expression was detected in synchronized third instar larvae. In addition, we clarified the evolutionary history of the three most common Dca promoter alleles (Dca237-1, Dca237-2, and Dca247) and showed that the insertion allele (Dca247), whose frequency increases with latitude, is associated with larger wing centroid size and higher average cell number in male flies. Finally, we showed that the overall linkage disequilibrium (LD) was low in the Dca promoter and that the insertion/deletion polymorphism that defines the Dca alleles was in strong LD with two other upstream sites. Our results provide strong support that Dca is a candidate for climatic adaptation in D. melanogaster.

Genetic mapping of adaptive wing size variation in Drosophila simulans

Many ecologically important traits exhibit latitudinal variation. Body size clines have been described repeatedly in insects across multiple continents, suggesting that similar selective forces are shaping these geographical gradients. It is unknown whether these parallel clinal patterns are controlled by the same or different genetic mechanism(s). We present here, quantitative trait loci (QTL) analysis of wing size variation in Drosophila simulans. Our results show that much of the wing size variation is controlled by a QTL on Chr 3L with relatively minor contribution from other chromosome arms. Comparative analysis of the genomic positions of the QTL indicates that the major QTL on Chr 3 are distinct in D. simulans and D. melanogaster, whereas the QTL on Chr 2R might overlap between species. Our results suggest that parallel evolution of wing size clines could be driven by non-identical genetic mechanisms but in both cases involve a major QTL as well as smaller effects of other genomic regions.

Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura

In fewer than two decades after invading the Americas, the fly Drosophila subobscura evolved latitudinal clines for chromosomal inversion frequencies and wing size that are parallel to the long-standing ones in native Palearctic populations. By sharp contrast, wing shape clines also evolved in the New World, but the relationship with latitude was opposite to that in the Old World. Previous work has suggested that wing trait differences among individuals are partially due to the association between chromosomal inversions and particular alleles which influence the trait under consideration. Furthermore, it is well documented that a few number of effective individuals founded the New World populations, which might have modified the biometrical effect of inversions on quantitative traits. Here we evaluate the relative contribution of chromosomal inversion clines in shaping the parallel clines in wing size and contrasting clines in wing shape in native and colonizing populations of the species. Our results reveal that inversion-size and inversion-shape associations in native and colonizing (South America) populations are generally different, probably due to the bottleneck effect. Contingent, unpredictable evolution was suggested as an explanation for the different details involved in the otherwise parallel wing size clines between Old and New World populations of D. subobscura. We challenge this assertion and conclude that contrasting wing shape clines came out as a correlated response of inversion clines that might have been predicted considering the genetic background of colonizers.

Positive selection has driven the evolution of the Drosophila insulin-like receptor (InR) at different timescales

The highly conserved insulin-signaling pathway influences very diverse processes including intermediary metabolism, reproduction, aging, and growth. The first pathway component is the insulin receptor that upon insulin binding triggers the signal-transduction cascade. Its variation, like that of other pathway components, might therefore affect many organismal traits. Variation at the "Insulin-like" receptor (InR) gene was surveyed both within Drosophila melanogaster and between species across the Drosophila phylogeny. In D. melanogaster, the level and pattern of variation at the approximately 8-kb region surveyed did not provide any indication of a recent selective event in this region. Maximum likelihood (ML) analyses revealed the past action of purifying selection acting differentially both across the phylogeny and along the studied gene. Moreover, the ML analyses and the McDonald and Kreitman test revealed the footprint of positive selection driving amino acid changes to fixation in the branch separating the Sophophora and the Drosophila subgenera, and in the D. melanogaster lineage, respectively. The oldest selective events could have affected either the insulin binding or the signal-transduction capacities of the receptor, whereas mutations affecting signal transduction would seem to underlie the more recent events.

QTL for the thermotolerance effect of heat hardening, knockdown resistance to heat and chill-coma recovery in an intercontinental set of recombinant inbred lines of Drosophila melanogaster

The thermotolerance effect of heat hardening (also called short-term acclimation), knockdown resistance to high temperature (KRHT) with and without heat hardening and chill-coma recovery (CCR) are important phenotypes of thermal adaptation in insects and other organisms. Drosophila melanogaster from Denmark and Australia were previously selected for low and high KRHT, respectively. These flies were crossed to construct recombinant inbred lines (RIL). KRHT was higher in heat-hardened than in nonhardened RIL. We quantify the heat-hardening effect (HHE) as the ratio in KRHT between heat-hardened and nonhardened RIL. Composite interval mapping revealed a more complex genetic architecture for KRHT without heat-hardening than for KRHT in heat-hardened insects. Five quantitative trait loci (QTL) were found for KRHT, but only two of them were significant after heat hardening. KRHT and CCR showed trade-off associations for QTL both in the middle of chromosome 2 and the right arm of chromosome 3, which should be the result of either pleiotropy or linkage. The major QTL on chromosome 2 explained 18% and 27-33% of the phenotypic variance in CCR and KRHT in nonhardened flies, respectively, but its KRHT effects decreased by heat hardening. We discuss candidate loci for each QTL. One HHE-QTL was found in the region of small heat-shock protein genes. However, HHE-QTL explained only a small fraction of the phenotypic variance. Most heat-resistance QTL did not colocalize with CCR-QTL. Large-effect QTL for CCR and KRHT without hardening (basal thermotolerance) were consistent across continents, with apparent transgressive segregation for CCR. HHE (inducible thermotolerance) was not regulated by large-effect QTL.

Autosomal variation for male body size and sperm competition phenotypes is uncorrelated in Drosophila melanogaster

Correlations between male body size and phenotypes impacting post-copulatory sexual selection are commonly observed during the manipulation of male body size by environmental rearing conditions. Here, we control for environmental influences and test for genetic correlations between natural variation in male body size and phenotypes affecting post-copulatory sexual selection in Drosophila melanogaster. Dry weights of virgin males from 90 second-chromosome and 88 third-chromosome substitution lines were measured. Highly significant line effects (p<0.001) documented a genetic basis to variation in male body size. No significant correlations were identified between male body size and the components of sperm competitive ability. These results suggest that natural autosomal variation for male body size has little impact on post-copulatory sexual selection. If genetic correlations exist between male body size and post-copulatory sexual selection then variation in the sex chromosomes are likely candidates, as might be expected if sexually antagonistic coevolution was responsible.

A major QTL affects temperature sensitive adult lethality and inbreeding depression in life span in Drosophila melanogaster

BACKGROUND

The study of inbreeding depression has major relevance for many disciplines, including conservation genetics and evolutionary biology. Still, the molecular genetic basis of this phenomenon remains poorly characterised, as knowledge on the mechanistic causes of inbreeding depression and the molecular properties of genes that give rise to or modulate its deleterious effects is lacking. These questions warrant the detailed study of genetic loci giving rise to inbreeding depression. However, the complex and polygenic nature of general inbreeding depression makes this a daunting task. Study of inbreeding effects in specific traits, such as age-specific mortality and life span, provide a good starting point, as a limited set of genes is expected to be involved.

RESULTS

Here we report on a QTL mapping study on inbreeding related and temperature sensitive lethality in male Drosophila melanogaster. The inbreeding effect was expressed at moderately high temperature, and manifested itself as severe premature mortality in males, but not in females. We used a North Carolina crossing design 3 to estimate average dominance ratio and heritability. We found the genetic basis of the lethal effect to be relatively simple, being due mainly to a single recessive QTL on the left arm of chromosome 2. This locus colocalised with a QTL that conditioned variation in female life span, acting as an overdominant locus for this trait. Male life span was additionally affected by variation at the X-chromosome.

CONCLUSION

This demonstrates that analysis of large conditional lethal effects is a viable strategy for delineating genes which are sensitive to inbreeding depression.

A population genomic approach to map recent positive selection in model species

Based on nearly complete genome sequences from a variety of organisms data on naturally occurring genetic variation on the scale of hundreds of loci to entire genomes have been collected in recent years. In parallel, new statistical tests have been developed to infer evidence of recent positive selection from these data and to localize the target regions of selection in the genome. These methods have now been successfully applied to Drosophila melanogaster, humans, mice and a few plant species. In genomic regions of normal recombination rates, the targets of positive selection have been mapped down to the level of individual genes.

Investigating latitudinal clines for life history and stress resistance traits in Drosophila simulans from eastern Australia

Latitudinal clines have been demonstrated for many quantitative traits in Drosophila and are assumed to be due to climatic selection. However, clinal studies are often performed in species of Drosophila that contain common cosmopolitan inversion polymorphisms that also show clinal patterns. These inversion polymorphisms may be responsible for much of the observed clinal variation. Here, we consider latitudinal clines for quantitative traits in Drosophila simulans from eastern Australia. Drosophila simulans does not contain cosmopolitan inversion polymorphisms, so allows the study of clinal selection on quantitative traits that are not confounded by associations with inversions. Body size showed a strong linear cline for both females and males. Starvation resistance exhibited a weak linear cline in females, whereas chill-coma recovery exhibited a significant nonlinear cline in females only. No clinal pattern was evident for development time, male chill-coma recovery, desiccation or heat resistance. We discuss these results with reference to the role inversion polymorphisms play in generating clines in quantitative traits of Drosophila.

QTL mapping of inbreeding-related cold sensitivity and conditional lethality in Drosophila melanogaster

Inbreeding depression is a central theme within genetics, and is of specific interest for researchers within evolutionary and conservation genetics and animal and plant breeding. Inbreeding effects are thought to be caused by the joint expression of conditional and unconditional deleterious alleles. Whenever the expression of deleterious alleles is conditional, this can result in extreme environmental sensitivity in certain inbred lineages. Analysis of conditional lethal effects can reveal some of the loci that are sensitive to inbreeding. We performed a QTL (quantitative trait locus) mapping study of inbreeding-related and conditionally expressed lethality in Drosophila melanogaster. The lethal effect was triggered by exposure to a cold shock. We used a North Carolina crossing Design 3 to establish the mapping population, as well as to estimate the average dominance ratio and heritability. We found two QTL on the major autosomes carrying recessive lethals that caused male mortality, one of which also affected female mortality. More detailed study of these loci will provide information on the mechanistic basis and environmental sensitivity of inbreeding depression.

Genomic analysis of adaptive differentiation in Drosophila melanogaster

Drosophila melanogaster shows clinal variation along latitudinal transects on multiple continents for several phenotypes, allozyme variants, sequence variants, and chromosome inversions. Previous investigation suggests that many such clines are due to spatially varying selection rather than demographic history, but the genomic extent of such selection is unknown. To map differentiation throughout the genome, we hybridized DNA from temperate and subtropical populations to Affymetrix tiling arrays. The dense genomic sampling of variants and low level of linkage disequilibrium in D. melanogaster enabled identification of many small, differentiated regions. Many regions are differentiated in parallel in the United States and Australia, strongly supporting the idea that they are influenced by spatially varying selection. Genomic differentiation is distributed nonrandomly with respect to gene function, even in regions differentiated on only one continent, providing further evidence for the role of selection. These data provide candidate genes for phenotypes known to vary clinally and implicate interesting new processes in genotype-by-environment interactions, including chorion proteins, proteins regulating meiotic recombination and segregation, gustatory and olfactory receptors, and proteins affecting synaptic function and behavior. This portrait of differentiation provides a genomic perspective on adaptation and the maintenance of variation through spatially varying selection.

Patterns of sequence variability and divergence at the diminutive gene region of Drosophila melanogaster: complex patterns suggest an ancestral selective sweep

To identify putatively swept regions of the Drosophila melanogaster genome, we performed a microsatellite screen spanning a 260-kb region of the X chromosome in populations from Zimbabwe, Ecuador, the United States, and China. Among the regions identified by this screen as showing a complex pattern of reduced heterozygosity and a skewed frequency spectrum was the gene diminutive (dm). To investigate the microsatellite findings, nucleotide sequence polymorphism data were generated in populations from both China and Zimbabwe spanning a 25-kb region and encompassing dm. Analysis of the sequence data reveals strongly reduced nucleotide variation across the entire gene region in both the non-African and the African populations, an extended haplotype pattern, and structured linkage disequilibrium, as well as a rejection of neutrality in favor of selection using a composite likelihood-ratio test. Additionally, unusual patterns of synonymous site evolution were observed at the second exon of this locus. On the basis of simulation studies as well as recently proposed methods for distinguishing between selection and nonequilibrium demography, we find that this "footprint" is best explained by a selective sweep in the ancestral population, the signal of which has been somewhat blurred via founder effects in the non-African samples.

Candidate genes and thermal phenotypes: identifying ecologically important genetic variation for thermotolerance in the Australian Drosophila melanogaster cline

Clinal variation in traits often reflects climatic adaptation; in Drosophila melanogaster clinal variation provides an opportunity to link variation in chromosomal inversions, microsatellite loci and various candidate genes to adaptive variation in traits. We undertook association studies with crosses from a single population of D. melanogaster from eastern Australia to investigate the association between genetic markers and traits showing clinal variation. By genotyping parents and phenotyping offspring, we minimized genotyping costs but had the power to detect association between markers and quantitative traits. Consistent with prior studies, we found strong associations between the clinal chromosomal inversion In(3R)Payne and markers within it, as well as among these markers. We also found an association between In(3L)Payne and one marker located within this inversion. Of the five predicted associations between markers and traits, four were detected (increased heat, decreased cold resistance and body size with the heat shock gene hsr-omega S, increased cold resistance with the inversion In(3L)Payne), while one was not detected (heat resistance and the heat shock gene hsp68). In a set of eight exploratory tests, we detected one positive association (between hsp23a and heat resistance) but no associations of heat resistance with alleles at the hsp26, hsp83, Desat 2, alpha-Gpdh, hsp70 loci, while cold resistance was not associated with Frost and Dca loci. These results confirm interactions between hsr-omega and thermal resistance, as well as between In(3L)Payne and cold resistance, but do not provide evidence for associations between thermal responses and alleles at other clinically varying marker genes.

X-linked QTL for knockdown resistance to high temperature in Drosophila melanogaster

Knockdown Resistance to High Temperature (KRHT) is an adaptive trait of thermotolerance in insects. An interval mapping was performed on chromosome X of Drosophila melanogaster to search for quantitative trait loci (QTL) affecting KRHT. A backcross population was obtained from two lines that dramatically differ for KRHT. Microsatellites were used as markers. Composite interval mapping identified a large-effect QTL in the region of band 10 where putative candidate genes map. To further test for this QTL a set of recombinant (but non-inbred) lines was obtained from backcrosses between the parental lines used for the interval mapping. Recombinant line analysis confirmed that one QTL is targeted by band 10. We identify and discuss candidate loci contained within our QTL region.

Mapping regions within cosmopolitan inversion In(3R)Payne associated with natural variation in body size in Drosophila melanogaster

Associations between genotypes for inversions and quantitative traits have been reported in several organisms, but little has been done to localize regions within inversions controlling variation in these traits. Here, we use an association mapping technique to identify genomic regions controlling variation in wing size within the cosmopolitan inversion In(3R)Payne in Drosophila melanogaster. Previous studies have shown that this inversion strongly influences variation in wing size, a trait highly correlated with body size. We found three alleles from two separate regions within In(3R)Payne with significant additive effects on wing size after the additional effect of the inversion itself had been taken into account. There were also several alleles with significant genotype-by-inversion interaction effects on wing size. None of the alleles tested had a significant additive effect on development time, suggesting different genes control these traits and that clinal patterns in them have therefore arisen independently. The presence of multiple regions within In(3R)Payne controlling size is consistent with the idea that inversions persist in populations because they contain multiple sets of locally adapted alleles, but more work needs to be done to test if they are indeed coadapted.

Knockdown resistance to heat stress and slow recovery from chill coma are genetically associated in a quantitative trait locus region of chromosome 2 inDrosophila melanogaster

In insects, two ecologically relevant traits of thermal adaptation are knockdown resistance to high temperature (KRHT) and chill-coma recovery (CCR). Chromosome 2 of Drosophila melanogaster was tested for quantitative trait loci (QTL) affecting both CCR and KRHT in backcrosses between homosequential lines that are fixed for the standard (noninverted) sequence of this autosome. These lines were obtained by artificial selection on KRHT and subsequent inbreeding from a stock that was derived from a single wild population. Heat-induced expression of the 70KD heat-shock protein (Hsp70) was also examined for variation between the lines. Composite interval mapping was performed for each trait on each reciprocal backcross, identifying one QTL region in the middle of chromosome 2 for both KRHT and CCR. The largest estimates of additive effects were found in pericentromeric regions of chromosome 2, accounting for 10-14% (CCR) and 10-17% (KRHT) of the phenotypic variance in BC populations. No QTL was found in the region of the heat-shock factor (hsf) gene. However, the two parental lines have diverged in the heat-induced Hsp70 expression. Distribution of KRHT QTL on chromosome 2 was similar between this study based on crosses between lines selected from a single wild population and previous work based on crosses between selection lines from different continents. Colocalized QTL showed a trade-off association between CCR and KRHT, which should be the result of either multiple, tightly linked trait-specific genes or a single gene with pleiotropic effects on the traits. We discuss candidate loci contained within the QTL regions.

The association between inversion In(3R)Payne and clinally varying traits in Drosophila melanogaster

In Drosophila melanogaster, inversion In(3R)Payne increases in frequency towards low latitudes and has been putatively associated with variation in size and thermal resistance, traits that also vary clinally. To assess the association between size and inversion, we obtained isofemale lines of inverted and standard karyotype of In(3R)Payne from the ends of the Australian D. melanogaster east coast cline. In the northern population, there was a significant association between In(3R)Payne and body size, with standard lines from this population being relatively larger than inverted lines. In contrast, the inversion had no influence on development time or cold resistance. We strengthened our findings further in a separate study with flies from populations from the middle of the cline as well as from the cline ends. These flies were scored for wing size and the presence of In(3R)Payne using a molecular marker. In females, the inversion accounted for around 30% of the size difference between cline ends, while in males the equivalent figure was 60%. Adaptive shifts in size but not in the other traits are therefore likely to have involved genes closely associated with In(3R)Payne. Because the size difference between karyotypes was similar in different populations, there was no evidence for coadaptation within populations.

Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia

Drosophila melanogaster invaded Australia around 100 years ago, most likely through a northern invasion. The wide range of climatic conditions in eastern Australia across which D. melanogaster is now found provides an opportunity for researchers to identify traits and genes that are associated with climatic adaptation. Allozyme studies indicate clinal patterns for at least four loci including a strong linear cline in Adh and a non-linear cline in alpha-Gpdh. Inversion clines were initially established from cytological studies but have now been validated with larger sample sizes using molecular markers for breakpoints. Recent collections indicate that some genetic markers (Adh and In(3R)Payne) have changed over the last 20 years reflecting continuing evolution. Heritable clines have been established for quantitative traits including wing length/area, thorax length and cold and heat resistance. A cline in egg size independent of body size and a weak cline in competitive ability have also been established. Postulated clinal patterns for resistance to desiccation and starvation have not been supported by extensive sampling. Experiments under laboratory and semi-natural conditions have suggested selective factors generating clinal patterns, particularly for reproductive patterns over winter. Attempts are being made to link clinal variation in traits to specific genes using QTL analysis and the candidate locus approach, and to identify the genetic architecture of trait variation along the cline. This is proving difficult because of inversion polymorphisms that generate disequilibrium among genes. Substantial gaps still remain in linking clines to field selection and understanding the genetic and physiological basis of the adaptive shifts. However D. melanogaster populations in eastern Australia remain an excellent resource for understanding past and future evolutionary responses to climate change.

Altitudinal patterns for latitudinally varying traits and polymorphic markers in Drosophila melanogaster from eastern Australia

Altitudinal changes in traits and genetic markers can complement the studies on latitudinal patterns and provide evidence of natural selection because of climatic factors. In Drosophila melanogaster, latitudinal variation is well known but altitudinal patterns have rarely been investigated. Here, we examine five traits and five genetic markers on chromosome 3R in D. melanogaster collected at high and low altitudes from five latitudes along the eastern coast of Australia. Significant altitudinal differentiation was observed for cold tolerance, development time, ovariole number in unmated females, and the microsatellite marker DMU25686. Differences tended to match latitudinal patterns, in that trait values at high altitudes were also found at high latitudes, suggesting that factors linked to temperature are likely selective agents. Cold tolerance was closely associated with average temperature and other climatic factors, but no significant associations were detected for the other traits. Genes around DMU25686 represent good candidates for climatic adaptation.

Patterns of diversity and linkage disequilibrium within the cosmopolitan inversion In(3R)Payne in Drosophila melanogaster are indicative of coadaptation

The cosmopolitan inversion In(3R)Payne in Drosophila melanogaster decreases in frequency with increasing distance from the equator on three continents, indicating it is subject to strong natural selection. We investigated patterns of genetic variation and linkage disequilibrium (LD) in 24 molecular markers located within and near In(3R)Payne to determine if different parts of the inversion responded to selection the same way. We found reduced variation in the markers we used compared to others distributed throughout the genome, consistent with the inversion having a relatively recent origin (<N(e) generations). LD between markers and In(3R)Payne varied significantly among markers within the inversion, with regions of high association interspersed by regions of low association. Several factors indicate that these patterns were not due to demographic factors such as admixture and bottlenecks associated with colonization, but instead reflected strong epistatic selection. Furthermore, we found that nonadjacent regions with high association to the inversion contained markers with the strongest clinal patterns in allele frequency; in most cases, the level of clinal variation was beyond what could be explained by hitchhiking with In(3R)Payne, indicating that genes within these regions are targets of selection. Our results provide some support for the hypothesis that inversions persist in natural populations because they hold together favorable combinations of alleles that act together to facilitate adaptive shifts.

Vernalization sensitivity in Arabidopsis thaliana (Brassicaceae): the effects of latitude and FLC variation

Latitudinal variation in climate is predicted to select for latitudinal differentiation in sensitivity to the environmental cues that signal plants to flower at the appropriate time for a given climate. In Arabidopsis thaliana, flowering is promoted by exposure to cold temperatures (vernalization), and several vernalization pathway loci are known. To test whether natural variation in vernalization sensitivity could account for a previously observed latitudinal cline in flowering time in A. thaliana, we exposed 21 European accessions to 0, 10, 20, or 30 d of vernalization and observed leaf number at flowering under short days in a growth chamber. We observed a significant latitudinal cline in vernalization sensitivity: southern accessions were more sensitive to vernalization than northern accessions. In addition, accessions that were late flowering in the absence of vernalization were more sensitive to vernalization cues. Allelic variation at the flowering time regulatory gene FLC was not associated with mean vernalization sensitivity, but one allele class exhibited greater variance in vernalization sensitivity.

Inversion frequencies in Drosophila serrata along an eastern Australian transect

Clinal patterns over broad geographic regions provide a way of identifying characteristics of species under selection and are increasingly being used in quantitative trait locus mapping of adaptive genetic variation in Drosophila. However, interpretations of clinal patterns can be complicated by inversions that also vary clinally and reduce recombination in some parts of the genome. Drosophila serrata (Malloch) is an Australian endemic species being used to investigate the genetic basis of geographic variation in climatic adaptation and mate recognition. Here we describe inversions in D. serrata populations from the east coast of Australia, covering tropical and temperate regions. Seven autosomal paracentric inversions and 1 apparently complex X chromosome arrangement were identified from these populations. All inverted arrangements were relatively more common in tropical populations; 2 common inversions showed clinal patterns over part of the range of D. serrata. Inversion polymorphism was relatively higher in tropical populations and almost absent in populations near the cooler southern border, in agreement with findings on other Drosophila species. While these patterns will complicate mapping of adaptive variation in D. serrata, they suggest that this species will be useful in investigatingthe dynamics of inversion-trait associations in natural populations.

The latitudinal cline in the In(3R)Payne inversion polymorphism has shifted in the last 20 years in Australian Drosophila melanogaster populations

Clinal variation has been described in a number of inversions in Drosophila but these clines are often characterized by cytological techniques using small sample sizes, and associations with specific genes are rarely considered. Here we have developed a molecular assay for In(3R)Payne in Drosophila melanogaster from eastern Australia populations. It shows in repeated samples that the inversion cline is very tightly associated with latitude and is almost fixed in tropical populations while relatively rare in temperate populations. This steep cline has shifted in position in the last 20 years. The heat shock gene, hsr-omega, located centrally inside the inversion sequence, shows a different clinal pattern to In(3R)Payne. These results suggest strong ongoing selection on In(3R)Payne over the last 100 years since the colonization of Australia that is partly independent of hsr-omega.

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.

The genetic covariance among clinal environments after adaptation to an environmental gradient in Drosophila serrata

We examined the genetic basis of clinal adaptation by determining the evolutionary response of life-history traits to laboratory natural selection along a gradient of thermal stress in Drosophila serrata. A gradient of heat stress was created by exposing larvae to a heat stress of 36 degrees for 4 hr for 0, 1, 2, 3, 4, or 5 days of larval development, with the remainder of development taking place at 25 degrees. Replicated lines were exposed to each level of this stress every second generation for 30 generations. At the end of selection, we conducted a complete reciprocal transfer experiment where all populations were raised in all environments, to estimate the realized additive genetic covariance matrix among clinal environments in three life-history traits. Visualization of the genetic covariance functions of the life-history traits revealed that the genetic correlation between environments generally declined as environments became more different and even became negative between the most different environments in some cases. One exception to this general pattern was a life-history trait representing the classic trade-off between development time and body size, which responded to selection in a similar genetic fashion across all environments. Adaptation to clinal environments may involve a number of distinct genetic effects along the length of the cline, the complexity of which may not be fully revealed by focusing primarily on populations at the ends of the cline.

Temperature-related genetic changes in laboratory populations of Drosophila subobscura: evidence against simple climatic-based explanations for latitudinal clines

Parallel latitudinal clines to the long-standing ones in the original Palearctic populations have independently evolved at different rates for chromosomal polymorphism and body size in South and North American populations of Drosophila subobscura since colonization around 25 years ago. This strongly suggests that (micro) evolutionary changes are largely predictable, but the underlying mechanisms are unknown. The putative role of temperature per se was investigated by using three sets of populations at each of three temperatures (13 degrees , 18 degrees , and 22 degrees C) spanning much of the tolerable range for this species. We found a lower chromosomal diversity at the warmest temperature; a quick and consistent shift in gene arrangement frequencies in response to temperature; an evolutionary decrease in wing size, mediated by both cell area and cell number, at 18 degrees C; no relationship between wing size and those inversions involved in latitudinal clines; and a shortening of the basal length of longitudinal vein IV relative to its total length with increasing standard dose. The trends for chromosomal polymorphism and body size were generally inconsistent from simple climatic-based explanations of worldwide latitudinal patterns. The findings are discussed in the light of available information on D. subobscura and results from earlier thermal selection experiments with various Drosophila species.

Genetic analysis of benzoquinone production in Tribolium confusum

Many species of tenebrionid beetles produce and secrete benzoquinones from specialized prothoracic and postabdominal glands. Tribolium confusum produces two compounds methyl-1,4-benzoquinone (MBQ) and ethyl-1,4-benzoquinone (EBQ). These compounds are hypothesized to function as external defense compounds, killing microbes and deterring predators, and their ability to evolve by natural selection depends on both selection and the genetic vs. environmental contribution to phenotypic variation. We crossed a strain of T. confusum that produces high quantities of benzoquinones, b-Pakistan, with a low-producing strain, b-+, and measured both the internal and external quantities of MBQ and EBQ for the two extreme strains and their F1 progeny. Internal amounts show a clear pattern of inheritance, with at least 50% of the phenotypic variation attributed to genotype. Additive and dominance coefficients for internal amounts indicate that the trait is additive with no significant dominance. In contrast, external quantities show little pattern of inheritance. The role of genetics and environment in determining quantities of secretory defensive compounds is important to elucidating the ecology and evolutionary potential of chemical defenses.

Quantitative trait loci affecting knockdown resistance to high temperature in Drosophila melanogaster

Knockdown resistance to high temperature is an ecologically important trait in small insects. A composite interval mapping was performed on the two major autosomes of Drosophila melanogaster to search for quantitative trait loci (QTL) affecting knockdown resistance to high temperature (KRHT). Two dramatically divergent lines from geographically different thermal environments were artificially selected on KRHT. These lines were crossed to produce two backcross (BC) populations. Each BC was analysed for 200 males with 18 marker loci on chromosomes 2 and 3. Three X-linked markers were used to test for X-linked QTL in an exploratory way. The largest estimate of autosome additive effects was found in the pericentromeric region of chromosome 2, accounting for 19.26% (BC to the low line) and 29.15% (BC to the high line) of the phenotypic variance in BC populations, but it could represent multiple closely linked QTL. Complete dominance was apparent for three QTL on chromosome 3, where heat-shock genes are concentrated. Exploratory analysis of chromosome X indicated a substantial contribution of this chromosome to KRHT. The results show that a large-effect QTL with dominant gene action maps on the right arm of chromosome 3. Further, the results confirm that QTL for heat resistance are not limited to chromosome 3.