Developmental acclimation affects clinal variation in stress resistance traits in Drosophila buzzatii

Patterns of clinal genetic variation in Drosophila are often characterized after rearing at constant temperatures. However, clinal patterns might change after acclimation if populations differ in their plastic response to fluctuating environments. We studied longevity, starvation and heat knock‐down resistance after development at either constant or fluctuating temperatures in nine Drosophila buzzatii populations collected along an altitudinal gradient in Tenerife, Spain. Flies that developed at fluctuating temperatures had higher stress resistance despite experiencing a slightly lower average temperature than those at constant temperatures. Genetic variation along the gradient was found in both stress‐resistance traits. Because Q_ST values greatly exceeded F_ST values, genetic drift could not explain this diversification. In general, differences among populations were larger after rearing at fluctuating temperatures, especially in heat knock‐down, for which clinal patterns disappeared when flies were reared at constant temperatures. This result emphasizes the importance of determining whether populations originating from different environments differ in their plastic responses to stress.     

Genetic covariances promote climatic adaptation in Australian Drosophila*

Evolutionary potential for adaptation hinges upon the orientation of genetic variation for traits under selection, captured by the additive genetic variance-covariance matrix (G), as well as the evolutionary stability of G. Yet studies that assess both the stability of G and its alignment with selection are extraordinarily rare. We evaluated the stability of G in three Drosophila melanogaster populations that have adapted to local climatic conditions along a latitudinal cline. We estimated population- and sex-specific G matrices for wing size and three climatic stress-resistance traits that diverge adaptively along the cline. To determine how G affects evolutionary potential within these populations, we used simulations to quantify how well G aligns with the direction of trait divergence along the cline (as a proxy for the direction of local selection) and how genetic covariances between traits and sexes influence this alignment. We found that G was stable across the cline, showing no significant divergence overall, or in sex-specific subcomponents, among populations. G also aligned well with the direction of clinal divergence, with genetic covariances strongly elevating evolutionary potential for adaptation to climatic extremes. These results suggest that genetic covariances between both traits and sexes should significantly boost evolutionary responses to environmental change.     

The genetic variance but not the genetic covariance of life‐history traits changes towards the north in a time‐constrained insect

Seasonal time constraints are usually stronger at higher than lower latitudes and can exert strong selection on life‐history traits and the correlations among these traits. To predict the response of life‐history traits to environmental change along a latitudinal gradient, information must be obtained about genetic variance in traits and also genetic correlation between traits, that is the genetic variance‐covariance matrix, G . Here, we estimated G for key life‐history traits in an obligate univoltine damselfly that faces seasonal time constraints. We exposed populations to simulated native temperatures and photoperiods and common garden environmental conditions in a laboratory set‐up. Despite differences in genetic variance in these traits between populations (lower variance at northern latitudes), there was no evidence for latitude‐specific covariance of the life‐history traits. At simulated native conditions, all populations showed strong genetic and phenotypic correlations between traits that shaped growth and development. The variance–covariance matrix changed considerably when populations were exposed to common garden conditions compared with the simulated natural conditions, showing the importance of environmentally induced changes in multivariate genetic structure. Our results highlight the importance of estimating variance–covariance matrixes in environments that mimic selection pressures and not only trait variances or mean trait values in common garden conditions for understanding the trait evolution across populations and environments.     

Testing limits to adaptation along altitudinal gradients in rainforest Drosophila

Given that evolution can generate rapid and dramatic shifts in the ecological tolerance of a species, what prevents populations adapting to expand into new habitat at the edge of their distributions? Recent population genetic models have focused on the relative costs and benefits of migration between populations. On the one hand, migration may limit adaptive divergence by preventing local populations from matching their local selective optima. On the other hand, migration may also contribute to the genetic variance necessary to allow populations to track these changing optima. Empirical evidence for these contrasting effects of gene flow in natural situations are lacking, largely because it remains difficult to acquire. Here, we develop a way to explore theoretical models by estimating genetic divergence in traits that confer stress resistance along similar ecological gradients in rainforest Drosophila. This approach allows testing for the coupling of clinal divergence with local density, and the effects of genetic variance and the rate of change of the optimum on the response to selection. In support of a swamping effect of migration on phenotypic divergence, our data show no evidence for a cline in stress-related traits where the altitudinal gradient is steep, but significant clinal divergence where it is shallow. However, where clinal divergence is detected, sites showing trait means closer to the presumed local optimum have more genetic variation than sites with trait means distant from their local optimum. This pattern suggests that gene flow also aids a sustained response to selection.     

