Starvation resistance and adult body composition in a latitudinal cline of Drosophila melanogaster

Abstract Latitudinal geographic variation in Drosophila melanogaster is pervasive. Parallel clines in traits such as body size, egg size, ovariole number, and development time have been found on several continents throughout the world. However, a cline in starvation resistance and fat content in D. melanogaster has so far been found only in India. Here we investigate starvation resistance and fat content in 10 populations from South America, in which clines in body size, egg size, and development time have previously been found. We find no evidence for a cline in starvation resistance or fat content in South America. We therefore suggest that the cline in starvation resistance in India may have evolved in response to specific climatic variation found only in India.     

A comprehensive assessment of geographic variation in heat tolerance and hardening capacity in populations of Drosophila melanogaster from eastern Australia

We examined latitudinal variation in adult and larval heat tolerance in Drosophila melanogaster from eastern Australia. Adults were assessed using static and ramping assays. Basal and hardened static heat knockdown time showed significant linear clines; heat tolerance increased towards the tropics, particularly for hardened flies, suggesting that tropical populations have a greater hardening response. A similar pattern was evident for ramping heat knockdown time at 0.06 °C min^?1 increase. There was no cline for ramping heat knockdown temperature (CT_max) at 0.1 °C min^?1 increase. Acute (static) heat knockdown temperature increased towards temperate latitudes, probably reflecting a greater capacity of temperate flies to withstand sudden temperature increases during summer in temperate Australia. Larval viability showed a quadratic association with latitude under heat stress. Thus, patterns of heat resistance depend on assay methods. Genetic correlations in thermotolerance across life stages and evolutionary potential for critical thermal limits should be the focus of future studies.     

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

A large‐scale latitudinal pattern of life‐history traits in a strictly univoltine damselfly

1. Variation in thermal conditions and season length along latitudinal gradients affect body size‐related traits over different life stages. Selection is expected to optimise these size traits in response to the costs and benefits. 2. Egg, hatchling, larval and adult size in males and females were estimated along a latitudinal gradient of 2730 km across Europe in the univoltine damselfly Lestes sponsa, using a combination of field‐collection and laboratory‐rearing experiments. In the laboratory, individuals were grown in temperatures and photoperiod simulating those at the latitude of origin, and in common‐garden conditions. 3. The size of adults sampled in nature was negatively correlated with latitude. In all populations the females were larger than the males. Results from simulated and common‐garden rearing experiments supported this pattern of size difference across latitudes and between sexes, suggesting a genetic component for the latitudinal size trend and female‐biased size dimorphism. In contrast, hatchling size showed a positive relationship with latitude, but egg size, although differing between latitudes, showed no such relationship. 4. The results support a converse Bergmann cline, i.e. a negative body size cline towards the north. This negative cline in body size is probably driven by progressively stronger seasonal time and temperature constraints towards the higher latitudes and by the obligate univoltine life cycle of L. sponsa. As egg size showed no relationship with latitude, other environmental factors besides temperature, such as desiccation risk, probably affect this trait.     

Polymorphism in the couch potato gene clines in eastern Australia but is not associated with ovarian dormancy in Drosophila melanogaster

Natural selection can generate parallel latitudinal clines in traits and gene frequencies across continents, but these have rarely been linked. An amino acid (isoleucine to lysine, or I462K) polymorphism of the couch potato (cpo) gene in Drosophila melanogaster is thought to control female reproductive diapause cline in North America (Schmidt et al. 2008, Proc Natl Acad Sci USA, 105, 16207-16211). Here, we show that under standard diapause-inducing conditions (12 °C and short photoperiod) (Saunders et al. 1989, Proc Natl Acad Sci USA, 86, 3748-3752), egg maturation in Australian flies is delayed, but not arrested at previtellogenic stages. At 12 °C, the phenotypic distribution in egg development was bimodal at stages 8 and 14 and showed a strong nonlinear pattern on the east coast of Australia, with incidence of egg maturation delay (ovarian dormancy) increasing both toward tropical and temperate climates. Furthermore, we found no evidence for an association between the cpo I462K polymorphism and ovarian dormancy at either 12 or 10 °C (when egg maturation was often delayed at stage 7). Owing to strong linkage disequilibrium, the latitudinal cline in cpo allele frequencies was no longer evident once variation in the In(3R)P inversion polymorphism was taken into account. Our results suggest that the standard diapause-inducing conditions (12 °C and short photoperiod) were not sufficient to cause the typical previtellogenic developmental arrest in Australian flies and that the cpo I462K polymorphism does not explain the observed delay in egg development. In conclusion, ovarian dormancy does not show a simple latitudinal cline, and the lack of cpo-dormancy association suggests a different genetic basis to reproductive dormancy in North America and Australia.     

