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

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

Evolutionary pace of chromosomal polymorphism in colonizing populations of Drosophila subobscura: an evolutionary time series

Abstract. Biologists have long debated the speed, uniformity, and predictability of evolutionary change. However, evaluating such patterns on a geographic scale requires time-series data on replicate sets of natural populations. Drosophila subobscura has proven an ideal model system for such studies. This fly is broadly distributed in the Old World, but was introduced into both North and South America just over two decades ago and then spread rapidly. Rapid, uniform, and predictable evolution would be demonstrated if the invading flies evolved latitudinal clines that progressively converged on those of the native populations. Evolutionary geneticists quickly capitalized on this opportunity to monitor evolutionary dynamics. Just a few years after the introduction, they surveyed chromosomal inversion frequencies in both North and South America. On both continents they detected incipient latitudinal clines in chromosome inversion frequencies that almost always had the same sign with latitude as in the Old World. Thus the initial evolution of chromosomal polymorphisms on a continental scale was remarkably rapid and consistent. Here we report newer samples of inversion frequencies for the colonizing populations: the time series now spans almost one decade for North America and almost two decades for South America. Almost all inversions in the New World continue to show the same sign of frequency with latitude as in the Old World. Nevertheless, inversion clines have not consistently increased in steepness over time; nor have they consistently continued to converge on the Old World baseline. However, five arrangements in South America show directional, continentwide shifts in frequency. Overall, the initial consistency of clinal evolutionary trajectories seen in the first surveys seems not to have been maintained.     

Genetic loads and coadaptivity of chromosomal inversions I. Inversion polymorphism and genetic load on chromosome O in a Drosophila subobscura population from Petnica

Inversion polymorphism on chromosome O and polymorphism for the viability of determining genes have been studied in a natural population of Drosophila subobscura from Petnica (Serbia). The range of inversion polymorphism and the abundance of particular gene arrangements in the study population agree with a general pattern of inversion polymorphism of D. subobscura in Europe. The data obtained on the amount of genetic loads show that the D. subobscura population from Petnica displays a moderate degree of that polymorphism, compared to the other studied populations of these species. Therefore, the D. subobscura population from Petnica could be tentatively classified as an ecologically central population. Examination association of chromosomal, thus, inversion polymorphism with gene polymorphism, in the form of genetic loads show that differences exist in the mean viability among certain gene arrangements. The distribution of deleterious genes among chromosome O gene arrangements were non-random.     

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

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

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

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

Nucleotide Polymorphism at the rp49 region of Drosophila subobscura: Lack of geographic subdivision within chromosomal arrangements in Europe

Restriction-site polymorphism at the rp49 gene region has been studied in 234 lines of Drosophila subobscura representing different gene arrangements for the O chromosome. The population surveyed (El Pedroso, Spain) was sampled four times in each of two consecutive years. The data indicate that the two chromosomal classes studied, O_[ST] and O_{[3 + 4]}, are genetically differentiated in El Pedroso. Comparison of the present results with those previously obtained for two other populations further supports that, for a given chromosomal class, European populations are not genetically differentiated. This lack of differentiation at the rp49 region within O_[ST] and within O_{[3 + 4]} stands in contrast to the clear latitudinal clines found in Europe for these arrangements.     

Genetic variation in the spread of Drosophila subobscura from a nonequilibrium population

Abstract.— Drosophila subobscura was first identified in North America in the early 1980s, and a newer D. subobscura population in Utah appears to have been established more than 10 years later. In this study, we use nuclear microsatellite allele frequencies, mitochondrial restriction fragment length polymorphism (RFLP) allele frequencies, and computer simulations to investigate possible scenarios of how this species has spread across North America. Our method develops a 95% confidence interval for the maximum and minimum number of founders that could have colonized the new population given various scenarios for spread. Unlike many other methods, it may be applied to nonequilibrium source populations given certain conditions. We find that observed allele frequency differences between newer and older D. subobscura populations are consistent with very few inseminated females being transported east from the coast in a single step or with larger numbers of colonizers invading after several intermediate steps. They are not consistent with a large, panmictic population of D. subobscura colonizing Utah in a single step.     

Relationship between chromosomal polymorphism and wing size in a natural population of Drosophila subobscura

