Variation in wing length in Eurasian natural populations of Drosophila melanogaster

Selection and geography shape male reproductive tract transcriptomes in Drosophila melanogaster

Transcriptome analysis of several animal clades suggests that male reproductive tract gene expression evolves quickly. However, the factors influencing the abundance and distribution of within-species variation, the ultimate source of interspecific divergence, are poorly known. Drosophila melanogaster, an ancestrally African species that has recently spread throughout the world and colonized the Americas in the last roughly 100 years, exhibits phenotypic and genetic latitudinal clines on multiple continents, consistent with a role for spatially varying selection in shaping its biology. Nevertheless, geographic expression variation in the Americas is poorly described, as is its relationship to African expression variation. Here, we investigate these issues through the analysis of two male reproductive tissue transcriptomes [testis and accessory gland (AG)] in samples from Maine (USA), Panama, and Zambia. We find dramatic differences between these tissues in differential expression between Maine and Panama, with the accessory glands exhibiting abundant expression differentiation and the testis exhibiting very little. Latitudinal expression differentiation appears to be influenced by the selection of Panama expression phenotypes. While the testis shows little latitudinal expression differentiation, it exhibits much greater differentiation than the accessory gland in Zambia vs American population comparisons. Expression differentiation for both tissues is non-randomly distributed across the genome on a chromosome arm scale. Interspecific expression divergence between D. melanogaster and D. simulans is discordant with rates of differentiation between D. melanogaster populations. Strongly heterogeneous expression differentiation across tissues and timescales suggests a complex evolutionary process involving major temporal changes in the way selection influences expression evolution in these organs.

Seasonal variation in wing size and shape of Drosophila melanogaster reveals rapid adaptation to environmental changes

Populations in seasonal fluctuating environments receive multiple environmental cues and must deal with this heterogenic environment to survive and reproduce. An enlarged literature shows that this situation can be resolved through rapid adaptation in Drosophila melanogaster populations. Long-term monitoring of a population in its natural habitat and quantitative measurement of its responses to seasonal environmental changes are important for understanding the adaptive response of D. melanogaster to temporal variable selection. Here, we use inbred lines of a D. melanogaster population collected at monthly intervals between May to October over a temporal scale spanning three consecutive years to understand the variation in wing size and wing shape over these timepoints. The wing size and shape of this population changed significantly between months and a seasonal cycle of this traits is repeated for three years. Our results suggest that the effects of environmental variables that generated variation in body size between populations such as latitudinal clines, are a selective pressure in a different manner in terms of seasonal variation. Temperature related variable have a significant nonlinear relation to this fluctuating pattern in size and shape, whereas precipitation and humidity have a sex-specific effect which is more significant in males.

Conserved and Divergent Aspects of Plasticity and Sexual Dimorphism in Wing Size and Shape in Three Diptera

The ability of powered flight in insects facilitated their great evolutionary success allowing them to occupy various ecological niches. Beyond this primary task, wings are often involved in various premating behaviors, such as the generation of courtship songs and the initiation of mating in flight. These specific functions imply special adaptations of wing morphology, as well as sex-specific wing morphologies. Although wing morphology has been extensively studied in Drosophila melanogaster (Meigen, 1830), a comprehensive understanding of developmental plasticity and the impact of sex on wing size and shape plasticity is missing for other Diptera. Therefore, we raised flies of the three Diptera species Drosophila melanogaster, Ceratitis capitata (Wiedemann, 1824) and Musca domestica (Linnaeus, 1758) at different environmental conditions and applied geometric morphometrics to analyze wing shape. Our data showed extensive interspecific differences in wing shape, as well as a clear sexual wing shape dimorphism in all three species. We revealed an impact of different rearing temperatures on wing shape in all three species, which was mostly explained by plasticity in wing size in D. melanogaster. Rearing densities had significant effects on allometric wing shape in D. melanogaster, while no obvious effects were observed for the other two species. Additionally, we did not find evidence for sex-specific response to different rearing conditions in D. melanogaster and C. capitata, while a male-specific impact of different rearing conditions was observed on non-allometric wing shape in M. domestica. Overall, our data strongly suggests that many aspects of wing morphology underly species-specific adaptations and we discuss potential developmental and functional implications of our results.

Constraints and Opportunities for the Evolution of Metamorphic Organisms in a Changing Climate

