Colonization of America by Drosophila subobscura: Experiment in natural populations that supports the adaptive role of chromosomal-inversion polymorphism

When adaptation is slowed down: Genomic analysis of evolutionary stasis in thermal tolerance during biological invasion in a novel climate

Research conducted during the past two decades has demonstrated that biological invasions are excellent models of rapid evolution. Even so, characteristics of invasive populations such as a short time for recombination to assemble optimal combinations of alleles may occasionally limit adaptation to new environments. Here, we investigated such genetic constraints to adaptation in the invasive brown anole (Anolis sagrei)-a tropical ectotherm that was introduced to the southeastern United States, a region with a much colder climate than in its native Caribbean range. We examined thermal physiology for 30 invasive populations and tested for a climatic cline in cold tolerance. Also, we used genomics to identify mechanisms that may limit adaptation. We found no support for a climatic cline, indicating that thermal tolerance did not shift adaptively. Concomitantly, population genomic results were consistent with the occurrence of recombination cold spots that comprise more than half of the genome and maintain long-range associations among alleles in invasive populations. These genomic regions overlap with both candidate thermal tolerance loci that we identified using a standard genome-wide association test. Moreover, we found that recombination cold spots do not have a large contribution to population differentiation in the invasive range, contrary to observations in the native range. We suggest that limited recombination is constraining the contribution of large swaths of the genome to adaptation in invasive brown anoles. Our study provides an example of evolutionary stasis during invasion and highlights the possibility that reduced recombination occasionally slows down adaptation in invasive populations.

A lack of genetic diversity and minimal adaptive evolutionary divergence in introduced Mysis shrimp after 50 years

The successes of introduced populations in novel habitats often provide powerful examples of evolution and adaptation. In the 1950s, opossum shrimp (Mysis diluviana) individuals from Clearwater Lake in Minnesota, USA were transported and introduced to Twin Lakes in Colorado, USA by fisheries managers to supplement food sources for trout. Mysis were subsequently introduced from Twin Lakes into numerous lakes throughout Colorado. Because managers kept detailed records of the timing of the introductions, we had the opportunity to test for evolutionary divergence within a known time interval. Here, we used reduced representation genomic data to investigate patterns of genetic diversity, test for genetic divergence between populations, and for evidence of adaptive evolution within the introduced populations in Colorado. We found very low levels of genetic diversity across all populations, with evidence for some genetic divergence between the Minnesota source population and the introduced populations in Colorado. There was little differentiation among the Colorado populations, consistent with the known provenance of a single founding population, with the exception of the population from Gross Reservoir, Colorado. Demographic modeling suggests that at least one undocumented introduction from an unknown source population hybridized with the population in Gross Reservoir. Despite the overall low genetic diversity we observed, F ST outlier and environmental association analyses identified multiple loci exhibiting signatures of selection and adaptive variation related to elevation and lake depth. The success of introduced species is thought to be limited by genetic variation, but our results imply that populations with limited genetic variation can become established in a wide range of novel environments. From an applied perspective, the observed patterns of divergence between populations suggest that genetic analysis can be a useful forensic tool to determine likely sources of invasive species.

The Adaptive Value of Chromosomal Inversions and Climatic Change-Studies on the Natural Populations of Drosophila subobscura from the Balkans

