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.     

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.     

Clinal variation in Drosophila serrata for stress resistance and body size

Clines for size and stress resistance traits have been described for several Drosophila species and replicable clines across different species may indicate climatic selection. Here we consider clines in stress resistance traits in an Australian endemic species, D. serrata, by comparing levels of variation within and among isofemale lines initiated with flies collected from the eastern coast of Australia. We also consider clinical variation in chill coma recovery, a trait that has recently been shown to exhibit high levels of variation among Drosophila species. Patterns were compared with those in the cosmopolitan species D. melanogaster from the same area. Both desiccation and starvation resistance showed no clinical pattern despite heritable variation among isofemale lines. In contrast chill coma resistance exhibited a linear cline in the anticipated direction, resistance increasing with latitude. Body size was measured as wing length and body weight. Both traits showed geographic variation and strong non-linear clines with a sharp reduction in size in the tropics. These results are discussed in the context of climatic selection and evolutionary processes limiting species borders.     

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.     

Complex patterns of local adaptation in heat tolerance in Drosophila simulans from eastern Australia

Latitudinal clines are considered a powerful means of investigating evolutionary responses to climatic selection in nature. However, most clinal studies of climatic adaptation in Drosophila have involved species that contain cosmopolitan inversion polymorphisms that show clinal patterns themselves, making it difficult to determine whether the traits or inversions are under selection. Further, although climatic selection is unlikely to act on only one life stage in metamorphic organisms, a few studies have examined clinal patterns across life stages. Finally, clinal patterns of heat tolerance may also depend on the assay used. To unravel these potentially confounding effects on clinal patterns of thermal tolerance, we examined adult and larval heat tolerance traits in populations of Drosophila simulans from eastern Australia using static and dynamic (ramping 0.06 °C min^?1) assays. We also used microsatellites markers to clarify whether demographic factors or selection are responsible for population differentiation along clines. Significant cubic clinal patterns were observed for adult static basal, hardened and dynamic heat knockdown time and static basal heat survival in larvae. In contrast, static, hardened larval heat survival increased linearly with latitude whereas no clinal association was found for larval ramping survival. Significant associations between adult and larval traits and climatic variables, and low population differentiation at microsatellite loci, suggest a role for climatic selection, rather than demographic processes, in generating these clinal patterns. Our results suggest that adaptation to thermal stress may be species and life‐stage specific, complicating our efforts to understand the evolutionary responses to selection for increasing thermotolerance.     

Genome‐wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome‐wide next‐generation sequencing of DNA pools (‘pool‐seq’) from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome‐wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well‐known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.     

An independent non-linear latitudinal cline for the sn-glycerol-3-phosphate (alpha- Gpdh ) polymorphism of Drosophila melanogaster from eastern Australia

Latitudinal variation of the polymorphic sn-glycerol-3-phosphate (alpha-Gpdh) locus in Drosophila melanogaster has been characterized on several continents; however, apparent clinal patterns are potentially confounded by linkage with an inversion, close associations with other genetic markers that vary clinally, and a tandem alpha-Gpdh pseudogene. Here we compare clinal patterns in alpha-Gpdh with those of other linked markers by testing field flies from eastern Australian locations collected in two separate years. The alpha-Gpdh variation exhibited a consistent non-linear cline reflecting an increase in the alpha-GpdhF allele at extreme latitudes. This pattern was not influenced by the In(2L)t inversion wherein this locus is located, nor was it influenced by the presence of the alpha-Gpdh pseudogene, whose presence was ubiquitous and highly variable among populations. The alpha-Gpdh pattern was also independent of a cline in allozyme frequencies at the alcohol dehydrogenase (Adh) locus, and two length polymorphisms in the Adh gene. These results suggest clinal selection at the alpha-Gpdh locus that is partially or wholly unrelated to linear climatic gradients along the eastern coast of Australia.     

Inversions fail to account for allozyme clines

