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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Selective sweeps in a 2-locus model for sex-ratio meiotic drive in Drosophila simulans

A way to identify loci subject to positive selection is to detect the signature of selective sweeps in given chromosomal regions. It is revealed by the departure of DNA polymorphism patterns from the neutral equilibrium predicted by coalescent theory. We surveyed DNA sequence variation in a region formerly identified as causing "sex-ratio" meiotic drive in Drosophila simulans. We found evidence that this system evolved by positive selection at 2 neighboring loci, which thus appear to be required simultaneously for meiotic drive to occur. The 2 regions are approximately 150-kb distant, corresponding to a genetic distance of 0.1 cM. The presumably large transmission advantage of chromosomes carrying meiotic drive alleles at both loci has not erased the individual signature of selection at each locus. This chromosome fragment combines a high level of linkage disequilibrium between the 2 critical regions with a high recombination rate. As a result, 2 characteristic traits of selective sweeps--the reduction of variation and the departure from selective neutrality in haplotype tests--show a bimodal pattern. Linkage disequilibrium level indicates that, in the natural population from Madagascar used in this study, the selective sweep may be as recent as 100 years.     

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.     

Clinal variation and laboratory adaptation in the rainforest species Drosophila birchii for stress resistance, wing size, wing shape and development time

Clinal variation has been described in many invertebrates including drosophilids but usually over broad geographical gradients. Here we describe clinal variation in the rainforest species Drosophila birchii from Queensland, Australia, and potential confounding effects of laboratory adaptation. Clinal variation was detected for starvation and development time, but not for size or resistance to temperature extremes. Starvation resistance was higher at southern locations. Wing shape components were not associated with latitude although they did differ among populations. Time in laboratory culture did not influence wing size or heat knockdown resistance, but increased starvation resistance and decreased recovery time following a cold shock. Laboratory culture also increased development time and altered wing shape. The results indicate that clinal patterns can be detected in Drosophila over a relatively narrow geographical area. Laboratory adaptation is unlikely to have confounded the detection of geographical patterns.     

Desiccation tolerance and starvation resistance exhibit opposite latitudinal clines in Indian geographical populations of Drosophila kikkawai

Abstract. 1. Desiccation tolerance and starvation resistance demonstrated significant differentiation among seven Indian geographical populations of Drosophila kikkawai, collected along a latitudinal range of 12.6–32.7 °N. Lack of significant differences in two successive generations suggested that these physiological traits were genetically controlled. 2. North Indian populations of D. kikkawai displayed significantly higher desiccation tolerance than southern populations, whereas there was a reverse trend for starvation tolerance (r > 0.90). Regression slope values indicated an increase of 0.61 h for desiccation and a decrease of 1.71 h per degree latitude for starvation tolerance at 17 °C. The traits evidenced opposite latitudinal clines, and such data also matched thermal climatic conditions on the Indian subcontinent. The survival duration for such traits was significantly higher at 17 than at 25 °C. 3. Significantly higher starvation tolerance in south Indian populations might be due to large population size, species interactions, and higher metabolic rates in the humid tropical environments. In contrast, prolonged unfavourable colder climatic conditions are known to favour starvation tolerance in temperate regions. Thus, the causes of desiccation and starvation tolerance seem quite different under tropical and temperate conditions. 4. Starvation tolerance was correlated negatively with body weight and ovariole number, which might be due to a trade-off in favour of greater allocation to non-lipidic reserves for sustaining starvation tolerance in the tropics. Reduction in metabolic rate may not be applicable for observed higher starvation tolerance in the tropical populations. 5. Multiple regression analysis demonstrated a major effect of coefficient of variation of mean monthly temperature for both the traits of ecological significance. Thus, Indian geographical populations of D. kikkawai provided evidence of independent genetic divergence for starvation and desiccation tolerance under natural conditions.     

Resistance to fipronil in Drosophila simulans: influence of two point mutations in the RDL GABA receptor subunit