Wolbachia transmission in a naturally bi-infected Drosophila simulans strain from New-Caledonia

Wolbachia are maternally-transmitted endocellular bacteria infecting several arthropod species. In order to study the possibility of Wolbachia segregation in a naturally bi-infected host, isofemale lines from a bi-infected Drosophila simulans (Sturtevant) strain from Nouméa (New Caledonia) were backcrossed using uninfected males carrying the same nuclear background. Uninfected males were used to avoid the cytoplasmic incompatibility syndrome (CI) associated with the presence of Wolbachia in males. Each line was established using a female infected simultaneously by the two different Wolbachia variants wHa and wNo. The backcross led to some individuals carrying only one type of infection being recovered among the progeny of the bi-infected foundress females. Rarely, uninfected individuals were also recovered. Isolated for the first time in its natural host, wNo exhibited a significantly weaker CI phenotype than the isolated wHa variant. Infection fate when backcross conditions were relaxed varied depending on rearing conditions of the host. Under favourable conditions, the infection was generally maintained, while it was frequently lost under unfavourable conditions. This result probably reflects the direct fitness dependence of the symbiont on its host.     

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.     

Reproductive characteristics of the flower breeding Drosophila hibisci Bock (Drosophilidae) along a latitudinal gradient in eastern Australia: relation to flower and habitat features

For many insect species, egg and larval substrate characteristics are significantly correlated with interspecific differences in female reproductive allocation and egg size-number tradeoffs. We tested the hypothesis that a similar pattern occurred within the Australian drosophilid, Drosophila hibisci, that is restricted throughout its life cycle to flowers of species in the genus Hibiscus. These plants occur as small, isolated, normally monospecific stands that should facilitate differentiation of the fly populations in relation to specific oviposition and larval substrates. Data from 38 sites ranging from 20.8̀ to 34.4̀ S latitude in eastern Australia indicated no relationship between female body size, egg size, or ovariole numbers and floral size or mass among four species of Hibiscus. However, the flies did show a latitudinal cline in ovariole number that was independent of floral variation. Females averaged 15–20 ovarioles per female in the south (32–34̀ S latitude) and 10–12 ovarioles in the north (21–22̀ S latitude). The increase in ovariole number with latitude was due to divergence in the ovariole number of the largest females. In contrast, small females in the north and south had the same number of ovarioles. Reproductive allocation of female flies in the northern region was less than females in the southern region. The latitudinal divergence in ovariole number was not associated with habitat differences (density of trees, density of flies and beetles), nor with differences in floral characteristics (flower weight, petal length, yeast species present). Short term weather patterns in daily temperature and rainfall preceding collections pardy explain the variation in ovariole number. These observations in conjunction with preliminary genetic results suggest the cline is associated with genetic differences that interact with environmental determinants such as the temperature during larval development.     

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.     

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.     

Presence of an invasive species reverses latitudinal clines of multiple traits in a native species

Understanding the processes driving formation and maintenance of latitudinal clines has become increasingly important in light of accelerating global change. Many studies have focused on the role of abiotic factors, especially temperature, in generating clines, but biotic factors, including the introduction of non‐native species, may also drive clinal variation. We assessed the impact of invasion by predatory fire ants on latitudinal clines in multiple fitness‐relevant traits—morphology, physiological stress responsiveness, and antipredator behavior—in a native fence lizard. In areas invaded by fire ants, a latitudinal cline in morphology is opposite both the cline found in museum specimens from historical populations across the species’ full latitudinal range and that found in current populations uninvaded by fire ants. Similarly, clines in stress‐relevant hormone response to a stressor and in antipredator behavior differ significantly between the portions of the fence lizard range invaded and uninvaded by fire ants. Changes in these traits within fire ant‐invaded areas are adaptive and together support increased and more effective antipredator behavior that allows escape from attacks by this invasive predator. However, these changes may mismatch lizards to the environments under which they historically evolved. This research shows that novel biotic pressures can alter latitudinal clines in multiple traits within a single species on ecological timescales. As global change intensifies, a greater understanding of novel abiotic and biotic pressures and how affected organisms adapt to them across space and time will be central to predicting and managing our changing environment.     

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.     