The adaptive value of the Drosophila subobscura chromosomal inversion polymorphism with regard to environmental effects is well-known. However, the specific details of the inversion adaptations to the global warming scenario deserve to be analyzed. Toward this aim, polymorphism and karyotypes were studied in 574 individuals from Petnica (Serbia) in annual samples taken in June for the period 2019-2022. Comparing the results of Petnica (Cfa: humid subtropical climate) with those from Avala (Serbia: Cfb, temperate oceanic climate) and Font Groga (Barcelona, Spain; Csa: hot-summer Mediterranean climate), significant differences were observed for their chromosomal polymorphism. In Petnica, inversions from U and E chromosomes mainly reacted significantly with regard to temperature, humidity, and rainfall. Moreover, the inversion polymorphism from Petnica (2019-2022) was compared with that from 1995. In this period, a significant increase in mean and maximum temperature was observed. However, to properly explain the observed variations of inversions over time, it was necessary to carefully analyze annual seasonal changes and particular heat wave episodes. Interestingly, yearly fluctuations of U chromosome 'warm'-adapted inversions corresponded with opposite changes in 'non-thermal' inversions. Perhaps these types of inversions were not correctly defined with regard to thermal adaptation, or these fluctuations were also due to adaptations to other physical and/or biological variables. Finally, a joint study of chromosomal inversion polymorphism from many Balkan populations of D. subobscura indicated that different climatic regions presented distinct composition, including thermal-adapted inversions.

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.

Past history shapes evolution of reproductive success in a global warming scenario

Adaptive evolution is critical for animal populations to thrive in the fast-changing natural environments. Ectotherms are particularly vulnerable to global warming and, although their limited coping ability has been suggested, few real-time evolution experiments have directly accessed their evolutionary potential. Here, we report a long-term experimental evolution study addressing the evolution of Drosophila thermal reaction norms, after ∼30 generations under different dynamic thermal regimes: fluctuating (daily variation between 15 and 21 °C) or warming (daily fluctuation with increases in both thermal mean and variance across generations). We analyzed the evolutionary dynamics of Drosophila subobscura populations as a function of the thermally variable environments in which they evolved and their distinct background. Our results showed clear differences between the historically differentiated populations: high latitude D. subobscura populations responded to selection, improving their reproductive success at higher temperatures whereas their low latitude counterparts did not. This suggests population variation in the amount of genetic variation available for thermal adaptation, an aspect that needs to be considered to allow for better predictions of future climate change responses. Our results highlight the complex nature of thermal responses in face of environmental heterogeneity and emphasize the importance of considering inter-population variation in thermal evolution studies.

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.

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.

High developmental temperature leads to low reproduction despite adult temperature

Phenotypic plasticity can help organisms cope with changing thermal conditions and it may depend on which life-stage the thermal stress is imposed: for instance, exposure to stressful temperatures during development can trigger a positive plastic response in adults. Here, we analyze the thermal plastic response of laboratory populations of Drosophila subobscura, derived from two contrasting latitudes of the European cline. We measured reproductive performance through fecundity characters, after the experimental populations were exposed to five thermal treatments, with different combinations of developmental and adult temperatures (14 °C, 18 °C, or 26 °C). Our questions were whether (1) adult performance changes with exposure to higher (or lower) temperatures during development; (2) flies raised at lower temperatures outperform those developed at higher ones, supporting the "colder is better" hypothesis; (3) there is a cumulative effect on adult performance of exposing both juveniles and adults to higher (or lower) temperatures; (4) there is evidence for biogeographical effects on adult performance. Our main findings were that (1) higher developmental temperatures led to low reproductive performance regardless of adult temperature, while at lower temperatures reduced performance only occurred when colder conditions were persistent across juvenile and adult stages; (2) flies raised at lower temperatures did not always outperform those developed at other temperatures; (3) there were no harmful cumulative effects after exposing both juveniles and adults to higher temperatures; (4) both latitudinal populations showed similar thermal plasticity patterns. The negative effect of high developmental temperature on reproductive performance, regardless of adult temperature, highlights the developmental stage as very critical and most vulnerable to climate change and associated heat waves.

Detecting inversions with PCA in the presence of population structure

Chromosomal inversions can lead to reproductive isolation and adaptation in insects such as Drosophila melanogaster and the non-model malaria vector Anopheles gambiae. Inversions can be detected and characterized using principal component analysis (PCA) of single nucleotide polymorphisms (SNPs). To aid in developing such methods, we formed a new benchmark derived from three publicly-available insect data. We then used this benchmark to perform an extended validation of our software for inversion analysis (Asaph). Through that process, we identified and characterized several problematic test cases liable to misinterpretation that can help guide PCA-based inversion detection. Lastly, we re-analyzed the 2R chromosome arm of 150 An. gambiae and coluzzii samples and observed two inversions (2Rc and 2Rd) that were previously known but not annotated in these particular individuals. The resulting benchmark data set and methods will be useful for future inversion detection based solely on SNP data.

