Evolution of Chilean colonizing populations of <Emphasis Type="Italic">Drosophila subobscura</Emphasis>: lethal genes and chromosomal arrangements

Knowledge of the frequency, distribution, and fate of lethal genes in chromosomal inversions helps to illuminate the evolution of recently founded populations. We analyze the relationship between lethal genes and inversions in two colonizing populations of D. subobscura in Chile. In the ancestral Palearctic populations of this species, lethal genes seem distributed at random on chromosomes. But in colonizing American populations, some lethal genes are associated with specific chromosomal arrangements. Some of these associated lethals were detected only during the first stages of the colonization (O ₃₊₄₊₂ ), and never thereafter, whereas others have persisted (O ₃₊₄₊₇ and O₅). However, most lethal genes in American populations have been observed only once: they have arisen by novel mutation and soon disappear. Finally, recombination between different inversions has been observed in America. However, the persistence of lethal genes associated with the heterotic inversions O ₃₊₄₊₇ and O₅ could indicate that recombination inside these inversions is rare. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Allozyme divergence between different gene arrangements in Drosophila subobscura

Allozyme and chromosomal polymorphism were analysed simultaneously in three population samples of Drosophila subobscura from Scotland, Germany and Italy. The study revealed significant genic divergence between the various chromosomal gene arrangements for the loci Acph, Lap and Pep, which are located inside inversions. No difference was found with the loci Me and Odh, which are both on the same chromosome but proximal to the structurally polymorphic region. Samples of a given gene arrangement taken from the different localities are genetically uniform. The data indicate that the observed non-random associations between enzyme genes and chromosomal inversions are the result of inhibited recombination in heterokaryotypes and gradual differentiation of the separated gene pools, rather than selection and coadaptation.     

Genetic coadaptation in the chromosomal polymorphism of Drosophila subobscura II. Changes of gametic disequilibrium in experimental populations

Experimental populations were examined for temporal changes of gametic disequilibria between allozyme loci (Lap and Pept-1) and gene arrangements of the O chromosome of Drosophila subobscura (O _st and O ₃₊₄₊₇) under several environmental conditions. In the foundation of the experimental populations a genetic perturbation was carried out in order to test the relevance of the current hypotheses used to explain the allozyme-inversion associations observed in natural populations. Differential changes of gametic disequilibria were detected over generations under the different environmental conditions. Mere mechanical or stochastic factors cannot explain the results and natural selection is probably the major agent generating the detected gametic associations. The observations are interpreted as a proof of coadaptation of D. subobscura inversions.     

An experimental study of evolution in progress: Clines for quantitative traits in colonizing and Palearctic populations ofDrosophila

Summary Drosophila subobscura has recently colonized the American continent and is an excellent model for studying evolution in action. Previous analyses have shown that these colonizing populations have significant latitudinal clines for the frequencies of some chromosomal arrangements that parallel those clines found in the Old World. These results strongly suggest that this polymorphism is adaptive. In the present study, significant latitudinal clines for continuous morphometric variables (flies are larger in the north) have been detected in Old World populations ofD. subobscura. The adaptive nature of these clines is reinforced by the fact that parallel latitudinal clines for body size have also been detected inDrosophila obscura, a closely related sympatric species, as well as previously in otherDrosophila. On the other hand, no significant latitudinal clines for continuous morphometric traits, not even when using an overall size index, have been detected in colonizing populations ofD. subobscura. This is a rather surprising result given the number of generations that have elapsed since the species was detected in America and given that significant clines in chromosomal inversions are already established. Thus, the adaptive response of quantitative morphometric variables is not as rapid as that found for chromosomal inversions. Nevertheless, canonical correlation analysis suggests that significant latitudinal clines for body size might soon be detected in the American continent. The results obtained here are included in a projected time series with the aim of documenting size evolution in action.     

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 O₃₊₄ 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.     

THE DETECTION OF GAMETIC DISEQUILIBRIUM BETWEEN ALLOZYME LOCI IN NATURAL POPULATIONS OF DROSOPHILA

