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.     

INCIPIENT REPRODUCTIVE ISOLATION BETWEEN TWO SYMPATRIC MORPHS OF THE INTERTIDAL SNAIL LITTORINA SAXATILIS

The study of speciation in recent populations is essentially a study of the evolution of reproductive isolation mechanisms between sub-groups of a species. Prezygotic isolation can be of central importance to models of speciation, either being a consequence of reinforcement of assortative mating in hybrid zones, or a pleiotropic effect of morphological or behavioral adaptation to different environments. To suggest speciation by reinforcement between incipient species one must at least know that gene flow occurs, or have recently occurred, and that assortative mating has been established in the hybrid zone. In Galician populations of the marine snail Littorina saxatilis, two main morphs appear on the same shores, one on the upper-shore barnacle belt and the other in the lower-shore mussel belt. The two morphs overlap in distribution in the midshore where hybrids are found together with pure forms. Allozyme variation indicates that the two parental morphs share a common gene pool, although within shores, gene flow between morphs is less than gene flow within morphs. In this study, we observed mating behavior in the field, and we found that mating was not random in midshore sites, with a deficiency of heterotypic pairs. Habitat selection, assortative mating, and possibly sexual selection among females contributed to the partial reproductive isolation between the pure morphs. Sizes of mates were often positively correlated, in particular, in the upper shore, indicating size-assortative mating too. However, this seemed to be a consequence of nonrandom microdistributions of snails of different sizes. Because we also argue that the hybrid zone is of primary rather than secondary origin, this seems to be an example of sympatric reproductive isolation, either established by means of reinforcement or as a by-product to divergent selection acting on other characters.     

Disruptive selection on female reproductive characters in a hybrid zone of Littorina saxatilis

We studied natural selection in a hybrid zone of the intertidal marine snail Littorina saxatilis located in Galicia (NW Spain) by measuring the number and average size of the embryos carried by females. We related these characters with the females' position on the phenotypic and environmental gradients between the two pure morphs of the hybrid zone. In contrast with previous interpretations of studies made in this hybrid zone, we found a depression in both embryo number and size for phenotypically intermediate females. This disruptive natural selection could play an important role in the maintenance of the population's bimodal phenotypic distribution and in that of the assortative mating between the pure morphs. We found also that intermediate environments tended to be unfavorable for all phenotypes. Although the precise causes for the found depression in female reproductive characters remain to be determined, these results serve to emphasize the importance of studying whole fitness surfaces in hybrid zones, across the complete phenotypic and environmental ranges, instead of merely comparing the fitness averages of the two pure morphs and the intermediate individuals, taken as three discrete phenotypic classes.     

Recent speciation and secondary contact in endemic ants

Gene flow is the main force opposing divergent selection, and its effects are greater in populations in close proximity. Thus, complete reproductive isolation between parapatric populations is not expected, particularly in the absence of ecological adaptation and sharp environmental differences. Here, we explore the biogeographical patterns of an endemic ant species, Cataglyphis floricola, for which two colour morphs (black and bicolour) coexist in parapatry throughout continuous sandy habitat in southern Spain. Discriminant analyses of six biometric measurements of male genitalia and 27 cuticular hydrocarbons reveal high differentiation between morphs. Furthermore, the low number of shared alleles derived from nuclear markers (microsatellites) between the morphs at their contact zone suggests the absence of recent gene flow. Mitochondrial DNA (COI) phylogenetic analysis and median‐joining networks show that the black morph is basal to the bicolour morph, with unique haplotypes recovered for each morph. Mismatch distribution analysis and Bayesian skyline plots suggest that they are undergoing different demographic changes, with the bicolour and black morphs at demographic equilibrium and expansion, respectively. Thus, our results show complete reproductive isolation between the two colour morphs as evidenced from genetic, chemical and morphological data. We suggest that these divergence events could be explained by historical vicariance during the Pleistocene, in which reproductive traits experienced strong divergent selection between the morphs initiating or culminating speciation.     

