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

Simple Summary

Microbiota have a significant functional role in the life of the host, including immunity, lifespan and reproduction. Drosophila species are attractive model organisms for investigating microbiota diversity from different aspects due to their simple gut microbiota, short generation time and high fertility. Considering such an important role of the microbiota in the life of Drosophila, we investigated the extent to which lead (Pb), as one of the most abundant heavy metals in the environment, affects the microbiota and the fitness of this insect host. The results indicate that different factors, such as population origin and sex, may affect individual traits differently and this could be species-specific. In addition, there are members of microbiota that help the host to overcome environmental stress and they could play a key role in reducing the fitness cost in such situations. Studying the influence of microbiota on the adaptive response to heavy metals and the potential implications on overall host fitness is of great pertinence.

Abstract

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 O_{3 + 4 + 7} and O_ST 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 O_ST and O_{3 + 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 O_{3 + 4 + 7} lines, which was not found in any of the six O_ST lines.

Conclusions

The equivalent hsp70 mRNA amounts in O_ST and O_{3 + 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.

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 A₂ (the inverted arrangement considered warm adapted) or A_ST (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 A_ST 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.

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.

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.

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.

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.

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₃₊₄ and O₃₊₄₊₂₂ (derived from the O₃₊₄ arrangement) showed significant increases in 2008 and 2009 with regard to the 1994 sample. There was also a significant increase (∼50%) in the U₁₊₂ arrangement, while U₁₊₈₊₂ (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.