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

Adaptation of Drosophila subobscura chromosomal inversions to climatic variables: the Balkan natural population of Avala

The adaptive value of chromosomal inversions continues raising relevant questions in evolutionary biology. In many species of the Drosophila genus, different inversions have been recognized to be related to thermal adaptation, but it is necessary to determine to which specific climatic variables the inversions are adaptive. With this aim, the behavior of thermal adapted inversions of Drosophila subobscura regarding climatic variables was studied in the natural population of Avala (Serbia) during the 2014-2017 period. The results obtained were compared with those previously reported in the Font Groga (Barcelona, Spain) population, which presents different climatic and environmental conditions. In both populations, it was observed that most thermal adapted inversions were significantly associated with the first, second or both principal components, which were related with maximum, minimum and mean temperatures. Moreover, a significant increase over years (2004-2017) for the minimum temperature was detected. In parallel, a significant variation over time in Avala was only observed for the frequencies of 'warm' and 'non-thermal' adapted inversions of the U chromosome. However, stability in the chromosomal inversion polymorphism was observed for the 2014-2017 period which might result from the temporal span of the study and/or selective process acting on the population.

Rate of change for the thermal adapted inversions in Drosophila subobscura

The changes of chromosomal inversion polymorphism composition of Drosophila subobscura in samples from Apatin (Serbia) were studied in a 24-years interval (1994-2018). The variation was significant for all autosomes and directional, increasing the inversions considered as 'warm', whereas those reported as 'cold' decreased. Furthermore, the Chromosomal Thermal Index (CTI), which allows studying the thermal adaptation of the whole karyotype increased significantly in that period of time. These results were in agreement with the indicators of global warming in Apatin: a trend to increase of the mean, maximum and minimum (this latter even significant) temperatures, and an erratic pattern of rainfall (also usual in global warming). The deviations from the Wright-Fisher model of genetic drift were used to consider the possible effect of migration or selection as evolutionary factors responsible for the change in inversion frequencies. To quantify approximately the rate of change in the frequencies, for each kind of inversions ('cold', 'warm' and 'non-thermal adapted'), the difference in frequency between the Apatin samples obtained in 1994 and 2018 was computed and then it was divided by the number of years elapsed. This rate was always higher (from twice as many as thirty times more depending on the autosome) for thermal adapted inversions ('cold' or 'warm') than the 'non-thermal' adapted. From this study, it could be concluded that the chromosomal inversions of D. subobscura could change (in composition and frequencies) in a predictable direction and a rather 'rapid' rhythm to adapt to the global warming scenario.

Long-read based assembly and synteny analysis of a reference Drosophila subobscura genome reveals signatures of structural evolution driven by inversions recombination-suppression effects

BACKGROUND

Drosophila subobscura has long been a central model in evolutionary genetics. Presently, its use is hindered by the lack of a reference genome. To bridge this gap, here we used PacBio long-read technology, together with the available wealth of genetic marker information, to assemble and annotate a high-quality nuclear and complete mitochondrial genome for the species. With the obtained assembly, we performed the first synteny analysis of genome structure evolution in the subobscura subgroup.

RESULTS

We generated a highly-contiguous ~ 129 Mb-long nuclear genome, consisting of six pseudochromosomes corresponding to the six chromosomes of a female haploid set, and a complete 15,764 bp-long mitogenome, and provide an account of their numbers and distributions of codifying and repetitive content. All 12 identified paracentric inversion differences in the subobscura subgroup would have originated by chromosomal breakage and repair, with some associated duplications, but no evidence of direct gene disruptions by the breakpoints. Between lineages, inversion fixation rates were 10 times higher in continental D. subobscura than in the two small oceanic-island endemics D. guanche and D. madeirensis. Within D. subobscura, we found contrasting ratios of chromosomal divergence to polymorphism between the A sex chromosome and the autosomes.

CONCLUSIONS

We present the first high-quality, long-read sequencing of a D. subobscura genome. Our findings generally support genome structure evolution in this species being driven indirectly, through the inversions' recombination-suppression effects in maintaining sets of adaptive alleles together in the face of gene flow. The resources developed will serve to further establish the subobscura subgroup as model for comparative genomics and evolutionary indicator of global change.

Climatic adaptation of chromosomal inversions in Drosophila subobscura

Drosophila subobscura is a species with a rich chromosomal polymorphism which is adaptive to different climatic conditions. Five samples of the Font Groga population (Barcelona, Spain) were sampled in autumn during 5 consecutive years (2011-2015) to obtain their inversion chromosomal polymorphism, and climatic data of several meteorological variables were also collected. The aim was to analyze the adaptive potential of inversions with regard to climatic variables, being the most relevant: mean temperature (T_mean), maximum temperature (T_max), minimum temperature (T_min), humidity (Hm) and rainfall (Rf). As expected, no significant variation in inversion frequencies were detected over this short period of time. However, from a climatic point of view it was possible to differentiate 'warm' and 'dry' from 'cold' and 'humid' samples. The joint study of maximum (T_max) and minimum (T_min) temperatures was a key element to understand the effect on adaptation of many inversions. It was also observed that temperature had to be considered in conjunction with humidity and rainfall. All these factors would condition the biota of D. subobscura habitat, and chromosomal inversions could provide an adaptive response to it.

Chromosomal Thermal Index: a comprehensive way to integrate the thermal adaptation of Drosophila subobscura whole karyotype

Drosophila has demonstrated to be an excellent model to study the adaptation of organisms to global warming, with inversion chromosomal polymorphism having a key role in this adaptation. Here, we introduce a new index (Chromosomal Thermal Index or CTI) to quantify the thermal adaptation of a population according to its composition of "warm" and "cold" adapted inversions. This index is intuitive, has good statistical properties, and can be used to hypothesis on the effect of global warming on natural populations. We show the usefulness of CTI using data from European populations of D. subobscura, sampled in different years. Out of 15 comparisons over time, nine showed significant increase of CTI, in accordance with global warming expectations. Although large regions of the genome outside inversions contain thermal adaptation genes, our results show that the total amount of warm or cold inversions in populations seems to be directly involved in thermal adaptation, whereas the interactions between the inversions content of homologous and non-homologous chromosomes are not relevant.