We argue that developmental hormones facilitate the evolution of novel phenotypic innovations and timing of life history events by genetic accommodation. Within an individual’s life cycle, metamorphic hormones respond readily to environmental conditions and alter adult phenotypes. Across generations, the many effects of hormones can bias and at times constrain the evolution of traits during metamorphosis; yet, hormonal systems can overcome constraints through shifts in timing of, and acquisition of tissue specific responses to, endocrine regulation. Because of these actions of hormones, metamorphic hormones can shape the evolution of metamorphic organisms. We present a model called a developmental goblet, which provides a visual representation of how metamorphic organisms might evolve. In addition, because developmental hormones often respond to environmental changes, we discuss how endocrine regulation of postembryonic development may impact how organisms evolve in response to climate change. Thus, we propose that developmental hormones may provide a mechanistic link between climate change and organismal adaptation.

Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication

Background

The highly polyphagous Queensland fruit fly (Bactrocera tryoni Froggatt) expanded its range substantially during the twentieth century and is now the most economically important insect pest of Australian horticulture, prompting intensive efforts to develop a Sterile Insect Technique (SIT) control program. Using a “common garden” approach, we have screened for natural genetic variation in key environmental fitness traits among populations from across the geographic range of this species and monitored changes in those traits induced during domestication.

Results

Significant variation was detected between the populations for heat, desiccation and starvation resistance and wing length (as a measure of body size). Desiccation resistance was correlated with both starvation resistance and wing length. Bioassay data for three resampled populations indicate that much of the variation in desiccation resistance reflects persistent, inherited differences among the populations. No latitudinal cline was detected for any of the traits and only weak correlations were found with climatic variables for heat resistance and wing length. All three stress resistance phenotypes and wing length changed significantly in certain populations with ongoing domestication but there was also a strong population by domestication interaction effect for each trait.

Conclusions

Ecotypic variation in heat, starvation and desiccation resistance was detected in Australian Qfly populations, and these stress resistances diminished rapidly during domestication. Our results indicate a need to select source populations for SIT strains which have relatively high climatic stress resistance and to minimise loss of that resistance during domestication.

Life-History Evolution and the Genetics of Fitness Components in Drosophila melanogaster

Life-history traits or "fitness components"-such as age and size at maturity, fecundity and fertility, age-specific rates of survival, and life span-are the major phenotypic determinants of Darwinian fitness. Analyzing the evolution and genetics of these phenotypic targets of selection is central to our understanding of adaptation. Due to its simple and rapid life cycle, cosmopolitan distribution, ease of maintenance in the laboratory, well-understood evolutionary genetics, and its versatile genetic toolbox, the "vinegar fly" Drosophila melanogaster is one of the most powerful, experimentally tractable model systems for studying "life-history evolution." Here, I review what has been learned about the evolution and genetics of life-history variation in D. melanogaster by drawing on numerous sources spanning population and quantitative genetics, genomics, experimental evolution, evolutionary ecology, and physiology. This body of work has contributed greatly to our knowledge of several fundamental problems in evolutionary biology, including the amount and maintenance of genetic variation, the evolution of body size, clines and climate adaptation, the evolution of senescence, phenotypic plasticity, the nature of life-history trade-offs, and so forth. While major progress has been made, important facets of these and other questions remain open, and the D. melanogaster system will undoubtedly continue to deliver key insights into central issues of life-history evolution and the genetics of adaptation.

Testing the drivers of the temperature-size covariance using artificial selection

Body size often declines with increasing temperature. Although there is ample evidence for this effect to be adaptive, it remains unclear whether size shrinking at warmer temperatures is driven by specific properties of being smaller (e.g., surface to volume ratio) or by traits that are correlated with size (e.g., metabolism, growth). We used 290 generations (22 months) of artificial selection on a unicellular phytoplankton species to evolve a 13-fold difference in volume between small-selected and large-selected cells and tested their performance at 22°C (usual temperature), 18°C (-4), and 26°C (+4). Warmer temperatures increased fitness in small-selected individuals and reduced fitness in large-selected ones, indicating changes in size alone are sufficient to mediate temperature-dependent performance. Our results are incompatible with the often-cited geometric argument of warmer temperature intensifying resource limitation. Instead, we find evidence that is consistent with larger cells being more vulnerable to reactive oxygen species. By engineering cells of different sizes, our results suggest that smaller-celled species are pre-adapted for higher temperatures. We discuss the potential repercussions for global carbon cycles and the biological pump under climate warming.

No evidence that warmer temperatures are associated with selection for smaller body sizes

Reductions in animal body size over recent decades are often interpreted as an adaptive evolutionary response to climate warming. However, for reductions in size to reflect adaptive evolution, directional selection on body size within populations must have become negative, or where already negative, to have become more so, as temperatures increased. To test this hypothesis, we performed traditional and phylogenetic meta-analyses of the association between annual estimates of directional selection on body size from wild populations and annual mean temperatures from 39 longitudinal studies. We found no evidence that warmer environments were associated with selection for smaller size. Instead, selection consistently favoured larger individuals, and was invariant to temperature. These patterns were similar in ectotherms and endotherms. An analysis using year rather than temperature revealed similar patterns, suggesting no evidence that selection has changed over time, and also indicating that the lack of association with annual temperature was not an artefact of choosing an erroneous time window for aggregating the temperature data. Although phenotypic trends in size will be driven by a combination of genetic and environmental factors, our results suggest little evidence for a necessary ingredient-negative directional selection-for declines in body size to be considered an adaptive evolutionary response to changing selection pressures.