The Cyclically Seasonal Drosophila subobscura Inversion O7 Originated From Fragile Genomic Sites and Relocated Immunity and Metabolic Genes

Chromosome inversions are important contributors to standing genetic variation in Drosophila subobscura. Presently, the species is experiencing a rapid replacement of high-latitude by low-latitude inversions associated with global warming. Yet not all low-latitude inversions are correlated with the ongoing warming trend. This is particularly unexpected in the case of O7 because it shows a regular seasonal cycle that peaks in summer and rose with a heatwave. The inconsistent behavior of O7 across components of the ambient temperature suggests that is causally more complex than simply due to temperature alone. In order to understand the dynamics of O7, high-quality genomic data are needed to determine both the breakpoints and the genetic content. To fill this gap, here we generated a PacBio long read-based chromosome-scale genome assembly, from a highly homozygous line made isogenic for an O3 + 4 + 7 chromosome. Then we isolated the complete continuous sequence of O7 by conserved synteny analysis with the available reference genome. Main findings include the following: (i) the assembled O7 inversion stretches 9.936 Mb, containing > 1,000 annotated genes; (ii) O7 had a complex origin, involving multiple breaks associated with non-B DNA-forming motifs, formation of a microinversion, and ectopic repair in trans with the two homologous chromosomes; (iii) the O7 breakpoints carry a pre-inversion record of fragility, including a sequence insertion, and transposition with later inverted duplication of an Attacin immunity gene; and (iv) the O7 inversion relocated the major insulin signaling forkhead box subgroup O (foxo) gene in tight linkage with its antagonistic regulatory partner serine/threonine-protein kinase B (Akt1) and disrupted concerted evolution of the two inverted Attacin duplicates, reattaching them to dFOXO metabolic enhancers. Our findings suggest that O7 exerts antagonistic pleiotropic effects on reproduction and immunity, setting a framework to understand its relationship with climate change. Furthermore, they are relevant for fragility in genome rearrangement evolution and for current views on the contribution of breakage versus repair in shaping inversion-breakpoint junctions.

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.

Resolving the genetic paradox of invasions: Preadapted genomes and postintroduction hybridization of bigheaded carps in the Mississippi River Basin

The genetic paradox of biological invasions is complex and multifaceted. In particular, the relative role of disparate propagule sources and genetic adaptation through postintroduction hybridization has remained largely unexplored. To add resolution to this paradox, we investigate the genetic architecture responsible for the invasion of two invasive Asian carp species, bighead carp (Hypophthalmichthys nobilis) and silver carp (H. molitrix) (bigheaded carps) that experience extensive hybridization in the Mississippi River Basin (MRB). We sequenced the genomes of bighead and silver carps (~1.08G bp and ~1.15G bp, respectively) and their hybrids collected from the MRB. We found moderate-to-high heterozygosity in bighead (0.0021) and silver (0.0036) carps, detected significantly higher dN/dS ratios of single-copy orthologous genes in bigheaded carps versus 10 other species of fish, and identified genes in both species potentially associated with environmental adaptation and other invasion-related traits. Additionally, we observed a high genomic similarity (96.3% in all syntenic blocks) between bighead and silver carps and over 90% embryonic viability in their experimentally induced hybrids. Our results suggest intrinsic genomic features of bigheaded carps, likely associated with life history traits that presumably evolved within their native ranges, might have facilitated their initial establishment of invasion, whereas ex-situ interspecific hybridization between the carps might have promoted their range expansion. This study reveals an alternative mechanism that could resolve one of the genetic paradoxes in biological invasions and provides invaluable genomic resources for applied research involving bigheaded carps.

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.

