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.     

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.     

The phylogeny of nine species of the Drosophila obscura group inferred by the banding homologies of chromosomal regions. I. Element B

The phylogenetic relationships among nine Drosophila species belonging to the obscura group were investigated by establishing the segments displaying banding homologies in their element B (equivalent to the U element of D. subobscura). The phylogenetic ordering of the species was accomplished using overlapping inversions. Two African species, D. kitumensis and D. microlabis, were investigated. These species are homosequential for their element B gene arrangement but differ from that of D. obscura by several rearrangements. Drosophila obscura seems to be most closely related to D. subsilvestris, from which the respective element B gene arrangements differ at least by six inversions. Three species, D. obscura, D. ambigua, and D. tristis, are closely related and form a cluster. Drosophila obscura displays an element B polymorphism for a pericentric inversion for which D. ambigua is fixed for one gene arrangement and D. tristis for the other. Both D. ambigua and D. tristis share a short distal inversion in the small arm of the chromosome, and differ in this respect from D. obscura. Drosophila madeirensis, D. guanche, and D. subobscura all share the same element B gene arrangement, which is acrocentric, but metacentric in all the other species mentioned. It was found that the gene arrangements of the species from the obscura cluster seem to occupy an intermediate position between those of the species of the D. subobscura cluster and those of the African one. The data reported generally are in good agreement with information provided in the literature.     

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.     

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.'     

Changes in chromosomal polymorphism and global warming: The case of Drosophila subobscura from Apatin (Serbia)

In this study, chromosomal inversion polymorphism data for a natural population of Drosophila subobscura from a swampy region near the town of Apatin (Serbia) were compared with data for the same population collected approximately 15 years earlier. The pattern of chromosomal inversion polymorphism changed over time. There were significant increases in the frequency of characteristic southern latitude ("warm" adapted) chromosomal arrangements and significant decreases in the frequency of characteristic northern latitude ("cold" adapted) chromosomal arrangements in the O and U chromosomes. The chromosomal arrangements O(3+4) and O(3+4) (+) (22) (derived from the O(3+4) arrangement) showed significant increases in 2008 and 2009 with regard to the 1994 sample. There was also a significant increase (～50%) in the U(1) (+) (2) arrangement, while U(1+8) (+) (2) (a typical southern arrangement) was detected for the first time. Since the Apatin swampy population of D. subobscura has existed for a long time in a stable habitat with high humidity that has not been changed by man our results indicate that natural selection has produced chromosomal changes in response to the increase in temperature that has occurred in the Balkan Peninsula of central southeastern European.     

CLINAL PATTERNS OF CHROMOSOMAL INVERSION POLYMORPHISMS IN DROSOPHILA SUBOBSCURA ARE PARTLY ASSOCIATED WITH THERMAL PREFERENCES AND HEAT STRESS RESISTANCE

Latitudinal clines in the frequency of various chromosomal inversions are well documented in Drosophila subobscura . Because these clines are roughly parallel on three continents, they have undoubtedly evolved by natural selection. Here, we address whether individuals carrying different chromosomal arrangements also vary in their thermal preferences ( T _p) and heat stress tolerance ( T _ko). Our results show that although T _p and T _ko were uncorrelated, flies carrying "cold-adapted" gene arrangements tended to choose lower temperatures in the laboratory or had a lower heat stress tolerance, in line with what could be expected from the natural patterns. Different chromosomes were mainly responsible for the underlying genetic variation in both traits, which explains why they are linearly independent. Assuming T _p corresponds closely with temperatures that maximize fitness our results are consistent with previous laboratory natural selection experiments showing that thermal optimum diverged among thermal lines, and that chromosomes correlated with T _p differences responded to selection as predicted here. Also consistent with data from the regular tracking of the inversion polymorphism since the colonization of the Americas by D. subobscura , we tentatively conclude that selection on tolerance to thermal extremes is more important in the evolution and dynamics of clinal patterns than the relatively "minor" adjustments from behavioral thermoregulation.     

Evolutionary pace of chromosomal polymorphism in colonizing populations of Drosophila subobscura: an evolutionary time series

