Seasonal cycles of allozyme-by-chromosomal-inversion gametic disequilibrium in Drosophila subobscura

Allozyme loci are frequently found non randomly associated to the chromosomal inversions in which they are included in Drosophila. Two opposite views compete to explain strong allozyme-by-inversion gametic disequilibria: they result from natural selection or, conversely, merely represent remnants of associations accidentally established at the origin of inversions. Empirical efforts aimed at deciding between adaptive and historical scenarios have focused on the spatial distribution of disequilibria. Yet, the evolutionary significance of these associations remains uncertain. I report here the results of a time-series analysis of the seasonal variation of alleles at six allozyme loci (Acph, Lap, Pept-1, Ao, Mpi, and Xdh) in connection with the O chromosomal polymorphisms of D. subobscura. The findings were: (1) in the segment I of the O chromosome, Lap and Pept-1 allozymes changed seasonally in a cyclical fashion within the ST gene arrangement, but they changed erratically within the 3 + 4 gene configuration; (2) the frequencies of Lap1.11 and Pept-1(0.40) within ST dropped to their lowest values in early and late summer, respectively, when the seasonal level of the ST arrangement is lowest. Furthermore, Lap1.11 and Pept-1(0.40) covary with ST only within these seasons, yet in a fashion inconsistent with these alleles having a major influence on the dynamics of the inversion; (3) seasonal cycling of alleles within inversions were not detected at Acph, Ao, Mpi, and Xdh, yet these loci are nearly monomorphic at the study population, and/or their sampled series were shorter than those for Lap and Pept-1; and (4) simply monitoring allozyme frequencies separately for each inversion proved to be superior, for evidencing the seasonal cycles of the disequilibria, to the use of the D' coefficient of association. Observed seasonal cycles of allozymes within inversions likely reflect natural selection.     

Latitudinal differences in sex chromosome inversions,sex linked allozymes,and mitochondrial DNA variation in Drosophila subobscura

Mitochondrial DNA, sex linked allozymes, and chromosome A gene arrangement data, from eight European natural populations of Drosophila subobscura, were analyzed to determine the existence of latitudinal clines. Strong north-south correlations with latitude were found for gene arrangements and for the Hbdh and 6Pgdh allozymes. Gametic associations between the A₂ gene arrangement, the Hbdh ⁹⁶ and the 6Pgdh ⁹⁶ alleles, point out some kind of epistatic interaction. At mtDNA level, the Hae III, A variant did not show a previously found north-south clinal variation.     

The D' measure of overall gametic disequilibrium between pairs of multiallelic loci

Abstract The D ' coefficient is one of the most commonly used measures of the extent of gametic disequilibrium between multiallelic loci. It has been suggested that the range of the D ' measure of overall disequilibrium between pairs of multiallelic loci depends on allele frequencies, except under some very restricted conditions. Nevertheless, the problem of dependence of the range of D ' has not been characterized under a wide set of possible polymorphisms. Evaluation of the utility of D ' as a measure of the strength of overall disequilibrium between all possible pairs of alleles at two multiallelic loci requires better knowledge of its range than is currently available. In this work, the conditions of polymorphism under which the range of D ' is frequency independent are given. It is found that the range of D ' is more often independent of allelic frequencies than is commonly thought. Furthermore, the range of D ' undergoes only small fluctuations as a function of the polymorphisms at the loci. Numerical cases and microsatellite data from humans are used for illustration. These observations indicate that the D ' coefficient is a useful tool for the estimation and comparison of the extent of overall disequilibrium across pairs of multiallelic loci.     

Inbreeding reveals interpopulation differences in inversion polymorphism of Drosophila subobscura

The adaptive significance of inversion polymorphism of Drosophila subobscura is well established. However, little is known about gene combinations within inversions which are coadapted because of population-specific effects. We studied this aspect of Dobzhansky's coadaptation hypothesis, using the systematic inbreeding method. Differences in magnitude and quality of inversion polymorphism in two ecologically and topologically distinct habitats were compared after several generations of continuous full-sib inbreeding. Populations from the two habitats differ in frequency of homokaryotypes after third and fifth generation of inbreeding and in the levels of homozygosity of different gene arrangements. The effect of homozygosity appears population and chromosome specific.     

Contrasting population genetic structures using allozymes and the inversion polymorphism in Drosophila buzzatii

