Joint effects of natural selection and recombination on gene flow between Drosophila ananassae populations

We estimated DNA sequence variation in a 5.7-kb fragment of the furrowed (fw) gene region within and between four populations of Drosophila ananassae; fw is located in a chromosomal region of very low recombination. We analyzed gene flow between these four populations along a latitudinal transect on the Indian subcontinent: two populations from southern, subtropical areas (Hyderabad, India, and Sri Lanka) and two from more temperate zones in the north (Nepal and Burma). Furthermore, we compared the pattern of differentiation at fw with published data from Om(1D), a gene located in a region of normal recombination. While differentiation at Om(1D) shows an isolation-by-distance effect, at fw the pattern of differentiation is quite different such that the frequencies of single nucleotide polymorphisms are homogenized over extended geographic regions (i.e., among the two populations of the northern species range from Burma and Nepal as well as among the two southern populations from India and Sri Lanka), but strongly differentiated between the northern and southern populations. To examine these differences in the patterns of variation and differentiation between the Om(1D) and fw gene regions, we determine the critical values of our previously proposed test of the background selection hypothesis (henceforth called F(ST) test). Using these results, we show that the pattern of differentiation at fw may be inconsistent with the background selection model. The data depart from this model in a direction that is compatible with the occurrence of recent selective sweeps in the northern as well as southern populations.     

Molecular Genetic Variation in the Centromeric Region of the X Chromosome in Three Drosophila Ananassae Populations. I. Contrasts between the Vermilion and Forked Loci

We have surveyed three natural populations of Drosophila ananassae for restriction map variation at the forked (f) and vermilion (v) loci, using 6-cutter restriction enzymes. Both loci are located in the centromeric region of the X chromosome. Two major conclusions can be drawn from the data. First, we found strong evidence for population subdivision, i.e., significant differences in the frequency distributions of polymorphisms and/or haplotypes between the Burma, India, and Brazil populations. Secondly, the pattern of DNA sequence variation between the two loci is unexpectedly different. The level of nucleotide variation in the v locus region is reduced (relative to f), especially in the Burma population. Furthermore, in contrast to v, we found no insertions/deletions larger than 700 bp and no significant linkage disequilibrium at f. The genetic differentiation among subpopulations can readily be attributed to restricted migration as the predominant evolutionary force. According to population genetics theory, the differences in DNA polymorphisms between the two loci are in qualitative agreement with the hypothesis that recombination is reduced in the v locus region ("centromere effect") but not at f. In order to test this hypothesis directly, we determined the cytogenetic positions of several loci in the centromeric region by in situ hybridization and found by comparison with the genetic map that recombination at v is indeed very low, much lower than at f.     

Analysis of a genetic hitchhiking model, and its application to DNA polymorphism data from Drosophila melanogaster

Begun and Aquadro have demonstrated that levels of nucleotide variation correlate with recombination rate among 20 gene regions from across the genome of Drosophila melanogaster. It has been suggested that this correlation results from genetic hitchhiking associated with the fixation of strongly selected mutants. The hitchhiking process can be described as a series of two-step events. The first step consists of a strongly selected substitution wiping out linked variation in a population; this is followed by a recovery period in which polymorphism can build up via neutral mutations and random genetic drift. Genetic hitchhiking has previously been modeled as a steady-state process driven by recurring selected substitutions. We show here that the characteristic parameter of this steady-state model is alpha v, the product of selection intensity (alpha = 2Ns) and the frequency of beneficial mutations v (where N is population size and s is the selective advantage of the favored allele). We also demonstrate that the steady-state model describes the hitchhiking process adequately, unless the recombination rate is very low. To estimate alpha v, we use the data of DNA sequence variation from 17 D. melanogaster loci from regions of intermediate to high recombination rates. We find that alpha v is likely to be > 1.3 x 10(-8). Additional data are needed to estimate this parameter more precisely. The estimation of alpha v is important, as this parameter determines the shape of the frequency distribution of strongly selected substitutions.      

Linkage disequilibria and the site frequency spectra in the su(s) and su(w(a)) regions of the Drosophila melanogaster X chromosome