Directional selection reduces developmental canalization against genetic and environmental perturbations in Drosophila wings

Natural selection may enhance or weaken the robustness of phenotypes against genetic or environmental perturbations. However, important aspects of the relationship between adaptive evolution and canalization remain unclear. Recent work showed that the evolution of larger wing size in a high altitude natural population of Drosophila melanogaster was accompanied by decanalized wing development--specifically a loss of robustness to genetic perturbation. But this study did not address environmental robustness, and it compared populations that may have numerous biological differences. Here, we perform artificial selection on this same trait in D. melanogaster (larger wing length) and directly test whether this directional selection resulted in decanalization. We find that in general, size-selected replicates show greater frequencies of wing defects than control replicates both after mutagenesis (genetic perturbation) and when subjected to high temperature stress (environmental perturbation), although the increase in defect frequency varies importantly among replicates. These results support the hypothesis that directional selection may result in the loss of both genetic and environmental robustness-offering a rare window into the relationship between adaptation and canalization.

Response of life-history traits to artificial and natural selection for virulence and nonvirulence in a Drosophila parastitoid, Asobara tabida

Co-evolution of host-parasitoid interactions is determined by the costs of host resistance, which received empirical evidence, and the costs of parasitoid virulence, which have been mostly hypothesized. Asobara tabida is a parasitoid, which mainly parasitizes Drosophila melanogaster and D. subobscura, the first species being able to resist to the parasitoid development while the second species is not. To parasitize resistant hosts, including D. melanogaster, A. tabida develops sticky eggs, which prevent encapsulation, but this virulence mechanism may be costly. Interindividual and interpopulation variation in the proportion of sticky eggs respectively allowed us to (i) artificially select and compare life-history traits of a virulent and a nonvirulent laboratory strain, and (ii) compare a virulent and a nonvirulent field strain, to investigate the hypothetical costs of virulence. We observed strong differences between the 2 laboratory strains. The nonvirulent strain invested fewer resources in reproduction and walked less than the virulent one but lived longer. Concerning the field strains, we observed that the nonvirulent strain had larger wings while the virulent one walked more and faster. All together, our results suggest that virulence may not always be costly, but rather that different life histories associated with different levels of virulence may coexist at both intra- and interpopulation levels.

PHENOTYPIC PLASTICITY AND REACTION NORMS IN TEMPERATE AND TROPICAL POPULATIONS OF DROSOPHILA MELANOGASTER: OVARIAN SIZE AND DEVELOPMENTAL TEMPERATURE

The plasticity of ovariole number relative to developmental temperature was studied in three populations of Drosophila melanogaster at both ends of the cline: a temperate French population and two equatorial Congolese. Ovary size was much greater in the French flies, in agreement with an already known latitudinal cline. Among isofemale lines, significant differences in genetic variability were observed between populations with a maximum variability at intermediate temperatures. Parameters of phenotypic variability (CV and FA) were not statistically different among lines or populations, but a significant increase at low temperature was demonstrated for both. The shapes of the response curves (i.e., the norm of reaction) were analyzed by adjusting the data to a quadratic equation. The parameters of the equation were highly variable among lines. On the other hand, the temperature for maximum value of ovarioles (TMV) was much less variable and exhibited only a slightly significant difference between temperate and tropical flies (22.2°C vs. 22.7°C). During its geographic extension toward colder places, D. melanogaster underwent a large, presumably adaptative, increase in ovariole number but very little change in the norm of reaction of that trait.

Natural Genetic Variation and Candidate Genes for Morphological Traits in Drosophila melanogaster

