The evolutionary history of Drosophila buzzatii. XXXV. Inversion polymorphism and nucleotide variability in different regions of the second chromosome

Inversions are portions of a chromosome where the gene order is reversed relative to a standard reference orientation. Because of reduced levels of recombination in heterokaryotypes, inversions have a potentially important effect on patterns of nucleotide variability in those genomic regions close to, or included in, the inverted fragments. Here we report sequence variation at three anonymous regions (STSs) located at different positions in relation to second-chromosome inversion breakpoints in 29 isochromosomal lines derived from an Argentinean population of Drosophila buzzatii. In agreement with previous findings in Drosophila, gene flux (crossing over and/or gene conversion) between arrangements seems to appreciably increase as we approach the middle sections of inversion 2j, and patterns of nucleotide variability within, as well as genetic differentiation between chromosome arrangements, are comparable to those observed at the molecular marker outside the inverted fragments. On the other hand, nucleotide diversity near the proximal breakpoint of inversion 2j is reduced when contrasted with that found at the other regions, particularly in the case of derived inverted chromosomes. Using the data from the breakpoint, we estimate that the inversion polymorphism is approximately 1.63 N generations old, where N is the effective population size. An excess of low-frequency segregating polymorphisms is detected; mostly in the ancestral 2st arrangement and probably indicating a population expansion that predates the coalescent time of inversion 2j. Heterogeneity in mutation rates between the markers linked to the inversions may be sufficient to explain the different levels of nucleotide diversity observed. When considered in the context of other studies on patterns of variation relative to physical distance to inversion breakpoints, our data appear to be consistent with the conclusion that inversions are unlikely to be "long-lived" balanced polymorphisms.     

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

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

Contrasting Histories of Three Gene Regions Associated with In(3l)payne of Drosophila Melanogaster

In the present report, we studied nucleotide variation in three gene regions of Drosophila melanogaster, spanning >5 kb and showing different degrees of association with the cosmopolitan inversion In(3-L)Payne. The analysis of sequence variation in the regions surrounding the breakpoints and the heat shock 83 (Hsp83) gene locus, located close to the distal breakpoint, revealed the absence of shared polymorphisms and the presence of a number of fixed differences between arrangements, indicating absence of genetic exchange. In contrast, for the esterase-6 gene region, located in the center of the inversion, we observed the presence of shared polymorphisms between arrangements suggesting genetic exchange. In the regions close to the breakpoints, the common St arrangement is 10 times more polymorphic than inverted chromosomes. We propose that the lack of recombination between arrangements in these regions coupled with genetic hitchhiking is the best explanation for the low heterozygosity observed in inverted lines. Using the data for the breakpoints, we estimate that this inversion polymorphism is around 0.36 million yr old. Although it is widely accepted that inversions are examples of balanced polymorphisms, none of the current neutrality tests including our Monte Carlo simulations showed significant departure from neutral expectations.     

Interchromosomal exchange of genetic information between gene arrangements on the third chromosome of Drosophila pseudoobscura

During the last 60 years, the inversion polymorphism on the third chromosome of Drosophila pseudoobscura has become a case study of the evolution of linked blocks of genes, isolated from each other by the suppression of recombination in heterozygotes for different inversions. Due to its location within inverted regions in most gene arrangements, the amylase (Amy) gene region can be used to elucidate the molecular pattern of evolution in these inversions. We studied this region in the Tree Line phylad of gene arrangements, with regard to both restriction site polymorphisms (RSP) and nucleotide sequences. The analysis of restriction maps, encompassing 26 kb, corroborates the cytogenetic phylogeny established on the basis of inversion breakpoints. However, we found that the 2.7 kb of nucleotide sequences of the AmyI gene are identical in both Estes Park and Hidalgo arrangements, despite the fact that these inversions arose independently from Tree Line. These contrasting results suggest that a homogenizing force, most likely gene conversion, is able to bring about localized exchanges between otherwise isolated gene arrangements.      

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.     

Haplotype-based genomic sequencing of a chromosomal polymorphism in the white-throated sparrow (Zonotrichia albicollis)

