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).     

Stark sexual display divergence among jumping spider populations in the face of gene flow

Gene flow can inhibit evolutionary divergence by eroding genetic differences between populations. A current aim in speciation research is to identify conditions in which selection overcomes this process. We focused on a state of limited differentiation, asking whether selection enables divergence with gene flow in a set of Habronattus americanus jumping spider populations that exhibit three distinct male sexual display morphs. We found that each population is at high frequency or fixed for a single morph. These strong phenotypic differences contrast with low divergence at 210 AFLP markers, suggesting selection has driven or maintains morph divergence. Coinciding patterns of isolation by distance and ‘isolation by phenotype’ (i.e. increased genetic divergence among phenotypically contrasting populations) across the study area support several alternative demographic hypotheses for display divergence, each of which entails gene flow. Display‐associated structure appears broadly distributed across the genome and the markers producing this pattern do not stand out from background levels of differentiation. Overall, the results suggest selection can promote stark sexual display divergence in the face of gene flow among closely related populations.     

Variation in spectral reflectance among populations of Dendrobates pumilio, the strawberry poison frog, in the Bocas del Toro Archipelago, Panama

Aim The goal of this study was to quantify levels of variation in spectral reflectance within and among populations of Dendrobates pumilio, the strawberry poison frog from the Bocas del Toro Archipelago. Location This study was carried out in the Bocas del Toro Archipelago, in the Republic of Panama, Central America. Methods Spectral reflectance was measured for samples of individuals from fifteen distinct island and mainland populations, using an Ocean Optics 2000 spectrometer and a BiLink portable computer. Results Our results provide quantitative evidence for extreme polymorphism among populations, and more limited levels of polymorphism within some populations. No obvious signs of sexual dimorphism were found. All the colour morphs appear to have relatively little reflectance in the ultraviolet part of the spectrum. There is some evidence for clinal variation in colour and pattern across some mainland populations. There is also at least one area where distinctly different morphs occur in sympatry, suggesting the possibility of incipient reproductive isolation. We argue that variation in coloration may have been enhanced by sexual selection.     

Flickers of speciation: Sympatric colour morphs of the arc‐eye hawkfish,Paracirrhites arcatus,reveal key elements of divergence with gene flow

One of the primary challenges of evolutionary research is to identify ecological factors that favour reproductive isolation. Therefore, studying partially isolated taxa has the potential to provide novel insight into the mechanisms of evolutionary divergence. Our study utilizes an adaptive colour polymorphism in the arc‐eye hawkfish (Paracirrhites arcatus) to explore the evolution of reproductive barriers in the absence of geographic isolation. Dark and light morphs are ecologically partitioned into basaltic and coral microhabitats a few metres apart. To test whether ecological barriers have reduced gene flow among dark and light phenotypes, we evaluated genetic variation at 30 microsatellite loci and a nuclear exon (Mc1r) associated with melanistic coloration. We report low, but significant microsatellite differentiation among colour morphs and stronger divergence in the coding region of Mc1r indicating signatures of selection. Critically, we observed greater genetic divergence between colour morphs on the same reefs than that between the same morphs in different geographic locations. We hypothesize that adaptation to the contrasting microhabitats is overriding gene flow and is responsible for the partial reproductive isolation observed between sympatric colour morphs. Combined with complementary studies of hawkfish ecology and behaviour, these genetic results indicate an ecological barrier to gene flow initiated by habitat selection and enhanced by assortative mating. Hence, the arc‐eye hawkfish fulfil theoretical expectations for the earliest phase of speciation with gene flow.     

Multivariate phenotypes and the potential for alternative phenotypic optima in wall lizard (Podarcis muralis) ventral colour morphs

A major goal in evolutionary biology is to determine how phenotypic variation arises and is maintained in natural populations. Recent studies examining the morphological, physiological and behavioural differences among discrete colour morphotypes (morphs) have revealed several mechanisms that maintain discrete variation within populations, including frequency‐dependence, density‐dependence and correlational selection. For example, trade‐offs over resource allocation to morphological, physiological and behavioural traits can drive correlational selection for morph‐specific phenotypic optima. Here, we describe a ventral colour polymorphism in the wall lizard (Podarcis muralis) and test the hypothesis that morphs differ along multivariate axes defined by trade‐offs in morphological, physiological, and immunological traits. We show that ventral colour is a discrete trait and that morphs differ in body size, prevalence of infection by parasites and infection intensity. We also find that morphs differ along multivariate phenotypic axes and experience different multivariate selection pressures. Our results suggest that multivariate selection pressures may favour alternative optimal morph‐specific phenotypes in P. muralis.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

Evidence for selection on coloration in a Panamanian poison frog: a coalescent‐based approach

