The Genetic Content of Chromosomal Inversions across a Wide Latitudinal Gradient

There is increasing evidence regarding the role of chromosomal inversions in relevant biological processes such as local adaptation and speciation. A classic example of the adaptive role of chromosomal polymorphisms is given by the clines of inversion frequencies in Drosophila subobscura, repeatable across continents. Nevertheless, not much is known about the molecular variation associated with these polymorphisms. We characterized the genetic content of ca. 600 individuals from nine European populations following a latitudinal gradient by analysing 19 microsatellite loci from two autosomes (J and U) and the sex chromosome (A), taking into account their chromosomal inversions. Our results clearly demonstrate the molecular genetic uniformity within a given chromosomal inversion across a large latitudinal gradient, particularly from Groningen (Netherlands) in the north to Málaga (Spain) in the south, experiencing highly diverse environmental conditions. This low genetic differentiation within the same gene arrangement across the nine European populations is consistent with the local adaptation hypothesis for th evolutionof chromosomal polymorphisms. We also show the effective role of chromosomal inversions in maintaining different genetic pools within these inverted genomic regions even in the presence of high gene flow. Inversions represent thus an important barrier to gene flux and can help maintain specific allelic combinations with positive effects on fitness. Consistent patterns of microsatellite allele-inversion linkage disequilibrium particularly in loci within inversions were also observed. Finally, we identified areas within inversions presenting clinal variation that might be under selection.     

Investigating latitudinal clines for life history and stress resistance traits in Drosophila simulans from eastern Australia

Latitudinal clines have been demonstrated for many quantitative traits in Drosophila and are assumed to be due to climatic selection. However, clinal studies are often performed in species of Drosophila that contain common cosmopolitan inversion polymorphisms that also show clinal patterns. These inversion polymorphisms may be responsible for much of the observed clinal variation. Here, we consider latitudinal clines for quantitative traits in Drosophila simulans from eastern Australia. Drosophila simulans does not contain cosmopolitan inversion polymorphisms, so allows the study of clinal selection on quantitative traits that are not confounded by associations with inversions. Body size showed a strong linear cline for both females and males. Starvation resistance exhibited a weak linear cline in females, whereas chill-coma recovery exhibited a significant nonlinear cline in females only. No clinal pattern was evident for development time, male chill-coma recovery, desiccation or heat resistance. We discuss these results with reference to the role inversion polymorphisms play in generating clines in quantitative traits of Drosophila.     

The effect of urbanisation on avian morphology and latitudinal gradients in body size

Urban areas occupy a large and growing proportion of the earth. Such sites exhibit distinctive characteristics relative to adjacent rural habitats, and many species have colonised and now successfully exploit urban habitats. The change in selection pressures as a result of urbanisation has led to trait divergence in some urban populations relative to their rural counterparts, but studies have generally been local in scale and the generality of differentiation thus remains unknown. The European blackbird Turdus merula is one of the commonest urban bird species in the Western Palearctic, but populations vary substantially in the length of time they have been urbanised. Here we investigate patterns of morphological variation in European blackbirds occupying 11 paired urban and rural habitats across much of the urbanised range of this species and spanning 25° of latitude. First, we assessed the extent to which urban and rural blackbirds are differentiated morphologically and the consistency of any differentiation across the range. Paired urban and rural Blackbird populations frequently exhibited significant morphological differences, but the magnitude and direction of differentiation was site dependent. We then investigated whether the nature of latitudinal gradients in body‐size differed between urban and rural populations, as predicted by differences in the climatic regimes of urban and rural areas. Blackbird body‐size exhibited strong latitudinal gradients, but their form did not differ significantly between urban and rural habitats. The latitudinal gradient in body size may be a consequence of Seebohm's rule, that more migratory populations occurring at high latitudes have longer wings. We conclude that while there can be substantial morphological variation between adjacent urban and rural bird populations, such differentiation may not apply across a species’ range. Locality specific differentiation of urban and rural blackbirds may arise if the selection pressures acting on blackbird morphology vary in an inconsistent manner between urban and rural habitats. Alternatively, phenotypic divergence could arise in a stochastic manner depending on the morphological traits of colonists, through founder effects.     