Body size is a complex character associated to several fitness related traits that vary within and between species as a consequence of environmental and genetic factors. Latitudinal and altitudinal clines for different morphological traits have been described in several species of Drosophila and previous work identified genomic regions associated with such variation in D. melanogaster. However, the genetic factors that orchestrate morphological variation have been barely studied. Here, our main objective was to investigate genetic variation for different morphological traits associated to the second chromosome in natural populations of D. melanogaster along latitudinal and altitudinal gradients in Argentina. Our results revealed weak clinal signals and a strong population effect on morphological variation. Moreover, most pairwise comparisons between populations were significant. Our study also showed important within-population genetic variation, which must be associated to the second chromosome, as the lines are otherwise genetically identical. Next, we examined the contribution of different candidate genes to natural variation for these traits. We performed quantitative complementation tests using a battery of lines bearing mutated alleles at candidate genes located in the second chromosome and six second chromosome substitution lines derived from natural populations which exhibited divergent phenotypes. Results of complementation tests revealed that natural variation at all candidate genes studied, invected, Fasciclin 3, toucan, Reticulon-like1, jing and CG14478, affects the studied characters, suggesting that they are Quantitative Trait Genes for morphological traits. Finally, the phenotypic patterns observed suggest that different alleles of each gene might contribute to natural variation for morphological traits. However, non-additive effects cannot be ruled out, as wild-derived strains differ at myriads of second chromosome loci that may interact epistatically with mutant alleles.

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.

Parallel Gene Expression Differences between Low and High Latitude Populations of Drosophila melanogaster and D. simulans

Gene expression variation within species is relatively common, however, the role of natural selection in the maintenance of this variation is poorly understood. Here we investigate low and high latitude populations of Drosophila melanogaster and its sister species, D. simulans, to determine whether the two species show similar patterns of population differentiation, consistent with a role for spatially varying selection in maintaining gene expression variation. We compared at two temperatures the whole male transcriptome of D. melanogaster and D. simulans sampled from Panama City (Panama) and Maine (USA). We observed a significant excess of genes exhibiting differential expression in both species, consistent with parallel adaptation to heterogeneous environments. Moreover, the majority of genes showing parallel expression differentiation showed the same direction of differential expression in the two species and the magnitudes of expression differences between high and low latitude populations were correlated across species, further bolstering the conclusion that parallelism for expression phenotypes results from spatially varying selection. However, the species also exhibited important differences in expression phenotypes. For example, the genomic extent of genotype × environment interaction was much more common in D. melanogaster. Highly differentiated SNPs between low and high latitudes were enriched in the 3’ UTRs and CDS of the geographically differently expressed genes in both species, consistent with an important role for cis-acting variants in driving local adaptation for expression-related phenotypes.

While gene expression variation in natural populations is common, the population genetic processes responsible for the maintenance of this variation remain obscure. Here we study geographic differences in gene expression in recently established low and high latitude populations of two closely related species of Drosophila. We observe substantial parallelism in expression differences and expression plasticity between populations, which supports the idea that spatially varying selection correlated with latitude contributes to the maintenance of gene expression variation in these species. Comparison of inter-population sequence differentiation and expression differentiation suggests that cis-acting variants play a role in geographic expression differentiation.

Effects of body-size variation on flight-related traits in latitudinal populations of Drosophila melanogaster

In the present study, we tested the hypothesis whether flight-related traits such as wing area, flight-muscle ratio, wing loading and dispersal yield evidence of geographical variation in nine wild-collected as well as laboratory-reared (at 21°C) latitudinal populations of Drosophila melanogaster from the Indian subcontinent. We observed positive clinal variation in the wing-thorax ratio, wing aspect ratio and wing area, along a latitudinal gradient for both the sexes. In contrast, geographical changes in three parameters of flight ability, i.e. flight-muscle ratio, wing loading and dispersal, showed negative correlation with latitude. On the basis of isofemale line variability, we observed positive correlation of wing loading with flight-muscle ratio as well as dispersal behaviour in both the sexes. We also found positive correlation between duration of development and wing area. Interestingly, southern populations of D. melanogaster from warm and humid habitats exhibited higher flight-muscle ratio as well as the higher wing loading than northern populations which occur in cooler and drier climatic conditions. Laboratory tests for dispersal-related walking behaviour showed significantly higher values for southern populations compared with northern populations of D. melanogaster. Multiple regression analysis of geographical changes in flight-muscle ratio, wing loading as well as walking behaviour as a function of average temperature and relative humidity of the origin of populations in wild-collected flies have suggested adaptive changes in flight-related traits in response to steeper gradients of climatic factors in the Indian subcontinent. Finally, adaptive latitudinal variations in flight-related traits in D. melanogaster are consistent with results of other studies from different continents despite differences due to specific climatic conditions in the Indian subcontinent.

Complexity of the genetic basis of ageing in nature revealed by a clinal study of lifespan and methuselah, a gene for ageing, in Drosophila from eastern Australia

