Population History Shapes Responses to Different Temperature Regimes in Drosophila subobscura.

Drosophila subobscura is considered a good model species for investigation of a population's ability to adapt and cope with climate changes. Decade long research has shown that inversion frequencies change in response to environmental factors indicating their role in adaptation to novel environments. The mechanisms behind organisms' responses to temperature are complex, involving changes in physiology, behavior, gene expression and regulation. On the other hand, a population's ability to respond to suboptimal conditions depends on standing genetic variation and population history. In order to elucidate the role of local adaptation in population response to the changing temperature, we investigated the response to temperature in D. subobscura individuals originating from two different altitudes by combining traditional cytogenetic techniques with assessing the levels of Hsp70 protein expression. Inversion polymorphism was assessed in the flies sampled from natural populations and in flies reared in laboratory conditions at three different temperatures after five and sixteen generations and Hsp70 protein expression profile in 12th generation flies at the basal level and after heat shock induction. Our results indicate that local adaptation and population history influence population response to the changing temperature.

Life History Traits in Two Drosophila Species Differently Affected by Microbiota Diversity under Lead Exposure.

Life history traits determine the persistence and reproduction of each species. Factors that can affect life history traits are numerous and can be of different origin. We investigated the influence of population origin and heavy metal exposure on microbiota diversity and two life history traits, egg-to-adult viability and developmental time, in Drosophila melanogaster and Drosophila subobscura, grown in the laboratory on a lead (II) acetate-saturated substrate. We used 24 samples, 8 larval and 16 adult samples (two species × two substrates × two populations × two sexes). The composition of microbiota was determined by sequencing (NGS) of the V3-V4 variable regions of the 16S rRNA gene. The population origin showed a significant influence on life history traits, though each trait in the two species was affected differentially. Reduced viability in D. melanogaster could be a cost of fast development, decrease in Lactobacillus abundance and the presence of Wolbachia. The heavy metal exposure in D. subobscura caused shifts in developmental time but maintained the egg-to-adult viability at a similar level. Microbiota diversity indicated that the Komagataeibacter could be a valuable member of D. subobscura microbiota in overcoming the environmental stress. Research on the impact of microbiota on the adaptive response to heavy metals and consequently the potential tradeoffs among different life history traits is of great importance in evolutionary research.

Basal hsp70 expression levels do not explain adaptive variation of the warm- and cold-climate O3 + 4 + 7 and OST gene arrangements of Drosophila subobscura.

BACKGROUND: Drosophila subobscura exhibits a rich inversion polymorphism, with some adaptive inversions showing repeatable spatiotemporal patterns in frequencies related to temperature. Previous studies reported increased basal HSP70 protein levels in homokaryotypic strains for a warm-climate arrangement compared to a cold-climate one. These findings do not match the similar hsp70 genomic organization between arrangements, where gene expression levels are expected to be similar. In order to test this hypothesis and understand the molecular basis for hsp70 expression, we compared basal hsp70 mRNA levels in males and females, and analysed the 5' and 3' regulatory regions of hsp70 genes in warm- and cold-climate isochromosomal O3 + 4 + 7 and OST lines of D. subobscura. RESULTS: We observed comparable mRNA levels between the two arrangements and a sex-biased hsp70 gene expression. The number of heat-shock elements (HSEs) and GAGA sites on the promoters were identical amongst the OST and O3 + 4 + 7 lines analysed. This is also true for 3' AU-rich elements where most A and B copies of hsp70 have, respectively, two and one element in both arrangements. Beyond the regulatory elements, the only notable difference between both arrangements is the presence in 3' UTR of a 14 bp additional fragment after the stop codon in the hsp70A copy in five O3 + 4 + 7 lines, which was not found in any of the six OST lines. CONCLUSIONS: The equivalent hsp70 mRNA amounts in OST and O3 + 4 + 7 arrangements provide the first evidence of a parallelism between gene expression and genetic organization in D. subobscura lines having these arrangements. This is reinforced by the lack of important differential features in the number and structure of regulatory elements between both arrangements, despite the genetic differentiation observed when the complete 5' and 3' regulatory regions were considered. Therefore, the basal levels of hsp70 mRNA cannot account, in principle, for the adaptive variation of the two arrangements studied. Consequently, further studies are necessary to understand the intricate molecular mechanisms of hsp70 gene regulation in D. subobscura.

