Tracking changes in chromosomal arrangements and their genetic content during adaptation

Heat-induced female biased sex ratio during development is not mitigated after prolonged thermal selection

BACKGROUND

The negative impacts of climate change on biodiversity are consistently increasing. Developmental stages are particularly sensitive in many ectotherms. Moreover, sex-specific differences in how organisms cope with thermal stress can produce biased sex ratios upon emergence, with potentially major impacts on population persistence. This is an issue that needs investigation, particularly testing whether thermal selection can alleviate sex ratio distortions in the long-term is a critical but neglected issue. Here, we report an experiment analyzing the sex ratio patterns at different developmental temperatures in Drosophila subobscura populations subjected to long-term experimental evolution (~ 30 generations) under a warming environment.

RESULTS

We show that exposure to high developmental temperatures consistently promotes sex ratio imbalance upon emergence, with a higher number of female than male offspring. Furthermore, we found that thermal selection resulting from evolution in a warming environment did not alleviate such sex ratio distortions generated by heat stress.

CONCLUSIONS

We demonstrate that heat stress during development can lead to clear sex ratio deviations upon emergence likely because of differential survival between sexes. In face of these findings, it is likely that sex ratio deviations of this sort occur in natural populations when facing environmental perturbation. The inability of many insects to avoid thermal shifts during their (more) sessile developmental stages makes this finding particularly troublesome for population subsistence in face of climate warming events.

No evidence for short-term evolutionary response to a warming environment in Drosophila

Adaptive evolution is key in mediating responses to global warming and may sometimes be the only solution for species to survive. Such evolution will expectedly lead to changes in the populations' thermal reaction norm and improve their ability to cope with stressful conditions. Conversely, evolutionary constraints might limit the adaptive response. Here, we test these expectations by performing a real-time evolution experiment in historically differentiated Drosophila subobscura populations. We address the phenotypic change after nine generations of evolution in a daily fluctuating environment with average constant temperature, or in a warming environment with increasing average and amplitude temperature across generations. Our results showed that (1) evolution under a global warming scenario does not lead to a noticeable change in the thermal response; (2) historical background appears to be affecting responses under the warming environment, particularly at higher temperatures; and (3) thermal reaction norms are trait dependent: although lifelong exposure to low temperature decreases fecundity and productivity but not viability, high temperature causes negative transgenerational effects on productivity and viability, even with high fecundity. These findings in such an emblematic organism for thermal adaptation studies raise concerns about the short-term efficiency of adaptive responses to the current rising temperatures.

Deleterious mutation accumulation and the long-term fate of chromosomal inversions

Chromosomal inversions contribute widely to adaptation and speciation, yet they present a unique evolutionary puzzle as both their allelic content and frequency evolve in a feedback loop. In this simulation study, we quantified the role of the allelic content in determining the long-term fate of the inversion. Recessive deleterious mutations accumulated on both arrangements with most of them being private to a given arrangement. This led to increasing overdominance, allowing for the maintenance of the inversion polymorphism and generating strong non-adaptive divergence between arrangements. The accumulation of mutations was mitigated by gene conversion but nevertheless led to the fitness decline of at least one homokaryotype under all considered conditions. Surprisingly, this fitness degradation could be permanently halted by the branching of an arrangement into multiple highly divergent haplotypes. Our results highlight the dynamic features of inversions by showing how the non-adaptive evolution of allelic content can play a major role in the fate of the inversion.

Distribution modelling of an introduced species: do adaptive genetic markers affect potential range?

Biological invasions have increased in the last few decades mostly due to anthropogenic causes such as globalization of trade. Because invaders sometimes cause large economic losses and ecological disturbances, estimating their origin and potential geographical ranges is useful. Drosophila subobscura is native to the Old World but was introduced in the New World in the late 1970s and spread widely. We incorporate information on adaptive genetic markers into ecological niche modelling and then estimate the most probable geographical source of colonizers; evaluate whether the genetic bottleneck experienced by founders affects their potential distribution; and finally test whether this species has spread to all its potential suitable habitats worldwide. We find the environmental space occupied by this species in its native and introduced distributions are notably the same, although the introduced niche has shifted slightly towards higher temperature and lower precipitation. The genetic bottleneck of founding individuals was a key factor limiting the spread of this introduced species. We also find that regions in the Mediterranean and north-central Portugal show the highest probability of being the origin of the colonizers. Using genetically informed environmental niche modelling can enhance our understanding of the initial colonization and spread of invasive species, and also elucidate potential areas of future expansions worldwide.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