Effects of photoperiod on life-history and thermal stress resistance traits across populations of Drosophila subobscura

INTRODUCTION

Organisms use environmental cues to match their phenotype with the future availability of resources and environmental conditions. Changes in the magnitude and frequency of environmental cues such as photoperiod and temperature along latitudes can be used by organisms to predict seasonal changes. While the role of temperature variation on the induction of plastic and seasonal responses is well established, the importance of photoperiod for predicting seasonal changes is less explored.

MATERIALS AND METHODS

Here we studied changes in life-history and thermal stress resistance traits in Drosophila subobscura in response to variation in photoperiod (6:18, 12:12 and 18:6 light:dark cycles) mimicking seasonal variations in day length. The populations of D. subobscura were collected from five locations along a latitudinal gradient (from North Africa and Europe). These populations were exposed to different photoperiods for two generations, whereafter egg-to-adult viability, productivity, dry body weight, thermal tolerance, and starvation resistance were assessed.

RESULTS

We found strong effects of photoperiod, origin of populations, and their interactions on life-history and stress resistance traits. Thermal resistance varied between the populations and the effect of photoperiod depended on the trait and the method applied for the assessment of thermal resistance.

PERSPECTIVES

Our results show a strong effect of the origin of population and photoperiod on a range of fitness-related traits and provide evidence for local adaptation to environmental cues (photoperiod by population interaction). The findings emphasize an important and often neglected role of photoperiod in studies on thermal resistance and suggest that cues induced by photoperiod may provide some buffer enabling populations to cope with a more variable and unpredictable future climate.

The molecular characterization of fixed inversions breakpoints unveils the ancestral character of the Drosophila guanche chromosomal arrangements

Cytological studies revealed that the number of chromosomes and their organization varies across species. The increasing availability of whole genome sequences of multiple species across specific phylogenies has confirmed and greatly extended these cytological observations. In the Drosophila genus, the ancestral karyotype consists of five rod-like acrocentric chromosomes (Muller elements A to E) and one dot-like chromosome (element F), each exhibiting a generally conserved gene content. Chromosomal fusions and paracentric inversions are thus the major contributors, respectively, to chromosome number variation among species and to gene order variation within chromosomal element. The subobscura cluster of Drosophila consists in three species that retain the genus ancestral karyotype and differ by a reduced number of fixed inversions. Here, we have used cytological information and the D. guanche genome sequence to identify and molecularly characterize the breakpoints of inversions that became fixed since the D. guanche-D. subobscura split. Our results have led us to propose a modified version of the D. guanche cytological map of its X chromosome, and to establish that (i) most inversions became fixed in the D. subobscura lineage and (ii) the order in which the four X chromosome overlapping inversions occurred and became fixed.

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.

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.

Inversion evolutionary rates might limit the experimental identification of inversion breakpoints in non-model species

Chromosomal inversions are structural changes that alter gene order but generally not gene content in the affected region. In Drosophila, extensive cytological studies revealed the widespread character of inversion polymorphism, with evidence for its adaptive character. In Drosophila subobscura, polymorphism affects both its four large autosomal elements and its X (A) chromosome. The characterization of eight of these autosomal inversions breakpoints revealed that most of them originated through the staggered-breaks mechanism. Here, we have performed chromosomal walks to identify the breakpoints of two X-chromosome widely distributed inversions -A₂ and A₁- of D. subobscura. Inversion A₂ is considered a warm-adapted arrangement that exhibits parallel latitudinal clines in the species ancestral distribution area and in both American subcontinents, whereas inversion A₁ is only present in the Palearctic region where it presents an east-west cline. The duplication detected at the A₂ inversion breakpoints is consistent with its origin by the staggered-breaks mechanism. Inversion A₁ breakpoints could not be molecularly identified even though they could be narrowly delimited. This result points to chromosome walking limitations when using as a guide the genome of other species. Limitations stem from the rate of evolution by paracentric inversions, which in Drosophila is highest for the X chromosome.