Abstract. Biologists have long debated the speed, uniformity, and predictability of evolutionary change. However, evaluating such patterns on a geographic scale requires time-series data on replicate sets of natural populations. Drosophila subobscura has proven an ideal model system for such studies. This fly is broadly distributed in the Old World, but was introduced into both North and South America just over two decades ago and then spread rapidly. Rapid, uniform, and predictable evolution would be demonstrated if the invading flies evolved latitudinal clines that progressively converged on those of the native populations. Evolutionary geneticists quickly capitalized on this opportunity to monitor evolutionary dynamics. Just a few years after the introduction, they surveyed chromosomal inversion frequencies in both North and South America. On both continents they detected incipient latitudinal clines in chromosome inversion frequencies that almost always had the same sign with latitude as in the Old World. Thus the initial evolution of chromosomal polymorphisms on a continental scale was remarkably rapid and consistent. Here we report newer samples of inversion frequencies for the colonizing populations: the time series now spans almost one decade for North America and almost two decades for South America. Almost all inversions in the New World continue to show the same sign of frequency with latitude as in the Old World. Nevertheless, inversion clines have not consistently increased in steepness over time; nor have they consistently continued to converge on the Old World baseline. However, five arrangements in South America show directional, continentwide shifts in frequency. Overall, the initial consistency of clinal evolutionary trajectories seen in the first surveys seems not to have been maintained.     

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.     

Evolutionary history of the third chromosome gene arrangements of Drosophila pseudoobscura inferred from inversion breakpoints

The third chromosome of Drosophila pseudoobscura is polymorphic for numerous gene arrangements that form classical clines in North America. The polytene salivary chromosomes isolated from natural populations revealed changes in gene order that allowed the different gene arrangements to be linked together by paracentric inversions representing one of the first cases where genetic data were used to construct a phylogeny. Although the inversion phylogeny can be used to determine the relationships among the gene arrangements, the cytogenetic data are unable to infer the ancestral arrangement or the age of the different chromosome types. These are both important properties if one is to infer the evolutionary forces responsible for the spread and maintenance of the chromosomes. Here, we employ the nucleotide sequences of 18 regions distributed across the third chromosome in 80-100 D. pseudoobscura strains to test whether five gene arrangements are of unique or multiple origin, what the ancestral arrangement was, and what are the ages of the different arrangements. Each strain carried one of six commonly found gene arrangements and the sequences were used to infer their evolutionary relationships. Breakpoint regions in the center of the chromosome supported monophyly of the gene arrangements, whereas regions at the ends of the chromosome gave phylogenies that provided less support for monophyly of the chromosomes either because the individual markers did not have enough phylogenetically informative sites or genetic exchange scrambled information among the gene arrangements. A data set where the genetic markers were concatenated strongly supported a unique origin of the different gene arrangements. The inversion polymorphism of D. pseudoobscura is estimated to be about a million years old. We have also shown that the generated phylogeny is consistent with the cytological phylogeny of this species. In addition, the data presented here support hypothetical as the ancestral arrangement. One of the youngest arrangements, Arrowhead, has one of the highest population frequencies suggesting that selection has been responsible for its rapid increase.     

A photographic map of Drosophila madeirensis polytene chromosomes.

A photographic map of salivary gland polytene chromosomes of Drosophila madeirensis has been constructed showing homologies and differences with respect to the standard gene arrangement of D. subobscura. Only two paracentric inversions in the X chromosome and some slight minor dissimilarities of one or two bands in the autosomes differentiate the chromosomes of these species.     

Chromosomal polymorphism in natural populations of Drosophila subobscura near Zürich,Switzerland: a contribution to long-term comparisons

The chromosomal inversion polymorphism in natural populations of D. subobscura was investigated near Zürich, Seitzerland, in 1986 and 1987. These data are compared with earlier data from 1963, 1964 and 1984, collected in the same region. In all five acrocentric chromosomes, significant differences in the proportions of the dominant gene arrangements occurred; in the chromosomes A, I, O and U the standard order decreased from the sixties to the eighties, whereas the arrangements I₁, O₃₊₄, O₃₊₄₊₈ and U₁₊₂ became more frequent.     

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.     

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.     

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.     

Genetic load and coadaptation of chromosomal inversions. II. O-chromosomes in Drosophila subobscura populations

We have analysed the inversion polymorphism and genetic load of O-chromosomes in three populations of D. subobscura from southeastern Europe. As expected for a central populations the inversion polymorphism was extensive. In a like fashion, the genetic load, in particular the frequency of lethals, was heavy in all three populations. There were significant differences in the frequency of moderately deleterious genes. These differences in viability can be attributed to balancing selection. A comparison of these two kinds of genetic polymorphism indicates that there are differences in mean viability among different gene arrangements of O-chromosomes in the three populations. The differences observed are due to an unequal distribution of various viability classes among O-chromosome gene arrangements. We here show for the first time a specific distribution of lethal genes among these arrangements within the Palearctic distribution area of D. subobscura. The lethal allelism test showed lethals are non-randomly associated with the Ost gene arrangement. The amount of genetic load is heavy in gene arrangements with a high frequency, in comparison with the ones with a low frequency. Lethal genes may be protected in combinations of low and moderate frequency gene arrangements that harbor more lethal genes, as the Ost in the one population. Some arrangements that are less protected against recombination have a higher load than ones that are more protected against recombination. This can be taken as evidence for coadaptation.     

