Searching for the bull's eye: agents and targets of selection vary among geographically disparate cyanogenesis clines in white clover (Trifolium repens L.)

The recurrent evolution of adaptive clines within a species can be used to elucidate the selective factors and genetic responses that underlie adaptation. White clover is polymorphic for cyanogenesis (HCN release with tissue damage), and climate-associated cyanogenesis clines have evolved throughout the native and introduced species range. This polymorphism arises through two independently segregating Mendelian polymorphisms for the presence/absence of two required components: cyanogenic glucosides and their hydrolyzing enzyme linamarase. Cyanogenesis is commonly thought to function in herbivore defense; however, the individual cyanogenic components may also serve other physiological functions. To test whether cyanogenesis clines have evolved in response to the same selective pressures acting on the same genetic targets, we examined cyanogenesis cline shape and its environmental correlates in three world regions: southern New Zealand, the central United States and the US Pacific Northwest. For some regional comparisons, cline shapes are remarkably similar despite large differences in the spatial scales over which clines occur (40–1600 km). However, we also find evidence for major differences in both the agents and targets of selection among the sampled clines. Variation in cyanogenesis frequency is best predicted using a combination of minimum winter temperature and aridity variables. Together, our results provide evidence that recurrent adaptive clines do not necessarily reflect shared adaptive processes.     

A clinally varying promoter polymorphism associated with adaptive variation in wing size in Drosophila

Body size often shows adaptive clines in many ectotherms across altitude and latitude, but little is known about the genetic basis of these adaptive clines. Here we identify a polymorphism in the Dca (Drosophila cold acclimation) gene in Drosophila melanogaster that influences wing size, affects wing:thorax allometry and also controls a substantial proportion of the clinal wing‐size variation. A polymorphism in the promoter region of Dca had two common alleles showing strong reciprocal clinal variation in frequency with latitude along the east coast of Australia. The Dca‐237 allele increased towards the tropics where wing size is smaller. A within‐population association study highlighted that an increase in the frequency of this allele decreased wing size but did not influence thorax size. A manipulated increase in the level of expression of Dca achieved through UAS‐GAL4 was associated with a decrease in wing size but had no effect on thorax size. This was consistent with higher Dca expression levels in family lines with higher frequency of the Dca‐237 allele. Genetic variation in the promoter region of the Dca gene appears to influence adaptive size variation in the eastern Australian cline of Drosophila melanogaster and accounts for more than 10% of the genetic variation in size within and between populations.     

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping‐by‐sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.     

Rapid evolution of an adaptive cyanogenesis cline in introduced North American white clover (Trifolium repens L.)

White clover is polymorphic for cyanogenesis (HCN production after tissue damage), and this herbivore defence polymorphism has served as a classic model for studying adaptive variation. The cyanogenic phenotype requires two interacting biochemical components; the presence/absence of each component is controlled by a simple Mendelian gene (Ac/ac and Li/li). Climate-associated cyanogenesis clines occur in both native (Eurasian) and introduced populations worldwide, with cyanogenic plants predominating in warmer locations. Moreover, previous studies have suggested that epistatic selection may act within populations to maintain cyanogenic (AcLi) plants and acyanogenic plants that lack both components (acli plants) at the expense of plants possessing a single component (Acli and acLi plants). Here, we examine the roles of selection, gene flow and demography in the evolution of a latitudinal cyanogenesis cline in introduced North American populations. Using 1145 plants sampled across a 1650 km transect, we determine the distribution of cyanogenesis variation across the central United States and investigate whether clinal variation is adaptive or an artefact of population introduction history. We also test for the evidence of epistatic selection. We detect a clear latitudinal cline, with cyanogenesis frequencies increasing from 11% to 86% across the transect. Population structure analysis using nine microsatellite loci indicates that the cline is adaptive and not a by-product of demographic history. However, we find no evidence for epistatic selection within populations. Our results provide strong evidence for rapid adaptive evolution in these introduced populations, and they further suggest that the mechanisms maintaining adaptive variation may vary among populations of a species.     

Aridity shapes cyanogenesis cline evolution in white clover (Trifolium repens L.)