The Eyguieres 42 strain of Drosophila simulans, obtained by laboratory selection, displayed approximately 20,000-fold resistance to the insecticide fipronil. Molecular cloning of the cDNA encoding the RDL GABA receptor subunit of this strain revealed the presence of two mutations: the Rdl mutation (A301G) and an additional mutation in the third transmembrane domain (T350M). In order to assess the individual and combined roles of the two mutations in fipronil resistance, the functional properties of wild-type, A301G, T350M and A301G/T350M homomultimeric RDL receptors were compared by expression in Xenopus oocytes. In wild-type receptors, the inhibition of GABA (EC(30))-induced currents by fipronil and picrotoxin was enhanced by repeated GABA applications. The A301G mutation nearly abolished this effect, decreased the sensitivity to fipronil and picrotoxin and increased the reversibility of inhibition. The T350M mutation also reduced the sensitivity to both antagonists. Of the four receptor variants tested, the double mutant showed the highest resistance to fipronil, following repeated GABA applications. In conclusion, the present study emphasizes new aspects of the pharmacological alterations induced by the Rdl mutation and shows that resistance to GABA receptor-directed insecticides may implicate a mutation distinct from Rdl.     

Selection for adult desiccation resistance in Drosophila melanogaster: fitness components, larval resistance and stress correlations

In previous experiments we found that Drosophila melanogaster lines selected for increased adult desiccation resistance had increased resistance to other environmental stresses at the adult stage including starvation, intense ⁶⁰Co-γ radiation and a toxic ethanol level. In further studies on these lines, we now show that selection did not alter resistance to desiccation and ethanol at the larval stage. As well as having a lower early fecundity, selected lines showed increased adult male longevity and increased viability at high larval densities compared with control lines. There were no changes in development time or mating success. The increased male longevity is consistent with the reduced metabolic rate of the selected lines.
A genetic correlation between resistance to different stresses was confirmed by an analysis of isofemale lines derived from a population founded by flies from a stress-resistant line and an unselected line. The results are consistent with the existence of genes segregating in natural populations conferring increased general stress resistance.     

Altitudinal patterns for latitudinally varying traits and polymorphic markers in Drosophila melanogaster from eastern Australia

Altitudinal changes in traits and genetic markers can complement the studies on latitudinal patterns and provide evidence of natural selection because of climatic factors. In Drosophila melanogaster, latitudinal variation is well known but altitudinal patterns have rarely been investigated. Here, we examine five traits and five genetic markers on chromosome 3R in D. melanogaster collected at high and low altitudes from five latitudes along the eastern coast of Australia. Significant altitudinal differentiation was observed for cold tolerance, development time, ovariole number in unmated females, and the microsatellite marker DMU25686. Differences tended to match latitudinal patterns, in that trait values at high altitudes were also found at high latitudes, suggesting that factors linked to temperature are likely selective agents. Cold tolerance was closely associated with average temperature and other climatic factors, but no significant associations were detected for the other traits. Genes around DMU25686 represent good candidates for climatic adaptation.     

Identification of a candidate adaptive polymorphism for Drosophila life history by parallel independent clines on two continents

Life history traits are critical components of fitness and frequently reflect adaptive responses to environmental pressures. However, few genes that contribute to natural life history variation have been identified. Insulin signalling mediates the determination of life history traits in many organisms, and single gene manipulation in Drosophila melanogaster suggests that individual genes in the pathway have the potential to produce major effects on these quantitative traits. We evaluated allelic variation at two insulin signalling genes, the Insulin‐like Receptor (InR) and its substrate, chico, in natural populations of D. melanogaster. We found different patterns of variation: InR shows evidence of positive selection and clines in allele frequency across latitude; chico exhibits neutral patterns of evolution. The clinal patterns at InR are replicated between North America and Australia, showing striking similarity in the distribution of specific alleles and the rate at which allele frequencies change across latitude. Moreover, we identified a polymorphism at InR that appears to be functionally significant and consistent with hypothetical patterns of selection across geography. This polymorphism provides new characterization of genic regions of functionality within InR, and is likely a component in a suite of genes and traits that respond adaptively to climatic variation.     

