Testing alternative mechanisms of evolutionary divergence in an African rain forest passerine bird

Abstract Models of speciation in African rain forests have stressed either the role of isolation or ecological gradients. Here we contrast patterns of morphological and genetic divergence in parapatric and allopatric populations of the Little Greenbul, Andropadus virens, within different and similar habitats. We sampled 263 individuals from 18 sites and four different habitat types in Upper and Lower Guinea. We show that despite relatively high rates of gene flow among populations, A. virens has undergone significant morphological divergence across the savanna-forest ecotone and mountain-forest boundaries. These data support a central component of the divergence-with-gene-flow model of speciation by suggesting that despite large amounts of gene flow, selection is sufficiently intense to cause morphological divergence. Despite evidence of isolation based on neutral genetic markers, we find little evidence of morphological divergence in fitness-related traits between hypothesized refugial areas. Although genetic evidence suggests populations in Upper and Lower Guinea have been isolated for over 2 million years, morphological divergence appears to be driven more by habitat differences than geographic isolation and suggests that selection in parapatry may be more important than geographic isolation in causing adaptive divergence in morphology.     

Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.     

Parallel and nonparallel ecological,morphological and genetic divergence in lake–stream stickleback from a single catchment

Parallel phenotypic evolution in similar environments has been well studied in evolutionary biology; however, comparatively little is known about the influence of determinism and historical contingency on the nature, extent and generality of this divergence. Taking advantage of a novel system containing multiple lake–stream stickleback populations, we examined the extent of ecological, morphological and genetic divergence between three‐spined stickleback present in parapatric environments. Consistent with other lake–stream studies, we found a shift towards a deeper body and shorter gill rakers in stream fish. Morphological shifts were concurrent with changes in diet, indicated by both stable isotope and stomach contents analysis. Performing a multivariate test for shared and unique components of evolutionary response to the distance gradient from the lake, we found a strong signature of parallel adaptation. Nonparallel divergence was also present, attributable mainly to differences between river locations. We additionally found evidence of genetic substructuring across five lake–stream transitions, indicating that some level of reproductive isolation occurs between populations in these habitats. Strong correlations between pairwise measures of morphological, ecological and genetic distance between lake and stream populations supports the hypothesis that divergent natural selection between habitats drives adaptive divergence and reproductive isolation. Lake–stream stickleback divergence in Lough Neagh provides evidence for the deterministic role of selection and supports the hypothesis that parallel selection in similar environments may initiate parallel speciation.     

Morphological divergence within and between the Mexican sailfin mollies, Poecilia velifera and Poecilia petenensis

