Genetic loads and coadaptivity of chromosomal inversions I. Inversion polymorphism and genetic load on chromosome O in a Drosophila subobscura population from Petnica

Inversion polymorphism on chromosome O and polymorphism for the viability of determining genes have been studied in a natural population of Drosophila subobscura from Petnica (Serbia). The range of inversion polymorphism and the abundance of particular gene arrangements in the study population agree with a general pattern of inversion polymorphism of D. subobscura in Europe. The data obtained on the amount of genetic loads show that the D. subobscura population from Petnica displays a moderate degree of that polymorphism, compared to the other studied populations of these species. Therefore, the D. subobscura population from Petnica could be tentatively classified as an ecologically central population. Examination association of chromosomal, thus, inversion polymorphism with gene polymorphism, in the form of genetic loads show that differences exist in the mean viability among certain gene arrangements. The distribution of deleterious genes among chromosome O gene arrangements were non-random.     

Viabilities of originally natural O-chromosomal inversion homo- and heterokaryotypes in Drosophila subobscura

The influence of epistatic interactions of lethal and non-lethal genes upon viability of Drosophila inversion karyotypes is poorly known. In this paper we present comparative results for viabilities of 21 originally natural O-inversion homo- and 38 heterokaryotypes in a D. subobscura population. We observed strong heterotic effect in viability of O-lethal heterozygotes irrespective of different inversion backgrounds, which indicates a mechanism for protection of a great number of lethal genes. In O-non-lethal heterozygotes the heterotic effect in viability was absent irrespective of different inversion backgrounds. Our results showed a great extent of genetic load and high abundance of O-chromosomal arrangements in the population analyzed. It belongs to the set of central European populations. An epistasis of lethal genes present in O-inversion hetero- and not present in O-inversion homokaryotypes of moderate or low frequencies could be good example for co-adaptation of chromosomal inversions with regard to the genetic load. This represents a more efficient mechanism for limitation of genetic load than alternative mechanisms for protection of lethals. Except for lethal genes, possible epistatic interactions of mildly deleterious (subvital) genes, could also be responsible for limiting the extent of genetic load in natural D. subobscura populations. We demonstrated a non-random distribution of several combinations of viability classes among different O-inversion homo- and heterokaryotypes. As a consequence of that, the viabilities of the O-inversion homokaryotypes compared to heterokaryotypes were significantly higher at low frequencies than in moderate or high frequencies. This shows frequency-dependence as a mechanism of balancing selection for protection of O-chromosomal inversions in natural D. subobscura populations. In addition, the viabilities of the O-inversion homokaryotypes of lower frequency, compared to homokaryotypes of moderate or high frequency, were significantly higher. This again indicates the existence of supergene selection as another mechanism for protection of chromosomal inversions, as co-adapted complexes in natural D. subobscura populations.     

Molecular population genetics of the rp49 gene region in different chromosomal inversions of Drosophila subobscura.

Nucleotide variation at the ribosomal protein 49 (rp49) gene region has been studied in 75 lines of Drosophila subobscura belonging to four chromosomal arrangements (Ost, O3+4, O3+4+8, and O3+4+23). The location of the rp49 gene region within the inversion loop differs among heterokaryotypes: it is very close to one of the breakpoints in heterozygotes involving Ost chromosomes, while it is in a more central position in all other heterokaryotypes. The distribution of nucleotide polymorphism in the different arrangements is consistent with a monophyletic origin of the inversions. The data also provide evidence that gene conversion and possibly double crossover are involved in shuffling nucleotide variation among gene arrangements. The analyses reveal that the level of genetic exchange is higher when the region is located in a more central position of the inverted fragment than when it is close to the breakpoints. The pairwise difference distributions as well as the negative values of Tajima''s and Fu and Li''s statistics further support the hypothesis that nucleotide variation within chromosomal arrangements still reflects expansion after the origin of the inversions. Under the expansion model, we have estimated the time of origin of the studied inversions.     

Evidence of Extensive Genetic Exchange in the Rp49 Region among Polymorphic Chromosome Inversions in Drosophila Subobscura