Inversion polymorphisms have been linked to a variety of fundamental biological and evolutionary processes. Yet few studies have used large-scale genomic sequencing to directly compare the haplotypes associated with the standard and inverted chromosome arrangements. Here we describe the targeted genomic sequencing and comparison of haplotypes representing alternative arrangements of a common inversion polymorphism linked to a suite of phenotypes in the white-throated sparrow (Zonotrichia albicollis). More than 7.4 Mb of genomic sequence was generated and assembled from both the standard (ZAL2) and inverted (ZAL2(m)) arrangements. Sequencing of a pair of inversion breakpoints led to the identification of a ZAL2-specific segmental duplication, as well as evidence of breakpoint reusage. Comparison of the haplotype-based sequence assemblies revealed low genetic differentiation outside versus inside the inversion indicative of historical patterns of gene flow and suppressed recombination between ZAL2 and ZAL2(m). Finally, despite ZAL2(m) being maintained in a near constant state of heterozygosity, no signatures of genetic degeneration were detected on this chromosome. Overall, these results provide important insights into the genomic attributes of an inversion polymorphism linked to mate choice and variation in social behavior.     

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.     

Assaying chromosomal inversions by single-molecule haplotyping

Inversions are an important form of structural variation, but they are difficult to characterize, as their breakpoints often fall within inverted repeats. We have developed a method called 'haplotype fusion' in which an inversion breakpoint is genotyped by performing fusion PCR on single molecules of human genomic DNA. Fusing single-copy sequences bracketing an inversion breakpoint generates orientation-specific PCR products, exemplified by a genotyping assay for the int22 hemophilia A inversion on Xq28. Furthermore, we demonstrated that inversion events with breakpoints embedded within long (>100 kb) inverted repeats can be genotyped by haplotype-fusion PCR followed by bead-based single-molecule haplotyping on repeat-specific markers bracketing the inversion breakpoint. We illustrate this method by genotyping a Yp paracentric inversion sponsored by >300-kb-long inverted repeats. The generality of our methods to survey for, and genotype chromosomal inversions should help our understanding of the contribution of inversions to genomic variation, inherited diseases and cancer.     

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.     

Genetic flux between h1 and h2 haplotypes of the 17q21.31 inversion in European population

The chromosome 17q21.31 inversion is a 900-kb common structural polymorphism found primarily in European population.Although the genetic flux within inversion region was assumed to be considerable suppressed, it is still unclear about the details of genetic exchange between the H1(non-inverted sequence)and H2(inverted sequence) haplotypes of this inversion.Here we describe a refined map of genetic exchanges between pairs of gene arrangements within the 17q21.31 region.Using HapMap phase II data of 1,546 single nucleotide polymorphisms, we successfully deduced 96 H1 and 24 H2 haplotypes in European samples by neighbor-joining tree reconstruction. Furthermore, we identified 15 and 26 candidate tracts with reciprocal and non-reciprocal genetic exchanges, respectively.In all 15 regions harboring reciprocal exchange, haplotypes reconstructed by clone sequencing did not support these exchange events, suggesting that such signals of exchange between two sister chromosomes in certain heterozygous individual were caused by phasing error regions.On the other hand, the finished clone sequencing across 4 of 26 tracts with non-reciprocal genetic flux confirmed that this kind of genetic exchange was caused by gene conversion.In summary, as crossover between pairs of gene arrangements had been considerably suppressed, gene conversion might be the most important mechanism for genetic exchange at 17q21.31.     

Coalescent patterns for chromosomal inversions in divergent populations

Chromosomal inversions allow genetic divergence of locally adapted populations by reducing recombination between chromosomes with different arrangements. Divergence between populations (or hybridization between species) is expected to leave signatures in the neutral genetic diversity of the inverted region. Quantitative expectations for these patterns, however, have not been obtained. Here, we develop coalescent models of neutral sites linked to an inversion polymorphism in two locally adapted populations. We consider two scenarios of local adaptation: selection on the inversion breakpoints and selection on alleles inside the inversion. We find that ancient inversion polymorphisms cause genetic diversity to depart dramatically from neutral expectations. Other situations, however, lead to patterns that may be difficult to detect; important determinants are the age of the inversion and the rate of gene flux between arrangements. We also study inversions under genetic drift, finding that they produce patterns similar to locally adapted inversions of intermediate age. Our results are consistent with empirical observations, and provide the foundation for quantitative analyses of the roles that inversions have played in speciation.     

