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.     

Population genetics of Drosophila ananassae: genetic differentiation among Indian natural populations at the level of inversion polymorphism

The present study, which is one of the longest temporal (two decades) and largest spatial (different parts of India covered) investigations on inversion polymorphism in natural populations of D. ananassae, was undertaken to understand the dynamics of inversion polymorphism in a broad and comprehensive manner. Forty-five natural populations from different ecogeographic regions of the country (covering the regions from Kashmir to Kanniyakumari and Gujarat to Nagaland) were analysed for chromosomal inversions. All the populations show the presence of the three cosmopolitan inversions, frequencies of which vary among the populations analysed. Simple correlations between frequencies of different inversions and regression analysis of inversion frequencies with latitude, longitude and altitude were insignificant. This reinforces the concept of rigid polymorphism in D. ananassae. Genetic divergence (spatial and temporal) at the level of chromosomal polymorphism among natural populations was calculated. Results show spatial divergence but no temporal divergence. Rigid polymorphic systems of D. ananassae did not show long-term directional trends. On the basis of the present study, and after including comparisons with the studies conducted more than two decades ago, the most important conclusion to be drawn is that the three cosmopolitan inversions in D. ananassae segregate within populations at fairly similar frequencies, and the general geographic pattern has remained constant.     

Refugial isolation versus ecological gradients

Hypotheses for divergence and speciation in rainforests generally fall into two categories: those emphasizing the role of geographic isolation and those emphasizing the role of divergent selection along gradients. While a majority of studies have attempted to infer mechanisms based on the pattern of species richness and congruence of geographic boundaries, relatively few have tried to simultaneously test alternative hypotheses for diversification. Here we discuss four examples, taken from our work on diversification of tropical rainforest vertebrates, in which we examine patterns of genetic and morphological variation within and between biogeographic regions to address two alternative hypotheses. By estimating morphological divergence between geographically contiguous and isolated populations under similar and different ecological conditions, we attempt to evaluate the relative roles of geographic isolation and natural selection in population divergence. Results suggest that natural selection, even in the presence of appreciable gene flow, can result in morphological divergence that is greater than that found between populations isolated for millions of years and, in some cases, even greater than that found between congeneric, but distinct, species. The relatively small phenotypic divergence that occurs among long-term geographic isolates in similar habitats suggests that morphological divergence via drift may be negligible and/or that selection is acting to produce similar phenotypes in populations occupying similar habitats. Our results demonstrate that significant phenotypic divergence: (1) is not necessarily coupled with divergence in neutral molecular markers; and (2) can occur without geographic isolation in the presence of gene flow.     

Population genetic structure,differentiation, and diversity in Tetrix subulata pygmy grasshoppers: roles of population size and immigration

Genetic diversity within and among populations and species is influenced by complex demographic and evolutionary processes. Despite extensive research, there is no consensus regarding how landscape structure, spatial distribution, gene flow, and population dynamics impact genetic composition of natural populations. Here, we used amplified fragment length polymorphisms (AFLPs) to investigate effects of population size, geographic isolation, immigration, and gene flow on genetic structure, divergence, and diversity in populations of Tetrix subulata pygmy grasshoppers (Orthoptera: Tetrigidae) from 20 sampling locations in southern Sweden. Analyses of 1564 AFLP markers revealed low to moderate levels of genetic diversity (PPL = 59.5–90.1; Hj = 0.23–0.32) within and significant divergence among sampling localities. This suggests that evolution of functional traits in response to divergent selection is possible and that gene flow is restricted. Genetic diversity increased with population size and with increasing proportion of long‐winged phenotypes (a proxy of recent immigration) across populations on the island of Öland, but not on the mainland. Our data further suggested that the open water separating Öland from the mainland acts as a dispersal barrier that restricts migration and leads to genetic divergence among regions. Isolation by distance was evident for short interpopulation distances on the mainland, but gradually disappeared as populations separated by longer distances were included. Results illustrate that integrating ecological and molecular data is key to identifying drivers of population genetic structure in natural populations. Our findings also underscore the importance of landscape structure and spatial sampling scheme for conclusions regarding the role of gene flow and isolation by distance.     