Aim The strawberry poison frog, Oophaga pumilio, has undergone a remarkable radiation of colour morphs in the Bocas del Toro archipelago in Panama. This species shows extreme variation in colour and pattern between populations that have been geographically isolated for < 10,000 years. While previous research has suggested the involvement of divergent selection, to date no quantitative test has examined this hypothesis. Location Bocas del Toro archipelago, Panama. Methods We use a combination of population genetics, phylogeography and phenotypic analyses to test for divergent selection in coloration in O. pumilio. Tissue samples of 88 individuals from 15 distinct populations were collected. Using these data, we developed a gene tree using the mitochondrial DNA (mtDNA) d‐loop region. Using parameters derived from our mtDNA phylogeny, we predicted the coalescence of a hypothetical nuclear gene underlying coloration. We collected spectral reflectance and body size measurements on 94 individuals from four of the populations and performed a quantitative analysis of phenotypic divergence. Results The mtDNA d‐loop tree revealed considerable polyphyly across populations. Coalescent reconstructions of gene trees within population trees revealed incomplete genotypic sorting among populations. The quantitative analysis of phenotypic divergence revealed complete lineage sorting by colour, but not by body size: populations showed non‐overlapping variation in spectral reflectance measures of body coloration, while variation in body size did not separate populations. Simulations of the coalescent using parameter values derived from our empirical analyses demonstrated that the level of sorting among populations seen in colour cannot reasonably be attributed to drift. Main conclusions These results imply that divergence in colour, but not body size, is occurring at a faster rate than expected under neutral processes. Our study provides the first quantitative support for the claim that strong diversifying selection underlies colour variation in the strawberry poison frog.     

Genetic divergence and assortative mating between colour morphs of the sea urchin Paracentrotus gaimardi

Some species of sea urchins feature large variation in pigmentation. This variability may be the result of phenotypic plasticity or it may be associated with genetic divergence between morphs. Paracentrotus gaimardi exhibits five colour morphs (pink, brown, green, grey and black), which often occur side by side on the same rock. We studied genetic divergence between these morphs in three populations on the coast of Brazil. A fragment of the region encoding the mitochondrial ATPase 8 and 6 mitochondrial genes, a fragment of the intron of a nuclear histone and the entire nuclear gene coding for the sperm protein bindin were analysed. Mitochondrial DNA was differentiated between the pink and all other morphs, but the histone intron was similar in all colour morphs. In bindin, nine codons were found to be under positive selection and significant differences of allelic frequencies were observed in almost all pairwise comparisons between colour morphs. Although the molecular differentiation in bindin is not large enough to suggest reproductive isolation, some degree of assortative mating within morphs seems to be occurring in this species.     

Correlated morphological and colour differences among females of the damselfly Ischnura elegans

Abstract 1. The female‐limited colour polymorphic damselfly Ischnura elegans has proven to be an interesting study organism both as an example of female sexual polymorphism, and in the context of the evolution of colour polymorphism, as a model of speciation processes. 2. Previous research suggests the existence of correlations between colour morph and other phenotypic traits, and the different female morphs in I. elegans may be pursuing alternative phenotypically integrated strategies. However, previous research on morphological differences in southern Swedish individuals of this species was only carried out on laboratory‐raised offspring from a single population, leaving open the question of how widespread such differences are. 3. The present study therefore analysed multi‐generational data from 12 populations, investigating morphological differences between the female morphs in the field, differences in the pattern of phenotypic integration between morphs, and quantified selection on morphological traits. 4. It was found that consistent morphological differences indeed existed between the morphs across populations, confirming that the previously observed differences were not simply a laboratory artefact. It was also found, somewhat surprisingly, that despite the existence of sexual dimorphism in body size and shape, patterns of phenotypic integration differed most between the morphs and not between the sexes. Finally, linear selection gradients showed that female morphology affected fecundity differently between the morphs. 5. We discuss the relevance of these results to the male mimicry hypothesis and to the existence of potential ecological differences between the morphs.     

Male courtship preferences demonstrate discrimination against allopatric colour morphs in a cichlid fish

Whether premating isolation is achieved by male‐specific, female‐specific or sex‐independent assortative preferences often depends on the underlying evolutionary processes. Here we test mate preferences of males presented with females of different allopatric colour variants of the cichlid fish Tropheus sp., a Lake Tanganyika endemic with rich geographical colour pattern variation, in which the strength of sexual isolation varies between populations. We conducted two‐way mate choice experiments to compare behaviour of males of a red‐bodied morph (population Moliro) towards females from their own population with behaviour towards females from four allopatric populations at different stages of phylogenetic and phenotypic divergence. Males courted same‐population females significantly more intensely than females of other populations, and reduced their heteromorphic courtship efforts both with increasing genetic and increasing phenotypic distinctness of the females. In particular, females of a closely related red‐bodied population received significantly more courtship than either genetically distinct, similarly coloured females (‘Kirschfleck’ morph) or genetically related, differently coloured females (‘yellow‐blotch’ morph), both of which were courted similarly. Genetically and phenotypically distinct females (Tropheus polli) were not courted at all. Consistent with previous female‐choice experiments, female courtship activity also decreased with increasing genetic distance from the males’ population. Given successful experimental and natural introgression between colour morphs and the pervasive allopatry of related variants, we consider it unlikely that assortative preferences of both sexes were driven by direct selection during periods of secondary contact or, in turn, drove colour pattern differentiation in allopatry. Rather, we suggest that sexual isolation evolved as by‐product of allopatric divergence.     

