Footprints of selection in wild populations of Bicyclus anynana along a latitudinal cline

One of the major questions in ecology and evolutionary biology is how variation in the genome enables species to adapt to divergent environments. Here, we study footprints of thermal selection in candidate genes in six wild populations of the afrotropical butterfly Bicyclus anynana sampled along a c. 3000 km latitudinal cline. We sequenced coding regions of 31 selected genes with known functions in metabolism, pigment production, development and heat shock responses. These include genes for which we expect a priori a role in thermal adaptation and, thus, varying selection pressures along a latitudinal cline, and genes we do not expect to vary clinally and can be used as controls. We identified amino acid substitution polymorphisms in 13 genes and tested these for clinal variation by correlation analysis of allele frequencies with latitude. In addition, we used two F_ST‐based outlier methods to identify loci with higher population differentiation than expected under neutral evolution, while accounting for potentially confounding effects of population structure and demographic history. Two metabolic enzymes of the glycolytic pathway, UGP and Treh, showed clinal variation. The same loci showed elevated population differentiation and were identified as significant outliers. We found no evidence of clines in the pigmentation genes, heat shock proteins and developmental genes. However, we identified outlier loci in more localized parts of the range in the pigmentation genes yellow and black. We discuss that the observed clinal variation and elevated population divergence in UGP and Treh may reflect adaptation to a geographic thermal gradient.     

Contrasting molecular population genetics of four hexokinases in Drosophila melanogaster, D. simulans and D. yakuba

As part of a larger study contrasting patterns of variation in regulatory and nonregulatory enzymes of the central metabolic pathways we have examined the molecular variation in four uncharacterized hexokinase genes unique to muscle, fat body, and testis in Drosophila melanogaster, D. simulans, and D. yakuba. Earlier isoenzyme studies had designated these genes as Hex-A, Hex-C, and Hex-t. There are two tightly linked testes-specific genes designated here as Hex-t1 and Hex-t2. Substantial and concordant differences across species are seen in levels of both amino acid and silent polymorphism. The flight muscle form Hex-A is the most conserved followed by the fat body hexokinase Hex-C and testis-specific hexokinases Hex-t1 and Hex-t2. While constraints acting at the amino acid level are expected, the silent polymorphisms follow this pattern as well. All genes are in regions of normal recombination, therefore hitchhiking and background selection are not likely causes of interlocus differences. In D. melanogaster latitudinal clines are seen for amino acid polymorphisms at the Hex-C and Hex-t2 loci. There is evidence for accelerated amino acid substitution in Hex-t1 that has lost residues known to be associated with glucose and glucose-6-phosphate binding. D. simulans shows substantial linkage phase structuring that suggests historical population subdivision.     

The functional impact of Pgm amino acid polymorphism on glycogen content in Drosophila melanogaster

Earlier studies of the common PGM allozymes in Drosophila melanogaster reported no in vitro activity differences. However, our study of nucleotide variation observed that PGM allozymes are a heterogeneous mixture of amino acid polymorphisms. In this study, we analyze 10 PGM protein haplotypes with respect to PGM activity, thermostability, and adult glycogen content. We find a twofold difference in activity among PGM protein haplotypes that is associated with a threefold difference in glycogen content. The latitudinal clines for several Pgm amino acid polymorphisms show that high PGM activity, and apparently higher flux to glycogen synthesis, parallel the low activity clines at G6PD for reduced pentose shunt flux in northern latitudes. This suggests that amino acid polymorphism is under selection at this branch point and may be favored for increased metabolic storage associated with stress resistance and adaptation to temperate regions.     

Latitudinal clines: an evolutionary view on biological rhythms

Properties of the circadian and annual timing systems are expected to vary systematically with latitude on the basis of different annual light and temperature patterns at higher latitudes, creating specific selection pressures. We review literature with respect to latitudinal clines in circadian phenotypes as well as in polymorphisms of circadian clock genes and their possible association with annual timing. The use of latitudinal (and altitudinal) clines in identifying selective forces acting on biological rhythms is discussed, and we evaluate how these studies can reveal novel molecular and physiological components of these rhythms.     