Clinal studies are a powerful tool for understanding the genetic basis of climatic adaptation. However, while clines in quantitative traits and genetic polymorphisms have been observed within and across continents, few studies have attempted to demonstrate direct links between them. The gene methuselah in Drosophila has been shown to have a major effect on stress response and longevity phenotypes based largely on laboratory studies of induced mutations in the mth gene. Clinal patterns in the most common mth haplotype and for lifespan (both increasing with latitude) have been observed in North American populations of D. melanogaster, implicating climatic selection. While these clinal patterns have led some to suggest that mth influences ageing in natural populations, limited evidence on the association between the two has so far been collected. Here, we describe a significant cline in the mth haplotype in eastern Australian D. melanogaster populations that parallel the cline in North America. We also describe a cline in mth gene expression. These findings further support the idea that mth is itself under selection. In contrast, we show that lifespan has a strong nonlinear clinal pattern, increasing southwards from the tropics, but then decreasing again from mid-latitudes. Furthermore, in association studies, we find no evidence for a direct link between mth haplotype and lifespan. Thus, while our data support a role for mth variation being under natural selection, we found no link to naturally occurring variation in lifespan and ageing in Australian populations of D. melanogaster. Our results indicate that the mth locus likely has genetic background and environment-specific effects.

Altitudinal clinal variation in wing size and shape in African Drosophila melanogaster: one cline or many?

Geographical patterns of morphological variation have been useful in addressing hypotheses about environmental adaptation. In particular, latitudinal clines in phenotypes have been studied in a number of Drosophila species. Some environmental conditions along latitudinal clines-for example, temperature-also vary along altitudinal clines, but these have been studied infrequently and it remains unclear whether these environmental factors are similar enough for convergence or parallel evolution. Most clinal studies in Drosophila have dealt exclusively with univariate phenotypes, allowing for the detection of clinal relationships, but not for estimating the directions of covariation between them. We measured variation in wing shape and size in D. melanogaster derived from populations at varying altitudes and latitudes across sub-Saharan Africa. Geometric morphometrics allows us to compare shape changes associated with latitude and altitude, and manipulating rearing temperature allows us to quantify the extent to which thermal plasticity recapitulates clinal effects. Comparing effect vectors demonstrates that altitude, latitude, and temperature are only partly associated, and that the altitudinal shape effect may differ between Eastern and Western Africa. Our results suggest that selection responsible for these phenotypic clines may be more complex than just thermal adaptation.

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.

Molecular basis of adaptive shift in body size in Drosophila melanogaster: functional and sequence analyses of the Dca gene

Latitudinal body size clines in animals conforming to Bergmann's rule occur on many continents but isolating their underlying genetic basis remains a challenge. In Drosophila melanogaster, the gene Dca accounts for approximately 5-10% of the natural wing size variation (McKechnie SW, Blacket MJ, Song SV, Rako L, Carroll X, Johnson TK, Jensen LT, Lee SF, Wee CW, Hoffmann AA. 2010. A clinally varying promoter polymorphism associated with adaptive variation in wing size in Drosophila. Mol Ecol. 19:775-784). We present here functional evidence that Dca is a negative regulator of wing size. A significant negative latitudinal cline of Dca gene expression was detected in synchronized third instar larvae. In addition, we clarified the evolutionary history of the three most common Dca promoter alleles (Dca237-1, Dca237-2, and Dca247) and showed that the insertion allele (Dca247), whose frequency increases with latitude, is associated with larger wing centroid size and higher average cell number in male flies. Finally, we showed that the overall linkage disequilibrium (LD) was low in the Dca promoter and that the insertion/deletion polymorphism that defines the Dca alleles was in strong LD with two other upstream sites. Our results provide strong support that Dca is a candidate for climatic adaptation in D. melanogaster.

Heterozygosity maintains developmental stability of sternopleural bristles in Drosophila subobscura interpopulation hybrids

Interpopulation hybridization can lead to outbreeding depression within affected populations due to breakdown of coadapted gene complexes or heterosis in hybrid populations. One of the principal methods commonly used to estimate the level of developmental instability (DI) is fluctuating asymmetry (FA). We used three genetically differentiated Drosophila subobscura populations according to inversion polymorphism analysis and measured the variability of sternopleural bristle number and change in FA across generations P, F1, and F2 between intra- and interpopulation hybrids of D. subobscura. The mean variability of sternopleural bristle number in intra- and interpopulation hybrids of D. subobscura across generations cannot determine whether the changes at the level of developmental homeostasis are due exclusively to genomic coadaptation or to heterozygosity. Phenotypic variance (V(p)) and FA of sternopleural bristle number was higher in interpopulation than in intrapopulation hybrids across generations. F1 hybrids were more developmentally stable compared to each parental population in both intra- and interpopulation hybrids. The most probable mechanism providing developmental homeostasis is heterozygote or hybrid superiority, also called overdominace. However, V(p) was higher and FA lower in the F2 generation when compared to F1, due mainly to crossing-over in the formation of F2.

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.

Are latitudinal clines in body size adaptive?