Wing trait–inversion associations in Drosophila subobscura can be generalized within continents,but may change through time

Clinal variation is one of the most emblematic examples of the action of natural selection at a wide geographical range. In Drosophila subobscura, parallel clines in body size and inversions, but not in wing shape, were found in Europe and South and North America. Previous work has shown that a bottleneck effect might be largely responsible for differences in wing trait–inversion association between one European and one South American population. One question still unaddressed is whether the associations found before are present across other populations of the European and South American clines. Another open question is whether evolutionary dynamics in a new environment can lead to relevant changes in wing traits–inversion association. To analyse geographical variation in these associations, we characterized three recently laboratory founded D. subobscura populations from both the European and South American latitudinal clines. To address temporal variation, we also characterized the association at a later generation in the European populations. We found that wing size and shape associations can be generalized across populations of the same continent, but may change through time for wing size. The observed temporal changes are probably due to changes in the genetic content of inversions, derived from adaptation to the new, laboratory environment. Finally, we show that it is not possible to predict clinal variation from intrapopulation associations. All in all this suggests that, at least in the present, wing traits–inversion associations are not responsible for the maintenance of the latitudinal clines in wing shape and size.     

Trait variation in response to varying winter temperatures,diversity patterns and signatures of selection along the latitudinal distribution of the widespread grassland plant Arrhenatherum elatius

Across Europe, genetic diversity can be expected to decline toward the North because of stochastic and selective effects which may imply diminished phenotypic variation and less potential for future genetic adaptations to environmental change. Understanding such latitudinal patterns can aid provenance selection for breeding or assisted migration approaches. In an experiment simulating different winter temperatures, we assessed quantitative trait variation, genetic diversity, and differentiation for natural populations of the grass Arrhenatherum elatius originating from a large latitudinal gradient. In general, populations from the North grew smaller and had a lower flowering probability. Toward the North, the absolute plastic response to the different winter conditions as well as heritability for biomass production significantly declined. Genetic differentiation in plant height and probability of flowering were very strong and significantly higher than under neutral expectations derived from SNP data, suggesting adaptive differentiation. Differentiation in biomass production did not exceed but mirrored patterns for neutral genetic differentiation, suggesting that migration‐related processes caused the observed clinal trait variation. Our results demonstrate that genetic diversity and trait differentiation patterns for A. elatius along a latitudinal gradient are likely shaped by both local selection and genetic drift.     

Ecological genomics of Anopheles gambiae along a latitudinal cline: a population-resequencing approach

The association between fitness-related phenotypic traits and an environmental gradient offers one of the best opportunities to study the interplay between natural selection and migration. In cases in which specific genetic variants also show such clinal patterns, it may be possible to uncover the mutations responsible for local adaptation. The malaria vector, Anopheles gambiae, is associated with a latitudinal cline in aridity in Cameroon; a large inversion on chromosome 2L of this mosquito shows large differences in frequency along this cline, with high frequencies of the inverted karyotype present in northern, more arid populations and an almost complete absence of the inverted arrangement in southern populations. Here we use a genome resequencing approach to investigate patterns of population divergence along the cline. By sequencing pools of individuals from both ends of the cline as well as in the center of the cline- where the inversion is present in intermediate frequency- we demonstrate almost complete panmixia across collinear parts of the genome and high levels of differentiation in inverted parts of the genome. Sequencing of separate pools of each inversion arrangement in the center of the cline reveals large amounts of gene flux (i.e., gene conversion and double crossovers) even within inverted regions, especially away from the inversion breakpoints. The interplay between natural selection, migration, and gene flux allows us to identify several candidate genes responsible for the match between inversion frequency and environmental variables. These results, coupled with similar conclusions from studies of clinal variation in Drosophila, point to a number of important biological functions associated with local environmental adaptation.     

Microgeographic Differentiation of Chromosomal and Enzyme Polymorphisms in DROSOPHILA PERSIMILIS

We studied microgeographic and temporal genetic differentiation in natural populations of Drosophila persimilis with respect to chromosome inversion and enzyme polymorphisms. Both inversion frequencies and allozyme frequencies varied significantly over short distances. Neither differed significantly between morning and evening collections. Because several studies of the dispersal behavior of this species have been performed, we attempt to fit the observed data to mathematical models which relate dispersion to random genetic drift and to spatially varying selection coefficients. We conclude that the observations are due at least partly to behavioral differences among genotypes. i.e., habitat preferences. These results have implications for genetic load theory and models of selection in heterogeneous environments.     