Lack of mitochondrial genetic structure in hamlets (Hypoplectrus spp.): recent speciation or ongoing hybridization?

Species in the genus Hypoplectrus (hamlet fish) have been recognized primarily on the basis of colour morphology, which varies substantially. Limited differentiation in other morphological characters, however, has led to ongoing debate about their taxonomy. Our analysis of mtDNA sequences demonstrates neither reciprocal monophyly nor significant genetic differentiation among hamlet colour morphs. These data are potentially consistent with a model of recent speciation due to sexual selection on colour morphology. The presence within hamlets of two divergent mtDNA clades, however, suggests a longer history during which hybridization and gene flow have prevented the differentiation of hamlet colour morphs, at least in neutral genetic markers.     

Selection on hybrids of ecologically divergent ecotypes of a marine snail: the relative importance of exogenous and endogenous barriers

Unravelling the form of selection acting on hybrids of ecotypes undergoing ecological speciation is essential to understand the mechanisms behind the evolution of reproductive isolation in the face of gene flow. Shell phenotype is known to be affected by natural selection and is involved in the fitness of the marine snail Littorina saxatilis. Here, we studied the association between shell traits and fitness in hybrids in order to determine the relative role of exogenous and endogenous selection in this hybrid zone of L. saxatilis. We show that directional selection is the predominant mode of selection among hybrids. We also show its heterogeneity, affecting different shell traits, within populations at the level of the microhabitat. Therefore, endogenous selection mechanisms are most probably lacking in this hybrid zone and exogenous barriers (pre‐ and post‐zygotic) are possibly one of the main forces behind the evolution of barriers to gene flow between these ecologically divergent ecotypes. This study shows how this barrier might represent an important type of reproductive isolation within ecological speciation, and this should be taken into account in future studies of speciation in hybrid zones. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 391–400.     

Morphological variation in a secondary contact between divergent lineages of brown trout (Salmo trutta) from the Iberian Peninsula

The aim of this study was to analyze the morphological variation of brown trout (Salmo trutta) in the Duero basin, an Atlantic river basin in the Iberian Peninsula, where a spatial segregation of two divergent lineages was previously reported, based on isozyme, microsatellite and mtDNA data. In these studies, two divergent pure regions (Pisuerga and Lower-course) and several hybrid populations between them were identified. Morphological variation was evaluated in 11 populations representative of the genetic differentiation previously observed in the Duero basin, using multivariate analysis on 12 morphometric and 4 meristic traits. A large differentiation between populations was observed (interpopulation component of variance: 41.8%), similar to that previously detected with allozymes and microsatellites. Morphometric differentiation was also reflected by the high classification success of pure and hybrid individuals to their respective populations, using multivariate discriminant functions (94.1% and 79.0%, respectively). All multivariate and clustering analyses performed demonstrated a strong differentiation between the pure regions. The hybrid populations, though showing large differentiation among them, evidenced an intermediate position between the pure samples. Head and body shape traits were the most discriminant among the morphometric characters, while pectoral rays and gillrakers were the most discriminant among the meristic traits. These results confirmed the high divergence of the brown trout from the Duero basin and suggest some traits on which selection could be acting to explain the spatial segregation observed.     

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.     

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata)

Plant–pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator‐mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (F_ST) for nuclear and chloroplast microsatellite loci and used a half‐sib design to calculate quantitative trait divergence (Q_ST) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean F_ST to be 0.032. Q_ST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean F_ST. We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non‐neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.     

Evidence for divergent natural selection of a Lake Tanganyika cichlid inferred from repeated radiations in body size