Expression of thermal tolerance genes in two Drosophila species with different acclimation capacities.

Heat tolerance increases at higher acclimation temperatures in D. melanogaster, but not in D. subobscura. The two species represent separate lineages of the subgenus Sophophora of Drosophila with contrasting tropical African and temperate Palearctic evolutionary histories. D. melanogaster has five copies of the inducible hsp70 gene distributed in two clusters, named A (with two copies) and B (three copies), while D. subobscura has only two copies arranged similarly to cluster A of D. melanogaster. The hsp70s of the two species also differ in their cis-regulatory regions, with D. melanogaster exhibiting features of a faster and more productive promoter. We predicted that the interspecific variation in acclimation capacity of heat tolerance is explained by evolved variation in expression of the major group of heat shock proteins. To test this prediction, we compared basal levels of gene expression at different developmental temperatures within each of the two species. Furthermore, we explored the heat hardening dynamics by measuring the induction of gene expression during a ramping assay. The prediction of a stronger heat shock protein response in D. melanogaster as compared to D. subobscura was confirmed for both long-term acclimation and short-term hardening. For D. melanogaster the upregulation with temperature ramping ranged from less than two fold (hsp26) to 2500 fold (hsp70A) increase. In all cases induction in D. melanogaster exceeded that of D. subobscura homologs. These differences correlate with structural differences in the regulatory regions of hsp70, and might explain differences in acclimation capacity among species. Finally, in D. melanogaster we found an indication of an inverse relationship between basal and induced levels of hsp70A and hsp83 expression, suggesting a divergent role for thermal adaptation of these genes at benign and stressful temperatures, respectively.

How phenotypic convergence arises in experimental evolution.

Evolutionary convergence is a core issue in the study of adaptive evolution, as well as a highly debated topic at present. Few studies have analyzed this issue using a "real-time" or evolutionary trajectory approach. Do populations that are initially differentiated converge to a similar adaptive state when experiencing a common novel environment? Drosophila subobscura populations founded from different locations and years showed initial differences and variation in evolutionary rates in several traits during short-term (∼20 generations) laboratory adaptation. Here, we extend that analysis to 40 more generations to analyze (1) how differences in evolutionary dynamics among populations change between shorter and longer time spans, and (2) whether evolutionary convergence occurs after 60 generations of evolution in a common environment. We found substantial variation in longer term evolutionary trajectories and differences between short- and longer term evolutionary dynamics. Although we observed pervasive patterns of convergence toward the character values of long-established populations, populations still remain differentiated for several traits at the final generations analyzed. This pattern might involve transient divergence, as we report in some cases, indicating that more generations should lead to final convergence. These findings highlight the importance of longer term studies for understanding convergent evolution.

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 O3+4 exhibits non-HS levels of Hsp70 protein like those of "cold-climate" OST 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 OST , O3+4 , O3+4+8 and O3+4+16 , 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 OST, 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 OST than in O3+4 ; 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.

Strong hybrid male incompatibilities impede the spread of a selfish chromosome between populations of a fly.

Meiotically driving sex chromosomes manipulate gametogenesis to increase their transmission at a cost to the rest of the genome. The intragenomic conflicts they produce have major impacts on the ecology and evolution of their host species. However, their ecological dynamics remain poorly understood. Simple population genetic models predict meiotic drivers will rapidly reach fixation in populations and spread across landscapes. In contrast, natural populations commonly show spatial variation in the frequency of drivers, with drive present in clines or mosaics across species ranges. For example, Drosophila subobscura harbors a sex ratio distorting drive chromosome (SRs) at 15-25% frequency in North Africa, present at less than 2% frequency in adjacent southern Spain, and absent in other European populations. Here, we investigate the forces preventing the spread of the driver northward. We show that SRs has remained at a constant frequency in North Africa, and failed to spread in Spain. We find strong evidence that spread is impeded by genetic incompatibility between SRs and Spanish autosomal backgrounds. When we cross SRs from North Africa onto Spanish genetic backgrounds we observe strong incompatibilities specific to hybrids bearing SRs. The incompatibilities increase in severity in F2 male hybrids, leading to almost complete infertility. We find no evidence supporting an alternative hypothesis, that there is resistance to drive in Spanish populations. We conclude that the source of the stepped frequency variation is genetic incompatibility between the SRs chromosome and the genetic backgrounds of the adjacent population, preventing SRs spreading northward. The low frequency of SRs in South Spain is consistent with recurrent gene flow across the Strait of Gibraltar combined with selection against the SRs element through genetic incompatibility. This demonstrates that incompatibilities between drive chromosomes and naïve populations can prevent the spread of drive between populations, at a continental scale.