Evolution of additive and nonadditive genetic variance in development time along a cline in Drosophila serrata

Abstract. Latitudinal clines provide natural systems that may allow the effect of natural selection on the genetic variance to be determined. Ten clinal populations of Drosophila serrata collected from the eastern coast of Australia were used to examine clinal patterns in the trait mean and genetic variance of the life-history trait egg-to-adult development time. Development time significantly lengthened from tropical areas to temperate areas. The additive genetic variance for development time in each population was not associated with latitude but was associated with the population mean development time. Additive genetic variance tended to be larger in populations with more extreme development times and appeared to be consistent with allele frequency change. In contrast, the nonadditive genetic variance was not associated with the population mean but was associated with latitude. Levels of nonadditive genetic variance were greatest in the region of the cline where the gradient in the change in mean was greatest, consistent with Barton's (1999) conjecture that the generation of linkage disequilibrium may become an important component of the genetic variance in systems with a spatially varying optimum.     

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

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

Geographic variation in diapause incidence, life-history traits, and climatic adaptation in Drosophila melanogaster

In Drosophila melanogaster, exposure of females to low temperature and shortened photoperiod can induce the expression of reproductive quiescence or diapause. Diapause expression is highly variable within and among natural populations and has significant effects on life-history profiles, including patterns of longevity, fecundity, and stress resistance. We hypothesized that if diapause expression is associated with overwintering mechanisms and adaptation to temperate environments, the frequency of diapause incidence would exhibit a latitudinal cline among natural populations. Because stress resistance and reproductive traits are also clinal in this species, we also examined how patterns of fecundity and longevity varied with geography and how stress resistance and associated traits differed constitutively between diapause and nondiapause lines. Diapause incidence was shown to vary predictably with latitude, ranging from 35% to 90% among natural populations in the eastern United States Survivorship under starvation stress differed between diapause and nondiapause lines; diapause phenotypes were also distinct for total body triglyceride content and the developmental distribution of oocytes in the ovary following stress exposure. Patterns of longevity, fecundity, and ovariole number also varied with geography. The data suggest that, for North American populations, diapause expression is functionally associated with overwintering mechanisms and may be an integral life-history component in natural populations.     

The effect of developmental temperature on the genetic architecture underlying size and thermal clines in Drosophila melanogaster and D. simulans from the east coast of Australia

Body size and thermal tolerance clines in Drosophila melanogaster occur along the east coast of Australia. However the extent to which temperature affects the genetic architecture underlying the observed clinal divergence remains unknown. Clinal variation in these traits is associated with cosmopolitan chromosome inversions that cline in D. melanogaster. Whether this association influences the genetic architecture for these traits in D. melanogaster is unclear. Drosophila simulans shows linear clines in body size, but nonlinear clines in cold resistance. Clinally varying inversions are absent in D. simulans. Line-cross and clinal analyses were performed between tropical and temperate populations of D. melanogaster and D. simulans from the east coast of Australia to investigate whether clinal patterns and genetic effects contributing to clinal divergence in wing centroid size, thorax length, wing-to-thorax ratio, cold and heat resistance differed under different developmental temperatures (18 °C, 25 °C, and 29 °C). Developmental temperature influenced the genetic architecture in both species. Similarities between D. melanogaster and D. simulans suggest clinally varying inversion polymorphisms have little influence on the genetic architecture underlying clinal divergence in size in D. melanogaster. Differing genetic architectures across different temperatures highlight the need to consider different environments in future evolutionary and molecular studies of phenotypic divergence.     

Ecological genomics of Anopheles gambiae along a latitudinal cline: a population-resequencing approach

The association between fitness-related phenotypic traits and an environmental gradient offers one of the best opportunities to study the interplay between natural selection and migration. In cases in which specific genetic variants also show such clinal patterns, it may be possible to uncover the mutations responsible for local adaptation. The malaria vector, Anopheles gambiae, is associated with a latitudinal cline in aridity in Cameroon; a large inversion on chromosome 2L of this mosquito shows large differences in frequency along this cline, with high frequencies of the inverted karyotype present in northern, more arid populations and an almost complete absence of the inverted arrangement in southern populations. Here we use a genome resequencing approach to investigate patterns of population divergence along the cline. By sequencing pools of individuals from both ends of the cline as well as in the center of the cline- where the inversion is present in intermediate frequency- we demonstrate almost complete panmixia across collinear parts of the genome and high levels of differentiation in inverted parts of the genome. Sequencing of separate pools of each inversion arrangement in the center of the cline reveals large amounts of gene flux (i.e., gene conversion and double crossovers) even within inverted regions, especially away from the inversion breakpoints. The interplay between natural selection, migration, and gene flux allows us to identify several candidate genes responsible for the match between inversion frequency and environmental variables. These results, coupled with similar conclusions from studies of clinal variation in Drosophila, point to a number of important biological functions associated with local environmental adaptation.     

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

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