Adaptive differentiation between populations is often proposed to be the product of multiple interacting selective pressures, although empirical support for this is scarce. In white clover, populations show adaptive differentiation in frequencies of cyanogenesis, the ability to produce hydrogen cyanide after tissue damage. This polymorphism arises through independently segregating polymorphisms for the presence/absence of two required cyanogenic components, cyanogenic glucosides and their hydrolysing enzyme. White clover populations worldwide have evolved a series of recurrent, climate‐associated clines, with higher frequencies of cyanogenic plants in warmer locations. These clines have traditionally been hypothesized to reflect a fitness trade‐off between chemical defence in herbivore‐rich areas (warmer climates) and energetic costs of producing cyanogenic components in areas of low herbivore pressure (cooler climates). Recent observational studies suggest that cyanogenic components may also be beneficial in water‐stressed environments. We investigated fitness trade‐offs associated with temperature‐induced water stress in the cyanogenesis system using manipulative experiments in growth chambers and population surveys across a longitudinal precipitation gradient in the central United States. We find that plants producing cyanogenic glucosides have higher relative fitness in treatments simulating a moderate, persistent drought stress. In water‐neutral treatments, there are energetic costs to producing cyanogenic components, but only in treatments with nutrient stress. These fitness trade‐offs are consistent with cyanogenesis frequencies in natural populations, where we find clinal variation in the proportion of plants producing cyanogenic glucosides along the precipitation gradient. These results suggest that multiple selective pressures interact to maintain this adaptive polymorphism and that modelling adaptation will require knowledge of environment‐specific fitness effects.     

Drift,selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long‐ and short‐term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre‐date the effects of recent drift, and that functional variation in these loci persists despite small short‐term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of “drift debt,” a nonequilibrium state where present‐day levels of variation overestimate the amount of functional genetic diversity present in future populations.     

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.     

Phenotypic differentiation of the Solidago virgaurea complex along an elevational gradient: Insights from a common garden experiment and population genetics

Plant species distributed along wide elevational or latitudinal gradients show phenotypic variation due to their heterogeneous habitats. This study investigated whether phenotypic variation in populations of the Solidago virgaurea complex along an elevational gradient is caused by genetic differentiation. A common garden experiment was based on seeds collected from nine populations of the S. virgaurea complex growing at elevations from 1,597 m to 2,779 m a.s.l. on Mt. Norikura in central Japan. Population genetic analyses with microsatellite markers were used to infer the genetic structure and levels of gene flow between populations. Leaf mass per area was lower, while leaf nitrogen and chlorophyll concentrations were greater for higher elevations at which seeds were originally collected. For reproductive traits, plants derived from higher elevations had larger flower heads on shorter stems and flowering started earlier. These elevational changes in morphology were consistent with the clines in the field, indicating that phenotypic variation along the elevational gradient would have been caused by genetic differentiation. However, population genetic analysis using 16 microsatellite loci suggested an extremely low level of genetic differentiation of neutral genes among the nine populations. Analysis of molecular variance also indicated that most genetic variation was partitioned into individuals within a population, and the genetic differentiation among the populations was not significant. This study suggests that genome regions responsible for adaptive traits may differ among the populations despite the existence of gene flow and that phenotypic variation of the S. virgaurea complex along the elevational gradient is maintained by strong selection pressure.     

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X‐chromosome, Y‐chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95–4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co‐localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X‐ and Y‐chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.     

Selection and sex‐biased dispersal in a coastal shark: the influence of philopatry on adaptive variation

Sex‐biased dispersal is expected to homogenize nuclear genetic variation relative to variation in genetic material inherited through the philopatric sex. When site fidelity occurs across a heterogeneous environment, local selective regimes may alter this pattern. We assessed spatial patterns of variation in nuclear‐encoded, single nucleotide polymorphisms (SNPs) and sequences of the mitochondrial control region in bonnethead sharks (Sphyrna tiburo), a species thought to exhibit female philopatry, collected from summer habitats used for gestation. Geographic patterns of mtDNA haplotypes and putatively neutral SNPs confirmed female philopatry and male‐mediated gene flow along the northeastern coast of the Gulf of Mexico. A total of 30 outlier SNP loci were identified; alleles at over half of these loci exhibited signatures of latitude‐associated selection. Our results indicate that in species with sex‐biased dispersal, philopatry can facilitate sorting of locally adaptive variation, with the dispersing sex facilitating movement of potentially adaptive variation among locations and environments.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

The molecular genetics of clinal variation: a case study of ebony and thoracic trident pigmentation in Drosophila melanogaster from eastern Australia

Widespread pigmentation diversity coupled with a well‐defined genetic system of melanin synthesis and patterning in Drosophila provides an excellent opportunity to study phenotypes undergoing evolutionary change. Pigmentation variation is highly correlated with different ecological variables and is thought to reflect adaptations to different environments. Several studies have linked candidate genes from Drosophila melanogaster to intra‐population variation and interspecific morphological divergence, but less clearly to variation among populations forming pigmentation clines. We characterized a new thoracic trident pigmentation cline in D. melanogaster populations from eastern Australia, and applied a candidate gene approach to explain the majority of the geographically structured phenotypic variation. More melanized populations from higher latitudes tended to express less ebony than their tropical counterparts, and an independent artificial selection experiment confirmed this association. By partitioning temperature dependent effects, we showed that the genetic differences underlying clinal patterns for trident variation at 25 °C do not explain the patterns observed at 16 °C. Changes in thoracic trident pigmentation could be a common evolutionary response to climatically mediated environmental pressures. On the Australian east coast most of the changes appear to be associated with regulatory divergence of the ebony gene but this depends on temperature.     