Patterns of geographic variation of thermal adapted candidate genes in Drosophila subobscura sex chromosome arrangements.

BACKGROUND: The role of chromosomal arrangements in adaptation is supported by the repeatable clinal variation in inversion frequencies across continents in colonizing species such as Drosophila subobscura. However, there is a lack of knowledge on the genetic variation in genes within inversions, possibly targets of climatic selection, across a geographic latitudinal gradient. In the present study we analysed four candidate loci for thermal adaptation, located close to the breakpoints, in two chromosomal arrangements of the sex (A) chromosome of Drosophila subobscura with different thermal preferences. Individual chromosomes with A2 (the inverted arrangement considered warm adapted) or AST (the standard ancestral arrangement considered cold adapted) were sequenced across four European localities at varying latitudes, up to ~ 2500 Kms apart. RESULTS: Importantly, we found very low differentiation for each specific arrangement across populations as well as no clinal patterns of genomic variation. This suggests wide gene exchange along the cline. Differentiation between the sex chromosome arrangements was significant in the two more proximal regions relative to the AST orientation but not in the distal ones, independently of their location inside or outside the inversion. This can be possibly due to variation in the levels of gene flux and/or selection acting in these regions. CONCLUSIONS: Gene flow appears to have homogenized the genetic content within-arrangement at a wide geographical scale, despite the expected diverse selective pressures in the specific natural environments of the different populations sampled. It is thus likely that the inversion frequency clines in this species are being maintained by local adaptation in face of gene flow. The differences between arrangements at non-coding regions might be associated with the previously observed differential gene expression in different thermal regimes. Higher resolution genomic scans for individual chromosomal arrangements performed over a large environmental gradient are needed to find the targets of selection and further elucidate the adaptive mechanisms maintaining chromosomal inversion polymorphisms.

Reproductive interference by male Drosophila subobscura on female D. persimilis: A laboratory experiment.

While females often reject courtship attempts by heterospecific males, reproductive interference by harassment from such males can nonetheless occur, potentially reducing female fitness. Such effects may be profound following a range expansion, when males from a new species may suddenly encounter (and perhaps even become abundant relative to) females of related native species. Drosophila subobscura recently invaded North America and may impact native species through reproductive interference and other processes. We test for the potential for reproductive interference by D. subobscura males on D. persimilis females in the laboratory. D. subobscura males aggressively copulated with D. persimilis females, including many females that exhibit rejection behaviors. Despite females attempting to dismount the males, the heterospecific copulations are on average longer than conspecific copulations, and females exhibit some reluctance to remate with conspecific males following this harassment. Females confined with both conspecific and heterospecific males produce fewer adult progeny than those with either conspecific males only or with conspecific males and distantly related D. simulans males that do not engage in female harassment. Overall, our results illustrate how reproductive interference by an invasive species can have negative effects on resident natural populations.

Tracking changes in chromosomal arrangements and their genetic content during adaptation.

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life-history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness-related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O3+4 arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.

Individual inversions or their combinations: which is the main selective target in a natural population of Drosophila subobscura?