Over the last decade, surveys of DNA sequence variation in natural populations of several Drosophila species and other taxa have established that polymorphism is reduced in genomic regions characterized by low rates of crossing over per physical length. Parallel studies have also established that divergence between species is not reduced in these same genomic regions, thus eliminating explanations that rely on a correlation between the rates of mutation and crossing over. Several theoretical models (directional hitchhiking, background selection, and random environment) have been proposed as population genetic explanations. In this study samples from an African population (n = 50) and a European population (n = 51) were surveyed at the su(s) (1955 bp) and su(w(a)) (3213 bp) loci for DNA sequence polymorphism, utilizing a stratified SSCP/DNA sequencing protocol. These loci are located near the telomere of the X chromosome, in a region of reduced crossing over per physical length, and exhibit a significant reduction in DNA sequence polymorphism. Unlike most previously surveyed, these loci reveal substantial skews toward rare site frequencies, consistent with the predictions of directional hitchhiking and random environment models and inconsistent with the general predictions of the background selection model (or neutral theory). No evidence for excess geographic differentiation at these loci is observed. Although linkage disequilibrium is observed between closely linked sites within these loci, many recombination events in the genealogy of the sampled alleles can be inferred and the genomic scale of linkage disequilibrium, measured in base pairs between sites, is the same as that observed for loci in regions of normal crossing over. We conclude that gene conversion must be high in these regions of low crossing over.     

Selection at linked sites in the partial selfer Caenorhabditis elegans

Natural selection can produce a correlation between local recombination rates and levels of neutral DNA polymorphism as a consequence of genetic hitchhiking and background selection. Theory suggests that selection at linked sites should affect patterns of neutral variation in partially selfing populations more dramatically than in outcrossing populations. However, empirical investigations of selection at linked sites have focused primarily on outcrossing species. To assess the potential role of selection as a determinant of neutral polymorphism in the context of partial self-fertilization, we conducted a multivariate analysis of single-nucleotide polymorphism (SNP) density throughout the genome of the nematode Caenorhabditis elegans. We based the analysis on a published SNP data set and partitioned the genome into windows to calculate SNP densities, recombination rates, and gene densities across all six chromosomes. Our analyses identify a strong, positive correlation between recombination rate and neutral polymorphism (as estimated by noncoding SNP density) across the genome of C. elegans. Furthermore, we find that levels of neutral polymorphism are lower in gene-dense regions than in gene-poor regions in some analyses. Analyses incorporating local estimates of divergence between C. elegans and C. briggsae indicate that a mutational explanation alone is unlikely to explain the observed patterns. Consequently, we interpret these findings as evidence that natural selection shapes genome-wide patterns of neutral polymorphism in C. elegans. Our study provides the first demonstration of such an effect in a partially selfing animal. Explicit models of genetic hitchhiking and background selection can each adequately describe the relationship between recombination rate and SNP density, but only when they incorporate selfing rate. Clarification of the relative roles of genetic hitchhiking and background selection in C. elegans awaits the development of specific theoretical predictions that account for partial self-fertilization and biased sex ratios.     

DNA sequence polymorphism and divergence at the erect wing and suppressor of sable loci of Drosophila melanogaster and D. simulans

Several evolutionary models of linked selection (e.g., genetic hitchhiking, background selection, and random environment) predict a reduction in polymorphism relative to divergence in genomic regions where the rate of crossing over per physical distance is restricted. We tested this prediction near the telomere of the Drosophila melanogaster and D. simulans X chromosome at two loci, erect wing (ewg) and suppressor of sable [su(s)]. Consistent with this prediction, polymorphism is reduced at both loci, while divergence is normal. The reduction is greater at ewg, the more distal of the two regions. Two models can be discriminated by comparing the observed site frequency spectra with those predicted by the models. The hitchhiking model predicts a skew toward rare variants in a sample, while the spectra under the background-selection model are similar to those of the neutral model of molecular evolution. Statistical tests of the fit to the predictions of these models require many sampled alleles and segregating sites. Thus we used SSCP and stratified DNA sequencing to cover a large number of randomly sampled alleles (approximately 50) from each of three populations. The result is a clear trend toward negative values of Tajima's D, indicating an excess of rare variants at ewg, the more distal of the two loci. One fixed difference among the populations and high FST values indicate strong population subdivision among the three populations at ewg. These results indicate genetic hitchhiking at ewg, in particular, geographically localized hitchhiking events within Africa. The reduction of polymorphism at su(s) combined with the excess of high-frequency variants in D. simulans is inconsistent with the hitchhiking and background-selection models.     

Contrasting patterns of nucleotide polymorphism at the alcohol dehydrogenase locus in the outcrossing Arabidopsis lyrata and the selfing Arabidopsis thaliana