Countergradient variation in carotenoid use between sympatric morphs of sockeye salmon (Oncorhynchus nerka) exposes nonanadromous hybrids in the wild by their mismatched spawning colour

Understanding the mechanisms that decrease gene flow between diverging populations is critical to understanding speciation. Anadromous (sockeye) and nonanadromous (kokanee) morphs of the Pacific sockeye salmon Oncorhynchus nerka spawn sympatrically and interbreed, yet allele frequency differences at neutral loci indicate restricted gene flow. Disruptive natural selection associated with strong selective differences between anadromous and nonanadromous life histories is thought to maintain the genetic differentiation of the morphs. Recently, a putative third morph of O. nerka exhibiting green rather than red breeding colour has been found on the spawning grounds sympatric with sockeye and kokanee. We investigated the ancestry of these green fish in a 2‐year controlled breeding study by using previously documented heritable, countergradient variation in red breeding colour to distinguish pure and hybrid morphs. Stabilizing sexual selection for similar red breeding colour in sockeye and kokanee has led to adaptive differences in the efficiency of carotenoid uptake between the morphs given differences in carotenoid availability between marine and lacustrine habitats. On the same diet, offspring parented by the green fish were intermediate in colour and in the concentration of dietary carotenoid pigments in their flesh and skin to those parented by either sockeye or kokanee; they were most similar to those parented by known kokanee × sockeye hybrids. This countergradient variation in carotenoid use results in a genotype‐environment mismatch in nonanadromous hybrids that exposes them by their breeding colour on the spawning grounds. Given that red colour is important in mate choice, sexual selection will almost certainly reduce reproductive opportunities for these hybrids and promote sympatric divergence of these incipient species. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society, 2005, 84 , 287–305.     

Evolution of adaptive diversity and genetic connectivity in Arctic charr (Salvelinus alpinus) in Iceland

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.     

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.     

Isolation by distance,not incipient ecological speciation,explains genetic differentiation in an Andean songbird (Aves: Furnariidae: Cranioleuca antisiensis,Line‐cheeked Spinetail) despite near threefold body size change across an environmental gradient

During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a P_ST–F_ST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.     

Interacting selection diversifies warning signals in a polytypic frog: an examination with the strawberry poison frog

Aposematic signals represent one of the most accessible traits to evaluate the interaction of natural and sexual selection on signal evolution. Here we investigate the contributions of these two selective forces on the aposematic signal evolution of the highly polytypic strawberry poison frog, Oophaga pumilio, of Bocas del Toro, Panama. Previous research has shown that the brightness of O. pumilio warning coloration can inform predators of the toxicity levels associated with different populations of the archipelago. Other studies suggest that sexual selection may be influencing warning signal brightness within populations via female mate choice (Isla Solarte, Isla Bastimentos, and Aquacate Peninsula populations) and male–male competition (Isla Solarte). Here we present two non-exclusive scenarios for how natural and sexual selection interact to drive phenotypic variation across this archipelago: (1) predators impose a selective regime whereby populations above a toxicity-brightness threshold are at liberty to diversify via sexual selection and below which populations are constrained to maintain a stricter resemblance to a more cryptic population mean, and (2) synergistic/additive effects of inter- and intrasexual selection drive the evolution of brighter males within populations above this toxicity threshold. We investigate whether aposematic patterns of divergence across the archipelago relative to the common mainland phenotype meet these predictions using existing data on O. pumilio morph toxicity measures and overall conspicuousness estimates to an avian predator. Using standardized z-scores to evaluate the range of trait values, we find that indeed the population representative of the common mainland phenotype (Almirante) represents an intermediate level of both toxicity and conspicuousness, and that derived Bocas del Toro populations vary in each of those components in directions predicted by the proposed scenarios. Furthermore, we find greater divergence towards conspicuousness than crypsis, a pattern suggestive of sexual and natural selection acting synergistically in morphs with high toxicity.     

Sexual selection and non-random mating for shell colour in a natural population of the marine snailLittorina mariae (Gastropoda: Prosobranchia)

