Numts help to reconstruct the demographic history of the ocellated lizard (Lacerta lepida) in a secondary contact zone

In northwestern Iberia, two largely allopatric Lacerta lepida mitochondrial lineages occur, L5 occurring to the south of Douro River and L3 to the north, with a zone of putative secondary contact in the region of the Douro River valley. Cytochrome b sequence chromatograms with polymorphisms at nucleotide sites diagnostic for the two lineages were detected in individuals in the region of the Douro River and further north within the range of L3. We show that these polymorphisms are caused by the presence of four different numts (I-IV) co-occurring with the L3 genome, together with low levels of heteroplasmy. Two of the numts (I and II) are similar to the mitochondrial genome of L5 but are quite divergent from the mitochondrial genome of L3 where they occur. We show that these numts are derived from the mitochondrial genome of L5 and were incorporated in L3 through hybridization at the time of secondary contact between the lineages. The additional incidence of these numts to the north of the putative contact zone is consistent with an earlier postglacial northward range expansion of L5, preceding that of L3. We show that genetic exchange between the lineages responsible for the origin of these numts in L3 after secondary contact occurred prior to, or coincident with, the northward expansion of L3. This study shows that, in the context of phylogeographic analysis, numts can provide evidence for past demographic events and can be useful tools for the reconstruction of complex evolutionary histories.     

Strong selection against hybrids maintains a narrow contact zone between morphologically cryptic lineages in a rainforest lizard

Phenotypically cryptic lineages comprise an important yet understudied part of biodiversity; in particular, we have much to learn about how these lineages are formed and maintained. To better understand the evolutionary significance of such lineages, we studied a hybrid zone between two morphologically cryptic phylogeographic lineages in the rainforest lizard, Lampropholis coggeri. Analyzing a multilocus genetic dataset through cline inference, individual-based methods and population measures of disequilibrium and using simulations to explore our genetic results in context of theoretical expectations, we inferred the processes maintaining this hybrid zone. We find that these lineages meet in a hybrid zone that is narrow (≈400 m) relative to inferred dispersal rate. Further, the hybrid zone exhibits substantial genetic disequilibrium and sharply coincident and largely concordant clines. Based on our knowledge about the region's biogeography, the species' natural history, and our simulation results, we suggest that strong selection against hybrids structures this system. As all clines show a relatively narrow range of introgression, we posit that this hybrid zone might not yet be in equilibrium. Nonetheless, our results clearly show that phylogeographic lineages can evolve substantial reproductive isolation without concomitant morphological diversification, suggesting that such lineages can constitute a significant component of evolutionary diversity.     

Reproductive isolation between oviparous and viviparous lineages of the Eurasian common lizard Zootoca vivipara in a contact zone

Contact zones between two evolutionary lineages are often useful for understanding the process of speciation because the observed genetic pattern reflects the history of differentiation. The Eurasian lacertid lizard Zootoca vivipara is a potentially interesting model for studying the role of reproductive mode in the speciation of squamate reptiles because it has both oviparous (Zootoca vivipara carniolica) and viviparous (Zootoca vivipara vivipara) populations that have recently been shown to be genetically distinct. We studied a newly‐discovered syntopic area of these two Zootoca subspecies in the central Italian Alps using genetic markers to investigate the level of introgression between them. Patterns of genetic differentiation in a fragment of the mitochondrial DNA cytb gene and a set of nuclear microsatellites show that the speciation process is complete in this area, with no evidence of recent introgression. Phylogenetic and genotypic divergence suggests that the two subspecies have experienced long independent evolutionary histories, during which genetic and phenotypic differences evolved. The possible roles of biogeography, reproductive mode, and cytogenetic differentiation in this speciation process are discussed. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 566–573.     

Genetic structure, introgression, and a narrow hybrid zone between northern and California spotted owls (Strix occidentalis)

