On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre‐eminent system for the study of selective pressures that arise from host–pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population‐genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.     

Reduced MHC variation in a threatened tuatara species

The relationship between neutral and adaptive genetic diversity is important to understand in assessing the implications of a population bottleneck. Fitness-related genes, such as those of the major histocompatibility complex (MHC), may be influenced by selection, and so retain diversity even when it is lost at neutral markers. We measured MHC class I variation in an archaic reptile species Sphenodon guntheri [North Brother Island (NBI) tuatara], which naturally occurs on one 4 ha island in Cook Strait, New Zealand, and has low levels of microsatellite diversity. MHC variation in S. guntheri was compared with microsatellite DNA variation, and with MHC variation in a large population of Sphenodon punctatus (Cook Strait tuatara) on Stephens Island. The NBI population shows significantly decreased levels of genetic diversity compared with the Stephens Island population. Only three different MHC sequences and three genotypes were found on NBI, compared with 15 sequences and 21 genotypes in a similar sample size from Stephens Island. Two sequences appear to be unique to the NBI population. The assortment of sequence variants into genotypes suggests strong gametic disequilibrium between two MHC class I loci in S. guntheri , and only two haplotypes that were present in Hardy–Weinberg proportions were identified. MHC diversity in NBI tuatara appears to be largely influenced by genetic drift, consistent with a recent population bottleneck. This may compromise the ability of this population to respond to novel disease threats.     

Drift and selection influence geographic variation at immune loci of prairie-chickens

Previous studies of immunity in wild populations have focused primarily on genes of the major histocompatibility complex (MHC); however, studies of model species have identified additional immune-related genes that also affect fitness. In this study, we sequenced five non-MHC immune genes in six greater prairie-chicken (Tympanuchus cupido) populations that have experienced varying degrees of genetic drift as a consequence of population bottlenecks and fragmentation. We compared patterns of geographic variation at the immune genes with six neutral microsatellite markers to investigate the relative effects of selection and genetic drift. Global F(ST) outlier tests identified positive selection on just one of five immune genes (IAP-1) in one population. In contrast, at other immune genes, standardized G'(ST) values were lower than those at microsatellites for a majority of pairwise population comparisons, consistent with balancing selection or with species-wide positive or purifying selection resulting in similar haplotype frequencies across populations. The effects of genetic drift were also evident as summary statistics (e.g., Tajima's D) did not differ from neutrality for the majority of cases, and immune gene diversity (number of haplotypes per gene) was correlated positively with population size. In summary, we found that both genetic drift and selection shaped variation at the five immune genes, and the strength and type of selection varied among genes. Our results caution that neutral forces, such as drift, can make it difficult to detect current selection on genes.     

A comparison of variability and population structure for major histocompatibility complex and microsatellite loci in California coastal steelhead (Oncorhynchus mykiss Walbaum)

The major histocompatibility complex (MHC) contains genes integral to immune response in vertebrates. MHC genes have been shown to be under selection in a number of vertebrate taxa, making them intriguing for population genetic studies. We have conducted a survey of genetic variation in an MHC class II gene for steelhead trout from 24 sites in coastal California and compared this variation to that observed at 16 presumably neutral microsatellite loci. A high amount of allelic variation was observed at the MHC when compared to previously published studies on other Pacific salmonids. Elevated nonsynonymous substitutions, relative to synonymous substitutions, were detected at the MHC gene, indicating the signature of historical balancing selection. The MHC data were tested for correlations to and deviations from the patterns found with the microsatellite data. Estimates of allelic richness for the MHC gene and for the microsatellites were positively correlated, as were estimates of population differentiation (F(ST)). An analysis for F(ST) outliers indicates that the MHC locus has an elevated F(ST) relative to the neutral expectation, although a significant result was found for only one particular geographical subgroup. Relatively uniform allele frequency distributions were detected in four populations, although this finding may be partially due to recent population bottlenecks. These results indicate that, at the scale studied here, drift and migration play a major role in the observed geographical variability of MHC genes in steelhead, and that contemporary selection is relatively weak and difficult to detect.     

Adaptive and neutral genetic differentiation among Scottish and endangered Irish red grouse (<Emphasis Type="Italic">Lagopus lagopus scotica</Emphasis>)