Effect of inversion polymorphism on the neutral nucleotide variability of linked chromosomal regions in Drosophila

Recombination is a main factor determining nucleotide variability in different regions of the genome. Chromosomal inversions, which are ubiquitous in the genus Drosophila, are known to reduce and redistribute recombination, and thus their specific effect on nucleotide variation may be of major importance as an explanatory factor for levels of DNA variation. Here, we use the coalescent approach to study this effect. First, we develop analytical expressions to predict nucleotide variability in old inversion polymorphisms that have reached mutation-drift-flux equilibrium. The effects on nucleotide variability of a new arrangement appearing in the population and reaching a stable polymorphism are then studied by computer simulation. We show that inversions modulate nucleotide variability in a complex way. The establishment of an inversion polymorphism involves a partial selective sweep that eliminates part of the variability in the population. This is followed by a slow convergence to the equilibrium values. During this convergence, regions close to the breakpoints exhibit much lower variability than central regions. However, at equilibrium, regions close to the breakpoints have higher levels of variability and differentiation between arrangements than regions in the middle of the inverted segment. The implications of these findings for overall variability levels during the evolution of Drosophila species are discussed.     

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.     

The structure and population genetics of the breakpoints associated with the cosmopolitan chromosomal inversion In(3R)Payne in Drosophila melanogaster

We report here the breakpoint structure and sequences of the Drosophila melanogaster cosmopolitan chromosomal inversion In(3R)P. Combining in situ hybridization to polytene chromosomes and long-range PCR, we have identified and sequenced the distal and proximal breakpoints. The breakpoints are not simple cut-and-paste structures; gene fragments and small duplications of DNA are associated with both breaks. The distal breakpoint breaks the tolkin (tok) gene and the proximal breakpoint breaks CG31279 and the tolloid (tld) gene. Functional copies of all three genes are found at the opposite breakpoints. We sequenced a representative sample of standard (St) and In(3R)P karyotypes for a 2-kb portion of the tok gene, as well as the same 2 kb from the pseudogene tok fragment found at the distal breakpoint of In(3R)P chromosomes. The tok gene in St arrangements possesses levels of polymorphism typical of D. melanogaster genes. The functional tok gene associated with In(3R)P shows little polymorphism. Numerous single-base changes, as well as deletions and duplications, are associated with the truncated copy of tok. The overall pattern of polymorphism is consistent with a recent origin of In(3R)P, on the order of Ne generations. The identification of these breakpoint sequences permits a simple PCR-based screen for In(3R)P.     

Pervasive gene conversion in chromosomal inversion heterozygotes

Chromosomal inversions shape recombination landscapes, and species differing by inversions may exhibit reduced gene flow in these regions of the genome. Though single crossovers within inversions are not usually recovered from inversion heterozygotes, the recombination barrier imposed by inversions is nuanced by noncrossover gene conversion. Here, we provide a genomewide empirical analysis of gene conversion rates both within species and in species hybrids. We estimate that gene conversion occurs at a rate of 1 × 10^–5 to 2.5 × 10^–5 converted sites per bp per generation in experimental crosses within Drosophila pseudoobscura and between D. pseudoobscura and its naturally hybridizing sister species D. persimilis. This analysis is the first direct empirical assessment of gene conversion rates within inversions of a species hybrid. Our data show that gene conversion rates in interspecies hybrids are at least as high as within‐species estimates of gene conversion rates, and gene conversion occurs regularly within and around inverted regions of species hybrids, even near inversion breakpoints. We also found that several gene conversion events appeared to be mitotic rather than meiotic in origin. Finally, we observed that gene conversion rates are higher in regions of lower local sequence divergence, yet our observed gene conversion rates in more divergent inverted regions were at least as high as in less divergent collinear regions. Given our observed high rates of gene conversion despite the sequence differentiation between species, especially in inverted regions, gene conversion has the potential to reduce the efficacy of inversions as barriers to recombination over evolutionary time.     

The foldback-like transposon Galileo is involved in the generation of two different natural chromosomal inversions of Drosophila buzzatii