Clinal variation for amino acid polymorphisms at the Pgm locus in Drosophila melanogaster

Clinal variation is common for enzymes in the glycolytic pathway for Drosophila melanogaster and is generally accepted as an adaptive response to different climates. Although the enzyme phosphoglucomutase (PGM) possesses several allozyme polymorphisms, it is unique in that it had been reported to show no clinal variation. Our recent DNA sequence investigation of Pgm found extensive cryptic amino acid polymorphism segregating with the allozyme alleles. In this study, we characterize the geographic variation of Pgm amino acid polymorphisms at the nucleotide level along a latitudinal cline in the eastern United States. A survey of 15 SNPs across the Pgm gene finds significant clinal differentiation for the allozyme polymorphisms as well as for many of the cryptic amino acid polymorphisms. A test of independence shows that pervasive linkage disequilibrium across this gene region can explain many of the amino acid clines. A single Pgm haplotype defined by two amino acid polymorphisms shows the strongest correlation with latitude and the steepest change in allele frequency across the cline. We propose that clinal selection at Pgm may in part explain the extensive amino acid polymorphism at this locus and is consistent with a multilocus response to selection in the glycolytic pathway.     

Extensive amino acid polymorphism at the pgm locus is consistent with adaptive protein evolution in Drosophila melanogaster

PGM plays a central role in the glycolytic pathway at the branch point leading to glycogen metabolism and is highly polymorphic in allozyme studies of many species. We have characterized the nucleotide diversity across the Pgm gene in Drosophila melanogaster and D. simulans to investigate the role that protein polymorphism plays at this crucial metabolic branch point shared with several other enzymes. Although D. melanogaster and D. simulans share common allozyme mobility alleles, we find these allozymes are the result of many different amino acid changes at the nucleotide level. In addition, specific allozyme classes within species contain several amino acid changes, which may explain the absence of latitudinal clines for PGM allozyme alleles, the lack of association of PGM allozymes with the cosmopolitan In(3L)P inversion, and the failure to detect differences between PGM allozymes in functional studies. We find a significant excess of amino acid polymorphisms within D. melanogaster when compared to the complete absence of fixed replacements with D. simulans. There is also strong linkage disequilibrium across the 2354 bp of the Pgm locus, which may be explained by a specific amino acid haplotype that is high in frequency yet contains an excess of singleton polymorphisms. Like G6pd, Pgm shows strong evidence for a branch point enzyme that exhibits adaptive protein evolution.     

Opposing clines for high and low temperature resistance in Drosophila melanogaster

In insects, species comparisons suggest a weak association between upper thermal limits and latitude in contrast to a stronger association for lower limits. To compare this to latitudinal patterns of thermal responses within species, we considered latitudinal variation in heat and cold resistance in Drosophila melanogaster. We found opposing clines in resistance to these temperature extremes in comparisons of 17–24 populations from coastal eastern Australia. Knockdown time following heat shock increased towards the tropics, whereas recovery time following cold shock decreased towards temperate latitudes. Mortality following cold shock also showed a clinal pattern. Clinal associations with latitude were linear and related to minimum temperatures in the coldest month (for cold resistance) and maximum temperatures in the warmest month (for heat resistance). This suggests that within species both high and low temperature responses can vary with latitude as a consequence of direct or indirect effects of selection.     

Variation in Drosophila melanogaster central metabolic genes appears driven by natural selection both within and between populations

In this report, we examine the hypothesis that the drivers of latitudinal selection observed in the eastern US Drosophila melanogaster populations are reiterated within seasons in a temperate orchard population in Pennsylvania, USA. Specifically, we ask whether alleles that are apparently favoured in northern populations are also favoured early in the spring, and decrease in frequency from the spring to autumn with the population expansion. We use SNP data collected for 46 metabolic genes and 128 SNPs representing the central metabolic pathway and examine for the aggregate SNP allele frequencies whether the association of allele change with latitude and that with increasing days of spring–autumn season are reversed. Testing by random permutation, we observe a highly significant negative correlation between these associations that is consistent with this expectation. This correlation is stronger when we confine our analysis to only those alleles that show significant latitudinal changes. This pattern is not caused by association with chromosomal inversions. When data are resampled using SNPs for amino acid change the relationship is not significant but is supported when SNPs associated with cis-expression are only considered. Our results suggest that climate factors driving latitudinal molecular variation in a metabolic pathway are related to those operating on a seasonal level within populations.     