An easy route to the massive karyotyping of complex chromosomal arrangements in Drosophila

Inversion polymorphism is widespread in the Drosophila genus as well as in other dipteran genera. The presence of polytene chromosomes in some insect organs and, thus, the possibility to observe the different arrangements generated by inversions through a microscope enhanced the cytological study of this structural polymorphism. In several Drosophila species, these studies provided evidence for the adaptive character of this polymorphism, which together with the standing interest to uncover targets of natural selection has led to a renewed interest for inversion polymorphism. Our recent molecular characterization of the breakpoint regions of five inversions of the E chromosome of D. subobscura has allowed us to design a PCR-based strategy to molecularly identify the different chromosomal arrangements and, most importantly, to determine the E chromosome karyotype of medium- and large-sized samples from natural populations. Individuals of a test sample that were both cytologically and molecularly karyotyped were used to establish the strategy that was subsequently applied to karyotype a larger sample. Our strategy has proved to be robust and time efficient, and it lays therefore the groundwork for future studies of the E chromosome structural polymorphism through space and time, and of its putative contribution to adaptation.

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.

Genomics of Parallel Experimental Evolution in Drosophila

What are the genomic foundations of adaptation in sexual populations? We address this question using fitness–character and whole-genome sequence data from 30 Drosophila laboratory populations. These 30 populations are part of a nearly 40-year laboratory radiation featuring 3 selection regimes, each shared by 10 populations for up to 837 generations, with moderately large effective population sizes. Each of 3 sets of the 10 populations that shared a selection regime consists of 5 populations that have long been maintained under that selection regime, paired with 5 populations that had only recently been subjected to that selection regime. We find a high degree of evolutionary parallelism in fitness phenotypes when most-recent selection regimes are shared, as in previous studies from our laboratory. We also find genomic parallelism with respect to the frequencies of single-nucleotide polymorphisms, transposable elements, insertions, and structural variants, which was expected. Entirely unexpected was a high degree of parallelism for linkage disequilibrium. The evolutionary genetic changes among these sexual populations are rapid and genomically extensive. This pattern may be due to segregating functional genetic variation that is abundantly maintained genome-wide by selection, variation that responds immediately to changes of selection regime.

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.

COLONIZATION OF AMERICA BY DROSOPHILA SUBOBSCURA: LETHAL GENES AND THE PROBLEM OF THE O5 INVERSION

In this work, the process of colonization of North and South America by the species Drosophila subobscura has been studied by analyzing the variability of lethal genes. The genetic structures of a Palearctic natural central population (Bordils, Spain) and a colonizer population from America (Gilroy, California) have been compared. The frequencies of lethal chromosomes and their allelism are 29.007% and 0.0069 in the first population and 14.414% and 0.0526 in the American population. A founder effect is detected after the computation of some population parameters (N_e , h, he and the lethal load). Furthermore, the allelism of lethal chromosomes has revealed a strong association between a lethal gene and the O₅ inversion both in Gilroy and in the population of Puerto Montt (Chile). The interpopulation allelism shows that the O₅ arrangement from the USA and Chile is the same, confirming that the colonizing processes of North and South America are correlated. The O₅ arrangement can also be useful as a genetic marker to trace the origin of the colonization. The frequency of the O₅ arrangement in the original population of the colonization could be used to estimate the number of colonizers. This population is still unknown, but taking the extreme values of the frequency of the O₅ inversion in natural Palearctic populations (1-15%), the number of colonizers could vary between 9 and 149 individuals.

Dense gene physical maps of the non-model species Drosophila subobscura