It is generally accepted that chromosomal inversions have been key elements in adaptation and speciation processes. In this context, Drosophila subobscura has been, and still is, an excellent model species due to its rich chromosomal polymorphism. In this species, many analyses from natural populations have demonstrated the adaptive potential of individual inversions (and their overlapped combinations, the so-called arrangements). However, little information is available on the evolutionary role of combinations generated by inversions located in homologous and nonhomologous chromosomes. The aim of this research was to ascertain whether these combinations are also a target for natural selection. For this objective, we have studied the inversion composition of homologous and nonhomologous chromosomes from a D. subobscura sample collected in a well-studied population, Mount Avala (Serbia). No significant deviation from H-W expectations was detected, and when comparing particular karyotypic combinations, likelihood ratios close to 1 were obtained. Thus, it seems that for each pair of homologous chromosomes inversions no deviation from randomness was detected. Finally, no linkage disequilibrium was observed between inversions located in different chromosomes of the karyotype. For all these reasons, it can be assumed that, at the cytological level, the individual inversions rather than their combinations in different chromosomes are the main target of selection.

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.

Evolution of mating behavior between two populations adapting to common environmental conditions.

Populations from the same species may be differentiated across contrasting environments, potentially affecting reproductive isolation among them. When such populations meet in a novel common environment, this isolation may be modified by biotic or abiotic factors. Curiously, the latter have been overlooked. We filled this gap by performing experimental evolution of three replicates of two populations of Drosophila subobscura adapting to a common laboratorial environment, and simulated encounters at three time points during this process. Previous studies showed that these populations were highly differentiated for several life-history traits and chromosomal inversions. First, we show initial differentiation for some mating traits, such as assortative mating and male mating rate, but not others (e.g., female mating latency). Mating frequency increased during experimental evolution in both sets of populations. The assortative mating found in one population remained constant throughout the adaptation process, while disassortative mating of the other population diminished across generations. Additionally, differences in male mating rate were sustained across generations. This study shows that mating behavior evolves rapidly in response to adaptation to a common abiotic environment, although with a complex pattern that does not correspond to the quick convergence seen for life-history traits.

Drosophila subobscura flies adapted to low lead concentration carry no fitness cost.

As a response to the long-term presence of heavy metals in the environment, populations can evolve resistance. Its maintenance may have detrimental effect on population's fitness, causing a fitness cost. Lead is one of the widely distributed elements in the environment exhibiting high toxicity on organisms. By analyzing developmental stages viability and developmental time, we evaluated fitness cost in Drosophila subobscura flies adapted to low lead concentration and control flies derived from the same wild population, as well as their hybrids. Significant changes in specific developmental stages viability were detected in both lines, as well as their hybrids, suggesting complex response to low lead concentration. The results show that a long-term exposure to low lead concentration may have a significant impact on a population's survival, especially in a changing environment conditions.

How much can history constrain adaptive evolution? A real-time evolutionary approach of inversion polymorphisms in Drosophila subobscura.

Chromosomal inversions are present in a wide range of animals and plants, having an important role in adaptation and speciation. Although empirical evidence of their adaptive value is abundant, the role of different processes underlying evolution of chromosomal polymorphisms is not fully understood. History and selection are likely to shape inversion polymorphism variation to an extent yet largely unknown. Here, we perform a real-time evolution study addressing the role of historical constraints and selection in the evolution of these polymorphisms. We founded laboratory populations of Drosophila subobscura derived from three locations along the European cline and followed the evolutionary dynamics of inversion polymorphisms throughout the first 40 generations. At the beginning, populations were highly differentiated and remained so throughout generations. We report evidence of positive selection for some inversions, variable between foundations. Signs of negative selection were more frequent, in particular for most cold-climate standard inversions across the three foundations. We found that previously observed convergence at the phenotypic level in these populations was not associated with convergence in inversion frequencies. In conclusion, our study shows that selection has shaped the evolutionary dynamics of inversion frequencies, but doing so within the constraints imposed by previous history. Both history and selection are therefore fundamental to predict the evolutionary potential of different populations to respond to global environmental changes.

Inbreeding and thermal adaptation in Drosophila subobscura.