Divergent natural selection is thought to play a vital role in speciation, but clear, measurable examples from nature are still few. Among the many possible sources of divergent natural selection, predation pressure may be important because predators are ubiquitous in food webs. Here, we show evidence for divergent natural selection in a Lake Tanganyika cichlid, Telmatochromis temporalis , which uses burrows under stones or empty snail shells as shelters. This species contains normal and dwarf morphs at several localities. The normal morph inhabits rocky shorelines, whereas the dwarf morph invariably inhabits shell beds, where empty snail shells densely cover the lake bottom. Genetic evidence suggested that the dwarf morph evolved independently from the normal morph at two areas, and morphological analysis and evaluation of habitat structure revealed that the body sizes of morphs closely matched the available shelter sizes in their habitats. These findings suggest that the two morphs repeatedly evolved through divergent natural selection associated with the strategy for sheltering from predators.     

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.     

Frequency- and density-dependent sexual selection in natural populations of Galician Littorina saxatilis Olivi

Galician exposed shore populations of the direct developing periwinkle Littorina saxatilis are strikingly polymorphic, with an ornamented and banded upper shore form and a smooth and unbanded lower shore form. Intermediates between the two pure forms occur in a narrow mid shore zone together with the parental forms. We have previously shown that the two pure forms share the same gene pool but that mating between them is non-random. This is due to a non-random microdistribution in the zone of overlap, and also to assortative mating. In this study we present data which show that intermediate (hybrid) females mate less often than pure females in micropatches dominated by either of the pure forms, but not in micropatches in which the two pure forms are equally common. Thus, sexual fitness in intermediate females depends on the frequency of both pure morphs. Furthermore, sexual selection against intermediate females also varies with the densities of snails within each micro patch. The biological mechanisms which may explain this particular reduction of female hybrid fitness are discussed. Assortative mating between the pure morphs is sometimes almost complete, while both morphs do not mate the intermediates assortatively. In the light of this, sexual selection against intermediate females may contribute considerably to restrict gene flow between the pure forms.     

DEVELOPMENTAL DECOUPLING OF ALTERNATIVE PHENOTYPES: INSIGHTS FROM THE TRANSCRIPTOMES OF HORN‐POLYPHENIC BEETLES

Developmental mechanisms play an important role in determining the costs, limits, and evolutionary consequences of phenotypic plasticity. One issue central to these claims is the hypothesis of developmental decoupling, where alternate morphs result from evolutionarily independent developmental pathways. We address this assumption through a microarray study that tests whether differences in gene expression between alternate morphs are as divergent as those between sexes, a classic example of developmental decoupling. We then examine whether genes with morph‐biased expression are less conserved than genes with shared expression between morphs, as predicted if developmental decoupling relaxes pleiotropic constraints on divergence. We focus on the developing horns and brains of two species of horned beetles with impressive sexual‐ and morph‐dimorphism in the expression of horns and fighting behavior. We find that patterns of gene expression were as divergent between morphs as they were between sexes. However, overall patterns of gene expression were also highly correlated across morphs and sexes. Morph‐biased genes were more evolutionarily divergent, suggesting a role of relaxed pleiotropic constraints or relaxed selection. Together these results suggest that alternate morphs are to some extent developmentally decoupled, and that this decoupling has significant evolutionary consequences. However, alternative morphs may not be as developmentally decoupled as sometimes assumed and such hypotheses of development should be revisited and refined.     

Morphological evidence of correlational selection and ecological segregation between dextral and sinistral forms in a polymorphic flatfish, Platichthys stellatus

Phenotypic polymorphisms in natural systems are often maintained by ecological selection, but only if niche segregation between morphs exists. Polymorphism for eyed-side direction is rare among the approximately 700 species of flatfish (Pleuronectiformes), and the evolutionary mechanisms that maintain it are unknown. Platichthys stellatus (starry flounder) is a polymorphic pleuronectid flatfish exhibiting large, clinal variation in proportion of left-eyed (sinistral) morphs, from 50% in California to 100% in Japan. Here I examined multiple traits related to swimming and foraging performance between sinistral and dextral morphs of P. stellatus from 12 sites to investigate if the two morphs differ in ways that may affect function and ecology. Direction of body asymmetry was correlated with several other characters: on an average, dextral morphs had longer, wider caudal peduncles, shorter snouts and fewer gill rakers than sinistral morphs. Although the differences were small in magnitude, they were consistent in direction across samples, implying that dextral and sinistral starry flounder may be targeting different prey types. Morphological differences between morphs were greatest in samples where the chances of competitive interactions between them were the greatest. These results suggest that the two morphs are not ecologically identical, may represent a rare example of divergent selection maintaining polymorphism of asymmetric forms, and that correlational selection between body asymmetry and other characters may be driven by competitive interactions between sinistral and dextral flatfish. This study is one of very few that demonstrates the ecological significance of direction in a species with polymorphic asymmetric forms.     