Genome‐wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome‐wide next‐generation sequencing of DNA pools (‘pool‐seq’) from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome‐wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well‐known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.     

Variation in Drosophila melanogaster central metabolic genes appears driven by natural selection both within and between populations

In this report, we examine the hypothesis that the drivers of latitudinal selection observed in the eastern US Drosophila melanogaster populations are reiterated within seasons in a temperate orchard population in Pennsylvania, USA. Specifically, we ask whether alleles that are apparently favoured in northern populations are also favoured early in the spring, and decrease in frequency from the spring to autumn with the population expansion. We use SNP data collected for 46 metabolic genes and 128 SNPs representing the central metabolic pathway and examine for the aggregate SNP allele frequencies whether the association of allele change with latitude and that with increasing days of spring–autumn season are reversed. Testing by random permutation, we observe a highly significant negative correlation between these associations that is consistent with this expectation. This correlation is stronger when we confine our analysis to only those alleles that show significant latitudinal changes. This pattern is not caused by association with chromosomal inversions. When data are resampled using SNPs for amino acid change the relationship is not significant but is supported when SNPs associated with cis-expression are only considered. Our results suggest that climate factors driving latitudinal molecular variation in a metabolic pathway are related to those operating on a seasonal level within populations.     

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.     

Does thermal plasticity align with local adaptation? An interspecific comparison of wing morphology in sepsid flies

Although genetic and plastic responses are sometimes considered as unrelated processes, their phenotypic effects may often align because genetic adaptation is expected to mirror phenotypic plasticity if adaptive, but run counter to it when maladaptive. Because the magnitude and direction of this alignment has further consequences for both the tempo and mode of adaptation, they are relevant for predicting an organisms’ reaction to environmental change. To better understand the interplay between phenotypic plasticity and genetic change in mediating adaptive phenotypic variation to climate variability, we here quantified genetic latitudinal variation and thermal plasticity in wing loading and wing shape in two closely related and widespread sepsid flies. Common garden rearing of 16 geographical populations reared across multiple temperatures revealed that wing loading decreases with latitude in both species. This pattern could be driven by selection for increased dispersal capacity in the cold. However, although allometry, sexual dimorphism, thermal plasticity and latitudinal differentiation in wing shape all show similar patterns in the two species, the relationship between the plastic and genetic responses differed between them. Although latitudinal differentiation (south to north) mirrored thermal plasticity (hot to cold) in Sepsis punctum, there was no relationship in Sepsis fulgens. While this suggests that thermal plasticity may have helped to mediate local adaptation in S. punctum, it also demonstrates that genetic wing shape differentiation and its relation to thermal plasticity may be complex and idiosyncratic, even among ecologically similar and closely related species. Hence, genetic responses can, but do not necessarily, align with phenotypic plasticity induced by changing environmental selection pressures.     

Latitudinal variation in plant chemical defences drives latitudinal patterns of leaf herbivory

A long‐standing paradigm in ecology holds that herbivore pressure and thus plant defences increase towards lower latitudes. However, recent work has challenged this prediction where studies have found no relationship or opposite trends where herbivory or plant defences increase at higher latitudes. Here we tested for latitudinal variation in herbivory, chemical defences (phenolic compounds), and nutritional traits (phosphorus and nitrogen) in leaves of a long‐lived tree species, the English oak Quercus robur. We further investigated the underlying climatic and soil factors associated with such variation. Across 38 populations of Q. robur distributed along an 18° latitudinal gradient, covering almost the entire latitudinal and climatic range of this species, we observed strong but divergent latitudinal gradients in leaf herbivory and leaf chemical defences and nutrients. As expected, there was a negative relationship between latitude and leaf herbivory where oak populations from lower latitudes exhibited higher levels of leaf herbivory. However, counter to predictions there was a positive relationship between leaf chemical defences and latitude where populations at higher latitudes were better defended. Similarly, leaf phosphorus and nitrogen increased with latitude. Path analysis indicated a significant (negative) effect of plant chemical defences (condensed tannins) on leaf herbivory, suggesting that the latitudinal gradient in leaf herbivory was driven by an inverse gradient in defensive investment. Leaf nutrients had no independent influence on herbivory. Further, we found significant indirect effects of precipitation and soil porosity on leaf herbivory, which were mediated by plant chemical defences. These findings suggest that abiotic factors shape latitudinal variation in plant defences and that these defences in turn underlie latitudinal variation in leaf herbivory. Overall, this study contributes to a better understanding of latitudinal variation in plant–herbivore interactions by determining the identity and modus operandi of abiotic factors concurrently shaping plant defences and herbivory.     