Temporal and spatial variation of inversion polymorphism in two natural populations of Drosophila buzzatii

The inversion polymorphism of the cactophilic fly Drosophila buzzatii was studied in two natural populations. We assessed the temporal changes and microspatial population structure. We observed a significant increase in the frequency of arrangement 2J at the expense of 2ST in both populations. These gene arrangements appear to affect the life-history of flies differently. Environmental heterogeneity explains the karyotype coexistence in nature. The analysis of population structure showed that differentiation of inversion frequencies among individual breeding sites, the rotting clacodes of Opuntia vulgaris, was highly significant. The karyotypic frequencies did not depart significantly from Hardy-Weinberg expectations, neither in individual rots nor in the total population. These results suggest that the observed population structure can be easily accounted by random genetic drift.     

Diurnal variability of gene arrangement frequencies in Drosophila subobscura populations from two habitats*

The rich inversion polymorphism of chromosomes of Drosophila subobscura varies in association with environmental changes at spatial and temporal level. Due to random factors this might become less regular and this might be the reason that there is little evidence of altitudinal and seasonal, microgeographic and habitat‐related variability. The variability of gene arrangement frequencies over daytime period was investigated in populations of D. subobscura from two ecologically different habitats. According to gene arrangement frequencies and inversion polymorphism parameters populations fit into the existing patterns of regional polymorphism. Even though populations display daytime variability for the frequencies of arrangements of chromosomes U and J, the nonsignificant differences in other arrangement frequencies, as well as data obtained from genetic distances between samples from different times of day, rule out the existence of subpopulations distinguished from each other by diurnal activity within two habitats. Some particular arrangement frequencies vary according to temperature and humidity and some changes appear for different karyotypic combinations between beech and oak habitat, which suggests different adaptive advantages and selective mechanisms. Cluster analysis of gene arrangement frequencies indicates some importance of A chromosome arrangements for diurnal activity and show in dependence to other ecological data that internally fixed daily activity rhythm might exists in D. subobscura. Analysis of light, temperature and humidity factors in relation to the frequency data obtained indicate a choice of different ecological niche in a temporally structured habitat due to behavioural genetic information.     

Rapid response to perturbation of chromosome frequencies in natural populations of Drosophila robusta

Perturbation of gene or chromosome frequencies in natural populations is one of the most powerful ways of demonstrating whether natural selection maintains genetic polymorphism or if other evolutionary forces are at work. Gene arrangement frequencies in two natural populations of Drosophila robusta were perturbed multiple times by releasing adult flies with contrasting karyotypes and carefully monitoring post-perturbation presence of hybrids and chromosome frequencies. In all cases, frequencies quickly returned to pre-perturbation levels, and in the following sampling periods, no evidence of the introduced chromosomes was apparent. Analysis of post-perturbation frequency changes included tests for heterogeneity among chromosome arrangements in rates of return to equilibrium values using population admixture analysis. In several cases, significant heterogeneity was detected indicating some form of natural selection was operating. Technical challenges to carrying out perturbation experiments in the wild are also discussed.     

The evolutionary history of Drosophila buzzatti. XXVI. Macrogeographic patterns of inversion polymorphism in New World populations

Inversion polymorphisms in the second and fourth chromosomes of the cactophilic Drosophila buzzatti in the native distribution range of the species are described. Over 5,000 flies from 26 localities were scored revealing interesting geographic structuring of arrangement frequencies. Multiple regression and partial correlation approaches showed that the frequencies of second and fourth chromosome arrangements vary clinically along latitudinal and altitudinal gradients and to a lesser extent with longitude. Although many non selective explanations can account for this pattern, its resemblance to the clinal pattern described in recently established Australian populations of Drosophila buzzatii, strongly suggests a selective explanation. Additionally, the correlated variation observed between the frequencies of arrangements 2St on the second chromosome and 4St on the fourth suggests a pattern of interchromosomal association, which, when considering the vast area surveyed, might be explained as the result of epistatic interactions. The analysis of population structure revealed a significant regional pattern, concordant with previously described phytogeographic regions. F-statistics showed that the patterns of variation were different not only between the second and fourth chromosomes, but also between second chromosome arrangements, suggesting that selective differentiation might have contributed to population structure. Since D. buzzatii breeds and feeds on the decaying tissues of diverse cactus species present in different phytogeographic regions, it is difficult to distinguish the underlying causes of the geographic patterns observed. However, inversion heterozygosity is not correlated with the diversity of potential cactus hosts. The evidence presented suggests that differential selection may be the main cause for the population structure. It is also possible to conclude that the inversion polymorphism of D. buzzatti is flexible rather than rigid.