Wing trait–inversion associations in Drosophila subobscura can be generalized within continents,but may change through time

Clinal variation is one of the most emblematic examples of the action of natural selection at a wide geographical range. In Drosophila subobscura, parallel clines in body size and inversions, but not in wing shape, were found in Europe and South and North America. Previous work has shown that a bottleneck effect might be largely responsible for differences in wing trait–inversion association between one European and one South American population. One question still unaddressed is whether the associations found before are present across other populations of the European and South American clines. Another open question is whether evolutionary dynamics in a new environment can lead to relevant changes in wing traits–inversion association. To analyse geographical variation in these associations, we characterized three recently laboratory founded D. subobscura populations from both the European and South American latitudinal clines. To address temporal variation, we also characterized the association at a later generation in the European populations. We found that wing size and shape associations can be generalized across populations of the same continent, but may change through time for wing size. The observed temporal changes are probably due to changes in the genetic content of inversions, derived from adaptation to the new, laboratory environment. Finally, we show that it is not possible to predict clinal variation from intrapopulation associations. All in all this suggests that, at least in the present, wing traits–inversion associations are not responsible for the maintenance of the latitudinal clines in wing shape and size.     

Rapid differentiation of plasticity in life history and morphology during invasive range expansion and concurrent local adaptation in the horned beetle Onthophagus taurus

Understanding the interplay between genetic differentiation, ancestral plasticity, and the evolution of plasticity during adaptation to environmental variation is critical to predict populations’ responses to environmental change. However, the role of plasticity in rapid adaptation in nature remains poorly understood. We here use the invasion of the horned beetle Onthophagus taurus in the United States during the last half century to study the contribution of ancestral plasticity and post-invasion evolution of plastic responses in rapid population differentiation. We document latitudinal variation in life history and morphology, including genetic compensation in development time and body size, likely adaptive responses to seasonal constraints in the North. However, clinal variation in development time and size was strongly dependent on rearing temperature, suggesting that population differentiation in plasticity played a critical role in successful adaptation on ecological timescales. Clinal variation in wing shape was independent of ancestral plasticity, but correlated with derived plasticity, consistent with evolutionary interdependence. In contrast, clinal variation in tibia shape aligned poorly with thermal plasticity. Overall, this study suggests that post-invasion evolution of plasticity contributed to range expansions and concurrent adaptation to novel climatic conditions.     

Swift laboratory thermal evolution of wing shape (but not size) in Drosophila subobscura and its relationship with chromosomal inversion polymorphism

Latitudinal clinal variation in wing size and shape has evolved in North American populations of Drosophila subobscura within about 20 years since colonization. While the size cline is consistent to that found in original European populations (and globally in other Drosophila species), different parts of the wing have evolved on the two continents. This clearly suggests that 'chance and necessity' are simultaneously playing their roles in the process of adaptation. We report here rapid and consistent thermal evolution of wing shape (but not size) that apparently is at odds with that suggestion. Three replicated populations of D. subobscura derived from an outbred stock at Puerto Montt (Chile) were kept at each of three temperatures (13, 18 and 22 degrees C) for 1 year and have diverged for 27 generations at most. We used the methods of geometric morphometrics to study wing shape variation in both females and males from the thermal stocks, and rates of genetic divergence for wing shape were found to be as fast or even faster than those previously estimated for wing size on a continental scale. These shape changes did not follow a neat linear trend with temperature, and are associated with localized shifts of particular landmarks with some differences between sexes. Wing shape variables were found to differ in response to male genetic constitution for polymorphic chromosomal inversions, which strongly suggests that changes in gene arrangement frequencies as a response to temperature underlie the correlated changes in wing shape because of gene-inversion linkage disequilibria. In fact, we also suggest that the shape cline in North America likely predated the size cline and is consistent with the quite different evolutionary rates between inversion and size clines. These findings cast strong doubts on the supposed 'unpredictability' of the geographical cline for wing traits in D. subobscura North American colonizing populations.     

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

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

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.     

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.     

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.     