The comparative analysis of genetic and physical maps as well as of whole genome sequences had revealed that in the Drosophila genus, most structural rearrangements occurred within chromosomal elements as a result of paracentric inversions. Genome sequence comparison would seem the best method to estimate rates of chromosomal evolution, but the high-quality reference genomes required for this endeavor are still scanty. Here, we have obtained dense physical maps for Muller elements A, C, and E of Drosophila subobscura, a species with an extensively studied rich and adaptive chromosomal polymorphism. These maps are based on 462 markers: 115, 236, and 111 markers for elements A, C, and E, respectively. The availability of these dense maps will facilitate genome assembly and will thus greatly contribute to obtaining a good reference genome, which is a required step for D. subobscura to attain the model species status. The comparative analysis of these physical maps and those obtained from the D. pseudoobscura and D. melanogaster genomes allowed us to infer the number of fixed inversions and chromosomal evolutionary rates for each pairwise comparison. For all three elements, rates inferred from the more closely related species were higher than those inferred from the more distantly related species, which together with results of relative-rate tests point to an acceleration in the D. subobscura lineage at least for elements A and E.

Multiple and diverse structural changes affect the breakpoint regions of polymorphic inversions across the Drosophila genus

Chromosomal polymorphism is widespread in the Drosophila genus, with extensive evidence supporting its adaptive character in diverse species. Moreover, inversions are the major contributors to the genus chromosomal evolution. The molecular characterization of a reduced number of polymorphic inversion breakpoints in Drosophila melanogaster and Drosophila subobscura supports that their inversions would have mostly originated through a mechanism that generates duplications —staggered double-strand breaks— and has thus the potential to contribute to their adaptive character. There is also evidence for inversion breakpoint reuse at different time scales. Here, we have characterized the breakpoints of two inversions of D. subobscura —O₄ and O₈— involved in complex arrangements that are frequent in the warm parts of the species distribution area. The duplications detected at their breakpoints are consistent with their origin through the staggered-break mechanism, which further supports it as the prevalent mechanism in D. subobscura. The comparative analysis of inversions breakpoint regions across the Drosophila genus has revealed several genes affected by multiple disruptions due not only to inversions but also to single-gene transpositions and duplications.

The origin of chromosomal inversions as a source of segmental duplications in the Sophophora subgenus of Drosophila

Chromosomal inversions can contribute to the adaptation of organisms to their environment by capturing particular advantageous allelic combinations of a set of genes included in the inverted fragment and also by advantageous functional changes due to the inversion process itself that might affect not only the expression of flanking genes but also their dose and structure. Of the two mechanisms originating inversions —ectopic recombination, and staggered double-strand breaks and subsequent repair— only the latter confers the inversion the potential to have dosage effects and/or to generate advantageous chimeric genes. In Drosophila subobscura, there is ample evidence for the adaptive character of its chromosomal polymorphism, with an important contribution of some warm-climate arrangements such as E_1+2+9+12. Here, we have characterized the breakpoints of inversion E₁₂ and established that it originated through the staggered-break mechanism like four of the five inversions of D. subobscura previously studied. This mechanism that also predominates in the D. melanogaster lineage might be prevalent in the Sophophora subgenus and contribute to the adaptive character of the polymorphic and fixed inversions of its species. Finally, we have shown that the D. subobscura inversion breakpoint regions have generally been disrupted by additional structural changes occurred at different time scales.

The Evolution of Genome Structure by Natural and Sexual Selection

Progress on understanding how genome structure evolves is accelerating with the arrival of new genomic, comparative, and theoretical approaches. This article reviews progress in understanding how chromosome inversions and sex chromosomes evolve, and how their evolution affects species' ecology. Analyses of clines in inversion frequencies in flies and mosquitoes imply strong local adaptation, and roles for both over- and under dominant selection. Those results are consistent with the hypothesis that inversions become established when they capture locally adapted alleles. Inversions can carry alleles that are beneficial to closely related species, causing them to introgress following hybridization. Models show that this "adaptive cassette" scenario can trigger large range expansions, as recently happened in malaria mosquitoes. Sex chromosomes are the most rapidly evolving genome regions of some taxa. Sexually antagonistic selection may be the key force driving transitions of sex determination between different pairs of chromosomes and between XY and ZW systems. Fusions between sex-chromosomes and autosomes most often involve the Y chromosome, a pattern that can be explained if fusions are mildly deleterious and fix by drift. Sexually antagonistic selection is one of several hypotheses to explain the recent discovery that the sex determination system has strong effects on the adult sex ratios of tetrapods. The emerging view of how genome structure evolves invokes a much richer constellation of forces than was envisioned during the Golden Age of research on Drosophila karyotypes.

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.