Latitudinal clines in<Emphasis Type="Italic">Drosophila melanogaster</Emphasis>: Body size,allozyme frequencies,inversion frequencies,and the insulin-signalling pathway

Many latitudinal clines exist inDrosophila melanogaster: in adult body size, in allele frequency at allozyme loci, and in frequencies of common cosmopolitan inversions. The question is raised whether these latitudinal clines are causally related. This review aims to connect data from two very different fields of study, evolutionary biology and cell biology, in explaining such natural genetic variation inD. melanogaster body size and development time. It is argued that adult body size clines, inversion frequency clines, and clines in allele frequency at loci involved in glycolysis and glycogen storage are part of the same adaptive strategy. Selection pressure is expected to differ at opposite ends of the clines. At high latitudes, selection onD. melanogaster would favour high larval growth rate at low temperatures, and resource storage in adults to survive winter. At low latitudes selection would favour lower larval critical size to survive crowding, and increased male activity leading to high male reproductive success. Studies of the insulin-signalling pathway inD. melanogaster point to the involvement of this pathway in metabolism and adult body size. The genes involved in the insulin-signalling pathway are associated with common cosmopolitan inversions that show latitudinal clines. Each chromosome region connected with a large common cosmopolitan inversion possesses a gene of the insulin transmembrane complex, a gene of the intermediate pathway and a gene of the TOR branch. The hypothesis is presented that temperateD. melanogaster populations have a higher frequency of a ’thrifty’ genotype corresponding to high insulin level or high signal level, while tropical populations possess a more ’spendthrift’ genotype corresponding to low insulin or low signal level.     

Latitudinal variation in wild populations of Drosophila melanogaster: heritabilities and reaction norms

Large amounts of genetic variation for wing length and wing area were demonstrated both within and between Drosophila melanogaster populations along a latitudinal gradient in South America. Wing length and wing area showed a strong positive correlation with latitude in both wild flies and laboratory-raised descendants. Large population differences were observed for heritability and coefficient of variation of these two traits, whereas relatively small population differences were found for development time, viability, pupal mortality, sex ratio and their norms of reaction to four developmental temperatures. No clear-cut latitudinal clines were established for these life-history characters. These results are discussed in the light of Bergmann's Rule and the relation between larval development and adult body size.     

Natural adaptation to environmental stress via physiological clock-regulation of stress resistance in Drosophila

In the present study we describe a novel mechanism by which populations, in an interplay between physiology and behaviour, can evolve adaptation to environmental extremes. Comparing Drosophila populations from different climate zones, we found diel shifts in high temperature resistance after maintenance at 25 °C for several generations. Resistance changes that co-occurred with field and laboratory activity of the populations are controlled by the physiological clock and appear to be a consequence of local adaptation to the thermal profiles of the environment of origin.     

Origin and decay of the P element-associated latitudinal cline in Australian Drosophila melanogaster

The latitudinal cline in P transposable element-associated characteristics in eastern Australian populations of Drosophila melanogaster has changed between 1986 and 1991–1994. New collections were made in 1991–1994 from localities along the eastern coast of Australia. P element-associated properties of 256 isofemale lines from 43 localities were evaluated using gonadal dysgenesis and/or singed-weak hypermutability assays. The overall results indicate that both P activity and P susceptibility have declined, with all populations showing a tendency towards a state with little P activity potential but with P repressor function (neutral or ‘Q’). P repressor function is strong in all populations except some of the most southerly. P activity potential peaks at about 27° SLat, and drops off to the south (as in 1983–1986 collections) and to the north (in contrast to 1983–1986 collections); thus the cline is no longer a simple P-to-Q-to-M pattern from north to south, but is now Q-P-Q-M. A mtDNA RFLP that putatively distinguishes North American and European populations varies in frequency among the populations but the frequency does not vary clinally with latitude, ruling out massive introductions from North America and Europe as causing the cline. This revised version was published online in July 2006 with corrections to the Cover Date.