Restriction map variation in 107 lines extracted from two natural populations of Drosophila subobscura was investigated with seven four-nucleotide-recognizing enzymes in a 1.6-kb region including the rp49 gene, that is located very close to the proximal breakpoint of inversion O3. Fourteen restriction site and 8 length polymorphisms, resulting in 73 haplotypes, have been identified. Estimated heterozygosity per nucleotide, pi = 0.0045, is comparable to the average nucleotide variation observed in Drosophila melanogaster. Because of the location of the rp49 region in D. subobscura, variation in three different gene arrangements-Ost, O3 + 4 and O3 + 4 + 8-has been compared. Out of 14 restriction site polymorphisms, 3 are shared by Ost, O3 + 4 and O3 + 4 + 8, and 3 additional ones are shared by Ost and O3 + 4, evidencing extensive genetic exchange among these polymorphic inversions. In agreement with previous data, the higher level of variation of O3 + 4 (as measured by haplotype diversity and nucleotide variation) suggests that O3 + 4 may be ancestral in relationship to extant gene arrangements.     

Lethal genes in O5, chromosomes of Drosophila subobscura from Europe and America

In populations of D. subobscura , a species that is know for its high chromosomal polymorphism, the O₅ inversion has a rather erratic frequency distribution in the Palearctic region. An O₅ lethal chromosomal line obtained from a Balkan population near Zanjic (South Adriatic, Montenegro, Yugoslavia) was tested for lethal allelism with other O₅ lethal chromosomal lines derived from American (USA and Chile) colonizing populations, and from the French population of Taulé. No allelism was found between the Balkan lethal gene and those from America and France. Thus, the lethal genes of the O₅ inversions are not of the same origin and it is most probable that the American colonizations did not start from the Zanjic population. The general difference in the chromosomal inversion polymorphism corroborates this conclusion. The cytological analysis confirms the assumption that all O₅ chromosomes studied are identical with respect to breakage points.     

Tracking the origin of the American colonization by Drosophila subobscura: genetic comparison between Eastern and Western Mediterranean populations

Several pieces of evidence indicate a Mediterranean origin of the colonization of America by Drosophila subobscura . To ascertain whether the origin was from the Eastern or the Western Mediterranean region, samples from Barcelona (Spain) and Mt Parnes (Greece) were collected and O chromosomal inversion polymorphism and lethal genes were analysed. The frequencies of lethal chromosomes were 0.244 ± 0.039 in Barcelona and 0.336 ± 0.043 in Mt Parnes, consistent with the expectations for large populations located in the central area of the species distribution. Lethal genes seem to be distributed at random along the O chromosome in both populations. The intra-populational allelism frequencies of Barcelona and Mt Parnes were 0.016 ± 0.007 and 0.012 ± 0.005 respectively. Thus, the estimates of the effective population size were high in both populations (between 6964 and 13 004 in Barcelona and 11 874 to 26 828 in Mt Parnes). The cases of allelism in Mt Parnes were observed only between individual lethal genes, but in Barcelona some concatenated clusters of allelism were detected. This pattern of allelism can be explained by synthetic lethality, hybrid dysgenesis, the induction of recurrent lethal mutations by different factors or an effect of microdifferentiation in subpopulations. In both populations, a reduction in fitness in the heterozygotes for lethal genes has been detected. Furthermore, the estimates of the migration coefficient (between 0.0085 and 0.0120 in Barcelona, and 0.0057 and 0.0087 in Mt Parnes) confirm the existence of gene flow between Palearctic populations of D. subobscura . Our lethal genes and chromosomal inversion results are consistent with a Mediterranean origin of the colonization of America by D. subobscura , but are inconclusive with regard to the identification of the population from which these colonizers came.     

Colonization of America by Drosophila subobscura: association between Odh gene haplotypes, lethal genes and chromosomal arrangements