Inversions and adaptation to the plant toxin ouabain shape DNA sequence variation within and between chromosomal inversions of Drosophila subobscura

Adaptation is defined as an evolutionary process allowing organisms to succeed in certain habitats or conditions. Chromosomal inversions have the potential to be key in the adaptation processes, since they can contribute to the maintenance of favoured combinations of adaptive alleles through reduced recombination between individuals carrying different inversions. We have analysed six genes (Pif1A, Abi, Sqd, Yrt, Atpα and Fmr1), located inside and outside three inversions of the O chromosome in European populations of Drosophila subobscura. Genetic differentiation was significant between inversions despite extensive recombination inside inverted regions, irrespective of gene distance to the inversion breakpoints. Surprisingly, the highest level of genetic differentiation between arrangements was found for the Atpα gene, which is located outside the O1 and O7 inversions. Two derived unrelated arrangements (O3+4+1 and O3+4+7) are nearly fixed for several amino acid substitutions at the Atpα gene that have been described to confer resistance in other species to the cardenolide ouabain, a plant toxin capable of blocking ATPases. Similarities in the Atpα variants, conferring ouabain resistance in both arrangements, may be the result of convergent substitution and be favoured in response to selective pressures presumably related to the presence of plants containing ouabain in the geographic locations where both inversions are present.

Introduced Drosophila subobscura populations perform better than native populations during an oviposition choice task due to increased fecundity but similar learning ability

The success of invasive species is tightly linked to their fitness in a putatively novel environment. While quantitative components of fitness have been studied extensively in the context of invasive species, fewer studies have looked at qualitative components of fitness, such as behavioral plasticity, and their interaction with quantitative components, despite intuitive benefits over the course of an invasion. In particular, learning is a form of behavioral plasticity that makes it possible to finely tune behavior according to environmental conditions. Learning can be crucial for survival and reproduction of introduced organisms in novel areas, for example, for detecting new predators, or finding mates or oviposition sites. Here we explored how oviposition performance evolved in relation to both fecundity and learning during an invasion, using native and introduced Drosophila subobscura populations performing an ecologically relevant task. Our results indicated that, under comparable conditions, invasive populations performed better during our oviposition task than did native populations. This was because invasive populations had higher fecundity, together with similar cognitive performance when compared to native populations, and that there was no interaction between learning and fecundity. Unexpectedly, our study did not reveal an allocation trade‐off (i.e., a negative relationship) between learning and fecundity. On the contrary, the pattern we observed was more consistent with an acquisition trade‐off, meaning that fecundity could be limited by availability of resources, unlike cognitive ability. This pattern might be the consequence of escaping natural enemies and/or competitors during the introduction. The apparent lack of evolution of learning may indicate that the introduced population did not face novel cognitive challenges in the new environment (i.e., cognitive “pre‐adaptation”). Alternatively, the evolution of learning may have been transient and therefore not detected.

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.

Revisiting classic clines in Drosophila melanogaster in the age of genomics

Adaptation to spatially varying environments has been studied for decades, but advances in sequencing technology are now enabling researchers to investigate the landscape of genetic variation underlying this adaptation genome wide. In this review we highlight some of the decades-long research on local adaptation in Drosophila melanogaster from well-studied clines in North America and Australia. We explore the evidence for parallel adaptation and identify commonalities in the genes responding to clinal selection across continents as well as discussing instances where patterns differ among clines. We also investigate recent studies utilizing whole-genome data to identify clines in D. melanogaster and several other systems. Although connecting segregating genomic variation to variation in phenotypes and fitness remains challenging, clinal genomics is poised to increase our understanding of local adaptation and the selective pressures that drive the extensive phenotypic diversity observed in nature.