Second chromosome inversion and genotypic frequencies at seven allozyme loci, differentially associated with inversions, were determined in seven natural populations of Drosophila buzzatii. The patterns of variation of allozymes and the inversion polymorphisms were significantly different, indicating the role of adaptive differentiation for the latter. Moreover, the patterns of population structure varied among allozyme loci, suggesting the operation of diversifying selection for certain loci. Differentiation was negligible for Leucyl‐amino peptidase (Lap) and Peptidase‐2 (Pep‐2), low to moderate for Aldehyde oxidase (Aldox), Peptidase‐1 (Pep‐1) and Esterase‐1 (Est‐1) and high for Esterase‐2 (Est‐2) and Xanthine dehydrogenase (Xdh). Significant linkage disequilibria were detected between inversions and Aldox, Est‐1, Est‐2 and Xdh. Multiple regression analyses of inversion and allele frequencies on environmental variables revealed the existence of clines for inversions, Est‐1, Est‐2, Xdh and Aldox along altitudinal, latitudinal and/or climatic gradients. Tests using conditional allele frequencies showed that Est‐1 and Aldox clines could be accounted for by hitchhiking with inversions, whereas natural selection should be invoked to explain the clines observed for Est‐2 and Xdh.     

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.     

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.     

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.     

Chromosomal inversion polymorphism of Drosophila subobscura from south‐eastern part of Europe

Five Drosophila subobscura populations from south‐eastern Europe were investigated with respect to chromosomal inversion polymorphism. They belong to central European populations and their gene arrangement frequencies are in accordance with latitudinal variation observed in this species. The new complex gene arrangement A₂₊₈₊₉, found in one of the populations analysed, may be specific for the southern Balkan refugium. A general tendency of gene arrangement frequencies to exhibit some southern characteristics was found for some of the autosomal standard arrangements.     

Preadult competition between Drosophila subobscura and Drosophila pseudoobscura

The Palaearctic species Drosophila subobscura has recently colonized a large area of North America where it coexists with Drosophila pseudoobscura. The viability and developmental rate of these species were studied at 13 d?C, 18d?C and 23 d?C and at densities of 10, 50, 100 and 200 eggs per vial. The two species were differently affected by density and temperature in the ranges studied. Both intra- and interspecific cultures showed that D. pseudoobscura was best adapted to 23 d?C, where it was clearly the dominant species. On the other hand, at 18 d?C and especially at 13 d?C, although D. subobscura was less viable than D. pseudoobscura, its developmental time was shorter, which may give advantage to this species. Results reported here agree with the observed distribution of these species in North America.     

Chromosomal inversions promote genomic islands of concerted evolution of Hsp70 genes in the Drosophila subobscura species subgroup

Heat‐shock (HS) assays to understand the connection between standing inversion variation and evolutionary response to climate change in Drosophila subobscura found that “warm‐climate” inversion O₃₊₄ exhibits non‐HS levels of Hsp70 protein like those of “cold‐climate” O_ST after HS induction. This was unexpected, as overexpression of Hsp70 can incur multiple fitness costs. To understand the genetic basis of this finding, we have determined the genomic sequence organization of the Hsp70 family in four different inversions, including O_ST, O₃₊₄ , O₃₊₄₊₈ and O₃₊₄₊₁₆ , using as outgroups the remainder of the subobscura species subgroup, namely Drosophila madeirensis and Drosophila guanche. We found (i) in all the assayed lines, the Hsp70 family resides in cytological locus 94A and consists of only two genes, each with four HS elements (HSEs) and three GAGA sites on its promoter. Yet, in O_ST, the family is comparatively more compact; (ii) the two Hsp70 copies evolve in concert through gene conversion, except in D. guanche; (iii) within D. subobscura, the rate of concerted evolution is strongly structured by inversion, being higher in O_ST than in O₃₊₄ ; and (iv) in D. guanche, the two copies accumulated multiple differences, including a newly evolved “gap‐type” HSE2. The absence of concerted evolution in this species may be related to a long‐gone‐unnoticed observation that it lacks Hsp70 HS response, perhaps because it has evolved within a narrow thermal range in an oceanic island. Our results point to a previously unrealized link between inversions and concerted evolution, with potentially major implications for understanding genome evolution.     

Chromosomal inversion polymorphism of Drosophila subobscura populations from Jastrebac Mountain shows temporal and habitat-related changes

Inversion polymorphism in populations of D. subobscura from a beech forest on Jastrebac mountain was studied in June 1990, 1993, and 1994, respectively. The same analysis was performed in 1990 for D. subobscura populations in a beech forest and an oak forest in the same region. Statistically significant differences in the frequencies of the gene arrangements of A₁, J and U chromosome were observed during the period of investigation. A tendency towards a decrease in the frequency of the standard gene arrangements was found for all chromosomes, but was particularly evident with chromosomes A and J. The frequency of the gene arrangements A₁ A₂, J₁ and U_1–2+6 increased at the same time. Differences in the frequency of the gene arrangements of A, J and U chromosomes were also observed when the populations from two ecologically different habitats (beech and oak forest) were compared in 1990.     

Relationship between chromosomal polymorphism and wing size in a natural population of Drosophila subobscura