Chromosomal polymorphism and wing size (as a measure of body size) were analysed simultaneously in two samples of Drosophila subobscura from Barcelona, Spain. The very rich chromosomal polymorphism of this species makes it difficult to detect the relationship of this polymorphism with any phenotypical character. However, a positive significant regression of wing size on the percentage of the autosome length with standard arrangement was found. Furthermore, for each polymorphic chromosome, except for the J chromosome, an association between the most frequent arrangements and wing size was observed. This trend, which was the same in the two samples, was that expected according to the latitudinal clines of both characters.     

Dissecting adaptive clinal variation: markers, inversions and size/stress associations in Drosophila melanogaster from a central field population

Many organisms show latitudinal variation for quantitative traits that is assumed to be due to climatic adaptation. These clines provide an opportunity to study the genetics of the adaptive process both at the phenotypic and the underlying molecular levels. Yet researchers rarely try to link variation in quantitative traits to their underlying molecular genetic basis. We describe a novel approach for exploring the genetic basis for clinal variation in size and stress traits in Drosophila melanogaster. We look for associations between genetic markers and traits that exhibit clinal patterns on the east coast of Australia using a single, geographically central population. There are strong associations between markers found within In(3R)Payne and variation in size, suggesting that this inversion explains much of the clinal variation in this trait. We also find that development time is associated with the Adh allozyme locus, cold resistance is negatively associated with the In(3L)Payne inversion and a genetic marker for Hsp70, a heat‐shock protein, is associated with heat resistance. Finally we discuss the importance of inversions in clinal variation for quantitative traits and for identifying quantitative trait loci.     

Tracking the origin of the American colonization by Drosophila subobscura: genetic comparison between Eastern and Western Mediterranean populations

Several pieces of evidence indicate a Mediterranean origin of the colonization of America by Drosophila subobscura . To ascertain whether the origin was from the Eastern or the Western Mediterranean region, samples from Barcelona (Spain) and Mt Parnes (Greece) were collected and O chromosomal inversion polymorphism and lethal genes were analysed. The frequencies of lethal chromosomes were 0.244 ± 0.039 in Barcelona and 0.336 ± 0.043 in Mt Parnes, consistent with the expectations for large populations located in the central area of the species distribution. Lethal genes seem to be distributed at random along the O chromosome in both populations. The intra-populational allelism frequencies of Barcelona and Mt Parnes were 0.016 ± 0.007 and 0.012 ± 0.005 respectively. Thus, the estimates of the effective population size were high in both populations (between 6964 and 13 004 in Barcelona and 11 874 to 26 828 in Mt Parnes). The cases of allelism in Mt Parnes were observed only between individual lethal genes, but in Barcelona some concatenated clusters of allelism were detected. This pattern of allelism can be explained by synthetic lethality, hybrid dysgenesis, the induction of recurrent lethal mutations by different factors or an effect of microdifferentiation in subpopulations. In both populations, a reduction in fitness in the heterozygotes for lethal genes has been detected. Furthermore, the estimates of the migration coefficient (between 0.0085 and 0.0120 in Barcelona, and 0.0057 and 0.0087 in Mt Parnes) confirm the existence of gene flow between Palearctic populations of D. subobscura . Our lethal genes and chromosomal inversion results are consistent with a Mediterranean origin of the colonization of America by D. subobscura , but are inconclusive with regard to the identification of the population from which these colonizers came.     

Chromosomal inversions promote genomic islands of concerted evolution of Hsp70 genes in the Drosophila subobscura species subgroup

Heat‐shock (HS) assays to understand the connection between standing inversion variation and evolutionary response to climate change in Drosophila subobscura found that “warm‐climate” inversion O₃₊₄ exhibits non‐HS levels of Hsp70 protein like those of “cold‐climate” O_ST after HS induction. This was unexpected, as overexpression of Hsp70 can incur multiple fitness costs. To understand the genetic basis of this finding, we have determined the genomic sequence organization of the Hsp70 family in four different inversions, including O_ST, O₃₊₄ , O₃₊₄₊₈ and O₃₊₄₊₁₆ , using as outgroups the remainder of the subobscura species subgroup, namely Drosophila madeirensis and Drosophila guanche. We found (i) in all the assayed lines, the Hsp70 family resides in cytological locus 94A and consists of only two genes, each with four HS elements (HSEs) and three GAGA sites on its promoter. Yet, in O_ST, the family is comparatively more compact; (ii) the two Hsp70 copies evolve in concert through gene conversion, except in D. guanche; (iii) within D. subobscura, the rate of concerted evolution is strongly structured by inversion, being higher in O_ST than in O₃₊₄ ; and (iv) in D. guanche, the two copies accumulated multiple differences, including a newly evolved “gap‐type” HSE2. The absence of concerted evolution in this species may be related to a long‐gone‐unnoticed observation that it lacks Hsp70 HS response, perhaps because it has evolved within a narrow thermal range in an oceanic island. Our results point to a previously unrealized link between inversions and concerted evolution, with potentially major implications for understanding genome evolution.     