Latitudinal clines for nucleotide polymorphisms in the Esterase 6 gene of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Previous studies have found non-neutral patterns of nucleotide polymorphism in the promoter and coding regions of Est6 in D. melanogaster. Coding region polymorphism peaks around two closely linked replacement differences associated with the EST6-F/EST6-S allozyme polymorphism. The promoter contains two common, highly diverged haplotype groups, P1 and P7, that differentially affect Est6 expression. Allozyme studies have also revealed latitudinal clines in EST6-F and EST6-S frequencies that recur across continents. Here we analyse nucleotide polymorphisms across the promoter and the region of peak coding sequence polymorphism in 10 Australian populations along a 25° latitudinal gradient in order to examine the basis for the allozyme clines. As with the earlier studies, we find an excess of intermediate to high frequency variants in both the P1/P7 region and around the two EST6-F/EST6-S replacements in some populations. The two EST6-F/EST6-S replacement polymorphisms show latitudinal clines whereas the P1 and P7 groups of promoter haplotypes do not. However the strongest clines are for three co-segregating silent site polymorphisms in a 4 bp stretch at the 3′ end of the sequenced region. Monte Carlo simulations show that the clines for those three sites can explain all others in the data but none of the others can explain those three. Thus the allozyme clines may not reflect selection on either the P1/P7 polymorphism or the two replacements previously associated with the EST6-F/EST-S difference.     

Latitudinal clines in the timing and temperature-sensitivity of photoperiodic reproductive diapause in Drosophila montana

Reproductive diapause is a primary mechanism used by arthropods to synchronize their life cycle with seasonal changes in temperate regions. Our study species, Drosophila montana, represents the northern insect species where flies enter reproductive diapause under short day conditions and where the precise timing of diapause is crucial for both survival and offspring production. We have studied clinal variation in the critical day length for female diapause induction (CDL) and their overall susceptibility to enter diapause (diapause incidence), as well as the temperature sensitivity of these traits. The study was performed using multiple strains from four latitudinal clines of the species – short clines in Finland and Alaska and long clines in the Rocky Mountains and the western coast of North America – and from one population in Kamchatka, Russia. CDL showed strong latitudinal clines on both continents, decreasing by one hour per five degrees decline in latitude, on average. CDL also decreased in all populations along with an increase in fly rearing temperature postponing the diapause to later calendar time, the effects of temperature being stronger in southern than in northern population. Female diapause incidence was close to 100% under short day/low temperature conditions in all populations, but decreased below 50% even under short days in 19°C in the southern North American western coast populations and in 22°C in most populations. Comparing a diversity of climatic data for the studied populations showed that while CDL is under a tight photoperiodic regulation linked with latitude, its length depends also on climatic factors determining the growing season length. Overall, the study deepens our understanding of how spatial and environmental parameters affect the seasonal timing of an important biological event, reproductive diapause and helps to estimate the evolutionary potential of insect populations to survive in changing climatic conditions.     