The northern spotted owl (Strix occidentalis caurina) is a threatened subspecies and the California spotted owl (Strix occidentalis occidentalis) is a subspecies of special concern in the western United States. Concern for their continued viability has arisen because of habitat loss caused by timber harvesting. The taxonomic status of the northern subspecies has been the subject of continuing controversy. We investigated the phylogeographical and population genetic structure of northern and California spotted owls with special reference to their region of contact. Mitochondrial DNA (mtDNA) control region sequences confirmed the existence of two well-differentiated lineages connected by a narrow hybrid zone in a region of low population density in north central California. Maximum-likelihood estimates indicated bidirectional gene flow between the lineages but limited introgression outside the region of contact. The lengths of both the mtDNA hybrid zone and the reduced density patch were similar and slightly exceeded estimates of natal dispersal distances. This suggests that the two subspecies were in secondary contact in a hybrid zone trapped by a population density trough. Consequently, the zone of interaction is expected to be geographically stable. We discovered a third, rare clade of haplotypes, which we interpreted to be a result of incomplete lineage sorting; those haplotypes result in a paraphyletic northern spotted owl with respect to the California spotted owl. A congeneric species, the barred owl (Strix varia), occasionally hybridizes with spotted owls; our results indicated an upper bound for the frequency of barred owl mtDNA haplotypes in northern spotted owl populations of 3%.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

When vicars meet: a narrow contact zone between morphologically cryptic phylogeographic lineages of the rainforest skink, Carlia rubrigularis

Phylogeographic analyses of the fauna of the Australian wet tropics rainforest have provided strong evidence for long-term isolation of populations among allopatric refugia, yet typically there is no corresponding divergence in morphology. This system provides an opportunity to examine the consequences of geographic isolation, independent of morphological divergence, and thus to assess the broader significance of historical subdivisions revealed through mitochondrial DNA phylogeography. We have located and characterized a zone of secondary contact between two long isolated (mtDNA divergence > 15%) lineages of the skink Carlia rubrigularis using one mitochondrial and eight nuclear (two intron, six microsatellite) markers. This revealed a remarkably narrow (width < 3 km) hybrid zone with substantial linkage disequilibrium and strong deficits of heterozygotes at two of three nuclear loci with diagnostic alleles. Cline centers were coincident across loci. Using a novel form of likelihood analysis, we were unable to distinguish between sigmoidal and stepped cline shapes except at one nuclear locus for which the latter was inferred. Given estimated dispersal rates of 90-133 m x gen(-1/2) and assuming equilibrium, the observed cline widths suggest effective selection against heterozygotes of at least 22-49% and possibly as high as 70%. These observations reveal substantial postmating isolation, although the absence of consistent deviations from Hardy-Weinberg equilibrium at diagnostic loci suggests that there is little accompanying premating isolation. The tight geographic correspondence between transitions in mtDNA and those for nuclear genes and corresponding evidence for selection against hybrids indicates that these morphologically cryptic phylogroups could be considered as incipient species. Nonetheless, we caution against the use of mtDNA phylogeography as a sole criterion for defining species boundaries.     

Hybridization and asymmetric introgression across a narrow zone of contact between Neotoma fuscipes and N. macrotis (Rodentia: Cricetidae)

When lineages diverge in allopatry and come into secondary contact, we have a unique opportunity to examine the degree to which they have become reproductively isolated from one another and the mechanisms that contribute to rates of interspecific gene flow. If hybridization and introgression have occurred in the past or are ongoing, examining patterns of variation in differentially inherited genetic markers can provide insight into underlying mechanisms determining interspecific reproductive interactions. We investigated genetic variation in a recently discovered contact zone between two species of woodrats (Neotoma fuscipes and N. macrotis) in central coastal California. Previous studies have found evidence of historic hybridization between these species, but an active site of sympatry had yet to be discovered. Here, we describe the first known area of sympatry between the two species and present evidence of ongoing hybridization at this site. We intensively sampled throughout the narrow area of interspecific contact and continuously into the adjacent areas where each taxon existed largely in isolation of the other. We genotyped 851 individuals sampled in a single year at 15 nuclear microsatellite loci and sequenced a portion of the cytochrome b gene in a subset of these to examine patterns of introgression. We find that the area of sympatry is less than 1 km in length, but hybridization appeared to be fairly common, with 15% of the woodrats being of hybrid origin. At least some hybrids are fertile because introgression was evident. However, introgression appeared to be asymmetric with backcrosses toward N. macrotis being more common, a pattern we suspect is due in part to directional mate choice related to differential body size. This study adds to our growing understanding of the nature of species boundaries, especially between lineages that are far along the continuum towards speciation. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 115 , 162–172.     