This study examined the patterns of morphological variation both between species, and between sexes and among populations within each species of the Mexican sailfin mollies, Poecilia velifera and Poecilia petenensis , using geometric morphometrics and linear measurements of morphological traits. While sexes within each species differed in characteristics that may be important in sexual selection, such as length of the dorsal fin, species differed in traits, such as body depth, that may also be influenced by natural selection due to differences in habitats. Within each species, many morphological traits were similar among populations, but important differences, including caudal peduncle depth in P. petenensis (but not in P. velifera ), suggested that habitat differences may also be important in shaping population divergence independently within each species. Indeed, the evolutionary vectors of male morphological population divergence for each species differed by an angle of 98·5°, representing nearly orthogonal vectors and suggesting independent shape divergence between these two molly species. Finally, geographic isolation did not explain the morphological differentiation seen among populations, suggesting that natural and sexual selection were strong forces promoting morphological diversification within these two species, despite the potential for a high degree of population connectivity and gene flow.     

Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as ‘isolation by adaptation’ (IBA). IBA can be complementary to the well‐known pattern of ‘isolation by distance’ (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome‐wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole‐genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome‐wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population‐level differentiation in this group of closely related populations.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

Divergence and reproductive isolation in the early stages of speciation

To understand speciation we need to identify the factors causing divergence between natural populations. The traditional approach to gaining such insights has been to focus on a particular theory and ask whether observed patterns of reproductive isolation between populations or species are consistent with the hypothesis in question. However, such studies are few and they do not allow us to compare between hypotheses, so often we cannot determine the relative contribution to divergence of different potential factors. Here, I describe a study of patterns of phenotypic divergence and premating and postmating isolation between populations of the grasshopper Chorthippus parallelus. Information on the phylogeographic relationships of the populations means that a priori predictions from existing hypotheses for the evolution of reproductive isolation can be compared with observations. I assess the relative contributions to premating isolation, postmating isolation and phenotypic divergence of long periods of allopatry, adaptation to different environments and processes associated with colonisation (such as population bottlenecks). Likelihood analysis reveals that long periods of allopatry in glacial refugia are associated with postmating reproductive isolation, but not premating isolation, which is more strongly associated with colonisation. Neither premating nor postmating isolation is higher between populations differing in potential environmental selection pressures. There are only weak correlations between patterns of genetic divergence and phenotypic divergence and no correlation between premating and postmating isolation. This suggests that the potential of a taxon to exercise mate choice may affect the types of factor that promote speciation in that group. I discuss the advantages and disadvantages of the general approach of simultaneously testing competing hypotheses for the evolution of reproductive isolation.     

Ecological reproductive isolation of coast and inland races of Mimulus guttatus

Adaptive divergence due to habitat differences is thought to play a major role in formation of new species. However it is rarely clear the extent to which individual reproductive isolating barriers related to habitat differentiation contribute to total isolation. Furthermore, it is often difficult to determine the specific environmental variables that drive the evolution of those ecological barriers, and the geographic scale at which habitat-mediated speciation occurs. Here, we address these questions through an analysis of the population structure and reproductive isolation between coastal perennial and inland annual forms of the yellow monkeyflower, Mimulus guttatus. We found substantial morphological and molecular genetic divergence among populations derived from coast and inland habitats. Reciprocal transplant experiments revealed nearly complete reproductive isolation between coast and inland populations mediated by selection against immigrants and flowering time differences, but not postzygotic isolation. Our results suggest that selection against immigrants is a function of adaptations to seasonal drought in inland habitat and to year round soil moisture and salt spray in coastal habitat. We conclude that the coast and inland populations collectively comprise distinct ecological races. Overall, this study suggests that adaptations to widespread habitats can lead to the formation of reproductively isolated species.     

Genetic structure is correlated with phenotypic divergence rather than geographic isolation in the highly polymorphic strawberry poison‐dart frog

Phenotypic and genetic divergence can be influenced by a variety of factors, including sexual and natural selection, genetic drift and geographic isolation. Investigating the roles of these factors in natural systems can provide insight into the relative influences of allopatric and ecological modes of biological diversification in nature. The strawberry poison frog, Dendrobates pumilio, presents an excellent opportunity for this kind of research, displaying a diverse array of colour morphs and inhabiting a heterogeneous landscape that includes oceanic islands, fragmented rainforest patches and wide expanses of suitable habitat. In this study, we use 15 highly polymorphic microsatellite loci to estimate population structure and gene flow among populations from across the range of D. pumilio and a causal modelling framework to statistically test 12 hypotheses regarding the geographic and phenotypic variables that explain genetic differentiation within this system. Our results demonstrate that the genetic distance between populations is most strongly associated with differences in dorsal coloration. Previous experimental studies have shown that phenotypic differences can result in sexual and natural selection against non‐native phenotypes, and our results now show that these forces lead to genetic isolation between different colour morphs in the wild, presenting a potential case of incipient speciation through selection.     

Parallels,nonparallels, and plasticity in population differentiation of threespine stickleback within a lake

The frequent occurrence of parallel phenotypic divergence in similar habitats is often evoked when emphasizing the role of ecology in adaptive radiation and speciation. However, because phenotypic plasticity can contribute to the observed pattern of divergence, confirmation of divergence at loci underlying phenotypic traits is important for confirming adaptive divergence. In the present study, we examine parallel morphological, neutral, and potentially adaptive genetic divergence of threespine stickleback inhabiting different habitats within a lake. Three genetic clusters best explained the neutral genetic structure within the lake; however, morphological differences were only weakly connected to genetic clusters and there was considerable phenotypic variation within clusters. Among the factors that could contribute to the observed pattern of morphological and genetic divergence are phenotypic plasticity, selective mortality of hybrids, and habitat choice based on morphology. Several loci are identified as outliers indicating divergent selection between the morphs and some parallels in morphological and adaptive genetic divergence are found in stickleback spawning at two lava sites. However, neutral genetic structure indicates considerable genetic connectivity among the two lava sites, and the parallels in morphology may therefore represent selective distribution of phenotypes rather than parallel divergence. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 803–813.     

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi‐interconnected ponds within a small geographic region (~7.5 km²) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.     

Ecological niche dimensionality and the evolutionary diversification of stick insects

The degree of phenotypic divergence and reproductive isolation between taxon pairs can vary quantitatively, and often increases as evolutionary divergence proceeds through various stages, from polymorphism to population differentiation, ecotype and race formation, speciation, and post-speciational divergence. Although divergent natural selection promotes divergence, it does not always result in strong differentiation. For example, divergent selection can fail to complete speciation, and distinct species pairs sometimes collapse ('speciation in reverse'). Widely-discussed explanations for this variability concern genetic architecture, and the geographic arrangement of populations. A less-explored possibility is that the degree of phenotypic and reproductive divergence between taxon pairs is positively related to the number of ecological niche dimensions (i.e., traits) subject to divergent selection. Some data supporting this idea stem from laboratory experimental evolution studies using Drosophila, but tests from nature are lacking. Here we report results from manipulative field experiments in natural populations of herbivorous Timema stick insects that are consistent with this 'niche dimensionality' hypothesis. In such insects, divergent selection between host plants might occur for cryptic colouration (camouflage to evade visual predation), physiology (to detoxify plant chemicals), or both of these niche dimensions. We show that divergent selection on the single niche dimension of cryptic colouration can result in ecotype formation and intermediate levels of phenotypic and reproductive divergence between populations feeding on different hosts. However, greater divergence between a species pair involved divergent selection on both niche dimensions. Although further replication of the trends reported here is required, the results suggest that dimensionality of selection may complement genetic and geographic explanations for the degree of diversification in nature.     

Consequences of divergence and introgression for speciation in Andean cloud forest birds

Divergence with gene flow is well documented and reveals the influence of ecological adaptation on speciation. Yet, it remains intuitive that gene exchange inhibits speciation in many scenarios, particularly among ecologically similar populations. The influence of gene flow on the divergence of populations facing similar selection pressures has received less empirical attention than scenarios where differentiation is coupled with local environmental adaptation. I used a paired study design to test the influence of genomic divergence and introgression on plumage differentiation between ecologically similar allopatric replacements of Andean cloud forest birds. Through analyses of short‐read genome‐wide sequences from over 160 individuals in 16 codistributed lineages, I found that plumage divergence is associated with deep genetic divergence, implicating a prominent role of geographic isolation in speciation. By contrast, lineages that lack plumage divergence across the same geographic barrier are more recently isolated or exhibit a signature of secondary genetic introgression, indicating a negative relationship between gene flow and divergence in phenotypic traits important to speciation. My results suggest that the evolutionary outcomes of cycles of isolation and divergence in this important theatre of biotic diversification are sensitive to time spent in the absence of gene flow.     

Population genomics reveals that refugial isolation and habitat change lead to incipient speciation in the Ground tit

Complex geological environments and Pleistocene glaciations contribute to speciation. For example, the unique geomorphological configuration and distinct geological evolution of the Qinghai–Tibet Plateau has led to the evolution of many endemic species adapted to a life at the plateau. Most of these endemic birds show no genetic divergence between different populations across the vast plateau region, with one particularly interesting exception. The ground tit (Pseudopodoces humilis) exhibits considerable genetic differences between the populations distributed on the central (platform) and those in the eastern margin (edge) of the plateau. Glacial refugial isolation and geographic barrier are hypotheses that have been proposed to explain this unique pattern of divergence. Here we use genome‐wide resequencing data to test these two hypotheses as explanation to the genetic divergence of the ground tit. The neighbor‐joining tree, principal component and STRUCTURE analyses suggest the presence of two incompletely diverged genetic lineages, which are geographically congruent with the populations in the platform and edge areas of the plateau. Our demographic model and population demography analyses estimated that the two populations diverged between 57 and 167 kya, and both populations experienced similar bottlenecks and expansions. Congruently, we also observed retracted suitable habitats for the two populations during the glacial time, suggesting that isolation in different glacial refugia. Interestingly, the genetic divergence between the two populations is further maintained and accumulated by habitat differentiation and morphological divergence. In accordance with this, we observed that highly divergent genomic regions harbour genes enriched in biological functions involved in muscle related processes and metabolic activities. Overall, the separation into different glacial refugia, the habitat difference and the morphological divergence contribute together to the adaptive population diversification in situ and we suggest that the platform and edge populations of the ground tit constitutes an interesting example of incipient speciation.     

Genetic population structure and call variation in a passerine bird, the satin bowerbird, Ptilonorhynchus violaceus

Geographic variation in vocalizations is widespread in passerine birds, but its origins and maintenance remain unclear. One hypothesis to explain this variation is that it is associated with geographic isolation among populations and therefore should follow a vicariant pattern similar to that typically found in neutral genetic markers. Alternatively, if environmental selection strongly influences vocalizations, then genetic divergence and vocal divergence may be disassociated. This study compared genetic divergence derived from 11 microsatellite markers with a metric of phenotypic divergence derived from male bower advertisement calls. Data were obtained from 16 populations throughout the entire distribution of the satin bowerbird, an Australian wet-forest-restricted passerine. There was no relationship between call divergence and genetic divergence, similar to most other studies on birds with learned vocalizations. Genetic divergence followed a vicariant model of evolution, with the differentiation of isolated populations and isolation-by-distance among continuous populations. Previous work on Ptilonorhynchus violaceus has shown that advertisement call structure is strongly influenced by the acoustic environment of different habitats. Divergence in vocalizations among genetically related populations in different habitats indicates that satin bowerbirds match their vocalizations to the environment in which they live, despite the homogenizing influence of gene flow. In combination with convergence of vocalizations among genetically divergent populations occurring in the same habitat, this shows the overriding importance that habitat-related selection can have on the establishment and maintenance of variation in vocalizations.     

A novel resource polymorphism in fish, driven by differential bottom environments: an example from an ancient lake in Japan

Divergent natural selection rooted in differential resource use can generate and maintain intraspecific eco-morphological divergence (i.e., resource polymorphism), ultimately leading to population splitting and speciation. Differing bottom environments create lake habitats with different benthos communities, which may cause selection in benthivorous fishes. Here, we document the nature of eco-morphological and genetic divergence among local populations of the Japanese gudgeon Sarcocheilichthys (Cyprinidae), which inhabits contrasting habitats in the littoral zones (rocky vs. pebbly habitats) in Lake Biwa, a representative ancient lake in East Asia. Eco-morphological analyses revealed that Sarcocheilichthys variegatus microoculus from rocky and pebbly zones differed in morphology and diet, and that populations from rocky environments had longer heads and deeper bodies, which are expected to be advantageous for capturing cryptic and/or attached prey in structurally complex, rocky habitats. Sarcocheilichthys biwaensis, a rock-dwelling specialist, exhibited similar morphologies to the sympatric congener, S. v. microoculus, except for body/fin coloration. Genetic analyses based on mitochondrial and nuclear microsatellite DNA data revealed no clear genetic differentiation among local populations within/between the gudgeon species. Although the morphogenetic factors that contribute to morphological divergence remain unclear, our results suggest that the gudgeon populations in Lake Biwa show a state of resource polymorphism associated with differences in the bottom environment. This is a novel example of resource polymorphism in fish within an Asian ancient lake, emphasizing the importance and generality of feeding adaptation as an evolutionary mechanism that generates morphological diversification.     

Geographic variation in genital morphology of Ciulfina praying mantids

Geographic variation in morphological traits is widespread and important to our current understanding of evolutionary processes. Although male genitalia are perhaps the most divergent morphological traits in animals, geographic variation in genital traits has received little attention and the mechanism driving such variation is unclear. The species isolation hypothesis of genital evolution makes explicit predictions about geographic variation in genitalia predicting patterns of genital divergence that reflect the risk of mating with related but incompatible species. The sexual selection and pleiotropy hypotheses, however, predict general levels of geographic variation that reflect divergent sexual selection pressures or genetic drift. To test these predictions, we investigated geographic variation in genital morphology in the praying mantid genus Ciulfina (Mantodea: Liturgusidae) using elliptic Fourier analysis. We found significant levels of geographic variation in the genital morphology of four Ciulfina species irrespective of the relative proximity of different populations to contact zones with other species. These results reject the species isolation hypothesis, and instead support either the sexual selection or pleiotropy hypotheses to explain patterns of genital evolution in this genus.     

Introgressive hybridization erodes morphological divergence between lentic and lotic habitats in an endangered minnow

Introgressive hybridization may erode phenotypic divergence along environmental gradients, collapsing locally adapted populations into a hybrid swarm. Alternatively, introgression may promote phenotypic divergence by providing variation on which natural selection can act. In freshwater fishes, water flow often selects for divergent morphological traits in lake versus stream habitats. We tested the effects of introgression on lake–stream morphological divergence in the minnow Owens Tui Chub (Siphateles bicolor snyderi), which has been rendered endangered by introgession from the introduced Lahontan Tui Chub (Siphateles bicolor obesa). Using geometric morphometric analysis of 457 individual Tui Chub from thirteen populations, we found that both native and introgressing parent taxa exhibited divergent body and caudal fin shapes in lake versus stream habitats, but their trajectories of divergence were distinct. In contrast, introgressed populations exhibited intermediate body and caudal fin shapes that were not differentiated by habitat type, indicating that introgression has eroded phenotypic divergence along the lentic–lotic gradient throughout the historic range of the Owens Tui Chub. Individuals within hybrid populations were less morphologically variable than those within parent populations, suggesting hybrid adaptation to selective agents other than water flow or loss of variance by drift.     

Testing evolutionary hypotheses for phenotypic divergence using landscape genetics

Understanding the evolutionary causes of phenotypic variation among populations has long been a central theme in evolutionary biology. Several factors can influence phenotypic divergence, including geographic isolation, genetic drift, divergent natural or sexual selection, and phenotypic plasticity. But the relative importance of these factors in generating phenotypic divergence in nature is still a tantalizing and unresolved problem in evolutionary biology. The origin and maintenance of phenotypic divergence is also at the root of many ongoing debates in evolutionary biology, such as the extent to which gene flow constrains adaptive divergence ( Garant et al. 2007 ) and the relative importance of genetic drift, natural selection, and sexual selection in initiating reproductive isolation and speciation ( Coyne & Orr 2004 ). In this issue, Wang & Summers (2010) test the causes of one of the most fantastic examples of phenotypic divergence in nature: colour pattern divergence among populations of the strawberry poison frog (Dendrobates pumilio) in Panama and Costa Rica ( Fig. 1 ). This study provides a beautiful example of the use of the emerging field of landscape genetics to differentiate among hypotheses for phenotypic divergence. Using landscape genetic analyses, Wang & Summers were able to reject the hypotheses that colour pattern divergence is due to isolation‐by‐distance (IBD) or landscape resistance. Instead, the hypothesis left standing is that colour divergence is due to divergent selection, in turn driving reproductive isolation among populations with different colour morphs. More generally, this study provides a wonderful example of how the emerging field of landscape genetics, which has primarily been applied to questions in conservation and ecology, now plays an essential role in evolutionary research.

Figure 1 Open in figure viewer PowerPoint Divergent colour morphs observed among populations of the strawberry poison frog, Dendrobates pumilio. Frogs are from San Cristobal (upper left), Cerro Brujo (upper right), Bastimentos (lower right), and Agua (lower left).     

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.