LOCAL ADAPTATION MAINTAINS CLINAL VARIATION IN MELANIN‐BASED COLORATION OF EUROPEAN BARN OWLS (TYTO ALBA)

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin‐based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish‐brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.     

INTEGRATING LANDSCAPE GENOMICS AND SPATIALLY EXPLICIT APPROACHES TO DETECT LOCI UNDER SELECTION IN CLINAL POPULATIONS

Uncovering the genetic basis of adaptation hinges on the ability to detect loci under selection. However, population genomics outlier approaches to detect selected loci may be inappropriate for clinal populations or those with unclear population structure because they require that individuals be clustered into populations. An alternate approach, landscape genomics, uses individual‐based approaches to detect loci under selection and reveal potential environmental drivers of selection. We tested four landscape genomics methods on a simulated clinal population to determine their effectiveness at identifying a locus under varying selection strengths along an environmental gradient. We found all methods produced very low type I error rates across all selection strengths, but elevated type II error rates under “weak” selection. We then applied these methods to an AFLP genome scan of an alpine plant, Campanula barbata, and identified five highly supported candidate loci associated with precipitation variables. These loci also showed spatial autocorrelation and cline patterns indicative of selection along a precipitation gradient. Our results suggest that landscape genomics in combination with other spatial analyses provides a powerful approach for identifying loci potentially under selection and explaining spatially complex interactions between species and their environment.     

Genetic architecture and phenotypic plasticity of thermally-regulated traits in an eruptive species, <Emphasis Type="Italic">Dendroctonus ponderosae</Emphasis>

Phenotypic plasticity in thermally-regulated traits enables close tracking of changing environmental conditions, and can thereby enhance the potential for rapid population increase, a hallmark of outbreak insect species. In a changing climate, exposure to conditions that exceed the capacity of existing phenotypic plasticity may occur. Combining information on genetic architecture and trait plasticity among populations that are distributed along a latitudinal cline can provide insight into how thermally-regulated traits evolve in divergent environments and the potential for adaptation. Dendroctonus ponderosae feed on Pinus species in diverse climatic regimes throughout western North America, and show eruptive population dynamics. We describe geographical patterns of plasticity in D. ponderosae development time and adult size by examining reaction norms of populations from multiple latitudes. The relative influence of additive and non-additive genetic effects on population differences in the two phenotypic traits at a single temperature is quantified using line-cross experiments and joint-scaling tests. We found significant genetic and phenotypic variation among D. ponderosae populations. Simple additive genetic variance was not the primary source of the observed variation, and dominance and epistasis contributed greatly to the genetic divergence of the two thermally-regulated traits. Hybrid breakdown was also observed in F₂ hybrid crosses between northern and southern populations, further indication of substantial genetic differences among clinal populations and potential reproductive isolation within D. ponderosae. Although it is unclear what maintains variation in the life-history traits, observed plasticity in thermally-regulated traits that are directly linked to rapid numerical change may contribute to the outbreak nature of D. ponderosae, particularly in a changing climate.     

Genetic and phenotypic sources of life history variation along a cline in voltinism in the cricket Allonemobius socius

Clinal variation in life histories can be genetically based, resulting from selection imposed by different environments, or it may be due to the differential expression of phenotypically plastic traits. We examined the cline in voltinism in the egg-diapausing cricket Allonemobius socius, with populations spanning the switch from a univoltine to a bivoltine phenology. A common garden experiment was employed, using environments that mimicked photoperiod and temperature conditions found in the field. There were only small differences in development time among populations, and the difference in phenology observed in the field is likely due to clinal variation in the length of the growing season. We found large genetically-based differences in the reaction norm for egg diapause that were further magnified by environmental cues. The synergism of genetic and environmental effects was an example of cogradient selection. In the zone of transition between phenologies, voltinism appeared to be a conditional strategy, rather than a genetic polymorphism. First-generation females from this area can lay either direct-developing or diapause eggs depending on the likelihood that a second generation will have sufficient time to develop. For this species, the cline in voltinism is the result of a combination of environmental effects on development, and genetic and environmental influences on egg diapause propensity.     

Correlated changes in thermotolerance traits and body color phenotypes in montane populations of Drosophila melanogaster: analysis of within- and between-population variations