Allozyme and inversion data from natural populations of Drosophila melanogaster from the eastern United States were analyzed to determine whether the clines at allozyme loci are due to nonrandom associations with common cosmopolitan inversions. All inversions show strong clines. Clines were large and significant for half of the eight allozyme loci. An analysis of the contribution of inversions to clines of allozyme genes revealed three outcomes: the inversion cline (1) enhanced the allozyme cline, but was only partly responsible, (2) reduced the allozyme cline, and (3) had no effect. The allozyme clines were mainly determined by the pattern of allele frequencies within the chromosomal arrangements. Consequently, it was concluded that allozyme clines would exist in the absence of inversion clines.     

An inversion supergene in Drosophila underpins latitudinal clines in survival traits

Chromosomal inversions often contribute to local adaptation across latitudinal clines, but the underlying selective mechanisms remain poorly understood. We and others have previously shown that a clinal inversion polymorphism in Drosophila melanogaster, In(3R)Payne, underpins body size clines along the North American and Australian east coasts. Here, we ask whether this polymorphism also contributes to clinal variation in other fitness‐related traits, namely survival traits (lifespan, survival upon starvation and survival upon cold shock). We generated homokaryon lines, either carrying the inverted or standard chromosomal arrangement, isolated from populations approximating the endpoints of the North American cline (Florida, Maine) and phenotyped the flies at two growth temperatures (18 °C, 25 °C). Across both temperatures, high‐latitude flies from Maine lived longer and were more stress resistant than low‐latitude flies from Florida, as previously observed. Interestingly, we find that this latitudinal pattern is partly explained by the clinal distribution of the In(3R)P polymorphism, which is at ~ 50% frequency in Florida but absent in Maine: inverted karyotypes tended to be shorter‐lived and less stress resistant than uninverted karyotypes. We also detected an interaction between karyotype and temperature on survival traits. As In(3R)P influences body size and multiple survival traits, it can be viewed as a ‘supergene’, a cluster of tightly linked loci affecting multiple complex phenotypes. We conjecture that the inversion cline is maintained by fitness trade‐offs and balancing selection across geography; elucidating the mechanisms whereby this inversion affects alternative, locally adapted phenotypes across the cline is an important task for future work.     

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.     

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.     

Chromosome Inversion Polymorphisms in DROSOPHILA MELANOGASTER. I. Latitudinal Clines and Associations between Inversions in Australasian Populations

Nineteen Australasian populations of Drosophila melanogaster have been screened for chromosome inversion polymorphisms. All 15 of the inversion types found are paracentric and autosomal, but only four of these, one on each of the major autosome arms, are common and cosmopolitan. North-south clines occur, with the frequencies of all four of the common cosmopolitan inversions increasing toward the equator. These clines in the Southern Hemisphere mirror north-south clines in the Northern Hemisphere, where the frequencies of all four of the common cosmopolitan inversions again increase towards the equator.-While few of the Australasian populations show significant disequilibrium between linked common cosmopolitan inversions, those that do invariably have excesses of coupling gametes, which is consistent with other reports. We also find nonrandom associations between the two major autosomes, with the northern populations in Australasia (those with high inversion frequencies) tending to be deficient in gametes with common cosmopolitan inversions on both major autosomes, while the southern populations in Australasia (low inversion frequencies) tend to have an excess of this class of gametes.-The clines and the nonrandom associations between the two major autosomes are best interpreted in terms of selection operating to maintain the common cosmopolitan inversion polymorphisms in natural populations of D. melanogaster.     

The latitudinal cline in the In(3R)Payne inversion polymorphism has shifted in the last 20 years in Australian Drosophila melanogaster populations

Clinal variation has been described in a number of inversions in Drosophila but these clines are often characterized by cytological techniques using small sample sizes, and associations with specific genes are rarely considered. Here we have developed a molecular assay for In(3R)Payne in Drosophila melanogaster from eastern Australia populations. It shows in repeated samples that the inversion cline is very tightly associated with latitude and is almost fixed in tropical populations while relatively rare in temperate populations. This steep cline has shifted in position in the last 20 years. The heat shock gene, hsr-omega, located centrally inside the inversion sequence, shows a different clinal pattern to In(3R)Payne. These results suggest strong ongoing selection on In(3R)Payne over the last 100 years since the colonization of Australia that is partly independent of hsr-omega.     

Latitudinal variation for two enzyme loci and an inversion polymorphism in Drosophila melanogaster from Central and South America