Disentangling the roles of history and local selection in shaping clinal variation of allele frequencies and gene expression in Norway spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (F(ST) = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to F(ST). The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.     

Are latitudinal clines in body size adaptive?

Body size of animals often increases with increasing latitude. These latitudinal clines in body size have interested biologists for over 150 years. However, the mechanisms that generate these clines in size are still unclear, though latitudinal gradients in temperature appear to play an important role. More importantly, many studies that examine latitudinal clines in body size and the mechanisms responsible for these clines use phenotypic data, confounding genetic (adaptive) and non‐genetic (plasticity) sources of variation. Yet, most of these studies make adaptive conclusions based on phenotypic measures of size. Here I show the dangers of making adaptive inferences from phenotypic measures of size. In addition, I use a specific form of plasticity in body size of ectotherms, called the temperature–size rule, to illustrate how confusion about genetic and non‐genetic contributions to phenotypic variation has hampered progress in understanding the evolution of latitudinal clines in size. Field‐based measurements of body size can no doubt be influenced by plasticity, but demonstrating that latitudinal clines have a genetic basis is necessary to show that these patterns are adaptive.     

Latitudinal pattern in body size in a cockroach,Eupolyphaga sinensis

Body size of insects with flexible life cycles is expected to conform to the saw‐tooth model, a model in which the relationship between size and developmental time depends on length of the growing season. In species with high variability in terms of voltinism, however, more complex patterns can be expected. Few empirical studies have demonstrated the existence of such relationships, or whether climatic factors contribute to these relationships. In this study, we investigated the geographic variation in body size of the Chinese cockroach, Eupolyphaga sinensis Walker (Blattaria: Polyphagidae), which has a variable life cycle length. The sizes of adults – collected from 14 localities ranging from temperate to subtropical zones in China – were measured, using body length, body width, and pronotum width as parameters. The relationship between size, latitude, and climate factors (encompassing 10 variables) was then investigated. We found that the body size of E. sinensis varied considerably with latitude: cockroaches were larger at low and high latitudes, but smaller at intermediate latitudes. Thus, the relationship between climate and body size conformed to a saw‐tooth pattern. Results indicate that two factors were significantly associated with body size clines: season length and variability in life cycle length. Our results also demonstrated that climatic factors contribute to latitudinal clines in body size, which has important ecological and evolutionary implications. It can be expected that global climate change may alter latitudinal clines in body size of E. sinensis.     

Contrasting plant physiological adaptation to climate in the native and introduced range of Hypericum perforatum

How introduced plants, which may be locally adapted to specific climatic conditions in their native range, cope with the new abiotic conditions that they encounter as exotics is not well understood. In particular, it is unclear what role plasticity versus adaptive evolution plays in enabling exotics to persist under new environmental circumstances in the introduced range. We determined the extent to which native and introduced populations of St. John's Wort (Hypericum perforatum) are genetically differentiated with respect to leaf-level morphological and physiological traits that allow plants to tolerate different climatic conditions. In common gardens in Washington and Spain, and in a greenhouse, we examined clinal variation in percent leaf nitrogen and carbon, leaf delta(13)C values (as an integrative measure of water use efficiency), specific leaf area (SLA), root and shoot biomass, root/shoot ratio, total leaf area, and leaf area ratio (LAR). As well, we determined whether native European H. perforatum experienced directional selection on leaf-level traits in the introduced range and we compared, across gardens, levels of plasticity in these traits. In field gardens in both Washington and Spain, native populations formed latitudinal clines in percent leaf N. In the greenhouse, native populations formed latitudinal clines in root and shoot biomass and total leaf area, and in the Washington garden only, native populations also exhibited latitudinal clines in percent leaf C and leaf delta(13)C. Traits that failed to show consistent latitudinal clines instead exhibited significant phenotypic plasticity. Introduced St. John's Wort populations also formed significant or marginally significant latitudinal clines in percent leaf N in Washington and Spain, percent leaf C in Washington, and in root biomass and total leaf area in the greenhouse. In the Washington common garden, there was strong directional selection among European populations for higher percent leaf N and leaf delta(13)C, but no selection on any other measured trait. The presence of convergent, genetically based latitudinal clines between native and introduced H. perforatum, together with previously published molecular data, suggest that native and exotic genotypes have independently adapted to a broad-scale variation in climate that varies with latitude.     

Evidence for a link between local and seasonal cycles in gene frequencies and latitudinal gene clines in a cyclic parthenogen

In an earlier study (Rhombergh et al., Can. J. Genet. Cytol. 27: 224-232, 1985) of natural populations of the cyclic parthenogenetic Rose aphids, Macrosiphum rosae, 6 out of 31 loci were found to be polymorphic and one locus (Esterase-4) showed cyclic seasonal changes in gene and genotypic frequencies. Assuming that the Est-4 polymorphism was balanced and due to some climatic factor that varies seasonally, and realizing that most environmental factors that vary seasonally also vary latitudinally, we predicted existence of a latitudinal gene cline at this locus. In the present study we surveyed four polymorphic loci (chosen to be used as markers) in six geographic populations spanning over 1200 km between the United States and Canada and found all four loci to have latitudinal clines. We think that the gene clines are due to a latitudinal cline in the degree of advancement of local populations through the seasonal cycle, and have called such a pattern a 'seasonal phase cline'. The results are discussed in relation to the temporal instability of local patterns and persistence of genetic variability on the large scale in aphids. It is argued that population structure of aphids makes retention of selectively neutral or weakly selected polymorphisms difficult.     

What determines conformity to Bergmann's rule?

Aim  Bergmann's rule, the tendency of body size within species in bird and mammal populations to be positively correlated with latitude, is among the best known biogeographical generalizations. The factors behind such clines, however, are not well understood. Here we use a large data base of 79 mammalian carnivore species to examine the factors affecting latitudinal size clines.
Location  Worldwide.
Methods  We measured the skulls and teeth of carnivores in natural history museums, and calculated the amount of variation in size explained by latitude, supplementing our measurements with published data. We examined the effects of a number of variables on the tendency to show latitudinal clines.
Results  We found that geographical range and latitudinal extent are strongly related to size clines. Minimum temperatures across the range, net primary productivity and habitat diversity also have some, albeit much less, influence.
Main conclusions  We suggest that species with large geographical ranges are likely to encounter significant heterogeneity in those factors that influence body size, and are thus likely to exhibit size clines. However, the key factors that determine body size may not always operate along a latitudinal (or other geographical) cline, but be spatially linked to patches in the species range. One such important factor is likely to be food availability, which we show is a strong predictor of size in the brown bear ( Ursus arctos ) but is not associated with a latitudinal cline. We argue that the spatial distribution of key resources within the species range constitutes a significant predictor of carnivore body size.     

Flight morphology along a latitudinal gradient in a butterfly: do geographic clines differ between agricultural and woodland landscapes?

Bergman and converse Bergman rules, amongst others, describe latitudinal variation in size of organisms, including flying ectotherms like butterflies. However, geographic clines in morphological traits of functional significance for flight performance and thermoregulation may also exist, although they have received less attention within a biogeographical context. Variation in flight‐related morphology has often been studied relative to landscape structure. However, the extent to which landscape effects interact with latitudinal clines of phenotypic variation has rarely been tested. Here we address the effect of latitude, landscape type and the interaction effect on body size and flight‐related morphology in the speckled wood butterfly Pararge aegeria. Male adult butterflies were collected from two replicate populations in each agricultural and woodland landscape types along a 700 km cline in six latitudinal zones. Overall size, adult body mass and wing area increased with latitude in line with Bergmann's rule. Forewing length, however, decreased with latitude. As predicted from thermoregulatory needs in ectotherms, the basal wing part was darker to the north. Latitudinal trends for flight‐related morphological traits were opposite to predictions about flight endurance under cooler conditions that were observed in some non‐lepidopteran insects, i.e. wing loading increased and wing aspect ratio decreased with latitude. Opposite trends can, however, be explained by other aspects of butterfly flight performance (i.e. mate‐location behaviour). As predicted from differences in environmental buffering in woodland landscapes along the latitudinal gradient, significant landscape×latitude interaction effects indicated stronger latitudinal clines and stronger phenotypic variation for size and flight morphology in the agricultural landscape compared to the woodland landscape. In agreement with significant interaction effects, morphological differentiation increased with latitude and was higher between population pairs of agricultural landscape than between population pairs of woodland landscape. These results demonstrate that landscape, latitude and their interaction contribute to the understanding of the complex geographic variation in P. aegeria adult phenotypes across Europe.