The capacity of the usual tests (chi-square and related tests) to detect gametic disequilibrium between allozyme loci in natural populations of Drosophila has been investigated. We analyzed a large collection of previously reported gametic samples from natural populations involving a variety of loosely linked allozyme loci located along the O chromosome of Drosophila subobscura and the second chromosome of D. melanogaster. It is found that the statistical power of the individual tests to detect the sample disequilibria between allozyme loci is remarkably low, being the average (over pairs of loci) of power estimates close to 0.20 in both species. Moreover, the average minimum disequilibrium (D‘_min) that would be required to reject (90% probability) the hypothesis of gametic equilibrium is higher than 0.50 given the observed degree of polymorphism and sample sizes used. This means that statistically significant associations between allozyme loci would rarely be detected by single-sample tests even when much disequilibrium is present in natural populations of Drosophila. However, an alternative approach based on the analysis of disquilibrium for large sets of gametic samples, combining probabilities from single independent tests and assessing significance by a bootstrap procedure, reveals that most of the locus pairs within segment I and II of the O chromosome of D. subobscura and left arm of the second chromosome of D. melanogaster present significant nonrandom associations. Within these chromosomal sections, the observed average absolute value of disquilibrium (D‘) between loci is around 0.25 (under the more conservative estimation). Also, a positive relationship between the magnitude of disequilibrium and linkage was detected. These findings suggest that weak or moderate values of disequilibrium between loosely linked allozyme loci are more frequent in natural populations of Drosophila than is currently believed.     

Geographic pattern of allozyme and inversion polymorphism on chromosome O of Drosophila subobscura and its evolutionary origin

The chromosome O of Drosophila subobscura was studied with respect to genetic variability at three enzyme loci (Odh, Me, and Lap-4) and with respect to inversion polymorphism. Population samples were taken from seven localities along a north-south gradient from Sweden and Scotland to Tunisia.The chromosomal analysis revealed clinal frequency changes for gene arrangements from north to south. With the enzyme loci Odh and Me allele frequencies are similar throughout the distribution range. Both loci are located outside the common inversion complex O₃₊₄. On the other hand, frequency changes parallel to those of the gene arrangements were observed for the alleles of the Lap-locus. Nonrandom associations between Lap-alleles and the superimposed gene arrangements O_ST, O₃₊₄, O₃₊₄₊₈, and O₃₊₄₊₂₃ were found. These gene arrangements differ from each other with respect to allele frequencies at the Lap-locus but for a given gene arrangement the relative frequencies of Lap-alleles remain relatively constant along the north-south gradient. Thus allele frequencies at the Lap-locus can be predicted from inversion frequencies.These observations can be interpreted in such a way that the pattern of allozyme variation within gene arrangements is due to founder effects caused by the unique origin of inversions. The gene blocks in the different inversions seem to represent more or less separated gene pools. In polymorphic populations the coexistence of genetically differentiated inversions presumably gives rise to heterotic interaction.     

Inversions fail to account for allozyme clines

Allozyme and inversion data from natural populations of Drosophila melanogaster from the eastern United States were analyzed to determine whether the clines at allozyme loci are due to nonrandom associations with common cosmopolitan inversions. All inversions show strong clines. Clines were large and significant for half of the eight allozyme loci. An analysis of the contribution of inversions to clines of allozyme genes revealed three outcomes: the inversion cline (1) enhanced the allozyme cline, but was only partly responsible, (2) reduced the allozyme cline, and (3) had no effect. The allozyme clines were mainly determined by the pattern of allele frequencies within the chromosomal arrangements. Consequently, it was concluded that allozyme clines would exist in the absence of inversion clines.     

Tracking the origin of an invasive species: Drosophila subobscura in Argentina

Biological invasions are excellent opportunities to study the evolutionary forces leading to the adaptation of a species to a new habitat. Knowledge of the introduction history of colonizing species helps tracking colonizing routes and assists in defining management strategies for invasive species. The Palearctic species Drosophila subobscura is a good model organism for tracking colonizations since it was detected in Chile and western North America three decades ago and later on in the Atlantic coast of Argentina. To unravel the origin of the Argentinean colonizers two populations have been analysed with several genetic markers. Chromosomal arrangements and microsatellite alleles found in Argentina are almost similar to those observed in Chile and USA. The lethal allelism test demonstrates that the lethal gene associated with the O₅ inversions in Argentina is identical to that found in Chile and USA, strongly supporting the hypothesis that all the American colonizing populations originated from the same colonization event. A secondary bottleneck is detected in the Argentinean populations and the genetic markers suggest that these populations originated from the invasion of 80–150 founding individuals from Chile.     