The colonization of America by Drosophila subobscura has been a unique exper iment in nature that has allowed us to explore the effects of evolution on a continental scale. To analyze this evolutionary event, nucleotide sequences of the Odh (Octanol dehydrogenase) gene were obtained for 43 lethal chromosomal lines from colonizing populations of North America and 5 from South America, in addition to 5 chromosomal lines from Europe with different viabilities and 2 from laboratory marker stocks. Since 10 different Odh haplotypes were found in America, the minimum number of colonizers would be 5 (or 3 mated females). Only one Odh haplotype was found in American O(5) inversions confirming that only one copy of this inversion was included among the sample of colonizers. The same Odh haplotypes were detected in association with the same chromosomal arrangements and with identical lethal genes in both North and South America indicating that exactly the same chromosome types reached both hemispheres. These observations indicate that the two continental colonizations are not independent. They are derived from the same colonization event. The population from which the colonization started should contain the O(5) inversion, a non-negligible frequency of the O(3+4+7) arrangement and all other arrangements found in America. So far the only populations that fulfill all these requirements are those from Greece, indicating that these populations can be considered good candidates as a starting point for an in depth analysis of the origin of the American colonization by D. subobscura.     

Colonization of the Americas by Drosophila Subobscura: Lethal-Gene Allelism and Association with Chromosomal Arrangements

Drosophila subobscura is a Palearctic species that has recently colonized the Americas. It was first found in 1978 in Puerto Montt, Chile, and in 1982 in Port Townsend, WA. The colonization and rapid expansion of the species in western South and North America provides distinctive opportunities for investigating the process of evolution in action. The inversion polymorphism in the O chromosome from populations of central California and northern Washington, separated by 1300 km, corresponds to a previously observed latitudinal cline, also observed in Europe. Recessive lethal genes are not randomly distributed among the chromosomal arrangements. The incidence of lethal allelism is high, yielding unrealistically low estimates of the effective size of these populations (on the order of 1000 individuals). The high incidence of lethal allelism is likely to be a consequence of the low number of the American colonizers (on the order of 10-100 individuals), but the persistence of the allelism over several years suggests that some lethal-carrying chromosomes may be heterotic owing to shared associations between lethal and other genes.     

Conditions for Protected Inversion Polymorphism under Supergene Selection

Conditions for protected inversion polymorphism under the operation of both karyotype and supergene selection in a viability model have been analytically determined. When supergene selection (the effect of recombination in homokaryotypes lowering the mean fitness of their offspring) is acting on gene arrangements and there is no karyotype selection, it is demonstrated that a polymorphic stable equilibrium is reached by the population, which is a function of only the recombination effects in homokaryotypes. Under both supergene and karyotype selection the degree of dominance (h) of karyotype selection is critical to produce a protected inversion polymorphism. In general, the opportunity for protected polymorphism increases as the degree of dominance decreases. For small s values, the conditions for protected polymorphism are r > 2sh and c > 2s(h - 1), where r and c are the average loss of viability for offspring of ST/ST and IN/IN homokaryotypes, respectively. These findings suggest that supergene selection may be an important balancing mechanism contributing to the maintenance of inversion polymorphism.     

Chromosomal inversion polymorphism leads to extensive genetic structure: a multilocus survey in Drosophila subobscura

The adaptive character of inversion polymorphism in Drosophila subobscura is well established. The O(ST) and O(3+4) chromosomal arrangements of this species differ by two overlapping inversions that arose independently on O(3) chromosomes. Nucleotide variation in eight gene regions distributed along inversion O(3) was analyzed in 14 O(ST) and 14 O(3+4) lines. Levels of variation within arrangements were quite similar along the inversion. In addition, we detected (i) extensive genetic differentiation between arrangements in all regions, regardless of their distance to the inversion breakpoints; (ii) strong association between nucleotide variants and chromosomal arrangements; and (iii) high levels of linkage disequilibrium in intralocus and also in interlocus comparisons, extending over distances as great as approximately 4 Mb. These results are not consistent with the higher genetic exchange between chromosomal arrangements expected in the central part of an inversion from double-crossover events. Hence, double crossovers were not produced or, alternatively, recombinant chromosomes were eliminated by natural selection to maintain coadapted gene complexes. If the strong genetic differentiation detected along O(3) extends to other inversions, nucleotide variation would be highly structured not only in D. subobscura, but also in the genome of other species with a rich chromosomal polymorphism.     

Genetic exchange versus genetic differentiation in a medium-sized inversion of Drosophila: the A2/Ast arrangements of Drosophila subobscura