Comparative analysis of adaptive and neutral markers of Drosophila mediopunctata populations dispersed among forest fragments

Comparison of adaptive and neutral genetic markers is a valuable approach to characterize the evolutionary consequences of populations living in environments threatened by anthropogenic disturbances, such as forest fragmentation. Shifts in allele frequencies, low genetic variability, and a small effective population size can be considered clear signs of forest fragmentation effects (due to genetic drift) over natural populations, while adaptive responses correlate with environmental variables. Brazilian Atlantic Forest had its landscape drastically reduced and fragmented. Now, several forest remnants are isolated from each other by urban and crop areas. We sampled Drosophila mediopunctata populations from eight forest remnants dispersed on two adjacent geomorphological regions, which are physiognomic and climatically quite distinct. Microsatellite data of inversion-free chromosomes (neutral genetic marker) indicate low structuration among populations suggesting that they were panmictic and greatly influenced by gene flow. Moreover, significant differences in chromosomal inversion frequencies (adaptive genetic marker) among populations and their correlations with climatic and geographical variables indicate that genetic divergence among populations could be an adaptive response to their environment. Nonetheless, we observed a significant difference in inversion frequencies of a population in two consecutive years that may be associated with edge and demographic effects. Also, it may be reflecting seasonal changes of inversion frequencies influenced by great temperature variation due to edge effects. Moreover, the forest fragment size does not affect genetic variation of neutral markers. Our data indicate that despite oscillations in chromosomal inversion frequencies, D. mediopunctata populations from Brazilian Atlantic Forest and their divergence may be driven by adaptive factors to local differences, perhaps because it is a small flying insect easily carried by the wind increasing its migration rates.

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.

Different Genomic Changes Underlie Adaptive Evolution in Populations of Contrasting History

Experimental evolution is a powerful tool to understand the adaptive potential of populations under environmental change. Here, we study the importance of the historical genetic background in the outcome of evolution at the genome-wide level. Using the natural clinal variation of Drosophila subobscura, we sampled populations from two contrasting latitudes (Adraga, Portugal and Groningen, Netherlands) and introduced them in a new common environment in the laboratory. We characterized the genome-wide temporal changes underlying the evolutionary dynamics of these populations, which had previously shown fast convergence at the phenotypic level, but not at chromosomal inversion frequencies. We found that initially differentiated populations did not converge either at genome-wide level or at candidate SNPs with signs of selection. In contrast, populations from Portugal showed convergence to the control population that derived from the same geographical origin and had been long-established in the laboratory. Candidate SNPs showed a variety of different allele frequency change patterns across generations, indicative of an underlying polygenic basis. We did not detect strong linkage around candidate SNPs, but rather a small but long-ranging effect. In conclusion, we found that history played a major role in genomic variation and evolution, with initially differentiated populations reaching the same adaptive outcome through different genetic routes.

Predictable phenotypic, but not karyotypic, evolution of populations with contrasting initial history

The relative impact of selection, chance and history will determine the predictability of evolution. There is a lack of empirical research on this subject, particularly in sexual organisms. Here we use experimental evolution to test the predictability of evolution. We analyse the real-time evolution of Drosophila subobscura populations derived from contrasting European latitudes placed in a novel laboratory environment. Each natural population was sampled twice within a three-year interval. We study evolutionary responses at both phenotypic (life-history, morphological and physiological traits) and karyotypic levels for around 30 generations of laboratory culture. Our results show (1) repeatable historical effects between years in the initial state, at both phenotypic and karyotypic levels; (2) predictable phenotypic evolution with general convergence except for body size; and (3) unpredictable karyotypic evolution. We conclude that the predictability of evolution is contingent on the trait and level of organization, highlighting the importance of studying multiple biological levels with respect to evolutionary patterns.