INCIPIENT SPECIATION DESPITE LITTLE ASSORTATIVE MATING: THE YELLOW‐RUMPED WARBLER HYBRID ZONE

Hybrid zones between recently diverged taxa are natural laboratories for speciation research, allowing us to determine whether there is reproductive isolation between divergent forms and the causes of that isolation. We present a study of a classic avian hybrid zone in North America between two subspecies of the yellow‐rumped warbler (Dendroica coronata). Although previous work has shown very little differentiation in mitochondrial DNA across this hybrid zone, we identified two nuclear loci (one sex‐linked and one autosomal) that show fixed differences across the hybrid zone, in a close concordance with patterns of plumage variation. Temporal stability and limited width of the hybrid zone, along with substantial linkage disequilibrium between these two diagnostic markers in the center of the zone, indicate that there is moderate reproductive isolation between these populations, with an estimated strength of selection maintaining the zone of 18%. Pairing data indicate that assortative mating is either very weak or absent, suggesting that this reproductive isolation is largely due to postmating barriers. Thus, despite extensive hybridization the two forms are distinct evolutionary groups carrying genes for divergent adaptive peaks, and this situation appears relatively stable.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by‐product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced‐representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus‐specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by‐product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.     

Using environmental effects on fecundity to compare the adaptive characteristics of the morphs in a hybrid zone of Littorina saxatilis

We studied the role of ecological processes in the maintenance of a hybrid zone for the intertidal snail Littorina saxatilis by comparing the fecundity responses of the females from different morphs to local environmental variation. This procedure makes it possible to discern if the environments maximizing fecundity are the same in these morphs, and therefore to study adaptive differences between them. We found that the pure morphs were adapted to their respective shore levels, thus confirming that this hybrid zone is associated with the occurrence of an ecotone between two habitats. However, fecundity differences did not explain the precise spatial location and width of the hybrid zone. We also found that the adaptive characteristics of hybrids were very different from those of the pure morphs, in sharp contrast with any expectations based on their intermediate phenotype. The comparison of fitness responses to environmental variation might be an useful procedure in the field study of adaptive differences between groups of individuals.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation

Differential gene flow, reductions in diversity following linked selection and/or features of the genome can structure patterns of genomic differentiation during the process of speciation. Possible sources of reproductive isolation are well studied between coastal and inland subspecies groups of Swainson's thrushes, with differences in seasonal migratory behaviour likely playing a key role in reducing hybrid fitness. We assembled and annotated a draft reference genome for this species and generated whole‐genome shotgun sequence data for populations adjacent to the hybrid zone between these groups. We documented substantial genomewide heterogeneity in relative estimates of genetic differentiation between the groups. Within population diversity was lower in areas of high relative differentiation, supporting a role for selective sweeps in generating this pattern. Absolute genetic differentiation was reduced in these areas, further suggesting that recurrent selective sweeps in the ancestral population and/or between divergent populations following secondary contact likely occurred. Relative genetic differentiation was also higher near centromeres and on the Z chromosome, suggesting that features of the genome also contribute to genomewide heterogeneity. Genes linked to migratory traits were concentrated in islands of differentiation, supporting previous suggestions that seasonal migration is under divergent selection between Swainson's thrushes. Differences in migratory behaviour likely play a central role in the speciation of many taxa; we developed the infrastructure here to permit future investigations into the role several candidate genes play in reducing gene flow between not only Swainson's thrushes but other species as well.