Testing evolutionary hypotheses about species borders: patterns of genetic variation towards the southern borders of two rainforest Drosophila and a related habitat generalist

Several evolutionary hypotheses help explain why only some species adapt readily to new conditions and expand distributions beyond borders, but there is limited evidence testing these hypotheses. In this study, we consider patterns of neutral (microsatellite) and quantitative genetic variation in traits in three species of Drosophila from the montium species group in eastern Australia. We found little support for restricted or asymmetrical gene flow in any species. In rainforest-restricted Drosophila birchii, there was evidence of selection for increased desiccation and starvation resistance towards the southern border, and a reduction in genetic diversity in desiccation resistance at this border. No such patterns existed for Drosophila bunnanda, which has an even more restricted distribution. In the habitat generalist Drosophila serrata, there was evidence for geographic selection for wing size and development time, although clinal patterns for increased cold and starvation resistance towards the southern border could not be differentiated from neutral expectations. These findings suggest that borders in these species are not limited by low overall genetic variation but instead in two of the species reflect patterns of selection and genetic variability in key traits limiting borders.     

In search of clinal variation in the period and clock timing genes in Australian Drosophila melanogaster populations

Clinal variation for repeat number in the Thr-Gly region of the period circadian timing gene in Drosophila melanogaster was described in Europe and has subsequently been used as evidence of thermal selection on period alleles. To test for clinal variation in this gene along the east coast of Australia, the period polymorphism was scored on flies from multiple samples collected repeatedly over a 5-year interval, along with variation at another circadian rhythm locus, clock. For period, there was no consistent evidence of clinal variation in the 17 and/or 20 repeat alleles, although when average allele length was examined a weak consistent clinal pattern was detected. For clock there was no evidence of clinal variation in the two most common alleles or in average repeat size. These data are inconsistent with the reported patterns in Europe and suggest that clinal variation in timing genes needs to be re-examined in this region.     

Genetic mapping of adaptive wing size variation in Drosophila simulans

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

Thermal evolution of pre-adult life history traits, geometric size and shape, and developmental stability in Drosophila subobscura

Replicated lines of Drosophila subobscura originating from a large outbred stock collected at the estimated Chilean epicentre (Puerto Montt) of the original New World invasion were allowed to evolve under controlled conditions of larval crowding for 3.5 years at three temperature levels (13, 18 and 22 degrees C). Several pre-adult life history traits (development time, survival and competitive ability), adult life history related traits (wing size, wing shape and wing-aspect ratio), and wing size and shape asymmetries were measured at the three temperatures. Cold-adapted (13 degrees C) populations evolved longer development times and showed lower survival at the highest developmental temperature. No divergence for wing size was detected following adaptation to temperature extremes (13 and 22 degrees C), in agreement with earlier observations, but wing shape changes were obvious as a result of both thermal adaptation and development at different temperatures. However, the evolutionary trends observed for the wing-aspect ratio were inconsistent with an adaptive hypothesis. There was some indication that wing shape asymmetry has evolutionarily increased in warm-adapted populations, which suggests that there is additive genetic variation for fluctuating asymmetry and that it can evolve under rapid environmental changes caused by thermal stress. Overall, our results cast strong doubts on the hypothesis that body size itself is the target of selection, and suggest that pre-adult life history traits are more closely related to thermal adaptation.     

Causes of variation in wing loading among Drosophila species

We examined influences on wing and body size in 11 species (12 strains) of Drosophila. Six measures of wing length and width were closely correlated with wing area and suggested little variation in wing shape among the species. Among ten species wing loading, an important factor in flight costs and manoeuvrability, increased as body mass increased at a rate consistent with expectations from allometric scaling of wing area and body mass to body length. Intraspecific variation in wing loading showed similar relationships to body mass. Density and temperature during larval development influenced wing loading through general allometric relations of body size and wing area. Temperature during the pupal stage, but not during wing hardening after eclosion, influenced wing area independently of body size. Wing area increased as growth temperature decreased. Individuals reared at cooler temperatures thus compensated for a potential allometric increase in wing loading by differentially enlarging the wing area during pupal development.     