Chromosomal inversions are the most common type of genome rearrangement in the genus Drosophila. Although the potential of transposable elements (TEs) for generating inversions has been repeatedly demonstrated in the laboratory, little is known on their role in the generation of natural inversions, which are those effectively contributing to the adaptation and/or evolution of species. We have cloned and sequenced the two breakpoints of the polymorphic inversion 2q7 of D. buzzatii. The sequence analysis of the breakpoint regions revealed the presence in the inverted chromosomes of large insertions, formed by complex assemblies of transposons, that are absent from the chromosomes without the inversion. Among the transposons inserted, the Foldback-like element Galileo, that was previously found responsible of the generation of the widespread inversion 2j of D. buzzatii, is present at both 2q7 breakpoints and is the most likely inducer of the inversion. A detailed study of the nucleotide and structural variation in the breakpoint regions of six chromosomal lines with the 2q7 inversion detected no nucleotide differences between them, which suggests a monophyletic and recent origin. In contrast, a remarkable degree of structural variation was observed in the same six chromosomal lines. It thus appears that the two breakpoints of the inverted chromosomes have become genetically unstable hotspots, as was previously found for the 2j inversion breakpoints. The possibility that this instability is caused by structural properties of Foldback elements is discussed.     

Unusual haplotype structure at the proximal breakpoint of In(2L)t in a natural population of Drosophila melanogaster.

The existence of temporally stable frequency clines for In(2L)t in natural populations of Drosophila melanogaster suggests a role for selection in the maintenance of this polymorphism. We have collected nucleotide polymorphism data from the proximal breakpoint junction regions of In(2L)t to infer its evolutionary history. The finding of a novel LINE-like element near the In(2L)t breakpoint junction in sampled inverted chromosomes supports a transposable element-mediated origin for this inversion. An analysis of nucleotide variation in a Costa Rican population sample of standard and inverted chromosomes indicates a unique and relatively recent origin for In(2L)t. Additional In(2L)t alleles from three geographically diverse populations reveal no detectable geographic differentiation. Low levels of In(2L)t nucleotide polymorphism suggest a recent increase in the inversion's frequency in tropical populations. An unusual feature of our sample of standard alleles is a marked heterogeneity in levels of linkage disequilibrium among polymorphic sites across the breakpoint region. We introduce a test of neutral equilibrium haplotype structure that corrects both for multiple tests and for an arbitrarily chosen window size. It reveals that an approximately 1.4-kb region immediately spanning the breakpoint has fewer haplotypes than expected under the neutral model, given the expected level of recombination in this genomic region. Certain features of our data suggest that the unusual pattern in standard chromosomes is the product of selection rather than demography.     

Transpecific polymorphisms in an inversion linked esterase locus in Drosophila buzzatii

Nucleotide variation was studied in a 1.1 kb section of the coding region of an Esterase gene (Est-A) that maps in the center of the segments rearranged by polymorphic inversions in the cactophilic Drosophila buzzatii. We examine 30 homozygous second-chromosome lines differing in gene arrangement and three D. koepferae isofemale lines as outgroups. Our data show that Est-A is a highly polymorphic gene at both synonymous and replacement sites. Significant departures from homogeneity in the distribution of the ratio of silent polymorphism to divergence predicted by the neutral theory reveals a local excess of silent polymorphism. This is consistent with the presence of two apparent narrow peaks of elevated silent polymorphism surrounding nonconservative amino acid substitutions. These polymorphisms as well as others at synonymous and nonsynonymous sites are shared with D. koepferae. We suggest that the presence of shared nucleotide polymorphisms is probably due to interspecific gene flow and/or balancing selection acting on replacement variants and/or to a decreased probability of loss of ancestral polymorphisms caused by linkage to an adaptive inversion polymorphism. Recurrent mutation and persistence of neutral ancestral polymorphisms cannot, however, be ruled out. The analysis of the distribution of nucleotide variation among the three chromosomal arrangements sampled reveals that derived arrangements (J and JZ(3)) are less polymorphic than the ancestral ST, and that the widely distributed ST and J arrangements are genetically differentiated. However, a significant number of polymorphisms are shared between arrangements, suggesting frequent exchange either from gene conversion or from double crossovers in heterokaryotypes. Finally, our present results in combination with data of sequence variation at the breakpoints of inversion J suggest that this old gene arrangement has risen in frequency in relatively recent times.     

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.