Haplotype structure,adaptive history and associations with exploratory behaviour of the DRD4 gene region in four great tit (Parus major) populations

The assessment of genetic architecture and selection history in genes for behavioural traits is fundamental to our understanding of how these traits evolve. The dopamine receptor D4 (DRD4) gene is a prime candidate for explaining genetic variation in novelty seeking behaviour, a commonly assayed personality trait in animals. Previously, we showed that a single nucleotide polymorphism in exon 3 of this gene is associated with exploratory behaviour in at least one of four Western European great tit (Parus major) populations. These heterogeneous association results were explained by potential variable linkage disequilibrium (LD) patterns between this marker and the causal variant or by other genetic or environmental differences among the populations. Different adaptive histories are further hypothesized to have contributed to these population differences. Here, we genotyped 98 polymorphisms of the complete DRD4 gene including the flanking regions for 595 individuals of the four populations. We show that the LD structure, specifically around the original exon 3 SNP is conserved across the four populations and does not explain the heterogeneous association results. Study‐wide significant associations with exploratory behaviour were detected in more than one haplotype block around exon 2, 3 and 4 in two of the four tested populations with different allele effect models. This indicates genetic heterogeneity in the association between multiple DRD4 polymorphisms and exploratory behaviour across populations. The association signals were in or close to regions with signatures of positive selection. We therefore hypothesize that variation in exploratory and other dopamine‐related behaviour evolves locally by occasional adaptive shifts in the frequency of underlying genetic variants.     

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

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

Adaptive introgression in animals: examples and comparison to new mutation and standing variation as sources of adaptive variation

Adaptive genetic variation has been thought to originate primarily from either new mutation or standing variation. Another potential source of adaptive variation is adaptive variants from other (donor) species that are introgressed into the (recipient) species, termed adaptive introgression. Here, the various attributes of these three potential sources of adaptive variation are compared. For example, the rate of adaptive change is generally thought to be faster from standing variation, slower from mutation and potentially intermediate from adaptive introgression. Additionally, the higher initial frequency of adaptive variation from standing variation and lower initial frequency from mutation might result in a higher probability of fixation of the adaptive variants for standing variation. Adaptive variation from introgression might have higher initial frequency than new adaptive mutations but lower than that from standing variation, again making the impact of adaptive introgression variation potentially intermediate. Adaptive introgressive variants might have multiple changes within a gene and affect multiple loci, an advantage also potentially found for adaptive standing variation but not for new adaptive mutants. The processes that might produce a common variant in two taxa, convergence, trans‐species polymorphism from incomplete lineage sorting or from balancing selection and adaptive introgression, are also compared. Finally, potential examples of adaptive introgression in animals, including balancing selection for multiple alleles for major histocompatibility complex (MHC), S and csd genes, pesticide resistance in mice, black colour in wolves and white colour in coyotes, Neanderthal or Denisovan ancestry in humans, mimicry genes in Heliconius butterflies, beak traits in Darwin's finches, yellow skin in chickens and non‐native ancestry in an endangered native salamander, are examined.     