Chromosomal inversion polymorphism affects nucleotide variation at loci associated with inversions. In Drosophila subobscura, a species with a rich chromosomal inversion polymorphism and the largest recombinational map so far reported in the Drosophila genus, extensive genetic structure of nucleotide variation was detected in the segment affected by the O(3) inversion, a moderately sized inversion at Muller's element E. Indeed, a strong genetic differentiation all over O(3) and no evidence of a higher genetic exchange in the center of the inversion than at breakpoints were detected. In order to ascertain, whether other polymorphic and differently sized inversions of D. subobscura also exhibited a strong genetic structure, nucleotide variation in 5 gene regions (P236, P275, P150, Sxl, and P125) located along the A(2) inversion was analyzed in A(st) and A(2) chromosomes of D. subobscura. A(2) is a medium-sized inversion at Muller's element A and forms a single inversion loop in heterokaryotypes. The lower level of variation in A(2) relative to A(st) and the significant excess of low-frequency variants at polymorphic sites indicate that nucleotide variation at A(2) is not at mutation-drift equilibrium. The closest region to an inversion breakpoint, P236, exhibits the highest level of genetic differentiation (F(ST)) and of linkage disequilibrium (LD) between arrangements and variants at nucleotide polymorphic sites. The remaining 4 regions show a higher level of genetic exchange between A(2) and A(st) chromosomes than P236, as revealed by F(ST) and LD estimates. However, significant genetic differentiation between the A(st) and A(2) arrangements was detected not only at P236 but also in the other 4 regions separated from the nearest breakpoint by 1.2-2.9 Mb. Therefore, the extent of genetic exchange between arrangements has not been high enough to homogenize nucleotide variation in the center of the A(2) inversion. A(2) can be considered a typical successful inversion of D. subobscura according to its relative length. Chromosomal inversion polymorphism of D. subobscura might thus cause the genome of this species to be highly structured and to harbor different gene pools that might contribute to maintain adaptations to particular environments.     

The effect of different concentrations of lead on inversion polymorphism in Drosophila subobscura

Drosophila subobscura is a wild Drosophila species that is spread over almost all of Europe. It possesses an uniquely rich inversion polymorphism on all five long chromosomes. This polymorphism is to a certain degree associated with the variation and dynamics of ecological factors in space and time. We analyzed the changes of inversion polymorphism components of Drosophila subobscura flies maintained on media with different concentrations of lead in laboratory conditions. The effects of lead on inversion polymorphism were observed by cytological analysis of gene arrangements on all of the five acrocentric chromosomes, as well as by cytological analysis of karyotypes on all of the four autosomes. The frequencies of particular gene arrangements on the four autosomes changed significantly in the samples maintained on medium not supplemented with lead. The frequencies of some gene arrangements on all of the five acrocentric chromosomes changed significantly in the flies maintained on media supplemented with lead. The length of exposure to different lead concentrations results in a significant change in the frequency of a few gene arrangements on two autosomes. However, the results show that different concentrations of lead, as well as the length of exposure, do not affect major parameters of inversion polymorphism. The results suggest that some gene arrangements could be linked with adaptive processes in evolving heavy metal resistance.     

Nucleotide variation at the yellow gene region is not reduced in Drosophila subobscura: a study in relation to chromosomal polymorphism

In contrast to Drosophila melanogaster and Drosophila simulans, the yellow (y) gene region of Drosophila subobscura is not located in a region with a strong reduction in recombination. In addition, this gene maps very close to the breakpoints of different inversions that segregate as polymorphic in natural populations of D. subobscura. Therefore, levels of variation at the y gene region in this species relative to those found in D. melanogaster and D. simulans may be affected not only by the change in the recombinational environment, but also by the presence of inversion polymorphism. To further investigate these aspects, an approximately 5.4-kb region of the A (=X) chromosome including the y gene was sequenced in 25 lines of D. subobscura and in the closely related species Drosophila madeirensis and Drosophila guanche. The D. subobscura lines studied differed in their A-chromosomal arrangements, A(st), A(2), and A(1). Unlike in D. melanogaster and D. simulans, levels of variation at the y gene region of D. subobscura are not reduced relative to those found at other genomic regions in the same species (rp49, Acp70A, and Acph-1). This result supports the effect of the change in the recombinational environment of a particular gene on the level of neutral variation. In addition, nucleotide variation is affected by chromosomal polymorphism. A strong genetic differentiation is detected between the A(1) arrangement and either A(st) or A(2), but not between A(st) and A(2). This result is consistent with the location of the y gene relative to the breakpoints of inversions A(1) and A(2). In addition, the pattern of nucleotide polymorphism in A(st)+A(2) and A(1) seems to point out that variation at the y gene region within these chromosomal classes is in the phase transient to equilibrium. The estimated ages of these arrangements assuming a star genealogy indicate that their origin cannot predate the D. madeirensis split. Therefore, the present results are consistent with a chromosomal phylogeny where Am(1), which is an arrangement present in D. madeirensis but absent in current populations of D. subobscura, would be the ancestral arrangement.     