Long-term changes in the chromosomal inversion polymorphism of Drosophila subobscura. I. Mediterranean populations from southwestern Europe

The chromosomal polymorphism of seven Mediterranean populations of Drosophila subobscura has been compared with that of the same populations collected 26 to 35 years ago. Significant latitudinal clines for the frequencies of A(ST), E(ST), O(ST). and U(ST) chromosomal arrangements have been detected in the old and new samples. Standard gene arrangements are frequent in the north and decrease in frequency towards the south. Significant negative regression coefficients between latitude and transformed frequency have also been observed for the more frequent nonstandard gene arrangements. The pattern of the clines is practically the same in the old and new collections. Furthermore, the frequencies of gene arrangements of all chromosomes have changed significantly during this period in a systematic way: an increase in the frequency of those arrangements typical of southern latitudes and a decrease for those more common in northern latitudes is observed in all populations. These changes could be due to climatic factors that are correlated with latitude, making the chromosomal composition of this species more "southern.'     

Resistance to environmental stress in Drosophila ananassae: latitudinal variation and adaptation among populations

Geographical variation in traits related to fitness is often the result of adaptive evolution. Stress resistance traits in Drosophila often show clinal variation, suggesting that selection affects resistance traits either directly or indirectly. Multiple stress resistance traits were investigated in 45 natural populations of Drosophila ananassae collected from all over India. There was significant positive correlation between starvation resistance and lipid content. Significant negative correlations between desiccation and lipid content and between desiccation and heat resistance were also found. Flies from lower latitudes had higher starvation resistance, heat resistance and lipid content but the pattern was reversed for desiccation resistance. These results suggest that flies from different localities varied in their susceptibility to starvation because of difference in their propensity to store body lipid. Multiple regression analysis provided evidence of climatic selection driven by latitudinal variation in the seasonal amplitude of temperature and humidity changes within the Indian. Finally, our results suggest a high degree of variation in stress resistance at the population level in D. ananassae.     

Long‐term changes in the chromosomal inversion polymorphism of Drosophila subobscura. II. European populations*

The chromosomal polymorphism of 13 European populations of Drosophila subobscura has been compared with that of the same populations collected 15–35 years ago. The chromosomal polymorphism of the old populations differs significantly from that of the new populations, mainly for chromosomes U and O. There is a very good agreement between the geographical space and the genetic space as shown by a graphical representation of the 26 statistical populations (13 old and 13 new) obtained by a principal coordinate analysis. This reflects both the existence of significant latitudinal clines for the frequencies of some chromosomal arrangements in the old and new samples and systematic changes that have taken place in these populations during the period that elapsed between the two surveys. An increase in the frequency of those arrangements typical of southern latitudes and a decrease for those common in northern latitudes is observed in all populations – Mediterranean, Atlantic and Central European. Furthermore, the genetic distances of the new populations to a southern population of reference have decreased in comparison with those of the old populations. These changes could be the result of climatic factors that are correlated with latitude. In particular, the assumption that global warming is responsible for all the changes observed appears rather likely. Whether these systematic changes of the chromosomal polymorphism are a consequence of local adaptations or have been produced by migration from the south remains an open question.     

Rapid micro-evolution and loss of chromosomal diversity in Drosophila in response to climate warming

Concern regarding the ecological impact of rapid global warming has encouraged research on climate-induced changes in biological systems. Critical problems, still poorly understood, are the potential for rapid adaptive responses and their genetic costs to populations. The O chromosomal polymorphisms of Drosophila subobscura have been monitored at a southern Palearctic locality experiencing sustained climate warming since the mid-1970s. Observations suggest that the population is rapidly evolving in response to the new environmental conditions, and has lost a significant amount of chromosomal diversity (18.3% in 16 years). These findings are consistent with results from another population of D. subobscura, which is also undergoing climate warming, and are in accord with what would be expected from latitudinal and seasonal patterns of the various inversions. In addition, data on the O chromosomal polymorphisms from other localities throughout t he range of this species suggest that other populations vary similarly.     