A Widespread Chromosomal Inversion Polymorphism Contributes to a Major Life-History Transition,Local Adaptation,and Reproductive Isolation

The role of chromosomal inversions in adaptation and speciation is controversial. Historically, inversions were thought to contribute to these processes either by directly causing hybrid sterility or by facilitating the maintenance of co-adapted gene complexes. Because inversions suppress recombination when heterozygous, a recently proposed local adaptation mechanism predicts that they will spread if they capture alleles at multiple loci involved in divergent adaptation to contrasting environments. Many empirical studies have found inversion polymorphisms linked to putatively adaptive phenotypes or distributed along environmental clines. However, direct involvement of an inversion in local adaptation and consequent ecological reproductive isolation has not to our knowledge been demonstrated in nature. In this study, we discovered that a chromosomal inversion polymorphism is geographically widespread, and we test the extent to which it contributes to adaptation and reproductive isolation under natural field conditions. Replicated crosses between the prezygotically reproductively isolated annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus, revealed that alternative chromosomal inversion arrangements are associated with life-history divergence over thousands of kilometers across North America. The inversion polymorphism affected adaptive flowering time divergence and other morphological traits in all replicated crosses between four pairs of annual and perennial populations. To determine if the inversion contributes to adaptation and reproductive isolation in natural populations, we conducted a novel reciprocal transplant experiment involving outbred lines, where alternative arrangements of the inversion were reciprocally introgressed into the genetic backgrounds of each ecotype. Our results demonstrate for the first time in nature the contribution of an inversion to adaptation, an annual/perennial life-history shift, and multiple reproductive isolating barriers. These results are consistent with the local adaptation mechanism being responsible for the distribution of the two inversion arrangements across the geographic range of M. guttatus and that locally adaptive inversion effects contribute directly to reproductive isolation. Such a mechanism may be partially responsible for the observation that closely related species often differ by multiple chromosomal rearrangements.     

Hawthorn-infesting populations of Rhagoletis pomonella in Mexico and speciation mode plurality

Categorizing speciation into dichotomous allopatric versus nonallopatric modes may not always adequately describe the geographic context of divergence for taxa. If some of the genetic changes generating inherent barriers to gene flow between populations evolved in geographic isolation, whereas others arose in sympatry, then the mode of divergence would be mixed. The apple maggot fly, Rhagoletis pomonella, has contributed to this emerging concept of a mixed speciation mode "plurality." Genetic studies have implied that a source of diapause life-history variation associated with inversions and contributing to sympatric host race formation and speciation for R. pomonella in the United States may have introgressed from the Eje Volcanico Trans Mexicano (EVTM; a.k.a. the Altiplano) in the past. A critical unresolved issue concerning the introgression hypothesis is how past gene flow occurred given the current 1200-km disjunction in the ranges of hawthorn-infesting flies in the EVTM region of Mexico and the southern extreme of the U.S. population in Texas. Here, we report the discovery of a hawthorn-infesting population of R. pomonella in the Sierra Madre Oriental Mountains (SMO) of Mexico. Sequence data from 15 nuclear loci and mitochondrial DNA imply that the SMO flies are related to, but still different from, U.S. and EVTM flies. The host affiliations, diapause characteristics, and phylogeography of the SMO population are consistent with it having served as a conduit for gene flow between Mexico and the United States. We also present evidence suggesting greater permeability of collinear versus rearranged regions of the genome to introgression, in accord with recent models of chromosomal speciation. We discuss the implications of the results in the context of speciation mode plurality. We do not argue for abandoning the terms sympatry or allopatry, but caution that categorizing divergence into either/or geographic modes may not describe the genetic origins of all species. For R. pomonella in the United States, the proximate selection pressures triggering race formation and speciation stem from sympatric host shifts. However, some of the phenological variation contributing to host-related ecological adaptation and reproductive isolation in sympatry at the present time appears to have an older history, having originated and become packaged into inversion polymorphism in allopatry.     