A hybrid zone between two subspecies of Chorthippus parallelus. X-chromosome variation through a contact zone

Two grasshopper subspecies, Chorthippus parallelus parallelus (Cpp) and Chorthippus parallelus erythropus (Cpe), meet along the Pyrenees where they hybridize and produce a hybrid zone. A contact zone located in the Col de Portalet has been analyzed for the distribution of chromosome markers on the sex (X) chromosome in pure and hybrid populations. C-banding allowed us to distinguish both pure subspecific forms and recombinant forms, and to demonstrate their different frequencies through the contact zone. Interestingly, the distal C-band (P) on the X chromosome that characterizes Cpp individuals occurs at very high frequency through the zone and then drops sharply, close to where pure Cpe populations are found. A novel interstitial C-band (H) has been found, probably derived from that characterizing Cpe individuals (E). This marker band is only present in hybrid populations. These data are discussed in the light of the dynamics of the hybrid zone.     

Genetic panmixia within a narrow contact zone between chromosomally and ecologically distinct black fly sibling species (Diptera: Simuliidae)

Hybrid zones are windows into the speciation process, and their study can give clues into the maintenance and breakdown of species boundaries. Using both genetic and ecological tools, we investigate lineage diversification across a contact zone characterized by chromosome rearrangements. We show that black fly sibling species, Simulium arcticum sensu stricto (s.s.) and Simulium saxosum, lack genetic differentiation at both microsatellite and mtDNA loci in allopatry and sympatry, as well as exhibit high levels of gene flow and continuous chromosome variation in sympatry. Furthermore, hybrid frequencies at the contact zone are similar to those seen between races, rather than species. In contrast, S. arcticum s.s. and S. saxosum maintain ecological differences and distinct habitat associations ‐ the contact zone situated at the margin of suitable habitat for each sibling species. Moreover, gene flow occurs only in a narrow band along an ecological transition. Except for the contact zone, S. arcticum s.s. and S. saxosum hybrids do not occur elsewhere within the sibling species' ranges. Although S. arcticum s.s. and S. saxosum maintain the potential to interbreed freely, we conclude that habitat associations and, perhaps, chromosome systems prevent expansion of ranges and assimilation of lineages.     

Individual performance in relation to cytonuclear discordance in a northern contact zone between long‐toed salamander (Ambystoma macrodactylum) lineages

Cytonuclear discordance in contact zones between related lineages is common, with mitochondrial clines often being displaced from clines in nuclear allele frequency. Proposed explanations for such a pattern include adaptive introgression of mtDNA or a neutral wake of mtDNA being left behind following hybrid zone movement. However, studies investigating these hypotheses are rare. Our previous survey of genetic variation in the long‐toed salamander (Ambystoma macrodactylum) highlighted a potential case of cytonuclear discordance between two lineages in western Canada. Here, we use additional markers and samples to clarify the extent of this discordance. We simultaneously assess the feeding performance of individuals in a common environment to test for an association between mitotype and individual performance. The genetic results confirm a general pattern of cytonuclear discordance in the focal region. However, we also observed more limited introgression of a diagnostic nuclear marker. Intriguingly, although there were differences in individual performance associated with the transition between mitotypes, these differences were not fully explained by mitotype. Instead, the lowest performance was observed in individuals demonstrating the greatest mismatch between mtDNA and all nuclear markers, suggesting the potential for cytonuclear incompatibilities to be acting. These results highlight the complexity of understanding the causes and consequences of mtDNA introgression and cytonuclear discordance in contact zones.     

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.     

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X‐chromosome, Y‐chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95–4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co‐localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X‐ and Y‐chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.     

Asymmetric and differential gene introgression at a contact zone between two highly divergent lineages of field voles (Microtus agrestis)