Genomics of natural populations: Evolutionary forces that establish and maintain gene arrangements in Drosophila pseudoobscura

The evolution of complex traits in heterogeneous environments may shape the order of genes within chromosomes. Drosophila pseudoobscura has a rich gene arrangement polymorphism that allows one to test evolutionary genetic hypotheses about how chromosomal inversions are established in populations. D. pseudoobscura has >30 gene arrangements on a single chromosome that were generated through a series of overlapping inversion mutations with >10 inversions with appreciable frequencies and wide geographic distributions. This study analyses the genomic sequences of 54 strains of Drosophila pseudoobscura that carry one of six different chromosomal arrangements to test whether (i) genetic drift, (ii) hitchhiking with an adaptive allele, (iii) direct effects of inversions to create gene disruptions caused by breakpoints, or (iv) indirect effects of inversions in limiting the formation of recombinant gametes are responsible for the establishment of new gene arrangements. We found that the inversion events do not disrupt the structure of protein coding genes at the breakpoints. Population genetic analyses of 2,669 protein coding genes identified 277 outlier loci harbouring elevated frequencies of arrangement‐specific derived alleles. Significant linkage disequilibrium occurs among distant loci interspersed between regions with low levels of association indicating that distant allelic combinations are held together despite shared polymorphism among arrangements. Outlier genes showing evidence of genetic differentiation between arrangements are enriched for sensory perception and detoxification genes. The data presented here support the indirect effect of inversion hypothesis where chromosomal inversions are favoured because they maintain linked associations among multilocus allelic combinations among different arrangements.     

Meiotic transmission of Drosophila pseudoobscura chromosomal arrangements

Drosophila pseudoobscura harbors a rich gene arrangement polymorphism on the third chromosome generated by a series of overlapping paracentric inversions. The arrangements suppress recombination in heterokaryotypic individuals, which allows for the selective maintenance of coadapted gene complexes. Previous mapping experiments used to determine the degree to which recombination is suppressed in gene arrangement heterozygotes produced non-recombinant progeny in non-Mendelian ratios. The deviations from Mendelian expectations could be the result of viability differences between wild and mutant chromosomes, meiotic drive because of achiasmate pairing of homologues in heterokaryotypic females during meiosis, or a combination of both mechanisms. The possibility that the frequencies of the chromosomal arrangements in natural populations are affected by mechanisms other than adaptive selection led us to consider these hypotheses. We performed reciprocal crosses involving both heterozygous males and females to determine if the frequency of the non-recombinant progeny deviates significantly from Mendelian expectations and if the frequencies deviate between reciprocal crosses. We failed to observe non-Mendelian ratios in multiple crosses, and the frequency of the non-recombinant classes differed in only one of five pairs of reciprocal crosses despite sufficient power to detect these differences in all crosses. Our results indicate that deviations from Mendelian expectations in recombination experiments involving the D. pseudoobscura inversion system are most likely due to fitness differences of gene arrangement karyotypes in different environments.     

Excess of nonsynonymous polymorphism at Acph-1 in different gene arrangements of Drosophila subobscura