The relative influence of natural selection and geography on gene flow in guppies

Two general processes may influence gene flow among populations. One involves divergent selection, wherein the maladaptation of immigrants and hybrids impedes gene flow between ecological environments (i.e. ecological speciation). The other involves geographic features that limit dispersal. We determined the relative influence of these two processes in natural populations of Trinidadian guppies (Poecilia reticulata). If selection is important, gene flow should be reduced between different selective environments. If geography is important, gene flow should be impeded by geographic distance and physical barriers. We examined how genetic divergence, long-term gene flow, and contemporary dispersal within a watershed were influenced by waterfalls, geographic distance, predation, and habitat features. We found that waterfalls and geographic distance increased genetic divergence and reduced dispersal and long-term gene flow. Differences in predation or habitat features did not influence genetic divergence or gene flow. In contrast, differences in predation did appear to reduce contemporary dispersal. We suggest that the standard predictions of ecological speciation may be heavily nuanced by the mating behaviour and life history strategies of guppies.     

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.     

Multilevel analyses of genetic differentiation in Anopheles gambiae s.s. reveal patterns of gene flow important for malaria-fighting mosquito projects

Malaria control projects based on the introduction and spread of transgenes into mosquito populations depend on the extent of isolation between those populations. On the basis of the distribution of paracentric inversions, Anopheles gambiae has been subdivided into five subspecific chromosomal forms. Estimating gene flow between and within these forms of An. gambiae presents a number of challenges. We compared patterns of genetic divergence (F(ST)) between sympatric populations of the Bamako and Mopti forms at five sites. We used microsatellite loci within the j inversion on chromosome 2, which is fixed in the Bamako form but absent in the Mopti form, and microsatellites on chromosome 3, a region void of inversions. Estimates of genetic diversity and F(ST)'s suggest genetic exchanges between forms for the third chromosome but little for the j inversion. These results suggest a role for the inversion in speciation. Extensive gene flow within forms among sites resulted in populations clustering according to form despite substantial gene flow between forms. These patterns underscore the low levels of current gene flow between chromosomal forms in this area of sympatry. Introducing refractoriness genes in areas of the genome void of inversions may facilitate their spread within forms but their passage between forms may prove more difficult than previously thought.     

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.     

Gene flow and selection interact to promote adaptive divergence in regions of low recombination

Adaptation to new environments often occurs in the face of gene flow. Under these conditions, gene flow and recombination can impede adaptation by breaking down linkage disequilibrium between locally adapted alleles. Theory predicts that this decay can be halted or slowed if adaptive alleles are tightly linked in regions of low recombination, potentially favouring divergence and adaptive evolution in these regions over others. Here, we compiled a global genomic data set of over 1,300 individual threespine stickleback from 52 populations and compared the tendency for adaptive alleles to occur in regions of low recombination between populations that diverged with or without gene flow. In support of theory, we found that putatively adaptive alleles (F_ST and d_XY outliers) tend to occur more often in regions of low recombination in populations where divergent selection and gene flow have jointly occurred. This result remained significant when we employed different genomic window sizes, controlled for the effects of mutation rate and gene density, controlled for overall genetic differentiation, varied the genetic map used to estimate recombination and used a continuous (rather than discrete) measure of geographic distance as proxy for gene flow/shared ancestry. We argue that our study provides the first statistical evidence that the interaction of gene flow and selection biases divergence toward regions of low recombination.     

Divergence with gene flow in Anopheles funestus from the Sudan Savanna of Burkina Faso, West Africa