Climatic adaptation of chromosomal inversions in <Emphasis Type="Italic">Drosophila subobscura</Emphasis>

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.     

Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura

In fewer than two decades after invading the Americas, the fly Drosophila subobscura evolved latitudinal clines for chromosomal inversion frequencies and wing size that are parallel to the long‐standing ones in native Palearctic populations. By sharp contrast, wing shape clines also evolved in the New World, but the relationship with latitude was opposite to that in the Old World. Previous work has suggested that wing trait differences among individuals are partially due to the association between chromosomal inversions and particular alleles which influence the trait under consideration. Furthermore, it is well documented that a few number of effective individuals founded the New World populations, which might have modified the biometrical effect of inversions on quantitative traits. Here we evaluate the relative contribution of chromosomal inversion clines in shaping the parallel clines in wing size and contrasting clines in wing shape in native and colonizing populations of the species. Our results reveal that inversion‐size and inversion‐shape associations in native and colonizing (South America) populations are generally different, probably due to the bottleneck effect. Contingent, unpredictable evolution was suggested as an explanation for the different details involved in the otherwise parallel wing size clines between Old and New World populations of D. subobscura. We challenge this assertion and conclude that contrasting wing shape clines came out as a correlated response of inversion clines that might have been predicted considering the genetic background of colonizers.     

DNA variation at the rp49 gene region in Drosophila madeirensis and D. subobscura from Madeira: inferences about the origin of an insular endemic species

An ~1.6 ‐ kb fragment spanning the rp49 gene was sequenced in 16 lines of Drosophila subobscura from Madeira and in 22 lines of the endemic species D. madeirensis. Nucleotide diversity in D. subobscura from Madeira (π=0.0081) was similar to that in lines from Spain carrying the O₃₊₄ chromosomal arrangement (π=0.0080). No significant genetic differentiation was detected between insular and continental O₃₊₄ lines of D. subobscura. These results are compatible both with a rather recent and massive colonization, and with multiple colonization events from the continent. Nucleotide diversity in D. madeirensis (π=0.0076) was similar to that in D. subobscura, which deviates from the expectation, under strict neutrality, of a lower level of variation in an insular species with a small population size. The observed numbers of shared polymorphisms and of fixed differences between D. madeirensis and D. subobscura are compatible with the isolation model of speciation, where shared polymorphisms are due to common ancestry.     

Analysis of F ST outliers at allozyme loci in Pacific salmon: implications for natural selection

Natural selection has been invoked to explain the observed geographic distribution of allozyme allele frequencies for a number of teleost species. The effects of selection on allozyme loci in three species of Pacific salmon were tested. A simulation-based approach to estimate the null distribution of population differentiation (F _ST) and test for F _ST outliers was used. This approach showed that a majority of allozyme loci conform to neutral expectations predicted by the simulation model, with relatively few F _ST outliers found. No consistent F _ST outlier loci were found across species. Analysis of population sub-groups based on geography and genetic identity reduced the number of outlier loci for some species, indicating that large geographic groups may include genetically divergent populations and/or that there is geographic heterogeneity in selection pressure upon allozyme loci. Two outlier allozyme loci found in this analysis, lactate dehydrogenase-B and malic enzyme, have been shown to be influenced by selection in other teleost species. This approach is also useful in identifying allozyme loci (or other genetic markers) that meet assumptions for population genetic study.     

The role of genic selection in the establishment of inversion polymorphism in Drosophila subobscura

A model for explaining the establishment of newly arisen inversions in natural populations, in which the inverted segment may be selected for if its load of deleterious mutations is smaller than the average load of the noninverted segment in the population, is tested for Drosophila subobscura. The results show that for new inversions, originally with no deleterious alleles, the expected cumulative distribution of inversion lengths fits fairly well with the observed one. Therefore, genic selection may be an important cause of the establishment of newly arisen inversions in natural populations of D. subobscura. The applicability of the model to the maintenance of the inversion polymorphism present in this species; is discussed.     

Linkage Disequilibrium in Natural Populations of DROSOPHILA MELANOGASTER. Seasonal Variation

Linkage disequilibrium among ten polymorphic allozyme loci and polymorphic inversions on chromosomes 2 and 3 in a natural population of Drosophila melanogaster was examined early and late in the annual season. Similar to previous studies, little linkage disequilibrium was observed among allozymes. The two significant cases that were observed in the first sample behaved in a contradictory way. One declined much more rapidly than expected due simply to recombination; the other declined slowly as expected. There was little change in allozyme or inversion frequencies during the season.     

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.     