Genomic Patterns of Geographic Differentiation in Drosophila simulans

Geographic patterns of genetic differentiation have long been used to understand population history and to learn about the biological mechanisms of adaptation. Here we present an examination of genomic patterns of differentiation between northern and southern populations of Australian and North American Drosophila simulans, with an emphasis on characterizing signals of parallel differentiation. We report on the genomic scale of differentiation and functional enrichment of outlier SNPs. While, overall, signals of shared differentiation are modest, we find the strongest support for parallel differentiation in genomic regions that are associated with regulation. Comparisons to Drosophila melanogaster yield potential candidate genes involved in local adaptation in both species, providing insight into common selective pressures and responses. In contrast to D. melanogaster, in D. simulans we observe patterns of variation that are inconsistent with a model of temperate adaptation out of a tropical ancestral range, highlighting potential differences in demographic and colonization histories of this cosmopolitan species pair.     

Quantifying effects of environmental and geographical factors on patterns of genetic differentiation

Elucidating the factors influencing genetic differentiation is an important task in biology, and the relative contribution from natural selection and genetic drift has long been debated. In this study, we used a regression-based approach to simultaneously estimate the quantitative contributions of environmental adaptation and isolation by distance on genetic variation in Boechera stricta, a wild relative of Arabidopsis. Patterns of discrete and continuous genetic differentiation coexist within this species. For the discrete differentiation between two major genetic groups, environment has larger contribution than geography, and we also identified a significant environment-by-geography interaction effect. Elsewhere in the species range, we found a latitudinal cline of genetic variation reflecting only isolation by distance. To further confirm the effect of environmental selection on genetic divergence, we identified the specific environmental variables predicting local genotypes in allopatric and sympatric regions. Water availability was identified as the possible cause of differential local adaptation in both geographical regions, confirming the role of environmental adaptation in driving and maintaining genetic differentiation between the two major genetic groups. In addition, the environment-by-geography interaction is further confirmed by the finding that water availability is represented by different environmental factors in the allopatric and sympatric regions. In conclusion, this study shows that geographical and environmental factors together created stronger and more discrete genetic differentiation than isolation by distance alone, which only produced a gradual, clinal pattern of genetic variation. These findings emphasize the importance of environmental selection in shaping patterns of species-wide genetic variation in the natural environment.     

Genomics of clinal variation in Drosophila: disentangling the interactions of selection and demography

Clines in phenotypes and genotype frequencies across environmental gradients are commonly taken as evidence for spatially varying selection. Classical examples include the latitudinal clines in various species of Drosophila, which often occur in parallel fashion on multiple continents. Today, genomewide analysis of such clinal systems provides a fantastic opportunity for unravelling the genetics of adaptation, yet major challenges remain. A well‐known but often neglected problem is that demographic processes can also generate clinality, independent of or coincident with selection. A closely related issue is how to identify true genic targets of clinal selection. In this issue of Molecular Ecology, three studies illustrate these challenges and how they might be met. Bergland et al. report evidence suggesting that the well‐known parallel latitudinal clines in North American and Australian D. melanogaster are confounded by admixture from Africa and Europe, highlighting the importance of distinguishing demographic from adaptive clines. In a companion study, Machado et al. provide the first genomic comparison of latitudinal differentiation in D. melanogaster and its sister species D. simulans. While D. simulans is less clinal than D. melanogaster, a significant fraction of clinal genes is shared between both species, suggesting the existence of convergent adaptation to clinaly varying selection pressures. Finally, by drawing on several independent sources of evidence, Bo?i?evi? et al. identify a functional network of eight clinal genes that are likely involved in cold adaptation. Together, these studies remind us that clinality does not necessarily imply selection and that separating adaptive signal from demographic noise requires great effort and care.     