Thermotolerance traits vary across geographical gradients but there is a lack of clinal variation in some Drosophila species. Thus, it is not clear whether thermotolerance or other correlated traits are the target of natural selection. In order to test selection responses, we investigated body melanization and thermotolerance traits in six altitudinal populations of Drosophila melanogaster . Based on rearing different geographical populations under uniform growth conditions at 21 °C (common garden experiments), clinal variations for cold resistance are in the direction opposite to heat resistance along an altitudinal gradient, that is darker flies from highland populations evidenced higher levels of cold resistance while lowland populations showed higher heat resistance. Phenotypic plastic responses for body melanization at 17–28 °C showed significant correlations with thermotolerance traits. At 17 °C, regression coefficients as a function of altitude are highly significant and positive for cold resistance but negative for heat knockdown. However, for flies reared at 28 °C, there is no elevational change in melanization as well as thermotolerance traits. Thus, both genetic and plastic changes of body melanization and thermotolerance traits suggest a correlated selection response. Further, within-population analyses of body melanization (based on dark, intermediate and light color phenotypes) showed significant associations with thermotolerance traits. Correlated variations in body melanization and thermal tolerances are associated with climatic thermal variability ( T _cv) but not with T _min. or T _max. along an altitudinal gradient.     

Latitudinal clines in body size,but not in thermal tolerance or heat‐shock cognate 70 (HSC70), in the highly‐dispersing intertidal gastropod Littorina keenae (Gastropoda: Littorinidae)

Natural populations of widely‐distributed animals often exhibit clinal variation in phenotypic traits or in allele frequencies of a particular gene over their geographical range. A planktotrophic intertidal snail, Littorina keenae is broadly distributed along the north‐eastern Pacific coast through a large latitudinal range (24°50′N–43°18′N). We tested for latitudinal clines in two complex phenotypic traits – thermal tolerance and body size – and one single locus trait – heat shock cognate 70 (HSC70) – in L. keenae along almost its entire geographical range. We found only weak evidence for a latitudinal cline in the thermal tolerance and no evidence for a cline in allele frequencies at HSC70. However, as predicted by Bergmann's rule, we detected a strong latitudinal cline that accounted for 60% of the variance in body size (R² = 0.598; P < 0.001). In contrast, body size did not significantly affect thermal tolerance. HSC70 showed no genetic differentiation among the populations, supporting our previous mitochondrial gene‐based estimate of high gene flow during this snail's free‐swimming larval stage. Given that L. keenae experiences panmixia along its species range, the observed size cline may be partially or entirely caused by a phenotypically plastic response to local thermal environments rather than by genetic divergence in body size among populations in response to locally optimizing natural selection. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 494–505.     

Population structure in the barn swallow,Hirundo rustica: a comparison between neutral DNA markers and quantitative traits

Local adaptation to variable environments can generate clinal variation in morphological traits. Alternatively, similar patterns of clinal variation may be generated simply as a result of genetic drift/migration balance. Teasing apart these different processes is a continuing focus in evolutionary ecology. We compare genetic differentiation at molecular loci and quantitative traits to analyse the effect of these different processes in a morphological latitudinal cline of the barn swallow, Hirundo rustica, breeding across Europe. The results obtained show no structuring at neutral microsatellite loci, which contrasts with positive structuring at five quantitative morphometric traits. This supports the hypothesis that the observed morphometric cline in barn swallows is the result of selection acting in a spatially heterogeneous environment. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 306–314.     

Natural phenological variation in aspen (Populus tremula): the SwAsp collection

The genus Populus is currently the main model system for genetic, genomic, and physiological research in trees. Phenotypic variation in aspen (Populus tremula) populations growing in different environments across Sweden is expected to reflect genetic variation that is important for local adaptation. To analyze such natural phenotypic and genetic variation, the Swedish Aspen (SwAsp) Collection was established. Trees were taken from 12 different populations across Sweden, from 56° to 66° latitude north and planted in two common gardens in Ekebo (55.9°N) and Sävar (63.4°N). Data related to phenological and growth traits were collected during the second year of growth. Some traits like the date of bud set and leaf area duration showed strong clinal variation patterns with latitude in both field trials, but the date of bud flush did not change along a latitudinal cline. The phenological traits showed moderate within-populations heritabilities, although growth traits showed weaker clinal patterns and lower heritabilities than the phenological traits. This research forms the starting point for the development of the SwAsp collection, a resource facilitating analysis of the natural genetic variation in aspen, the elucidation of the structure and dynamics of aspen populations, and the future identification of the genes controlling adaptive traits using association mapping of selected candidate genes.     

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

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