Adaptation and constraint in the evolution of Drosophila melanogaster wing shape

SUMMARY We have taken advantage of parallel instances of natural selection on body size in Drosophila melanogaster to investigate constraints and adaptation affecting wing shape. Using recently developed techniques for statistical shape analysis, we have examined variation in wing shape in similar body size clines on three continents. Gender-related shape differences were constant among all populations, suggesting that gender differences represent a developmental constraint on wing shape. In contrast, the underlying shape varied significantly between continents and shape change within each cline (i.e., between small and large body size populations) also varied between continents. Therefore, variation at these two levels presumably results from either drift or natural selection. Functional considerations suggest that shape variation between the continents is unlikely to be adaptive. However, cline-related shape change, which we show has a significant allometric component, may be adaptive. The overall range of wing shape variation, across a large range of wing size, is extremely small, and the possibility that wing shape is subject to stabilizing selection (or canalization) is discussed.     

Levels of variation in stress resistance in drosophila among strains, local populations, and geographic regions: patterns for desiccation, starvation, cold resistance, and associated traits

Stress resistance traits in Drosophila often show clinal variation. Although these patterns suggest selection, there is generally no attempt to test how large differences at the geographical level are relative to levels of variation within and between local populations. Here we compare these levels in D. melanogaster from temperate Tasmania versus tropical northern Queensland by focusing on adult resistance to desiccation, cold and starvation stress, as well as associated traits (size, lipid content). For starvation and desiccation resistance, levels of variation were highest among strains from the same population. whereas there was little differentiation among local populations and a low level of differentiation at the geographic level. For adult cold resistance, there was local differentiation and strain variation but no geographic variation. For size (thorax length), geographic differentiation was higher despite some overlap among strains from the tropical and temperate locations. Finally, for lipid levels there was only evidence for variation among strains. The low level of differentiation among geographic locations for stress resistance was further verified with the characterization of isofemale strains from 18 locations along a coastal transect extending from Tasmania to northern Queensland. Crosses among some of the isofemale strains showed that results were not confounded by inbreeding effects. Strains derived from a cross between a tropical and temperate strain differed for all traits, and variation among strains for body size was higher than strain variation within the geographic regions. Unlike in previous studies, lipid content and starvation resistance were not correlated in any set of strains, but there was a correlation between cold resistance and lipid content. There was also a correlation between desiccation resistance and size but only in the geographic cross strains. These findings suggest a large amount of variation in stress resistance at the population level and inconsistent correlation patterns across experimental approaches.     

Genetic architecture of wing morphology in Drosophila simulans and an analysis of temperature effects on genetic parameter estimates

The Drosophila wing has been used as a model to investigate the mechanisms responsible for size and shape changes in nature, since such changes might underlie morphological evolution. To improve the understanding of wing morphological variation and the interpretation of genetic parameters estimates, we have established 59 lines from a Drosophila simulans laboratory population through single pair random matings. The offspring of each line were reared at three different temperatures, and the wing morphology of 12 individuals was analyzed by adjusting an ellipse to the wings' contour. Temperature, sex and line significantly affected wing trait variation, which was mainly characterized by longer wings having the second, fourth and fifth longitudinal veins closer together at the wing tip. As for the genetic parameter estimates, while the cross-environment heritability of some traits, such as wing size (SI), decreased with an increasing difference between the temperatures at which parents and offspring were reared, wing shape (SH) heritability did not seem to change. Since we found indications that neither an increase in the phenotypic variation nor the occurrence of genotype-environment interactions could fully explain the low heritabilities of SI estimated by cross-environment regressions, we discuss the importance of other effects for explaining this discrepancy between the SI and SH heritability estimates. In addition, although the genetic matrix was not entirely represented in the phenotypic matrix, several correspondences were identified, suggesting that the observed patterns of wing morphology variation are genetically controlled.     