Body size of animals often increases with increasing latitude. These latitudinal clines in body size have interested biologists for over 150 years. However, the mechanisms that generate these clines in size are still unclear, though latitudinal gradients in temperature appear to play an important role. More importantly, many studies that examine latitudinal clines in body size and the mechanisms responsible for these clines use phenotypic data, confounding genetic (adaptive) and non-genetic (plasticity) sources of variation. Yet, most of these studies make adaptive conclusions based on phenotypic measures of size. Here I show the dangers of making adaptive inferences from phenotypic measures of size. In addition, I use a specific form of plasticity in body size of ectotherms, called the temperature - size rule, to illustrate how confusion about genetic and non-genetic contributions to phenotypic variation has hampered progress in understanding the evolution of latitudinal clines in size. Field-based measurements of body size can no doubt be influenced by plasticity, but demonstrating that latitudinal clines have a genetic basis is necessary to show that these patterns are adaptive.

Wing morphometry as a tool for correct identification of primary and secondary New World screwworm fly

Cochliomyia hominivorax and Cochliomyia macellaria are endemic Neotropical Calliphoridae species. The former causes severe myiasis in hosts while the latter is Sarcosaprophagous, but commonly found as a second invader in wounds. Due to the morphological similarity between them and the potential losses that C. hominivorax represents for cattle breeders, the rapid and correct identification of these two species is very important. In addition to a correct identification of these species, a good knowledge of C. hominivorax biology can be helpful for designing control programs. We applied geometric morphometric methods to assess wing differences between C. hominivorax and C. macellaria and conduct a preliminary analysis of wing morphological variation in C. hominivorax populations. Canonical variate analysis, using wing shape data, correctly classified 100% of the individuals analyzed according to sex and species. This result demonstrates that wing morphometry is a simple and reliable method for identifying C. hominivorax and C. macellaria samples and can be used to monitor C. hominivorax. Both species show sexual dimorphism, but in C. hominivorax it is magnified. We suggest that this may reflect different histories of selection pressures operating on males and females. Significant differences in wing size and shape were obtained among C. hominivorax populations, with little correlation with latitude. This result suggests that wing variation is also a good morphological marker for studying population variation in C. hominivorax.

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.

Phenotypic variability of natural populations of an invasive drosophilid, Zaprionus indianus, on different continents: comparison of wild-living and laboratory-grown flies

Phenotypic variability in nature is the most important feature for Darwinian adaptation, yet it has been rarely investigated in invasive species. Zaprionus indianus is an Afrotropical drosophilid species that have recently invaded the Palearctic and the Neotropical regions. Here, we compared the variability of three size-related traits and one meristic trait the sternopleural (STP) bristle number, between wild-collected flies living under different conditions: a stressful Mediterranean environment in Egypt, and a benign tropical environment in Brazil. From each population, a F(1) generation was also grown under the stable conditions of the laboratory. Variability of size in nature had a variance 13 times greater than in the laboratory, but not affected by different climates. By contrast, STP variability was identical in nature and in the laboratory. Sexual dimorphism was also investigated with contrasting results between traits. It is suggested that the very high invasiveness of Z. indianus might be related to a better capacity to survive adverse conditions.

Impact of body melanisation on desiccation resistance in montane populations of D. melanogaster: Analysis of seasonal variation

In the montane localities of subtropical regions, winter is the dry season and ectothermic drosophilids are expected to evolve desiccation resistance to cope with drier climatic conditions. An analysis of six montane populations (600-2226m) of D. melanogaster showed variations for body melanisation (i.e. pigmentation) and desiccation resistance across seasons as well as along altitude. During winter season, plastic changes for melanisation of three posterior abdominal segments (5th+6th+7th) correspond with higher desiccation resistance. Thus, we analyzed genetic and plastic effects for these ecophysiological traits by comparing wild-caught and laboratory reared individuals of D. melanogaster for autumn as well as winter season. A ratio of slope values in wild vs. laboratory populations has shown a 1.64-fold plastic effect during autumn; and a two-fold effect during winter. For body melanisation and desiccation resistance, evolutionary response to altitudinal environmental gradient is similar to the phenotypic response across seasons. Thus, our observations are in agreement with the co-gradient hypothesis. Further, we tested the hypothesis whether a thicker cuticle (either due to melanisation or cuticular lipids) leads to lesser cuticular water loss and higher desiccation resistance across seasons as well as according to altitude. Based on within and between population analyses, body melanisation was found to be positively correlated with desiccation resistance but negatively with cuticular water loss. Interestingly, there were no changes in the amount of cuticular lipids per fly across seasons as well as along altitude; and therefore cuticular lipids did not account for desiccation resistance. Cuticular water loss exhibited negative correlation with body melanisation but not with cuticular lipids as well as with changes in body size across seasons. Thus, our data suggest that seasonal changes in body melanisation confer desiccation resistance in montane populations of D. melanogaster.