No patterns in thermal plasticity along a latitudinal gradient in Drosophila simulans from eastern Australia

Phenotypic plasticity may be an important initial mechanism to counter environmental change, yet we know relatively little about the evolution of plasticity in nature. Species with widespread distributions are expected to have evolved higher levels of plasticity compared with those with more restricted, tropical distributions. At the intraspecific level, temperate populations are expected to have evolved higher levels of plasticity than their tropical counterparts. However, empirical support for these expectations is limited. In addition, no studies have comprehensively examined the evolution of thermal plasticity across life stages. Using populations of Drosophila simulans collected from a latitudinal cline spanning the entire east coast of Australia, we assessed thermal plasticity, measured as hardening capacity (the difference between basal and hardened thermal tolerance) for multiple measures of heat and cold tolerance across both adult and larval stages of development. This allowed us to explicitly ask whether the evolution of thermal plasticity is favoured in more variable, temperate environments. We found no relationship between thermal plasticity and latitude, providing little support for the hypothesis that temperate populations have evolved higher levels of thermal plasticity than their tropical counterparts. With the exception of adult heat survival, we also found no association between plasticity and ten climatic variables, indicating that the evolution of thermal plasticity is not easily predicted by the type of environment that a particular population occupies. We discuss these results in the context of the role of plasticity in a warming climate.     

Sex-ratio distortion in Drosophila simulans: co-occurence of a meiotic drive and a suppressor of drive

A sex-ratio distortion factor was found at high frequency in D. simulans strains from Seychelles and New Caledonia. This factor is poorly or not expressed within those strains which are resistant to it. Its presence was detected by crossing females from New Caledonia or the Seychelles with males from a different geographic origin. Most of the F1 males obtained produced an excess of females (up to 99%) in their progeny. The two strains are infected with Wolbachia, but these micro-organisms are not involved in the sex-ratio distortion. The sex-ratio factor is shown to be an X-linked meiotic driver; nuclear resistance factor(s) act by suppressing the drive. It is likely that the same X-located driver invaded the two populations, which subsequently developed resistance factor(s) against it.     

Evidence for association of chromosomal form and development time from complex clines and geographic races in the grasshopper Caledia captiva (Orthoptera: Acridadae)

Evidence for an association between chromosomal form and development time in the grasshopper Caledia captiva (F.) was obtained through comparison of two geographic taxa and analysis of a complex latitudinal cline within one of the taxa. Northern populations of the Moreton taxon possess a metacentric genome and are slow-developing. In contrast, the Torresian taxon, distributed throughout northern, coastal Australia, a region of pronounced seasonality in rainfall, and southern populations of the Moreton taxon, which inhabit a region of pronounced seasonality in temperature, Soth have an acrocentric genome and are fast-developing. The convergence of chromosomal form and development time between Torresian and southern Moreton populations appears to be driven by convergence in life history. Seasonality limits grasshoppers to one generation per year and favours fast development. The transition between relatively acrocentric southern Moreton populations and relatively metacentric northern Moreton populations is gradual but not monotonic. Instead, a shift to a bivoltine life history in the middle of the transect occurs and is associated with shifts in both development time and chromosomal form. These results imply an adaptive role for chromosomal form, although the causative link between chromosomal variation and variation in development time remains to be established.     

Sequential evolution of a symbiont inferred from the host: Wolbachia and Drosophila simulans