Chromosomal polymorphism and wing size (as a measure of body size) were analysed simultaneously in two samples of Drosophila subobscura from Barcelona, Spain. The very rich chromosomal polymorphism of this species makes it difficult to detect the relationship of this polymorphism with any phenotypical character. However, a positive significant regression of wing size on the percentage of the autosome length with standard arrangement was found. Furthermore, for each polymorphic chromosome, except for the J chromosome, an association between the most frequent arrangements and wing size was observed. This trend, which was the same in the two samples, was that expected according to the latitudinal clines of both characters.     

Linkage disequilibrium in laboratory populations of Drosophila nasuta

Summary Natural populations of Drosophila nasuta are polymorphic for many gene arrangements. Two non-overlapping inversions of the third chromosome, III-2 and III-35, are most common and display extreme linkage disequilibrium. Six randomly mating laboratory stocks, each founded by one gravid female heterozygous in coupling for both III-2 and III-35, were observed after 32 generations. Significant linkage disequilibrium was observed in all stocks. Recombinants were found in only two stocks. The absence of effective recombination in some stocks and its presence in others might be due to different genotypic backgrounds. We suggest that natural selection, influencing recombination rates in several ways, and intrachromosomal epistasis between the two inversions were the main factors for the maintenance of linkage disequilibrium in D. nasuta.     

Long-term changes in the chromosomal inversion polymorphism of Drosophila subobscura. I. Mediterranean populations from southwestern Europe

The chromosomal polymorphism of seven Mediterranean populations of Drosophila subobscura has been compared with that of the same populations collected 26 to 35 years ago. Significant latitudinal clines for the frequencies of A(ST), E(ST), O(ST). and U(ST) chromosomal arrangements have been detected in the old and new samples. Standard gene arrangements are frequent in the north and decrease in frequency towards the south. Significant negative regression coefficients between latitude and transformed frequency have also been observed for the more frequent nonstandard gene arrangements. The pattern of the clines is practically the same in the old and new collections. Furthermore, the frequencies of gene arrangements of all chromosomes have changed significantly during this period in a systematic way: an increase in the frequency of those arrangements typical of southern latitudes and a decrease for those more common in northern latitudes is observed in all populations. These changes could be due to climatic factors that are correlated with latitude, making the chromosomal composition of this species more "southern.'     

Latitudinal variation for two enzyme loci and an inversion polymorphism in Drosophila melanogaster from Central and South America

Abstract.— Many organisms show latitudinal variation for various genetically determined traits. Such clines may involve neutral variation and originate from historical events or their maintenance may be explained by selection. For Drosophila melanogaster , latitudinal variation for allozymes, inversions, and quantitative traits has been found on several continents. We sampled D. melanogaster populations in Panama and along a transect of 40 latitudinal degrees on the west coast of South America. Negative correlations with latitude were found for Adh^s and _αGpdh^F allele frequencies and for the frequency of the cosmopolitan inversion In(2L)t in Adh^sαpdh^F chromosomes. A positive correlation existed between wing length and latitude. Significant correlations were found between these traits and climatic variables like temperature and rainfall. The observed clines show considerable resemblance to those found on other continents. Gametic disequilibrium between Adh^s and αGpdh^F occurred predominantly at higher latitudes and was caused by the presence of In(2L)t . The reasons for the clinal distributions are discussed and it is argued that selection is the most likely explanation. However, the exact nature of the selective force and the interactions of allozymes with each other and with In(2L)t are complex and not fully understood. In tropical regions In(2L)t -containing genotypes have higher fitness than ST/ST and Adh and αGpdh hitchhike with the inversion, but there is also evidence for balancing selection at the Adh locus.     

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.     

Long‐term changes in the chromosomal inversion polymorphism of Drosophila subobscura. II. European populations*

The chromosomal polymorphism of 13 European populations of Drosophila subobscura has been compared with that of the same populations collected 15–35 years ago. The chromosomal polymorphism of the old populations differs significantly from that of the new populations, mainly for chromosomes U and O. There is a very good agreement between the geographical space and the genetic space as shown by a graphical representation of the 26 statistical populations (13 old and 13 new) obtained by a principal coordinate analysis. This reflects both the existence of significant latitudinal clines for the frequencies of some chromosomal arrangements in the old and new samples and systematic changes that have taken place in these populations during the period that elapsed between the two surveys. An increase in the frequency of those arrangements typical of southern latitudes and a decrease for those common in northern latitudes is observed in all populations – Mediterranean, Atlantic and Central European. Furthermore, the genetic distances of the new populations to a southern population of reference have decreased in comparison with those of the old populations. These changes could be the result of climatic factors that are correlated with latitude. In particular, the assumption that global warming is responsible for all the changes observed appears rather likely. Whether these systematic changes of the chromosomal polymorphism are a consequence of local adaptations or have been produced by migration from the south remains an open question.