In order to estimate the three independent components of mating behaviour, sexual selection in females, sexual selection in males and mating pattern, we studied the distribution of shell colour morphs among mating pairs and between copulating and non-copulating snails in four subsamples of a natural population ofL. mariae. The colour of the shell, the sex and a qualitative estimate of age was recorded for every snail. We found sexual selection acting against one of the two commonest colours (yellow) among the young females. However, in males none of the eight shell colour morphs was favoured during matings. Male sexual choice or differences in female sexual activity may cause the sexual fitness disadvantage of yellow females. Moreover, individuals of different colour morphs did not mate at random, rather dissasortatively. A behavioural choice among shell colour morphs or a non-random microdistribution of the morphs may cause the departure from random mating in this population.     

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.     

Genetic and phenotypic variation in a colourful treefrog across five geographic barriers

Aim Understanding the patterns and processes underlying phenotype in a polytypic species provides key insights into microevolutionary mechanisms of diversification. The red‐eyed treefrog, Agalychnis callidryas, exhibits strong regional differentiation in colour pattern, corresponding to five admixed mitochondrial DNA clades. We evaluated spatial diversity patterns across multiple, putative barriers to examine the fine‐scale processes that mediate phenotypic divergence between some regions while maintaining homogeneity between others. Location We examined patterns of phenotypic diversification among 17 sites that span five putative biogeographic barriers in lower Central America (Costa Rica and Panama). Methods We tested the extent to which genetic distance (F_ST) derived from six multilocus nuclear genotypes covaried with measures of phenotypic distance (leg coloration) within and between biogeographic regions. We used linear regression analyses to determine the role of geographic and genetic factors in structuring spatial patterns of phenotypic diversity. Results The factors that best explained patterns of phenotypic diversity varied among biogeographic regions. We identified one geographic barrier that impeded gene exchange and resulted in concordant phenotypic divergence across the Continental Divide, isolating Caribbean and Pacific populations. Across Caribbean Costa Rican populations, one barrier structured phenotypic but not genetic diversity patterns, indicating a role for selection. In other regions, the putative barriers had no determining effect on either genetic or leg colour structure. Main conclusions The processes mediating the distribution and diversification of colour pattern in this polytypic, wide‐ranging treefrog varied among biogeographic regions. Spatially varying selection combined with the isolating effects of geographic factors probably resulted in the patchy distribution of colour diversity across Costa Rican and Panamanian populations.     

Female sexual polymorphism and fecundity consequences of male mating harassment in the wild

Genetic and phenotypic variation in female response towards male mating attempts has been found in several laboratory studies, demonstrating sexually antagonistic co-evolution driven by mating costs on female fitness. Theoretical models suggest that the type and degree of genetic variation in female resistance could affect the evolutionary outcome of sexually antagonistic mating interactions, resulting in either rapid development of reproductive isolation and speciation or genetic clustering and female sexual polymorphisms. However, evidence for genetic variation of this kind in natural populations of non-model organisms is very limited. Likewise, we lack knowledge on female fecundity-consequences of matings and the degree of male mating harassment in natural settings. Here we present such data from natural populations of a colour polymorphic damselfly. Using a novel experimental technique of colour dusting males in the field, we show that heritable female colour morphs differ in their propensity to accept male mating attempts. These morphs also differ in their degree of resistance towards male mating attempts, the number of realized matings and in their fecundity-tolerance to matings and mating attempts. These results show that there may be genetic variation in both resistance and tolerance to male mating attempts (fitness consequences of matings) in natural populations, similar to the situation in plant-pathogen resistance systems. Male mating harassment could promote the maintenance of a sexual mating polymorphism in females, one of few empirical examples of sympatric genetic clusters maintained by sexual conflict.     

Reinforcement shapes clines in female mate discrimination in Drosophila subquinaria

Reinforcement of species boundaries may alter mate recognition in a way that also affects patterns of mate preference among conspecific populations. In the fly Drosophila subquinaria, females sympatric with the closely related species D. recens reject mating with heterospecific males as well as with conspecific males from allopatric populations. Here, we assess geographic variation in behavioral isolation within and among populations of D. subquinaria and use cline theory to understand patterns of selection on reinforced discrimination and its consequences for sexual isolation within species. We find that selection has fixed rejection of D. recens males in sympatry, while significant genetic variation in this behavior occurs within allopatric populations. In conspecific matings sexual isolation is also asymmetric and stronger in populations that are sympatric with D. recens. The clines in behavioral discrimination within and between species are similar in shape and are maintained by strong selection in the face of gene flow, and we show that some of their genetic basis may be either shared or linked. Thus, while reinforcement can drive extremely strong phenotypic divergence, the long‐term consequences for incipient speciation depend on gene flow, genetic linkage of discrimination traits, and the cost of these behaviors in allopatry.