Selection does not favor larger body size at lower temperature in a seed-feeding beetle

Body size of many animals increases with increasing latitude, a phenomenon known as Bergmann's rule (Bergmann clines). Latitudinal gradients in mean temperature are frequently assumed to be the underlying cause of this pattern because temperature covaries systematically with latitude, but whether and how temperature mediates selection on body size is unclear. To test the hypothesis that the "relative" advantage of being larger is greatest at cooler temperatures we compare the fitness of replicate lines of the seed beetle, Stator limbatus, for which body size was manipulated via artificial selection ("Large,""Control," and "Small" lines), when raised at low (22 degrees C) and high (34 degrees C) temperatures. Large-bodied beetles (Large lines) took the longest to develop but had the highest lifetime fecundity, and highest fitness (r(C)), at both low and high temperatures. However, the relative difference between the Large and Small lines did not change with temperature (replicate 2) or was greatest at high temperature (replicate 1), contrary to the prediction that the fitness advantage of being large relative to being small will decline with increasing temperature. Our results are consistent with two previous studies of this seed beetle, but inconsistent with prior studies that suggest that temperature-mediated selection on body size is a major contributor to the production of Bergmann clines. We conclude that other environmental and ecological variables that covary with latitude are more likely to produce the gradient in natural selection responsible for generating Bergmann clines.

Patterns of sequence variability and divergence at the diminutive gene region of Drosophila melanogaster: complex patterns suggest an ancestral selective sweep

To identify putatively swept regions of the Drosophila melanogaster genome, we performed a microsatellite screen spanning a 260-kb region of the X chromosome in populations from Zimbabwe, Ecuador, the United States, and China. Among the regions identified by this screen as showing a complex pattern of reduced heterozygosity and a skewed frequency spectrum was the gene diminutive (dm). To investigate the microsatellite findings, nucleotide sequence polymorphism data were generated in populations from both China and Zimbabwe spanning a 25-kb region and encompassing dm. Analysis of the sequence data reveals strongly reduced nucleotide variation across the entire gene region in both the non-African and the African populations, an extended haplotype pattern, and structured linkage disequilibrium, as well as a rejection of neutrality in favor of selection using a composite likelihood-ratio test. Additionally, unusual patterns of synonymous site evolution were observed at the second exon of this locus. On the basis of simulation studies as well as recently proposed methods for distinguishing between selection and nonequilibrium demography, we find that this "footprint" is best explained by a selective sweep in the ancestral population, the signal of which has been somewhat blurred via founder effects in the non-African samples.

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.

Geographic variation in body size and sexual size dimorphism of a seed-feeding beetle

Body size of many animals varies with latitude: body size is either larger at higher latitudes (Bergmann's rule) or smaller at higher latitudes (converse Bergmann's rule). However, the causes underlying these patterns are poorly understood. Also, studies rarely explore how sexual size dimorphism varies with latitude. Here we investigate geographic variation in body size and sexual size dimorphism of the seed-feeding beetle Stator limbatus, collected from 95 locations along a 38 degrees range in latitude. We examine 14 variables to test whether clines in environmental factors are adequate to explain geographic patterns of body size. We found that body size and sexual size dimorphism of S. limbatus varied considerably with latitude; beetles were smaller but more dimorphic at lower latitudes. Body size was not correlated with a gradient in mean temperature, contrary to the commonly accepted hypothesis that clines are produced by latitudinal gradients in temperature. Instead, we found that three factors were adequate to explain the cline in body size: clinal variation in host plant seed size, moisture (humidity), and seasonality (variance in humidity, precipitation, and temperature). We also found that the cline in sexual size dimorphism was partially explainable by a gradient in moisture, though moisture alone was not sufficient to explain the cline. Other ecological or environmental variables must necessarily contribute to differences in selection on male versus female body size. The main implications of our study are that the sexes differ in the magnitude of clinal variation in body size, creating latitudinal variation in sexual size dimorphism, and that clines in body size of seed beetles are likely influenced by variation in host seed size, water availability, and seasonality.

Developmental time and size-related traits in Drosophila buzzatii along an altitudinal gradient from Argentina

Clinal analysis for fitness-related traits provides a well-known approach to investigate adaptive evolution. Several fitness-related traits (developmental time, thorax length, wing length and wing loading) were measured at two laboratory generations (G7 and G33) of D. buzzatii from an altitudinal gradient from northwestern Argentina, where significant thermal differences persist. Developmental time (DT) was positively correlated with altitude of origin of population. Further, DT was negatively correlated with maximal mean temperature at the site of origin of population, and this thermal variable decreases with altitude. Wing loading tended to be larger in highland than in lowland populations, suggesting that flight performance is subject to stronger selection pressure in highland populations. Developmental time showed a significant increase with laboratory generation number. There was no significant correlation between developmental time and body size across populations along the altitudinal cline of DT. This result illustrates that developmental time and body size do not always evolve in the same direction, even though both traits are often positively and genetically correlated in a well-known tradeoff in Drosophila.