Sex‐ratio in natural and experimental populations of Drosophila subobscura from North America

The sex‐ratio (SR), defined as the proportion of males, has been studied in three North American colonizing populations of Drosophila subobscura (Eureka, Davis and Gilroy). The proportion of sexes under laboratory conditions was studied using the one‐generation serial transfer technique in one‐ and two‐species populations, to infer whether biased SR affects the outcome when competing with Drosophila pseudoobscura, another member of the same group now in sympatry with D. subobscura in North America. The wild samples of D. subobscura yielded a significantly higher number of males than females during those months where the species is more abundant. However, there was no significant deviation in the 1 : 1 proportion of sexes in the descendants of D. subobscura at any of the experimental conditions. On the contrary, D. pseudoobscura produced a higher proportion of females which could be responsible for the exclusion of D. subobscura in laboratory competition experiments with overlapping generations. Thus, if sexes are equal at birth and survival is similar, the preponderance of males of D. subobscura in our wild collections could indicate greater activity and probably greater chance of dispersal of males versus females especially under favourable conditions.     

Lethal genes in O5, chromosomes of Drosophila subobscura from Europe and America

In populations of D. subobscura , a species that is know for its high chromosomal polymorphism, the O₅ inversion has a rather erratic frequency distribution in the Palearctic region. An O₅ lethal chromosomal line obtained from a Balkan population near Zanjic (South Adriatic, Montenegro, Yugoslavia) was tested for lethal allelism with other O₅ lethal chromosomal lines derived from American (USA and Chile) colonizing populations, and from the French population of Taulé. No allelism was found between the Balkan lethal gene and those from America and France. Thus, the lethal genes of the O₅ inversions are not of the same origin and it is most probable that the American colonizations did not start from the Zanjic population. The general difference in the chromosomal inversion polymorphism corroborates this conclusion. The cytological analysis confirms the assumption that all O₅ chromosomes studied are identical with respect to breakage points.     

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

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

Contrasting patterns of natural variation in global Drosophila melanogaster populations

Despite the popularity of Drosophila melanogaster in functional and evolutionary genetics, the global pattern of natural variation has not yet been comprehensively described in this species. For the first time, we report a combined survey using neutral microsatellites and mitochondrial sequence variation jointly. Thirty-five populations originating from five continents were compared. In agreement with previous microsatellite studies, sub-Saharan African populations were the most variable ones. Consistent with previous reports of a single 'out of Africa' habitat expansion, we found that non-African populations contained a subset of the African alleles. The pattern of variation detected for the mitochondrial sequences differed substantially. The most divergent haplotypes were detected in the Mediterranean region while Africa harboured most haplotypes, which were all closely related. In the light of the well-established African origin of D. melanogaster, our results cast severe doubts about the suitability of mtDNA for biogeographic inference in this model organism.     

Evolution of total net fitness in thermal lines: Drosophila subobscura likes it 'warm'

Fisher's fundamental theorem states that heritable variation for net fitness sets a limit to the rate of response to natural selection. How will temperate (i.e. cold‐tolerant) species cope with contemporary rapid global warming? Using three‐fold replicated lines of Drosophila subobscura that had been allowed to evolve for 4 years (between 32 and 59 generations) at 13 °C (cold), 18 °C (the supposed optimum temperature), and 22 °C (warm) I assess here how net fitness changes according to thermal environments. Net fitness was estimated following the classical approach in population genetics of competing over a number of generation in outbred experimental populations multiple wild‐type O chromosomes (homologous to arm 3R in D. melanogaster) independently derived from each base thermal stock in an otherwise homogeneous genetic background against a balancer chromosome. Warm‐adapted populations (‘warm‐adapted O chromosomes’) performed comparatively well at all tested temperatures. However, net fitness was severely reduced in cold‐adapted populations when transferred to warmer conditions. It seems, therefore, that thermal fitness breath for D. subobscura flies is positively associated to temperature. These findings are discussed in relation to the fast world‐wide clinal shifts in the frequency of genetic markers correlated with current climate change.