Abstract.— Many organisms show latitudinal variation for various genetically determined traits. Such clines may involve neutral variation and originate from historical events or their maintenance may be explained by selection. For Drosophila melanogaster , latitudinal variation for allozymes, inversions, and quantitative traits has been found on several continents. We sampled D. melanogaster populations in Panama and along a transect of 40 latitudinal degrees on the west coast of South America. Negative correlations with latitude were found for Adh^s and _αGpdh^F allele frequencies and for the frequency of the cosmopolitan inversion In(2L)t in Adh^sαpdh^F chromosomes. A positive correlation existed between wing length and latitude. Significant correlations were found between these traits and climatic variables like temperature and rainfall. The observed clines show considerable resemblance to those found on other continents. Gametic disequilibrium between Adh^s and αGpdh^F occurred predominantly at higher latitudes and was caused by the presence of In(2L)t . The reasons for the clinal distributions are discussed and it is argued that selection is the most likely explanation. However, the exact nature of the selective force and the interactions of allozymes with each other and with In(2L)t are complex and not fully understood. In tropical regions In(2L)t -containing genotypes have higher fitness than ST/ST and Adh and αGpdh hitchhike with the inversion, but there is also evidence for balancing selection at the Adh locus.     

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.     

Direct and correlated responses to selection for longevity in Drosophila buzzatii

The possible associations between longevity, early fecundity, and stress-resistance traits were explored using artificial selection on longevity in a laboratory population of Drosophila buzzatii . Three replicated lines were selected for increased lifespan (L lines) and compared with the respective unselected controls (C lines) after the 14th generation of selection. Mean longevity exhibited a significant response to selection. The baseline mortality tended to decrease in the L lines and a negative correlated response to longevity selection was found for early fecundity. Egg-to-adult developmental time increased in L lines. Longevity selection increased stress resistance for both high and low temperatures, as measured by heat knockdown resistance and chill-coma recovery. Starvation resistance also tended to be higher in L than in C lines. The results obtained are consistent with the hypothesis of trade-offs between longevity and early fecundity, and also suggest a trade-off association between adult longevity and developmental time. Correlated selection responses were generally consistent with correlations among the traits previously inferred from altitudinal clines for longevity and stress-resistance phenotypes.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 738–748.     

QTL mapping reveals a striking coincidence in the positions of genomic regions associated with adaptive variation in body size in parallel clines of Drosophila melanogaster on different continents

Latitudinal genetic clines in body size are common in many ectotherm species and are attributed to climatic adaptation. Here, we use Quantitative Trait Loci (QTL) mapping to identify genomic regions associated with adaptive variation in body size in natural populations of Drosophila melanogaster from extreme ends of a cline in South America. Our results show that there is a significant association between the positions of QTL with strong effects on wing area in South America and those previously reported in a QTL mapping study of Australian cline end populations (P < 0.05). In both continents, the right arm of the third chromosome is associated with QTL with the strongest effect on wing area. We also show that QTL peaks for wing area and thorax length are associated with the same genomic regions, indicating that the clinal variation in the body size traits may have a similar genetic basis. The consistency of the results found for the South American and Australian cline end populations indicate that the genetic basis of the two clines may be similar and future efforts to identify the genes producing the response to selection should be focused on the genomic regions highlighted by the present work.     

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.     

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.     

Chromosomal inversion polymorphisms and adaptation

Chromosomal inversion polymorphisms continue to be identified from an increasing number of populations of insects, plants, bacteria and humans. In the fruit fly Drosophila, chromosomal polymorphisms were used in classic studies of natural selection. Recent molecular genetic studies suggest that inversion polymorphisms are dynamical systems. These studies also indicate patterns of disequilibrium and variation that are consistent with co-adapted gene complexes. Although these complexes have yet to be identified, recent studies have identified traits, such as body size, that are linked to inversion polymorphisms. Selection acting on these polymorphisms is strong because latitudinal clines in inversion frequency become re-established rapidly after a new continent is colonized. A combined molecular and phenotypic approach is helping to identify the role of inversion polymorphisms in adaptive divergence, but the genes responsible for associations between traits and inversion polymorphisms remain to be identified.