Nucleotide variation in the Acph-1 gene region was analyzed in a natural population of Drosophila subobscura from Bizerte (Tunisia). The lines studied differed in their gene arrangement for segment I of the O chromosome: 21 lines were O3+4+8, 21 were O3+4+23, and 3 were O3+4. According to chromosomal phylogenies, O3+4 is a central arrangement from which O3+4+8 and O3+4+23 originated. Strong genetic differentiation at Acph-1 was detected among the different arrangements, which is reflected in strong linkage disequilibrium between the variants at informative polymorphic sites and the type of arrangement. Estimates of silent nucleotide diversity are slightly lower within O3+4+23 (pisilent = 0.0166) than within O3+4+8 (pisilent = 0.0228) or O3+4 (pisilent = 0.0234). In contrast, nonsynonymous nucleotide diversity estimates (around 0.1%) are similar in the three arrangements. Most nonsynonymous rare variants are singletons, which results in highly significant Tajima's neutrality tests within the young O3+4+8 and O3+4+23 arrangements. This test is not significant for nonsynonymous mutations within a large Spanish O3+4 sample. In addition, a significant and marginally significant excess of nonsynonymous polymorphism was detected by the McDonald and Kreitman test within O3+4+23 and O3+4+8, respectively. This excess results in a rather high neutrality index (NI = 5.25) when both arrangements are jointly analyzed, in contrast to its value within the old O3+4 arrangement (NI = 1.74). The pattern of variation at Acph-1 within the young arrangements is unusual for nuclear genes and has the same characteristics previously detected in most genes of the mitochondrial genome. Assuming that most nonsynonymous mutations at Acph-1 are under weak negative selection, a smaller effective size of the young arrangements relative to O3+4 might explain the observed results. The relatively low frequency of O3+4+8 and O3+4+23 in the distribution area of D. subobscura, the more recent origin of these arrangements relative to O3+4 and the suppression of recombination in heterokaryotypes might contribute to the relatively small effective size of the young arrangements. Therefore, present results indicate that the differences in effective size and recombination caused by chromosomal arrangements are modulating nonsynonymous variation at Acph-1.     

Colonization of America by Drosophila subobscura: spatial and temporal lethal-gene allelism

About twenty years ago Drosophila subobscura, a western Palearctic species, colonized both North and South America. Lethal genes in the O chromosome has been subject to much research. Lethal gene allelisms between American populations far away have been studied. These allelisms were not negligible, but all cases were due to the lethal gene completely associated to the O5 chromosomal inversion. Here we analyze the lethal genes in a new American population of D. subobscura (Centralia, Washington), located fairly close to a previously studied population (Bellingham, Washington) and separated in space and time with other American populations (Gilroy I and II in California and Santiago de Chile). The frequencies of lethal and semilethal genes of Centralia were 16.9+/-4.6 and 6.2+/-3.0, respectively. The intrapopulational allelism of Centralia was 0.122+/-0.036. Interpopulational allelisms were studied using the lethal genes from the populations separated in space and time from Centralia. The interpopulational allelisms between Centralia and Gilroy I (California) and between Centralia and Bellingham (Washington) were higher than the intrapopulational allelism (0.155+/-0.032 and 0.153+/-0.024, respectively). In all these cases allelism was due to a complete association between a lethal gene and the O5 chromosomal inversion. Accordingly, no other lethal genes are shared in these populations.     

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.     

Seasonal changes of chromosomal inversion polymorphism in a Drosophila subobscura natural population from southeastern European continental refugium of the last glaciation period

Chromosomal inversion polymorphism was studied in a natural population of Drosophila subobscura from Djerdap gorge of Danube river, that represents one of southeastern European continental refugiums during last glaciation period. The existence for changes of chromosomal inversion polymorphism during different periods of time was tested. Observed very high abundance of chromosomal arrangements and lower IFR value grouped this population to central European populations. It was characterized with O7 and O(3+4+5) arrangements, not found in other populations from the region analyzed. Except in this population from warmer habitat, observed E(1+2+9+12) chromosomal arrangement was found, in the region analyzed, only in one population from Mediterranean refugial area. Many similarities have been observed between population from Djerdap gorge and other populations from the region analyzed, indicating that process of post-glacial migration may have started from this southeastern European refugium area. Seasonal changes in the frequencies for Est and O(3+4+7) chromosomal arrangements where Est was more frequent during the spring period and O(3+4+7) during the summer period were found. The changes for J-chromosome (Jst/Jst), E-chromosome (Est/E(1+2+9); E(8)/E(1+2+9) and O-chromosome (Ost/O(3+4+22); O(3+4)/O(3+4); O(3+4+1)/O(3+4+7) whose karyotypic combination frequencies were associated with seasonal changes were also found. These changes most probably represent associations with seasonal climatic changes in the population analyzed, supporting adaptivity of chromosomal inversions in D. subobscura.