This study aims to unravel the biogeography of a model symbiont/host system by exploiting the prediction that a symbiont will leave a signature of infection on the host. Specifically, a global sample of 1,442 Drosophila simulans from 33 countries and 64 sampling localities was employed to infer the phylogeography of the maternally inherited alpha-proteobacteria Wolbachia. Phylogenetic analyses, from three symbiont genes and 24 mtDNA genomes (excluding the A + T-rich region), showed that each of four Wolbachia strains infected D. simulans once. The global distribution and abundance of the Wolbachia strains and the three mtDNA haplogroups (D. simulans siI, siII and siIII) was then determined. Finally, network analyses of variable regions within siI (584 bp from seven additional lines) and siII (1,701 bp from 383 lines) facilitated a detailed biogeographic discussion. There is little variation in siIII and the haplogroup is restricted in its distribution. These data show how the history of an infection can be mapped by combining data from the symbiont and the host. They say little about the organismal history of the host because the mtDNA genome is a biased representation of the whole genome.     

Oviposition preference and life history traits in cactophilic Drosophila koepferae and D. buzzatii in association with their natural hosts

Drosophila koepferae and D. buzzatii are two closely related cactophilic species inhabiting the arid lands of southern South America. Previous studies have shown that D. buzzatii breeds primarily on the necrotic cladodes of several Opuntia cacti and D. koepferae on the rotting stems of columnar cacti of the genera Trichocereus and Cereus. In this paper, we analyze the patterns of host plant utilization in a locality where both Drosophila species are sympatric. Field studies showed an absence of differential attraction of adult flies to the rots of two major host cacti: O. sulphurea and T. terschekii. However, the proportion of D. buzzatii flies emerged from the rotting cladodes of O. sulphurea was significantly higher than in T. terschekii. In laboratory experiments, egg to adult viability in single species cultures varied when both Drosophila species were reared in media prepared with O. sulphurea or T. terschekii. In addition, between-species comparisons of flies emerged from single species cultures showed that D. buzzatii adults were smaller and developed faster than D. koepferae. Furthermore, analysis of flies emerged in mixed species cultures showed differences in oviposition preference and oviposition behavior. We discuss the observed between-species differences and suggest that these traits are the result of adaptation to specific patterns of spatial and temporal predictability of their respective preferred host plants: columnar are less dense and less ephemeral resources, whereas the opuntias are more abundant, and fast rotting cacti. This revised version was published online in July 2006 with corrections to the Cover Date.     

Microsatellite variation and differentiation in African and non-African populations of Drosophila simulans

Drosophila simulans originated in sub-Saharan Africa or Madagascar and colonized the rest of the world after the last glaciation about 10 000 years ago. Consistent with this demographic history, sub-Saharan African populations have been shown to harbour higher levels of microsatellite and sequence variation than cosmopolitan populations. Nevertheless, only limited information is available on the population structure of D. simulans. Here, we analysed X-linked and autosomal microsatellite loci in four sub-Saharan African, one North African, one Israeli, and two European D. simulans populations. Bayesian clustering algorithms combined the North African, Israeli, and European populations into a single cosmopolitan group. The four sub-Saharan populations were split into two separate groups. Pairwise F(ST) analysis, however, indicated significant population differentiation between all eight populations surveyed. A significant signal for population reduction in cosmopolitan populations was found only for X-linked loci.     

Intraspecific variation in survival and mitochondrial oxidative phosphorylation in wild-caught Drosophila simulans

Lifespans of organisms vary greatly even among individuals of the same species. Under the framework of the free oxygen radical theory of aging, it is predicted that variation in individual lifespan within a species will correlate with variation in the accumulation of oxidative damage to cell components from reactive oxygen species. In this study we test the hypothesis that variation in survival of three wild-caught Drosophila simulans fly lines (HW09, NC48 and MD106) correlates with three key aspects of mitochondrial bioenergetics. The rank order of median survival was HW09 > MD106 > NC48. Young HW09 flies (11-18 days) had (i) highest ADP:O (quantity of oxygen consumed by mitochondria when provided with a quantity of ADP) when metabolizing both electron transport chain complex I and complex III substrates; (ii) lowest rate of mitochondrial hydrogen peroxide production from complex III; and (iii) highest cytochrome c oxidase activity from complex IV. Rate of hydrogen peroxide production increased and cytochrome c oxidase activity decreased in all lines in the age range 11-25 days. This is the first study to correlate natural variation in organism survival with natural variation in mitochondrial bioenergetics.