Studying patterns of intra-specific genetic variation among populations allows for a better understanding of population structure and local adaptation. However, those patterns may differ according to the genetic markers applied, as neutral genetic markers reflect demographic processes and random genetic drift, whereas adaptive markers also carry the footprint of selection. In combination, neutral and adaptive genetic markers permit to assess the relative roles of drift and selection in shaping population structure. Among the best understood adaptive genetic loci are the genes of the major histocompatibility complex (MHC). We here study variation and differentiation at neutral SNP markers and MHC class II genes in red grouse (Lagopus lagopus scotica) from Ireland and Scotland. Irish red grouse populations are fragmented and drastically declining, but red grouse are abundant in Scotland. We find evidence for positive selection acting on the MHC genes and variation in MHC gene copy numbers among Irish individuals. Furthermore, there was significant population differentiation among red grouse from Ireland and Scotland at the neutral SNP markers (F_ST = 0.084) and the MHC-BLB genes (F_ST: BLB1 = 0.116, BLB2 = 0.090, BLB3 = 0.104). Differentiation at the MHC-BLB1 was significantly higher than at the neutral SNP markers, suggesting that selection plays an important role in shaping MHC variation, in addition to genetic drift. We speculate that the observed differentiation pattern might be due to local adaptation to different parasite regimes. These findings have strong conservation implications and we advise against the introduction of Scottish red grouse to supplement Irish populations.     

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake,the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long‐term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.     

Can balancing selection on MHC loci counteract genetic drift in small fragmented populations of black grouse?

The ability of natural populations to adapt to new environmental conditions is crucial for their survival and partly determined by the standing genetic variation in each population. Populations with higher genetic diversity are more likely to contain individuals that are better adapted to new circumstances than populations with lower genetic diversity. Here, we use both neutral and major histocompatibility complex (MHC) markers to test whether small and highly fragmented populations hold lower genetic diversity than large ones. We use black grouse as it is distributed across Europe and found in populations with varying degrees of isolation and size. We sampled 11 different populations; five continuous, three isolated, and three small and isolated. We tested patterns of genetic variation in these populations using three different types of genetic markers: nine microsatellites and 21 single nucleotide polymorphisms (SNPs) which both were found to be neutral, and two functional MHC genes that are presumably under selection. The small isolated populations displayed significantly lower neutral genetic diversity compared to continuous populations. A similar trend, but not as pronounced, was found for genotypes at MHC class II loci. Populations were less divergent at MHC genes compared to neutral markers. Measures of genetic diversity and population genetic structure were positively correlated among microsatellites and SNPs, but none of them were correlated to MHC when comparing all populations. Our results suggest that balancing selection at MHC loci does not counteract the power of genetic drift when populations get small and fragmented.     

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus)

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA‐DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N_e < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.     

TEMPORAL VARIATION AT THE MHC CLASS IIB IN WILD POPULATIONS OF THE GUPPY (POECILIA RETICULATA)

Understanding genetic diversity in natural populations is a fundamental objective of evolutionary biology. The immune genes of the major histocompatibility complex (MHC) are excellent candidates to study such diversity because they are highly polymorphic in populations. Although balancing selection may be responsible for maintaining diversity at these functionally important loci, temporal variation in selection pressure has rarely been examined. We examine temporal variation in MHC class IIB diversity in nine guppy (Poecilia reticulata) populations over two years. We found that five of the populations changed significantly more at the MHC than at neutral (microsatellite) loci as measured by F_ST, which suggests that the change at the MHC was due to selection and not neutral processes. Additionally, pairwise population differentiation measures at the MHC were higher in 2007 than in 2006, with the signature of selection changing from homogenizing to diversifying selection or neutral evolution. Interestingly, within the populations the magnitude of the change at the MHC between years was related to the change in the proportion of individuals infected by a common parasite, indicating a link between genetic structure and the parasite. Our data thereby implicate temporal variation in selective pressure as an important mechanism maintaining diversity at the MHC in wild populations.     

Balancing selection on MHC class I in wild brown trout Salmo trutta

Evidence is reported for balancing selection acting on variation at major histocompatibility complex (MHC) in wild populations of brown trout Salmo trutta. First, variation at an MHC class I (satr-uba)-linked microsatellite locus (mhc1) is retained in small S. trutta populations isolated above waterfalls although variation is lost at neutral microsatellite markers. Second, populations across several catchments are less differentiated at mhc1 than at neutral markers, as predicted by theory. The population structure of these fish was also elucidated.     

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.     

Loss of MHC and neutral variation in Peary caribou: genetic drift is not mitigated by balancing selection or exacerbated by MHC allele distributions

Theory and empirical results suggest that the rate of loss of variation at Mhc and neutral microsatellite loci may differ because selection influences Mhc genes, and because a high proportion of rare alleles at Mhc loci may result in high rates of loss via drift. Most published studies compare Mhc and microsatellite variation in various contemporary populations to infer the effects of population size on genetic variation, even though different populations are likely to have different demographic histories that may also affect contemporary genetic variation. We directly compared loss of variation at Mhc and microsatellite loci in Peary caribou by comparing historical and contemporary samples. We observed that similar proportions of genetic variation were lost over time at each type of marker despite strong evidence for selection at Mhc genes. These results suggest that microsatellites can be used to estimate genome-wide levels of variation, but also that adaptive potential is likely to be lost following population bottlenecks. However, gene conversion and recombination at Mhc loci may act to increase variation following bottlenecks.     