Nucleotide variation at the alcohol dehydrogenase locus (Adh) was studied in the outcrossing Arabidopsis lyrata, a close relative of the selfing Arabidopsis thaliana. Overall, estimated nucleotide diversity in the North American ssp. lyrata and two European ssp. petraea populations was 0.0038, lower than the corresponding specieswide estimate for A. thaliana at the same set of nucleotide sites. The distribution of segregating sites across the gene differed between the two species. Estimated sequence diversity within an A. lyrata population with a large sample size (0.0023) was much higher than has previously been observed for A. thaliana. This North American population has an excess of sites at intermediate frequencies compared with neutral expectation (Tajima's D = 2.3, P < 0.005), suggestive of linked balancing selection or a recent population bottleneck. In contrast, an excess of rare polymorphisms has been found in A. thaliana. Polymorphism within A. lyrata and divergence from A. thaliana appear to be correlated across the Adh gene sequence. The geographic distribution of polymorphism was quite different from that of A. thaliana, for which earlier studies of several genes found low within-population nucleotide site polymorphism and no overall continental differentiation of variation despite large differences in site frequencies between local populations. Differences between the outcrossing A. lyrata and the selfing A. thaliana reflect the impact of differences in mating system and the influence of bottlenecks in A. thaliana during rapid colonization on DNA sequence polymorphism. The influence of additional variability-reducing mechanisms, such as background selection or hitchhiking, may not be discernible.     

A metapopulation perspective on genetic diversity and differentiation in partially self-fertilizing plants

Abstract.— Partial self-fertilization is common in higher plants. Mating system variation is known to have important consequences for how genetic variation is distributed within and among populations. Selfing is known to reduce effective population size, and inbreeding species are therefore expected to have lower levels of genetic variation than comparable out crossing taxa. However, several recent empirical studies have shown that reductions in genetic diversity within populations of inbreeding species are far greater than the expected reductions based on the reduced effective population size. Two different processes have been argued to cause these patterns, selective sweeps (or hitchhiking) and background selection. Both are expected to be most effective in reducing genetic variation in regions of low recombination rates. Selfing is known to reduce the effective recombination rate, and inbreeding taxa are thus thought to be particularly vulnerable to the effects of hitchhiking or background selection. Here I propose a third explanation for the lower-than-expected levels of genetic diversity within populations of selfing species; recurrent extinctions and recolonizations of local populations, also known as metapopulation dynamics. I show that selfing in a metapopulation setting can result in large reductions in genetic diversity within populations, far greater than expected based the lower effective population size inbreeding species is expected to have. The reason for this depends on an interaction between selfing and pollen migration.     

Molecular Population Genetics of the Distal Portion of the X Chromosome in Drosophila: Evidence for Genetic Hitchhiking of the Yellow-Achaete Region

We have estimated DNA sequence variation and differentiation within and between Drosophila melanogaster and its sibling species, Drosophila simulans, using six-cutter restriction site variation at yellow-achaete (y-ac), phosphogluconate dehydrogenase (Pgd), and period (per). These three gene regions are of varying distance from the telomere of the X chromosome and range from very low to moderate rates of recombination in D. melanogaster. According to Tajima's test of neutrality, the Pgd region has been influenced by balancing selection in D. melanogaster. This is consistent with previous data suggesting the allozyme polymorphism at this locus is visible to selection. The Hudson, Kreitman, Aguadé test of neutrality reveals a significant departure from neutrality for the y-ac region compared to the per or rosy regions in D. simulans. There is also a significant departure for the y-ac region compared to the Adh 5' flanking region in D. melanogaster. In both species the departure appears to be due to reduced variation at y-ac compared to that expected from divergence between D. simulans and D. melanogaster. We conclude that recent hitchhiking associated with the selective fixation of one or more advantageous mutants in the y-ac region is the best explanation for reduced variation at y-ac.     

Characterisation of Ascaris from human and pig hosts by nuclear ribosomal DNA sequences.

The sequences of the nuclear ribosomal DNA region spanning the first internal transcribed spacer, the 5.8S rRNA gene and the second internal transcribed spacer were determined for Ascaris samples from pigs and humans from different geographical regions. The sequences of the 5.8S gene and the second internal transcribed spacer were the same for all samples examined, whereas all Ascaris samples from humans had six (1.3%) nucleotide differences in the first internal transcribed spacer compared with those from pigs. These differences provided some support for the existence of separate species of Ascaris or population variation within this genus. Using a nucleotide difference within a site for the restriction enzyme HaeIII, a PCR-linked restriction fragment length polymorphism method was established which allowed the delineation of the Ascaris samples from pigs and humans used herein. Exploiting the sequence differences in the first internal transcribed spacer, a PCR-based single-strand conformation polymorphism method was established for future analysis of the genetic structure of pig and human Ascaris populations in sympatric and allopatric zones.     

Highly differentiated population structure of a Mangrove species, Bruguiera gymnorhiza (Rhizophoraceae) revealed by one nuclear GapCp and one chloroplast intergenic spacer trnF–trnL

To evaluate the genetic diversity of a mangrove species and clarify the genetic structure of its populations, we studied nucleotide polymorphism in two DNA regions of Bruguiera gymnorhiza collected from the southern islands of Japan, Thailand, Malaysia, Indonesia, Micronesia, and India. The two DNA sequences were the chloroplast (cp) intergenic spacer between trnL and trnF genes (ca. 300 bp), and a part (ca. 550 bp) of the nuclear gene coding for glyceraldehyde-3-phosphate dehydrogenase (GapCp). Little polymorphism was found within each of the three geographical regions, Pacific Ocean, Bay of Bengal and Arabian Sea. Throughout the vast regions east of the Malay peninsula including Indonesia, Thailand, Micronesia and the southern islands of Japan (Pacific Ocean), essentially only one haplotype (apart from variation in number of a T repeat) was present. A second haplotype was present on the western coast of Malay Peninsula and the eastern coast of India (Bay of Bengal). On the southwest of Malay Peninsula both of these haplotypes were present. Finally a third haplotype was found only on the western coast of India (Arabian Sea). When taken over all geographic populations, total nucleotide variation within the species was large (μ = 0.006, average of the two genes). Our results are consistent with the hypothesis that this low genetic diversity within any local population and differentiation between the different oceans or regions are caused by very low gene flow between each of the different oceans coupled with frequent fluctuation of population sizes due to the change in sea level. The significance of these results is discussed from evolutionary point of the mangrove forests.     

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.     

Genetic diversity and condition factor: a significant relationship in Flemish but not in German populations of the European bullhead (Cottus gobio L.)

The population genetics of the Antarctic neritic krill species Euphausia crystallorophias was examined by nucleotide sequence variation in its mitochondrial DNA. A 616 base pair region of the cytochrome c oxidase subunit I (COI) gene was screened for mutations by single-strand conformational polymorphism (SSCP) combined with restriction digestion. E. crystallorophias caught in three different regions of the Antarctic coastline were used--two samples from the Mertz Glacier Polynya and one sample each from the western side of the Antarctic Peninsula and from the Davis Sea. Significant genetic differences between krill samples were identified. However, the extent of these differences did not correlate with the degree of geographic separation between the sampling sites. This suggests that the genetic structuring may be the result of small-scale differentiation rather than differentiation between resident populations in separate parts of the Southern Ocean. The possibility that genetic differences between samples within a region are as important as differences between regions has implications for other studies of krill population genetics.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

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.     

Relationship between DNA Polymorphism and Fixation Time

When there is no recombination among nucleotide sites in DNA sequences, DNA polymorphism and fixation of mutants at nucleotide sites are mutually related. Using the method of gene genealogy, the relationship between the DNA polymorphism and the fixation of mutant nucleotide was quantitatively investigated under the assumption that mutants are selectively neutral, that there is no recombination among nucleotide sites, and that the population is a random mating population with N diploid individuals. The results obtained indicate that the expected number of nucleotide differences between two DNA sequences randomly sampled from the population is 42% less when a mutant at a particular nucleotide site reaches fixation than at a random time, and that heterozygosity is also expected to be less when fixation takes place than at a random time, but the amount of reduction depends on the value of 4Nv in this case, where v is the mutation rate per DNA sequence per generation. The formula for obtaining the expected number of nucleotide differences between the two DNA sequences for a given fixation time is also derived, and indicates that, even when it takes a large number of generations for a mutant to reach fixation, this number is 33% less than at a random time. The computer simulation conducted suggests that the expected number of nucleotide differences between the two DNA sequences at the time when an advantageous mutant becomes fixed is essentially the same as that of neutral mutant if the fixation time is the same. The effect of recombination on the amount of DNA polymorphism was also investigated by using computer simulation.     

Testing for effects of recombination rate on nucleotide diversity in natural populations of Arabidopsis lyrata

We investigated DNA sequence diversity for loci on chromosomes 1 and 2 in six natural populations of Arabidopsis lyrata and tested for the role of natural selection in structuring genomewide patterns of variability, specifically examining the effects of recombination rate on levels of silent polymorphism. In contrast with theoretical predictions from models of genetic hitchhiking, maximum-likelihood-based analyses of diversity and divergence do not suggest reduction of diversity in the region of suppressed recombination near the centromere of chromosome 1, except in a single population from Russia, in which the pericentromeric region may have undergone a local selective sweep or demographic process that reduced variability. We discuss various possibilities that might explain why nucleotide diversity in most A. lyrata populations is not related to recombination rate, including genic recombination hotspots, and low gene density in the low recombination rate region.     

Evolution at the tip and base of the X chromosome in an African population of Drosophila melanogaster

Hitchhiking effects of advantageous mutations have been invoked to explain reduced polymorphism in regions of low crossing-over in Drosophila. Besides reducing DNA heterozygosity, hitchhiking effects should produce strong linkage disequilibrium and a frequency spectrum skewed toward an excess of rare polymorphisms (compared to the neutral expectation). We measured DNA polymorphism in a Zimbabwe population of D. melanogaster at three loci, yellow, achaete, and suppressor of forked, located in regions of reduced crossing-over. Similar to previously published surveys of these genomic regions in other populations, we observed low levels of nucleotide variability. However, the frequency spectrum was compatible with a neutral model, and there was abundant evidence for recombination in the history of the yellow and ac genes. Thus, some aspects of the data cannot be accounted for by a simple hitchhiking model. An alternative hypothesis, background selection, might be compatible with the observed patterns of linkage disequilibrium and the frequency spectrum. However, this model cannot account for the observed reduction in nucleotide heterozygosity. Thus, there is currently no satisfactory theoretical model for the data from the tip and base of the X chromosome in D. melanogaster.      

Nucleotide variation at the dopa decarboxylase (<Emphasis Type="Italic">Ddc</Emphasis>) gene in natural populations of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

We studied nucleotide sequence variation at the gene coding for dopa decarboxylase (Ddc) in seven populations of Drosophila melanogaster. Strength and pattern of linkage disequilibrium are somewhat distinct in the extensively sampled Spanish and Raleigh populations. In the Spanish population, a few sites are in strong positive association, whereas a large number of sites in the Raleigh population are associated nonrandomly but the association is not strong. Linkage disequilibrium analysis shows presence of two groups of haplotypes in the populations, each of which is fairly diverged, suggesting epistasis or inversion polymorphism. There is evidence of two forms of natural selection acting on Ddc. The McDonald-Kreitman test indicates a deficit of fixed amino acid differences between D. melanogaster and D. simulans, which may be due to negative selection. An excess of derived alleles at high frequency, significant according to the H-test, is consistent with the effect of hitchhiking. The hitchhiking may have been caused by directional selection downstream of the locus studied, as suggested by a gradual decrease of the polymorphism-to-divergence ratio. Altogether, the Ddc locus exhibits a complicated pattern of variation apparently due to several evolutionary forces. Such a complex pattern may be a result of an unusually high density of functionally important genes.     

Divergence between the Drosophila pseudoobscura and D. persimilis genome sequences in relation to chromosomal inversions

As whole-genome sequence assemblies accumulate, a challenge is to determine how these can be used to address fundamental evolutionary questions, such as inferring the process of speciation. Here, we use the sequence assemblies of Drosophila pseudoobscura and D. persimilis to test hypotheses regarding divergence with gene flow. We observe low differentiation between the two genome sequences in pericentromeric and peritelomeric regions. We interpret this result as primarily a remnant of the correlation between levels of variation and local recombination rate observed within populations. However, we also observe lower differentiation far from the fixed chromosomal inversions distinguishing these species and greater differentiation within and near these inversions. This finding is consistent with models suggesting that chromosomal inversions facilitate species divergence despite interspecies gene flow. We also document heterogeneity among the inverted regions in their degree of differentiation, suggesting temporal differences in the origin of each inverted region consistent with the inversions arising during a process of divergence with gene flow. While this study provides insights into the speciation process using two single-genome sequences, it was informed by lower throughput but more rigorous examinations of polymorphism and divergence. This reliance highlights the need for complementary genomic and population genetic approaches for tackling fundamental evolutionary questions such as speciation.