Phenotypic plasticity of body size in a temperate population ofDrosophila melanogaster: When the temperature—size rule does not apply

A natural population of Drosophila melanogaster in southern France was sampled in three different years and 10 isofemale lines were investigated from each sample. Two size-related traits, wing and thorax length, were measured and the wing/thorax ratio was also calculated. Phenotypic plasticity was analysed after development at seven different constant temperatures, ranging from 12 degrees C to 31 degrees C. The three year samples exhibited similar reaction norms, suggesting a stable genetic architecture in the natural population. The whole sample (30 lines) was used to determine precisely the shape of each reaction norm, using a derivative analysis. The practical conclusion was that polynomial adjustments could be used in all cases, but with different degrees: linear for the wing/thorax ratio, quadratic for thorax length, and cubic for wing length. Both wing and thorax length exhibited concave reaction norms, with a maximum within the viable thermal range. The temperatures of the maxima were, however, quite different, around 15 degrees C for the wing and 19.5 degrees C for the thorax. Assuming that thorax length is a better estimate of body size, it is not possible to state that increasing the temperature results in monotonically decreasing size (the temperature-size rule), although this is often seen to be the case for genetic variations in latitudinal clines. The variability of the traits was investigated at two levels-within and between lines-and expressed as a coefficient of variation. The within-line (environmental) variability revealed a regular, quadratic convex reaction norm for the three traits, with a minimum around 21 degrees C. This temperature of minimum variability may be considered as a physiological optimum, while extreme temperatures are stressful. The between-line (genetic) variability could also be adjusted to quadratic polynomials, but the curvature parameters were not significant. Our results show that the mean values of the traits and their variance are both plastic, but react in different ways along a temperature gradient. Extreme low or high temperatures decrease the size but increase the variability. These effects may be considered as a functional response to environmental stress.

Isofemale lines in Drosophila: an empirical approach to quantitative trait analysis in natural populations

Founding isofemale lines from wild collected females is a basic tool for investigating the genetic architecture of Drosophila natural populations. The method permits the analysis of quantitative traits under laboratory conditions, with a much broader scope than the mere evidence of a significant genetic heterogeneity among lines. Genetic variability is generally demonstrated by a significant coefficient of intraclass correlation, but several experimental precautions are needed and explained here. The relationship between classical (additive) heritability and intraclass correlation is not straightforward, presumably because the genetic bottlenecks due to the initiation of the lines unravel a significant, nonadditive genetic variance due to dominance and epistatic effects. It is thus suggested to consider intraclass correlation as a specific genetic parameter that enables comparisons between different traits, different populations or different environments. The use of isofemale lines is, however, not restricted to the calculation of an intraclass correlation. It can be used to estimate genetic correlations among traits or environments. The method is also convenient for the analysis of phenotypic plasticity in relation to an environmental gradient. A precise description of the response curves (the reaction norms) is possible, distinguishing trait parameters and plasticity parameters. A fairly general conclusion is that, for a given trait, mean value and plasticity are genetically independent. It is also possible to analyze traits, which, like sexual dimorphism, must be measured on different individuals, and even to demonstrate their genetic variability. In many cases, further empirical and theoretical analyses are possible and needed. It is argued that, in the future, isofemale lines will have an increasing significance among the various techniques appropriate to the analysis of quantitative evolutionary genetics in a diversity of species.

Microclinal variation for ovariole number and body size in Drosophila melanogaster in ?Evolution Canyon?

Sites that display strong environmental contrasts in close proximity, such as 'Evolution Canyon' on Mt. Carmel, Israel, are natural theatres for investigating adaptive evolution in action. We reared Drosophila melanogaster from collection sites along altitudinal transects on the north- and south-facing canyon slopes in each of three temperature environments, and assessed genetic variation in ovariole number and body size between and within collection sites, and temperature plasticity. Both traits exhibited significant genetic variation within collection sites and phenotypic plasticity in response to temperature, but not genetic variation for plasticity. Between-site genetic variation in ovariole number was negatively correlated with altitude on both slopes of the canyon, and collections from the north- and south-facing slopes were genetically differentiated for male, but not female, body size. Genetic variation between sites within easy dispersal range is consistent with the action of strong natural selection, although neither the selective agent(s) nor the direct targets of selection are known.