Anopheles funestus is a major vector of malaria across Africa. Understanding its complex and nonequilibrium population genetic structure is an important challenge that must be overcome before vector populations can be successfully perturbed for malaria control. Here we examine the role of chromosomal inversions in structuring genetic variation and facilitating divergence in Burkina Faso, West Africa, where two incipient species (chromosomal forms) of A. funestus, defined principally by rearrangements of chromosome 3R, have been hypothesized. Sampling across an approximately 300-km east-west transect largely contained within the Sudan-Savanna ecoclimatic zone, we analyzed chromosomal inversions, 16 microsatellite loci distributed genomewide, and 834 bp of the mtDNA ND5 gene. Both molecular markers revealed high genetic diversity, nearly all of which was accounted for by within-population differences among individuals, owing to recent population expansion. Across the study area there was no correlation between genetic and geographic distance. Significant genetic differentiation found between chromosomal forms on the basis of microsatellites was not genomewide but could be explained by chromosome 3R alone on the basis of loci inside and near inversions. These data are not compatible with complete reproductive isolation but are consistent with differential introgression and sympatric divergence between the chromosomal forms, facilitated by chromosome 3R inversions.     

Chromosome arm‐specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector,Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub‐Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine‐scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of F_ST, absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.     

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.     

Multiple chromosomal inversions contribute to adaptive divergence of a dune sunflower ecotype

Both models and case studies suggest that chromosomal inversions can facilitate adaptation and speciation in the presence of gene flow by suppressing recombination between locally adapted alleles. Until recently, however, it has been laborious and time‐consuming to identify and genotype inversions in natural populations. Here we apply RAD sequencing data and newly developed population genomic approaches to identify putative inversions that differentiate a sand dune ecotype of the prairie sunflower (Helianthus petiolaris) from populations found on the adjacent sand sheet. We detected seven large genomic regions that exhibit a different population structure than the rest of the genome and that vary in frequency between dune and nondune populations. These regions also show high linkage disequilibrium and high heterozygosity between, but not within, arrangements, consistent with the behaviour of large inversions, an inference subsequently validated in part by comparative genetic mapping. Genome–environment association analyses show that key environmental variables, including vegetation cover and soil nitrogen, are significantly associated with inversions. The inversions colocate with previously described “islands of differentiation,” and appear to play an important role in adaptive divergence and incipient speciation within H. petiolaris.     

The evolution of host preference in allopatric vs. parapatric populations of Timema cristinae walking-sticks

Divergent habitat preferences can contribute to speciation, as has been observed for host-plant preferences in phytophagous insects. Geographic variation in host preference can provide insight into the causes of preference evolution. For example, selection against maladaptive host-switching occurs only when multiple hosts are available in the local environment and can result in greater divergence in regions with multiple vs. a single host. Conversely, costs of finding a suitable host can select for preference even in populations using a single host. Some populations of Timema cristinae occur in regions with only one host-plant species present (in allopatry, surrounded by unsuitable hosts) whereas others occur in regions with two host-plant species adjacent to one another (in parapatry). Here, we use host choice and reciprocal-rearing experiments to document genetic divergence in host preference among 33 populations of T. cristinae. Populations feeding on Ceanothus exhibited a stronger preference for Ceanothus than did populations feeding on Adenostoma. Both allopatric and parapatric pairs of populations using the different hosts exhibited divergent host preferences, but the degree of divergence tended to be greater between allopatric pairs. Thus, gene flow between parapatric populations apparently constrains divergence. Host preferences led to levels of premating isolation between populations using alternate hosts that were comparable in magnitude to previously documented premating isolation caused by natural and sexual selection against migrants between hosts. Our findings demonstrate how gene flow and different forms of selection interact to determine the magnitude of reproductive isolation observed in nature.     

Genetic differentiation and inversion clines in Indian natural populations of Drosophila melanogaster.