Adaptive and Neutral Genetic Variation and Colonization History of Atlantic Salmon, Salmo salar

Synopsis I combined neutral microsatellite markers with the major histocompatibility complex (MHC) class IIB to study genetic differentiation and colonization history in Atlantic salmon, Salmo salar, in the Baltic Sea and in the north-eastern Atlantic. Baltic salmon populations have lower levels of microsatellite genetic variation, in terms of heterozygosity and allelic richness than Atlantic populations, confirming earlier findings with other genetic markers, suggesting that the Baltic Sea populations have been exposed to genetic bottlenecks, most likely at a founding event. On the other hand, the level of MHC variation was similar in the Baltic and in the north-eastern Atlantic, indicating that positive balancing selection has increased the level of MHC-variation. Both microsatellite and MHC class IIB genetic variation give strong support to the hypothesis that the Baltic salmon are of a biphyletic origin, the southern population in this study is strongly differentiated from both the northern Baltic salmon populations and from the north-eastern Atlantic populations. Salmon may have colonized the northern Baltic Sea either from the south, via the so called “N?rke strait” or from the north, via a proposed historical connection between the White Sea and the northern Baltic. At microsatellites, no significant isolation-by distance was found at either colonization route. At the MHC, populations were significantly isolated by distance when assuming that colonization occurred via the “N?rke strait”.     

High allelic variation of MHC class II alpha antigen and the role of selection in wild and cultured Sparus aurata populations

Genetic polymorphism and differentiation in wild and cultured sea bream samples were studied after amplification, cloning, and partial sequence of the major histocompatibility complex (MHC) class II alpha antigen. Forty-one alleles were detected from 43 unrelated individuals and sequence alignment of the obtained alleles revealed 28 polymorphic sites. High heterozygosity values and allelic richness were unveiled for both wild and cultured populations. The substitution pattern (d_N /d_S = 0.7) is not consistent with the effect of diversifying selection, indicating lower selection pressure on the a2 domain, as well as that too few advantageous non-synonymous mutations have accumulated as substrate for the diversifying selection to act. Comparison with previously published results on microsatellite markers suggests that balancing selection in MHC genes reduces the genetic drift and bottleneck effects that are common in aquaculture and which are known to reduce genetic variation at neutral markers. The present study stresses that both coding and non-coding loci should be analyzed for designing proper management strategies.     

Variation in MHC genotypes in two populations of house sparrow (Passer domesticus) with different population histories

Small populations are likely to have a low genetic ability for disease resistance due to loss of genetic variation through inbreeding and genetic drift. In vertebrates, the highest genetic diversity of the immune system is located at genes within the major histocompatibility complex (MHC). Interestingly, parasite‐mediated selection is thought to potentially maintain variation at MHC loci even in populations that are monomorphic at other loci. Therefore, general loss of genetic variation in the genome may not necessarily be associated with low variation at MHC loci. We evaluated inter‐ and intrapopulation variation in MHC genotypes between an inbred (Aldra) and a relatively outbred population (Hestmannøy) of house sparrows (Passer domesticus) in a metapopulation at Helgeland, Norway. Genomic (gDNA) and transcribed (cDNA) alleles of functional MHC class I and IIB loci, along with neutral noncoding microsatellite markers, were analyzed to obtain relevant estimates of genetic variation. We found lower allelic richness in microsatellites in the inbred population, but high genetic variation in MHC class I and IIB loci in both populations. This suggests that also the inbred population could be under balancing selection to maintain genetic variation for pathogen resistance.     

The mode of evolution of molecular markers in populations of house mice under artificial selection for locomotor behavior

A complete understanding of the mode of evolution of molecular markers is important for making inferences about different population genetic parameters, especially because a number of studies have reported patterns of allelic variation at molecular markers that are not in agreement with neutral evolutionary expectations. In the present study, house mice (Mus domesticus) from the fourteenth generation of a selection experiment for increased voluntary wheel-running activity were used to test how selection on a complex behavior affects the distribution of allelic variation by examining patterns of variation at six microsatellite and four allozyme loci. This population had a hierarchical structure that allowed for simultaneous testing of the effects of selection and genetic drift on the distribution of allelic variation by comparing observed patterns of allele frequencies and estimates of genetic divergence at multiple hierarchical levels to expectations under models of neutral evolution. The levels of genetic divergence among replicate lines and between selection groups, estimated from microsatellite data or pooled microsatellite and allozyme data, were not significantly different from expectations under neutral evolution. Furthermore, the pattern of change of allele frequencies between the base population and generation 14 was largely in agreement with expectations under neutral evolution (although the PGM locus exhibited a pattern of change within populations that was difficult to explain under neutral evolution). Overall the results generally provide support for the neutral evolution of molecular markers.     