SIMILARITIES AND DIFFERENCES IN ALTITUDINAL VERSUS LATITUDINAL VARIATION FOR MORPHOLOGICAL TRAITS IN DROSOPHILA MELANOGASTER

Understanding how natural environments shape phenotypic variation is a major aim in evolutionary biology. Here, we have examined clinal, likely genetically based variation in morphology among 19 populations of the fruit fly (Drosophila melanogaster) from Africa and Europe, spanning a range from sea level to 3000 m altitude and including locations approximating the southern and northern range limit. We were interested in testing whether latitude and altitude have similar phenotypic effects, as has often been postulated. Both latitude and altitude were positively correlated with wing area, ovariole number, and cell number. In contrast, latitude and altitude had opposite effects on the ratio between ovariole number and body size, which was negatively correlated with egg production rate per ovariole. We also used transgenic manipulation to examine how increased cell number affects morphology and found that larger transgenic flies, due to a higher number of cells, had more ovarioles, larger wings, and, unlike flies from natural populations, increased wing loading. Clinal patterns in morphology are thus not a simple function of changes in body size; instead, each trait might be subject to different selection pressures. Together, our results provide compelling evidence for profound similarities as well as differences between phenotypic effects of latitude and altitude.     

Quantifying the relative contributions of the X chromosome,autosomes, and mitochondrial genome to local adaptation*

During local adaptation with gene flow, some regions of the genome are inherently more responsive to selection than others. Recent theory predicts that X‐linked genes should disproportionately contribute to local adaptation relative to other genomic regions, yet this prediction remains to be tested. We carried out a multigeneration crossing scheme, using two cline‐end populations of Drosophila melanogaster, to estimate the relative contributions of the X chromosome, autosomes, and mitochondrial genome to divergence in four traits involved in local adaptation (wing size, resistance to heat, desiccation, and starvation stresses). We found that the mitochondrial genome and autosomes contributed significantly to clinal divergence in three of the four traits. In contrast, the X made no significant contribution to divergence in these traits. Given the small size of the mitochondrial genome, our results indicate that it plays a surprisingly large role in clinal adaptation. In contrast, the X, which represents roughly 20% of the Drosophila genome, contributes negligibly—a pattern that conflicts with theoretical predictions. These patterns reinforce recent work implying a central role of mitochondria in climatic adaptation, and suggest that different genomic regions may play fundamentally different roles in processes of divergence with gene flow.     

Rapid evolution of wing size clines in Drosophila subobscura

Parallel latitudinal clines across species and continents provide dramatic evidence of the efficacy of natural selection, however little is known about the dynamics involved in cline formation. For example, several drosophilids and other ectotherms increase in body and wing size at higher latitudes. Here we compare evolution in an ancestral European and a recently introduced (North America) cline in wing size and shape in Drosophila subobscura. We show that clinal variation in wing size, spanning more than 15 degrees of latitude, has evolved in less than two decades. In females from Europe and North America, the clines are statistically indistinguishable however the cline for North American males is significantly shallower than that for European males. We document that while overall patterns of wing size are similar on two continents, the European cline is obtained largely through changing the proximal portion of the wing, whereas the North American cline is largely in the distal portion. We use data from sites collected in 1986/1988 (Pegueroles et al. 1995) and our 1997 collections to compare synchronic (divergence between contemporary populations that share a common ancestor) and allochronic (changes over time within a population) estimates of the rates of evolution. We find that, for these populations, allochronically estimated evolutionary rates within a single population are over 0.02 haldanes (2800 darwins), a value similar in magnitude to the synchronic estimates from the extremes of the cline. This paper represents an expanded analysis of data partially presented in Huey et al. (2000).