Using a well-adapted Drosophila subobscura population (Avala, Serbia), a drastic experiment of inbreeding was carried out to assess whether the expected level of homozygosity could be reached or if other evolutionary forces affected the process. In general, no significant changes of inversion (or arrangement) frequencies were detected after 12 brother-sister mating generations. Furthermore, no significant differences were obtained between observed and expected (under the inbreeding model) karyotypic frequencies. Thus, these results seemed to indicate that the main evolutionary factor in the experiment was inbreeding. However, in the G12 generation, complete chromosomal fixation was reached only in two out of the eight final inbred lines. In these lines, the chromosomal compositions were difficult to interpret, but they could be likely a consequence of adaptation to particular laboratory conditions (constant 18 °C, food, light period, etc.). Finally, in a second experiment, the inbred lines presented higher fertility at 18 °C than at 13 °C. Also, there was a significant line effect on fertility: inbred line number 6 (A1, J1, U1+2; U1+2+6, E8, and O3+4+7) presented the highest values, which maybe the result of an adaptation to laboratory conditions. Thus, the results obtained in our experiments reflect the adaptive potential of D. subobscura inversions.

Changes in chromosomal polymorphism and global warming: The case of Drosophila subobscura from Apatin (Serbia).

In this study, chromosomal inversion polymorphism data for a natural population of Drosophila subobscura from a swampy region near the town of Apatin (Serbia) were compared with data for the same population collected approximately 15 years earlier. The pattern of chromosomal inversion polymorphism changed over time. There were significant increases in the frequency of characteristic southern latitude ("warm" adapted) chromosomal arrangements and significant decreases in the frequency of characteristic northern latitude ("cold" adapted) chromosomal arrangements in the O and U chromosomes. The chromosomal arrangements O(3+4) and O(3+4) (+) (22) (derived from the O(3+4) arrangement) showed significant increases in 2008 and 2009 with regard to the 1994 sample. There was also a significant increase (∼50%) in the U(1) (+) (2) arrangement, while U(1+8) (+) (2) (a typical southern arrangement) was detected for the first time. Since the Apatin swampy population of D. subobscura has existed for a long time in a stable habitat with high humidity that has not been changed by man our results indicate that natural selection has produced chromosomal changes in response to the increase in temperature that has occurred in the Balkan Peninsula of central southeastern European.

Changes in chromosomal polymorphism and global warming: the case of Drosophila subobscura from Apatin (Serbia)

In this study, chromosomal inversion polymorphism data for a natural population of Drosophila subobscura from a swampy region near the town of Apatin (Serbia) were compared with data for the same population collected approximately 15 years earlier. The pattern of chromosomal inversion polymorphism changed over time. There were significant increases in the frequency of characteristic southern latitude ("warm" adapted) chromosomal arrangements and significant decreases in the frequency of characteristic northern latitude ("cold" adapted) chromosomal arrangements in the O and U chromosomes. The chromosomal arrangements O3+4 and O3+4+22 (derived from the O3+4 arrangement) showed significant increases in 2008 and 2009 with regard to the 1994 sample. There was also a significant increase (~50 percent) in the U1+2 arrangement, while U1+8+2 (a typical southern arrangement) was detected for the first time. Since the Apatin swampy population of D. subobscura has existed for a long time in a stable habitat with high humidity that has not been changed by man our results indicate that natural selection has produced chromosomal changes in response to the increase in temperature that has occurred in the Balkan Peninsula of central southeastern European.

INTERSPECIFIC LABORATORY COMPETITION OF THE RECENTLY SYMPATRIC SPECIES DROSOPHILA SUBOBSCURA AND DROSOPHILA PSEUDOOBSCURA.

-Drosophila subobscura and D. pseudoobscura are closely related species coexisting on the West Coast of North America, which was recently colonized by D. subobscura. In competition experiments with overlapping generations, D. subobscura is eliminated by D. pseudoobscura in a few generations at all four temperatures and two initial frequencies tested. Yet in one-species cultures, D. subobscura thrives at all experimental conditions. Single-generation competition experiments reveal lower survivorship and productivity of D. subobscura at all temperatures and frequencies. Productivity per female is dependent on the initial frequencies: greater for D. subobscura as its initial frequency becomes higher, but lower for D. pseudoobscura as its frequency becomes higher. Strains of D. subobscura from three disparate geographic origins yield similar results.