O-chromosome lethal frequencies in Serbian and Montenegrin <Emphasis Type="Italic">Drosophila subobscura</Emphasis> populations

Lethal chromosomal frequencies were obtained from three Drosophila subobscura samples from the Mt. Avala (Serbia) population in September 2003 (0.218), June 2004 (0.204) and September 2004 (0.250). These values and those from other Balkan populations studied previously (Petnica, Kamariste, Zanjic and Djerdap) were used to analyze the possible effect of population, year, month and altitude above sea level on lethal chromosomal frequencies. According to ANOVAS no effect were observed. Furthermore, the lethal frequencies of the Balkan populations did not vary according to latitude. This is probably due to the relative proximity and high gene flow between these populations. From a joint study of all the Palearctic D. subobscura populations so far analyzed, it can be deduced that the Balkan populations are located in the central area of the species distribution. Finally, it seems that lethal chromosomal frequencies are a consequence of the genetic structure of the populations.     

STRUCTURE AND POPULATION GENETICS OF THE BREAKPOINTS OF A POLYMORPHIC INVERSION IN DROSOPHILA SUBOBSCURA

Drosophila subobscura is a paleartic species of the obscura group with a rich chromosomal polymorphism. To further our understanding on the origin of inversions and on how they regain variation, we have identified and sequenced the two breakpoints of a polymorphic inversion of D. subobscura—inversion 3 of the O chromosome—in a population sample. The breakpoints could be identified as two rather short fragments (～300 bp and 60 bp long) with no similarity to any known transposable element family or repetitive sequence. The presence of the ～300‐bp fragment at the two breakpoints of inverted chromosomes implies its duplication, an indication of the inversion origin via staggered double‐strand breaks. Present results and previous findings support that the mode of origin of inversions is neither related to the inversion age nor species‐group specific. The breakpoint regions do not consistently exhibit the lower level of variation within and stronger genetic differentiation between arrangements than more internal regions that would be expected, even in moderately small inversions, if gene conversion were greatly restricted at inversion breakpoints. Comparison of the proximal breakpoint region in species of the obscura group shows that this breakpoint lies in a small high‐turnover fragment within a long collinear region (～300 kb).     

A linkage map of sugar beet (Beta vulgaris L.)

Summary We have established a first linkage map for beets based on RFLP, isozyme and morphological markers. The population studied consisted of 96 F₂ individuals derived from an intraspecific cross. As was expected for outbreeding species, a relatively high degree of polymorphism was found within sugar beet; 47% of the DNA markers were polymorphic for the chosen population. The map consists of 115 independent chromosomal loci designated by 108 genomic DNA probes, 6 isozyme and one morphological marker. The loci cover 789 cM with an average spacing of 6.9 cM. They are dispersed over nine linkage groups corresponding to the haploid chromosome number of Beta species. Eighteen markers (15.4%) showed distorted segregation which, in most instances, can be explained by gametic selection of linked lethal loci. The application of the linkage map in sugar beet breeding is discussed.     

Ultrastructure of regions containing homologous loci in polytene chromosomes of Drosophila melanogaster and Drosophila subobscura

We have used a new approach involving in situ hybridisation and electron microscopy to establish ultrastructural homologies between polytene chromosome regions of Drosophila melanogaster and Drosophila subobscura. Twelve probes were chosen to cover all the chromosomal elements: the myospheroid gene, the collagen type IV gene, the collagen-like gene, the w26 homeobox gene, the β3 tubulin gene, the kinesin heavy chain gene, the tryptophan hydrolase gene, the Hsp82, Hsp22–26 and Hsp23–28, Hsp68, Hsp70 genes and the β unit of the F₀–F₁ ATPase gene. Most of these loci were previously undescribed in D. subobscura and imprecisely located in D. melanogaster. We have demonstrated here, by an ultrastructural analysis of each chromosomal region, that homologous genetic loci tend to show a similar ultrastructure in the two species. With a few exceptions, the structural homology extends to the chromosomal regions surrounding the loci. In some cases, however, no structurally recognisable homology can be seen either in the locus or in its flanking regions. Received: 15 December 1996; in revised form: 15 October 1997 / Accepted: 28 January 1998