Latitudinal variation of morphological characteristics in the swallowtail Sericinus montelus Gray, 1798 (Lepidoptera: Papilionidae)

Patterns of latitudinal variation in the phenotype or genotype of an organism may provide evidence for natural selection. In this study, we investigated seven populations of swallowtail Sericinus montelus Gray, 1798 (Lepidoptera: Papilionidae), a non‐migratory species, to explore the latitudinal variation of morphological characteristics in adults. The results showed that body size and the development of dark pigmentation on wings in this species responded strongly to latitude. The body size of both male and female adult of S. montelus was negatively correlated with latitude. These findings provided solid evidence to support the converse Bergmann's rule. We considered that the observed variation in morphological characteristics was most likely mediated by the seasonal length and thermoperiod to adapt to different latitudinal environment (e.g. shortened developmental time of immature stages for smaller body size at higher latitude). Moreover, the tendency towards progressively darker colour patterns (only in adult males) at increasingly low latitudes was consistent with Gloger's rule. We suggested that the observed colour variation was most likely associated with thermoregulation. Slight variation in the morphology of the W‐shaped stripe on the forewing of adult females was also found, and we presumed that the functions of sexual preferences, mimicry and thermoregulation might be involved.     

Latitudinal clines in the grasshopper Dichroplus elongatus: Coevolution of the A genome and B chromosomes?

Argentine populations of Dichroplus elongatus (Orthoptera: Acrididae) are polymorphic for B chromosomes. Previous studies showed that B chromosomes affect body size and some fitness components in Northwestern populations. We studied phenotype and B′s variation patterns along a latitudinal cline as well as the relationship between karyotype and body size related traits in 17 populations from East. Body size related traits showed a ‘saw tooth’ pattern of variation being small at low and high latitudes and large at intermediate latitudes in most of the analysed populations. Analyses of variance and principal components demonstrated that in most analysed populations B carrier males are associated with a decrease in body size related traits with respect to individuals with standard karyotype. Accordingly with the relationship between karyotype and body size, an opposite pattern of latitudinal variation in the frequencies of B′s with respect to body size variation was observed in this area. i.e. smaller individuals tend to have a higher frequency of B chromosomes. The comparison of the differentiation of both karyotype and body size traits with molecular neutral markers demonstrated the relative importance of selection moulding chromosome and phenotype variation. The observed pattern of phenotypic variation is likely to be the result of local adaptation to season length along the latitudinal gradient. The observed contrary pattern of B′s clinal variation may reflect the population ability to maintain this chromosome in relation to the local adaptation. The available evidence indicates that the distribution of B chromosome frequency was shaped by selective factors.     

Low variation in the polymorphic Clock gene poly-Q region despite population genetic structure across barn swallow (Hirundo rustica) populations

Recent studies of several species have reported a latitudinal cline in the circadian clock gene, Clock, which influences rhythms in both physiology and behavior. Latitudinal variation in this gene may hence reflect local adaptation to seasonal variation. In some bird populations, there is also an among-individual association between Clock poly-Q genotype and clutch initiation date and incubation period. We examined Clock poly-Q allele variation in the Barn Swallow (Hirundo rustica), a species with a cosmopolitan geographic distribution and considerable variation in life-history traits that may be influenced by the circadian clock. We genotyped Barn Swallows from five populations (from three subspecies) and compared variation at the Clock locus to that at microsatellite loci and mitochondrial DNA (mtDNA). We found very low variation in the Clock poly-Q region, as >96% of individuals were homozygous, and the two other alleles at this locus were globally rare. Genetic differentiation based on the Clock poly-Q locus was not correlated with genetic differentiation based on either microsatellite loci or mtDNA sequences. Our results show that high diversity in Clock poly-Q is not general across avian species. The low Clock variation in the background of heterogeneity in microsatellite and mtDNA loci in Barn Swallows may be an outcome of stabilizing selection on the Clock locus.     

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.     

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.