Secondary contact zones have the potential to shed light on the mode and rate at which reproductive isolation accumulates during allopatric speciation. We investigated the population genetics of a contact zone between two highly divergent lineages of field voles (Microtus agrestis) in the Swiss Jura mountains. To shed light on the processes underlying introgression, we used maternally, paternally, and bi-parentally inherited markers. Though the two lineages maintained a strong genetic structure, we found some hybrids and evidence of gene flow. The extent of introgression varied with the mode of inheritance, being highest for mtDNA and absent for the Y chromosome. In addition, introgression was asymmetric, occurring only from the Northern to the Southern lineage. Both patterns seem parsimoniously explained by neutral processes linked to differences in effective sizes and sex-biased dispersal rates. The lineage with lower effective population size was also the more introgressed, and the mode-of-inheritance effect correlated with the male-biased dispersal rate of microtine rodents. We cannot exclude, however, that Haldane's effect contributed to the latter, as we found a marginally significant deficit in males (the heterogametic sex) among hybrids. We propose a possible demographic scenario to account for the patterns documented, and empirical extensions to further investigate this contact zone.     

Genetic analysis of a chromosomal hybrid zone in the Australian morabine grasshoppers (Vandiemenella, viatica species group)

Whether chromosomal rearrangements promote speciation by providing barriers to gene exchange between populations is one of the long-standing debates in evolutionary biology. This question can be addressed by studying patterns of gene flow and selection in hybrid zones between chromosomally diverse taxa. Here we present results of the first study of the genetic structure of a hybrid zone between chromosomal races of morabine grasshoppers Vandiemenella viatica , P24(XY) and viatica 17, on Kangaroo Island, Australia. Chromosomal and 11 nuclear markers revealed a narrow hybrid zone with strong linkage disequilibrium and heterozygote deficits, most likely maintained by a balance between dispersal and selection. Widths and positions of clines for these markers are concordant and coincident, suggesting that selection is unlikely to be concentrated on a few chromosomes. In contrast, a mitochondrial marker showed a significantly wider cline with centre offset toward the P24(XY) side. We argue that the discordance between the mitochondrial and nuclear/chromosomal clines and overall asymmetry of the clines suggest a secondary origin of the contact zone and potential movement of the zone after contact. Genome-wide scans using many genetic markers and chromosomal mapping of these markers are needed to investigate whether chromosomal differences directly reduce gene flow after secondary contact.     

Phylogeography of the Trans‐Volcanic bunchgrass lizard (Sceloporus bicanthalis) across the highlands of south‐eastern Mexico

We quantify the population divergence processes that shaped population genetic structure in the Trans‐Volcanic bunchgrass lizard (Sceloporus bicanthalis) across the highlands of south‐eastern Mexico. Multilocus genetic data from nine nuclear loci and mitochondrial (mt)DNA were used to estimate the population divergence history for 47 samples of S. bicanthalis. Bayesian clustering methods partitioned S. bicanthalis into three populations: (1) a southern population in Oaxaca and southern Puebla; (2) a population in western Puebla; and (3) a northern population with a broad distribution across Hidalgo, Puebla, and Veracruz. The multilocus nuclear data and mtDNA both supported a Late Pleistocene increase in effective population size, and the nuclear data revealed low levels of unidirectional gene flow from the widespread northern population into the southern and western populations. Populations of S. bicanthalis experienced different demographic histories during the Pleistocene, and phylogeographical patterns were similar to those observed in many co‐distributed highland taxa. Although we recommend continuing to recognize S. bicanthalis as a single species, future research on the evolution of viviparity could gain novel insights by contrasting physiological and genomic patterns among the different populations located across the highlands of south‐eastern Mexico. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 852–865.     

A three‐way contact zone between forms of Patella rustica (Mollusca: Patellidae) in the central Mediterranean Sea

Previous studies have reported the occurrence of three differentiated mtDNA lineages within Patella rustica in the Mediterranean Sea. Two hypotheses have been proposed to explain these observations: (1) the maintenance of ancestral polymorphism within a single species; (2) the occurrence of cryptic species not identified previously. To distinguish between these hypotheses, we screened the genetic variability at nine allozyme loci, an intron from the α‐amylase gene and a mitochondrial gene for 187 individuals of P. rustica sampled from seven Mediterranean localities. Eight additional localities were screened for the last two markers to place the differentiated lineages in a clear geographic context. Our results demonstrate that the three mtDNA lineages correspond to three distinct nuclear genotype clusters and provide further details on their distribution: the cluster corresponding to the mtDNA lineage from the Atlantic and western Mediterranean extends as far as the south coast of Italy, whereas the remaining two clusters occur in sympatry in the eastern Mediterranean. One of the eastern Mediterranean clusters is highly differentiated and seems to be reproductively isolated from the codistributed form; we therefore suggest that it corresponds to a new species. The remaining two clusters are less differentiated and form a contact zone across south Italian shores. This three‐way contact zone constitutes an interesting model for the study of speciation in the marine realm. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 154–169.     