The devil is in the details: genetic variation in introduced populations and its contributions to invasion

The influence of genetic variation on invasion success has captivated researchers since the start of the field of invasion genetics 50 years ago. We review the history of work on this question and conclude that genetic variation-as surveyed with molecular markers-appears to shape invasion rarely. Instead, there is a significant disconnect between marker assays and ecologically relevant genetic variation in introductions. We argue that the potential for adaptation to facilitate invasion will be shaped by the details of genotypes affecting phenotypes, and we highlight three areas in which we see opportunities to make powerful new insights. (i) The genetic architecture of adaptive variation. Traits shaped by large-effect alleles may be strongly impacted by founder events yet more likely to respond to selection when genetic drift is strong. Large-effect loci may be especially relevant for traits involved in biotic interactions. (ii) Cryptic genetic variation exposed during invasion. Introductions have strong potential to uncover masked variation due to alterations in genetic and ecological environments. (iii) Genetic interactions during admixture of multiple source populations. As divergence among sources increases, positive followed by increasingly negative effects of admixture should be expected. Although generally hypothesized to be beneficial during invasion, admixture is most often reported among sources of intermediate divergence, supporting the possibility that incompatibilities among divergent source populations might be limiting their introgression. Finally, we note that these details of invasion genetics can be coupled with comparative demographic analyses to link genetic changes to the evolution of invasiveness itself.

Inferring the demographic history of Drosophila subobscura from nucleotide variation at regions not affected by chromosomal inversions

Drosophila subobscura presents a rich and complex chromosomal inversion polymorphism. It can thus be considered a model system (i) to study the mechanisms originating inversions and how inversions affect the levels and patterns of variation in the inverted regions and (ii) to study adaptation at both the single-gene and chromosomal inversion levels. It is therefore important to infer its demographic history as previous information indicated that its nucleotide variation is not at mutation-drift equilibrium. For that purpose, we sequenced 16 noncoding regions distributed across those parts of the J chromosome not affected by inversions in the studied population and possibly either by other selective events. The pattern of variation detected in these 16 regions is similar to that previously reported within different chromosomal arrangements, suggesting that the latter results would, thus, mainly reflect recent demographic events rather than the partial selective sweep imposed by the origin and frequency increase of inversions. Among the simple demographic models considered in our Approximate Bayesian Computation analysis of variation at the 16 regions, the model best supported by the data implies a population size expansion soon after the penultimate glacial period. This model constitutes a better null model, and it is therefore an important resource for subsequent studies aiming among others to uncover selective events across the species genome. Our results also highlight the importance of introducing the possibility of multiple hits in the coalescent simulations with an outgroup.

Extreme cost of rivalry in a monandrous species: male-male interactions result in failure to acquire mates and reduced longevity

Mating system variation is profound in animals. In insects, female willingness to remate varies from mating with hundreds of males (extreme polyandry) to never remating (monandry). This variation in female behaviour is predicted to affect the pattern of selection on males, with intense pre-copulatory sexual selection under monandry compared to a mix of pre- and post-copulatory forces affecting fitness under polyandry. We tested the hypothesis that differences in female mating biology would be reflected in different costs of pre-copulatory competition between males. We observed that exposure to rival males early in life was highly costly for males of a monandrous species, but had lower costs in the polyandrous species. Males from the monandrous species housed with competitors showed reduced ability to obtain a mate and decreased longevity. These effects were specific to exposure to rivals compared with other types of social interactions (heterospecific male and mated female) and were either absent or weaker in males of the polyandrous species. We conclude that males in monandrous species suffer severe physiological costs from interactions with rivals and note the significance of male-male interactions as a source of stress in laboratory culture.