Genomics of adaptive divergence with chromosome‐scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome‐scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low‐crossover centres relative to the high‐crossover peripheries of chromosomes (“Chromosome Centre‐Biased Differentiation”, CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome‐scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome‐scale heterogeneity in crossover rate in evolutionary genomics.     

Parapatric diversification after post-glacial range expansion in the gall fly Urophora cardui (Tephritidae)

Aim Primary and secondary genetic clines in post-glacial colonized regions have different implications for biogeographic distributions and the origin of species. Primary clines arise in situ after colonization as adaptive responses to environmental gradients, while secondary clines are caused by contact between vicariant lineages. Here we analyse primary versus secondary origin of a genetic cline in the tephritid fly Urophora cardui in Jutland, Denmark, in a post-glacial landscape. Location Western Palaearctic. Methods Phylogeographic and demographic analyses of U. cardui based on mitochondrial DNA (mtDNA) genealogies, hierarchical genetic variance tests based on allozymes and distribution analysis of a rare allele from the Jutland cline. Results There was no phylogeographic divergence between the Jutland population of U. cardui north of the cline and neighbouring western European regional populations, which all shared the common western European mtDNA haplotype H1. At nuclear loci, by contrast, the North Jutland population was diverged above the mean level of divergence among regional populations and had no loss of genetic variation. A rare allozyme allele that was frequent in the cline area (up to 45%) and was missing north of the cline also occurred at low frequency (0–14%) elsewhere in the sampling range. Shallow phylogeographic divergence was observed between Russian and western European populations and between English and continental populations. Main conclusions The genetic variation patterns support primary cline evolution and parapatric divergence in Jutland following a demographic expansion of a western European ancestral source population of U. cardui, and suggest cryptic refugia and/or selection in other European population assemblages. The patterns of intra-specific regional divergence are discussed with respect to the interpretation of cryptic refugia in Europe after the most recent ice age.     

Adaptation in temporally variable environments: stickleback armor in periodically breaching bar‐built estuaries

The evolutionary consequences of temporal variation in selection remain hotly debated. We explored these consequences by studying threespine stickleback in a set of bar‐built estuaries along the central California coast. In most years, heavy rains induce water flow strong enough to break through isolating sand bars, connecting streams to the ocean. New sand bars typically re‐form within a few weeks or months, thereby re‐isolating populations within the estuaries. These breaching events cause severe and often extremely rapid changes in abiotic and biotic conditions, including shifts in predator abundance. We investigated whether this strong temporal environmental variation can maintain within‐population variation while eroding adaptive divergence among populations that would be caused by spatial variation in selection. We used neutral genetic markers to explore population structure and then analysed how stickleback armor traits, the associated genes Eda and Pitx1 and elemental composition (%P) varies within and among populations. Despite strong gene flow, we detected evidence for divergence in stickleback defensive traits and Eda genotypes associated with predation regime. However, this among‐population variation was lower than that observed among other stickleback populations exposed to divergent predator regimes. In addition, within‐population variation was very high as compared to populations from environmentally stable locations. Elemental composition was strongly associated with armor traits, Eda genotype and the presence of predators, thus suggesting that spatiotemporal variation in armor traits generates corresponding variation in elemental phenotypes. We conclude that gene flow, and especially temporal environmental variation, can maintain high levels of within‐population variation while reducing, but not eliminating, among‐population variation driven by spatial environmental variation.     

Latitudinal clines in<Emphasis Type="Italic">Drosophila melanogaster</Emphasis>: Body size,allozyme frequencies,inversion frequencies,and the insulin-signalling pathway

Many latitudinal clines exist inDrosophila melanogaster: in adult body size, in allele frequency at allozyme loci, and in frequencies of common cosmopolitan inversions. The question is raised whether these latitudinal clines are causally related. This review aims to connect data from two very different fields of study, evolutionary biology and cell biology, in explaining such natural genetic variation inD. melanogaster body size and development time. It is argued that adult body size clines, inversion frequency clines, and clines in allele frequency at loci involved in glycolysis and glycogen storage are part of the same adaptive strategy. Selection pressure is expected to differ at opposite ends of the clines. At high latitudes, selection onD. melanogaster would favour high larval growth rate at low temperatures, and resource storage in adults to survive winter. At low latitudes selection would favour lower larval critical size to survive crowding, and increased male activity leading to high male reproductive success. Studies of the insulin-signalling pathway inD. melanogaster point to the involvement of this pathway in metabolism and adult body size. The genes involved in the insulin-signalling pathway are associated with common cosmopolitan inversions that show latitudinal clines. Each chromosome region connected with a large common cosmopolitan inversion possesses a gene of the insulin transmembrane complex, a gene of the intermediate pathway and a gene of the TOR branch. The hypothesis is presented that temperateD. melanogaster populations have a higher frequency of a ’thrifty’ genotype corresponding to high insulin level or high signal level, while tropical populations possess a more ’spendthrift’ genotype corresponding to low insulin or low signal level.     

Genetic variation underlying local adaptation of diapause induction along a cline in a butterfly

Diapause is a life history strategy allowing individuals to arrest development until favourable conditions return, and it is commonly induced by shortened day length that is latitude specific for local populations. Although understanding the evolutionary dynamics of a threshold trait like diapause induction provides insights into the adaptive process and adaptive potential of populations, the genetic mechanism of variation in photoperiodic induction of diapause is not well understood. Here, we investigate genetic variation underlying latitudinal variation in diapause induction and the selection dynamics acting upon it. Using a genomewide scan for divergent regions between two populations of the butterfly Pararge aegeria that differ strongly in their induction thresholds, we identified and investigated the patterns of variation in those regions. We then tested the association of these regions with diapause induction using between‐population crosses, finding significant SNP associations in four genes present in two chromosomal regions, one with the gene period, and the other with the genes kinesin, carnitine O‐acetyltransferase and timeless. Patterns of allele frequencies in these two regions in population samples along a latitudinal cline suggest strong selection against heterozygotes at two genes within these loci (period, timeless). Evidence for additional loci modifying the diapause decision was found in patterns of allelic change in relation to induction thresholds over the cline, as well as in backcross analyses. Taken together, population‐specific adaptations of diapause induction appear to be due to a combination of alleles of larger and smaller effect size, consistent with an exponential distribution of effect sizes involved in local adaption.