Genetic structure of natural populations of Castanea sativa in Turkey: evidence of a hybrid zone

This study points out the evidence of a hybrid zone between two groups of genetically differentiated populations of chestnut (Castanea sativa Mill.) in Turkey. Genetic structure, gene flow and introgression levels, based on 16 allozyme loci, were investigated on 34 population samples spanning the entire C. sativa distribution area in this country. The occurrence of the hybrid zone, located in the Bithynian region, was inferred in a group of populations showing the following genetic characteristics: (i) enhanced genetic variability and intermediate allelic frequencies between those of the Western and Eastern groups of populations; (ii) sharp and concordant changes in allele frequencies; (iii) decreased gene flow with the Western and Eastern populations. Starting from the cline width estimated to be 324 km, strength of selection was evaluated from the gene flow distance, as indicated from the degree of genetic structuring outside the hybrid zone. Evolutionary processes shaping the observed genetic differentiation and introgression are discussed on the basis of palynological data, palaeoclimatic events and evidence of hybridization found in other plant and animal species in the same region.     

Evidence for low‐level hybridization between two allochronic populations of the pine processionary moth,Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

Divergence between populations sharing the same habitat can be initiated by different reproductive times, leading to allochronic differentiation. A spatially localized allochronic summer population (SP) of the pine processionary moth Thaumetopoea pityocampa, recently discovered in Portugal, occurs in sympatry with the local winter population (WP). We examined the level of genetic differentiation between the two populations and estimated the current gene flow within the spatial framework of their co‐occurrence. Mitochondrial data indicated that the two sympatric populations were genetically closer than other WP populations. Conversely, microsatellite genotyping uncovered greater differentiation between the two sympatric populations than between allopatric ones. While male trapping confirmed that reproduction of SP and WP occurred at distinct times, clustering approaches demonstrated the presence of a few LateSP individuals emerging within the WP flight period, although genetically identified as SP. We also identified rare recent hybridization events apparently occurring mainly in the margins of the current SP range. The ongoing gene flow detected between the ancestral and the emerging allochronic populations revealed an incomplete reproductive isolation, which must therefore be taken into account and integrated with studies focussed on ecological drivers, so that a complete understanding of the ongoing speciation process might be achieved.     

Genetic structure of a hybrid zone between two violets,Viola rossii Hemsl. and V. bissetii Maxim.: dominance of F1 individuals in a narrow contact range

The genetic composition of a hybrid zone can provide insight into the evolution of diversification in plants. We carried out morphological and amplified fragment length polymorphism analyses to investigate the genetic composition of a hybrid zone between two violets, Viola bissetii Hemsl. and Viola rossii Maxim. Our aim was to clarify the formation and maintenance of hybrids between these Viola species. We found that most hybrid individuals (V. bissetii × V. rossii) were of the F₁ generation, with a few of the F₂ generation. We found no backcrosses. The scarcity of post‐F₁ hybrids indicates that a species barrier is established between the parental species. The F₁‐dominated hybrid zone occupied only a narrow, intermediate ecotone between the parental habitats, suggesting that selection by environmental factors against hybrids may help to maintain the current conditions in this hybrid zone.     

Sexual selection on phenotypic traits in a hybrid zone of Littorina saxatilis (Olivi)

Step clinal transitions in inherited character(s) between genetically distinct populations are usually referred to as hybrid zones. An example is found in the population of the intertidal snail Littorina saxatilis in Galicia (NW Spain). We studied the shape of the overall fitness surface for sexual selection in this hybrid zone, and the position of hybrids and pure morphs on this surface. We found that sexual divergent selection acted on a combination of phenotypic traits separating the pure morphs, and therefore that sexual selection contributed to morph differentiation. The average fitness of hybrids as a group was not significantly different from that of the pure morphs, but they did show divergent sexual selection in some traits. These results are in agreement with a model of divergent selection favouring both the pure morph as well as those hybrids most resembling each morph. The finding of divergent selection is remarkable because quadratic selection gradients are usually weak in nature.