Balancing selection, random genetic drift, and genetic variation at the major histocompatibility complex in two wild populations of guppies (Poecilia reticulata)

Our understanding of the evolution of genes of the major histocompatibility complex (MHC) is rapidly increasing, but there are still enigmatic questions remaining, particularly regarding the maintenance of high levels of MHC polymorphisms in small, isolated populations. Here, we analyze the genetic variation at eight microsatellite loci and sequence variation at exon 2 of the MHC class IIB (DAB) genes in two wild populations of the Trinidadian guppy, Poecilia reticulata. We compare the genetic variation of a small (Ne, 100) and relatively isolated upland population to that of its much larger (Ne approximately 2400) downstream counterpart. As predicted, microsatellite diversity in the upland population is significantly lower and highly differentiated from the population further downstream. Surprisingly, however, these guppy populations are not differentiated by MHC genetic variation and show very similar levels of allelic richness. Computer simulations indicate that the observed level of genetic variation can be maintained with overdominant selection acting at three DAB loci. The selection coefficients differ dramatically between the upland (s > or = 0.2) and lowland (s < or = 0.01) populations. Parasitological analysis on wild-caught fish shows that parasite load is significantly higher on upland than on lowland fish, which suggests that large differences in selection intensity may indeed exist between populations. Based on the infection intensity, a substantial proportion of the upland fish would have suffered direct or indirect fitness consequences as a result of their high parasite loads. Selection by parasites plays a particularly important role in the evolution of guppies in the upland habitat, which has resulted in high levels of MHC diversity being maintained in this population despite considerable genetic drift.     

Extensive polymorphism and geographical variation at a positively selected MHC class II B gene of the lesser kestrel (Falco naumanni)

Understanding the selective forces that shape genetic variation in natural populations remains a high priority in evolutionary biology. Genes at the major histocompatibility complex (MHC) have become excellent models for the investigation of adaptive variation and natural selection because of their crucial role in fighting off pathogens. Here we present one of the first data sets examining patterns of MHC variation in wild populations of a bird of prey, the lesser kestrel, Falco naumanni . We report extensive polymorphism at the second exon of a putatively functional MHC class II gene, Fana- DAB*1. Overall, 103 alleles were isolated from 121 individuals sampled from Spain to Kazakhstan. Bayesian inference of diversifying selection suggests that several amino acid sites may have experienced strong positive selection (ω = 4.02 per codon). The analysis also suggests a prominent role of recombination in generating and maintaining MHC diversity (ρ = 4 Nc  = 0.389 per codon, θ = 0.017 per codon). Both the Fana -DAB*1 locus and a set of eight polymorphic microsatellite markers revealed an isolation-by-distance pattern across the Western Palaearctic ( r  = 0.67; P  = 0.01 and r  = 0.50; P  = 0.04, respectively). Nonetheless, geographical variation at the MHC contrasts with relatively uniform distributions in the frequencies of microsatellite alleles. In addition, we found lower fixation rates in the MHC than those predicted by genetic drift after controlling for neutral mitochondrial sequences. Our results therefore underscore the role of balancing selection as well as spatial variations in parasite-mediated selection regimes in shaping MHC diversity when gene flow is limited.     

Latitudinal divergence of common frog (Rana temporaria) life history traits by natural selection: evidence from a comparison of molecular and quantitative genetic data

The relative roles of natural selection and direct environmental induction, as well as of natural selection and genetic drift, in creating clinal latitudinal variation in quantitative traits have seldom been assessed in vertebrates. To address these issues, we compared molecular and quantitative genetic differentiation between six common frog (Rana temporaria) populations along an approximately 1600 km long latitudinal gradient across Scandinavia. The degree of population differentiation (QST approximately 0.81) in three heritable quantitative traits (age and size at metamorphosis, growth rate) exceeded that in eight (neutral) microsatellite loci (FST = 0.24). Isolation by distance was clear for both neutral markers and quantitative traits, but considerably stronger for one of the three quantitative traits than for neutral markers. QST estimates obtained using animals subjected to different rearing conditions (temperature and food treatments) revealed some environmental dependency in patterns of population divergence in quantitative traits, but in general, these effects were weak in comparison to overall patterns. Pairwise comparisons of FST and QST estimates across populations and treatments revealed that the degree of quantitative trait differentiation was not generally predictable from knowledge of that in molecular markers. In fact, both positive and negative correlations were observed depending on conditions where the quantitative genetic variability had been measured. All in all, the results suggest a very high degree of genetic subdivision both in neutral marker genes and genes coding quantitative traits across a relatively recently (< 9000 years) colonized environmental gradient. In particular, they give evidence for natural selection being the primary agent behind the observed latitudinal differentiation in quantitative traits.