To study the genetic differentiation and inversion clines in Indian natural populations of Drosophila melanogaster, 14 natural populations (6 from the north and 8 from the south) were screened for chromosome inversions. The chromosomal analysis revealed the presence of 23 paracentric inversions, which include 4 common cosmopolitan, 4 rare cosmopolitan, 2 recurrent endemic, and 13 unique endemic (new inversions detected for the first time) inversions. The difference in karyotype frequencies between populations from the north and south were highly significant and the level of inversion heterozygosity was higher in populations from the south. Statistically significant negative correlations were found between each of the four common cosmopolitan inversions and latitude. These findings are in accord with results from other worldwide geographic regions and show that Indian populations of D. melanogaster have undergone considerable genetic differentiation at the level of inversion polymorphism.     

Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow

Divergent natural selection, adaptive divergence and gene flow may interact in a number of ways. Recent studies have focused on the balance between selection and gene flow in natural populations, and empirical work has shown that gene flow can constrain adaptive divergence, and that divergent selection can constrain gene flow. A caveat is that phenotypic diversification may be under the direct influence of environmental factors (i.e. it may be due to phenotypic plasticity), in addition to partial genetic influence. In this case, phenotypic divergence may occur between populations despite high gene flow that imposes a constraint on genetic divergence. Plasticity may dampen the effects of natural selection by allowing individuals to rapidly adapt phenotypically to new conditions, thus slowing adaptive genetic divergence. On the other hand, plasticity may promote future adaptive divergence by allowing populations to persist in novel environments. Plasticity may promote gene flow between selective regimes by allowing dispersers to adapt to alternate conditions, or high gene flow may result in the selection for increased plasticity. Here I expand frameworks for understanding relationships among selection, adaptation and gene flow to include the effects of phenotypic plasticity in natural populations, and highlight its importance in evolutionary diversification.     

Population genetics of marine species: the interaction of natural selection and historically differentiated populations

High gene flow, particularly as mediated by larval dispersal, has usually been viewed as sufficient to limit geographic isolation as a major source of population differentiation among marine species. Despite the general observation of relatively little geographic variation among populations of high dispersal marine species many cases of divergence have been observed and natural selection has usually been invoked to explain geographic divergence. Detailed study of several allozyme polymorphisms provided additional evidence that selection may be the predominant force that determines genetic divergence in marine systems. There is, however, growing evidence that marine species with high dispersal are more subdivided than originally thought. The use of multi-locus approaches and the application of molecular techniques have provided new insight into the nature of population divergence in marine species. I argue that (1) many species, which were formerly thought to be unstructured, are in fact subdivided into genetically discrete groups, (2) it is often the case that genetically subdivided populations have distinct evolutionary histories, (3) in many cases, natural selection is the consequence of introgression between these groups, and (4) the combination of molecular assays of both nuclear and mitochondrial DNA and allozyme loci provides the best approach to understanding the evolutionary dynamics of these interacting populations.     

CHROMOSOMAL INVERSIONS AND SPECIES DIFFERENCES: WHEN ARE GENES AFFECTING ADAPTIVE DIVERGENCE AND REPRODUCTIVE ISOLATION EXPECTED TO RESIDE WITHIN INVERSIONS?

Many factors can promote speciation, and one which has received much attention is chromosomal inversions. A number of models propose that the recombination suppressing effects of inversions facilitate the maintenance of differences between interbreeding populations in genes affecting adaptive divergence and reproductive isolation. These models predict that such genes will disproportionately reside within inversions, rather than in collinear regions. This hypothesis has received some support, but exceptions exist. Additionally, the effects of known low levels of recombination within inversions on these models are uninvestigated. Here, simulations are used to compare the maintenance of genetic differences between populations following secondary contact and hybridization in different inversion models. We compare regions with no recombination within them to regions with low recombination and to collinear regions with free recombination. Our most general finding is that the low levels of recombination within an inversion often result in the loss of accentuated divergence in inverted regions compared to collinear ones. We conclude that inversions can facilitate the maintenance of species differences under some conditions, but that large or qualitative differences between inverted and collinear regions need not occur. We also find that strong selection facilitates maintenance of divergence in a manner analogous to inversions.