Mate Choice Contributes to the Maintenance of Shell Color Polymorphism in a Marine Snail via Frequency-Dependent Sexual Selection

Introduction

Negative frequency-dependent selection (NFDS) operates when the relative fitness of a morph within a natural population is an inverse function of its relative frequency, i.e., the more abundant a morph is, the lower its fitness (Wright, 1939). Although NFDS has been considered the most intuitive and powerful force maintaining polymorphisms (Ayala and Campbell, 1974; Fitzpatrick et al., 2007; Trotter and Spencer, 2008; Svensson and Connallon, 2019), relatively few examples exist. Thus, the biological causes for this pattern in natural populations remain largely unknown.

Patterns of NFDS can arise in several distinct, non-exclusive scenarios: driven by predation (apostatic selection), where less frequent morphs have a survival advantage (Tucker, 1991; Olendorf et al., 2006; Fitzpatrick et al., 2009; Weir, 2018); by sexual conflict, where less frequent morphs of one sex obtain fitness benefits in detriment or against fitness interests in the other sex (Svensson et al., 2005; Gosden and Svensson, 2007); by mate competition, where the less frequent morph within one sex obtains a greater number of matings than the common one when competing for them (Sinervo and Lively, 1996; Bleay et al., 2007; Rios-Cardenas et al., 2018); and by mate choice, in cases when rare or dissimilar morphs are preferred and therefore have relatively more matings (Hughes et al., 2013; Chouteau et al., 2017; Valvo et al., 2019). It becomes obvious that the main difference between these scenarios is the component of fitness that is exploited by the less frequent morph. Although the simultaneous study of all these components in one species would be a cumbersome task, both mate choice and mate competition scenarios can be easily studied together as they relate to the same component of fitness. In this study, we refer to these two scenarios as sexual fitness (following Hedrick, 1983, p.163; Bell, 2008, p.167), a component of fitness defined by the capacity of individuals to find mates.

Mate choice is the process by which certain traits in one sex lead to non-random mating (Kokko et al., 2003; Edward, 2015). This can lead to NFDS when there is a tendency by the choosy sex to mate with the less frequent or rare morphs, which consequently enjoy greater sexual fitness (Hughes et al., 2013; Valvo et al., 2019); or when there is a fixed preference to mate for dissimilar morphs (Chouteau et al., 2017). A usual proxy for inferring such a choice is disassortative mating. Whenever observed, a straightforward consequence of disassortative mating is that the polymorphism is protected from fixation or loss, independently of other findings (Takahashi and Hori, 2008; Hedrick et al., 2016, 2018; Chouteau et al., 2017). Disassortative mating is, however, a rare event in nature (Jiang et al., 2013). It has therefore not been extensively studied, though a recent review suggests it may be more frequent in marine invertebrates than previously assumed (Olsen et al., 2020). A common explanation for the evolution of disassortative mating is the avoidance of inbreeding and the consequent maladaptation due to the accumulation of deleterious recessive alleles (Pusey, 1987; Tregenza and Wedell, 2000; Duthie and Reid, 2016; Olsen et al., 2020), similar to the evolution of self-incompatibility in plants (reviewed by Delph and Kelly, 2014). Another common explanation comes from studies in vertebrates, where disassortative mating has been shown to increase MHC heterozygosity and overall offspring fitness (Consuegra and Garcia de Leaniz, 2008; Huchard et al., 2013).

Mate competition, the contest of polymorphic individuals within one sex to carry out a greater number of matings (Andersson, 1994), can also lead to patterns of NFDS. Perhaps the most well-known case is that of Uta stansburiana lizards. In this species, males of different color-related mating strategies increase their fitness when they are less frequent, following a rock-paper-scissors dynamic (Sinervo and Lively, 1996; Bleay et al., 2007). Each of the three existing colors is associated with a mating strategy that is in advantage or disadvantage in relation with the others. Another example of NFDS between different mating strategies comes from swordtail fish (Rios-Cardenas et al., 2018). Two strategies were observed in this species: males with large blue bodies and yellow fins that display a great courting behavior; and sneaker males, completely yellow or completely blue, that mate females without consent using coercive behavior. The higher predation faced by courting males results in an interaction between opposing mating success and survival, resulting in NFDS.

Littorina fabalis (formerly known as Littorina mariae; Sacchi and Rastelli, 1966; Reid, 1996) is a small intertidal gastropod distributed across NE Atlantic rocky shores ranging from Portugal to the Barents Sea (Reid, 1996; Rolán-Alvarez et al., 2015a). It can be found grazing upon epiphytes on brown algae of the genus Fucus (Williams, 1990). It is gonochoric and has direct development (oviparous), with females laying a mass of fertilized eggs on the algae where they develop until they emerge as crawl-aways (i.e., miniature adult snails). The species show a higher interpopulation genetic differentiation than the sibling species L. obtusata, presumably due to frequent bottlenecks during winter (more frequent storms in this season detach the thalluses of the Fucus sp. algae from the rock causing great mortality) and therefore pointing to a relatively low effective population size in the most exposed shores (Rolán-Alvarez et al., 1995). Other important selective agents in L. fabalis are crabs (e.g., Carcinus maenas), intertidal fish (e.g., Lipophrys pholis) and marine birds (Reimchen, 1979, 1989; Rolán-Alvarez et al., 1995). Among these, intertidal fishes are more likely to be visual predators than either crabs or birds (see Ekendahl, 1998 and Johannesson and Ekendahl, 2002 for experiments of predation). Mating occurs mostly during low tide when the snails are more active. The male actively follows the mucus trails of the females (Johannesson et al., 2010; Saltin et al., 2013) and upon encountering one, it crawls over the shell of the latter, by its right side, and inserts the penis into the pallial cavity.

Populations of L. fabalis, especially those with the highest densities, often show a striking shell color polymorphism (yellow, brown, orange, and olive among other colors, see Reid, 1996; Rolán-Alvarez et al., 2015a). Predation is a popular explanation for the maintenance of this polymorphism and has been thoroughly investigated in the past with attention to crypsis and apostatic selection. For example, Reimchen (1979) found that the blennid Lipophrys pholis (previously named as Blennius pholis) actively predates those colors that are conspicuous in relation to the background. In L. saxatilis and L. obtusata this was further supported (Johannesson and Ekendahl, 2002) while it was shown that crabs could not account for visual predation (Ekendahl, 1998). Crypsis and apostatic selection have also been thoroughly studied in mangrove littorinids with contrasting evidence of crabs being visual predators (Parsonage and Hughes, 2002; McKillup and McKillup, 2008).

One population from NW Spain, Abelleira, has been monitored for more than 20 years and shows a rather stable shell color polymorphism over this period (Rolán-Alvarez and Ekendahl, 1996; Rolán-Alvarez et al., 2012, 2015b). Furthermore, disassortative mating between yellow and brown individuals and sexual selection favoring brown males have been observed in this population. However, these studies did not explicitly test if a frequency-dependent selection process is behind this stable shell color polymorphism (see Figure 1), or if genetic drift alone could explain the fluctuations in color frequencies across years. Disassortative mating caused by mate choice is expected to generate frequency-dependent sexual selection (Pusey, 1987; Takahashi and Hori, 2008; Hedrick et al., 2016, 2018). In addition, disassortative mating can also be generated as a mechanism to avoid inbreeding and its consequent negative effects. Thus, disentangling between these two possible causes of disassortative mating is crucial to understand the role of frequency-dependent sexual selection in the observation of NFDS patterns.

FIGURE 1

Figure 1. Changes in the observed frequency of light (white) and dark (black) individuals of L. fabalis for a range of years in the population of Abelleira (NW Spain). Data of 1990 obtained from Rolán-Alvarez and Ekendahl (1996).

Here, we investigate the mechanisms behind the maintenance of color polymorphism in L. fabalis by analyzing data collected across 7 years from one Spanish population (Abelleira), complemented by 1-year data from two more septentrional populations located in the White Sea. More specifically, we pursue the following objectives: i) to assess the contribution of genetic drift to explain the observed frequency changes, as theoretical models indicate that the maintenance of polymorphisms could be due to stochastic processes; ii) to test for a NFDS pattern attending to the sexual component of fitness; iii) to use a multimodel inference approach based on information theory to investigate whether mate choice, mate competition, or an interaction of both could explain the observed mating frequencies; and iv) to evaluate whether inbreeding avoidance could explain the evolution of disassortative mating in this species.

Materials and Methods

Sampling Design

Sampling involved collecting random mating pairs along the four closest unmated individuals of L. fabalis (within 10 cm of diameter; Figure 2). Collections were performed within the Fucus sp. canopy area in the lower fringe of the intertidal, following the approach implemented in previous studies (Rolán-Alvarez et al., 2012, 2015b). The mating pair and the closest unmated individuals are defined in this study as a microarea (number of matings from Table 1). This design was carried out in different populations at the latitudinal limits of its distribution, Abelleira (NW Spain; 42° 48′ 0.30″N– 9° 1′ 14.87″W) and White Sea (Russia, in two locations: White Sea 1 66°20′10.73″N – 33°37′53.00″E and White Sea 2 66°20′9.80″N – 33°38′41.29″E). Although not conclusive, finding similar results would suggest that the same evolutionary forces are acting upon the color polymorphism across the species distribution range. The Spanish population was sampled once a year over a period of 7 years (2011–2017), while the Russian population was sampled in the two localities in 2015. All collections were performed in summer, since the number of adults and mating pairs increases during this season (Rolán-Alvarez and Ekendahl, 1996). Once in the laboratory, the snails were sexed and the shell color was determined by visual inspection. A printed color model was used to classify them into four colors (brown, olive, yellow and orange) following previous studies (Rolán-Alvarez et al., 2015b). In subsequent analyses we grouped these four colors into two categories, dark (brown, olive) and light (yellow, orange), as the raw spectral reflectance data from L. fabalis shells from Abelleira showed consistent similarities within these categories (Rolán-Alvarez et al., 2015b). Olive and orange were the least frequent in Abelleira, while orange and brown were the least frequent in White Sea.

FIGURE 2

Figure 2. Experimental design for the analysis of sexual selection in L. fabalis. (A) Representative shells of the different colors (brown, olive, yellow, orange) and color categories (dark, light) present in Abelleira and in the White Sea. (B) Example of heterotypic (dark/light) and homotypic (dark/dark) mating pairs on Fucus algae in the field. During the mating the male crawls on top of the female and positions on her right side. (C) Mating pairs sampling design carried out. Once a mating pair is found (center) and collected, then 4 unmated individuals are collected within a 10 cm diameter microarea (black circle).

TABLE 1

Table 1. Assortative mating (I_PSI) estimates (mean ± sd) of shell color for L. fabalis in seven consecutive years in Abelleira (AB) and in two subpopulations of the White Sea (WS-1, WS-2) in 2015.

Estimation of Effective Population Size and Inbreeding

The samples collected in Abelleira (NW Spain) in 2014 (N = 564; both the mating pair and unmated individuals) were genotyped for 8 microsatellite loci developed for L. fabalis (Carvalho et al., 2015, 2016). DNA extraction was performed following Galindo et al. (2013) and details are given in the Supplementary Material (see Supplementary Table 1).

The effective population size (N_e) was estimated from the microsatellite genotypes using the linkage disequilibrium method implemented by NeEstimator v.2.01 (Do et al., 2014) with a minor allele frequency of 0.05. This method estimates the correlation between each pair of allele frequencies using an unbiased estimation of linkage disequilibrium (Weir, 1979), and then the average correlation across loci to estimate N_e (see Waples and Do, 2008). The N_e value was subsequently used to test whether changes in allele frequencies for shell color could be explained by genetic drift alone. This assumed a model of one locus and two alleles with dark color being dominant, as suggested for the related species L. saxatilis by Johannesson and Butlin, 2017). The test was done by estimating the variance in the shift of allele frequencies in one generation,

$${\mathrm{\sigma }}_{\text{q}1}^2=\frac{{\text{p}}_0{\text{q}}_0}{2{\text{N}}_{\text{e}}}(1)$$

Wright (1922), where p₀ and q₀ are the allele frequencies at generation 0. Since L. fabalis is an annual species (see Williams, 1992) and collections in Abelleira were done across 7 years, this estimate was calculated six times, one for each annual interval. Notice that we refer here to the population of Abelleira where the species is annual, unlike snails from Russian populations, whose cycle is longer. The null hypothesis (i.e., frequency changes in a trait are caused by genetic drift) would not be rejected if the observed values of the allele frequencies across the studied years in Abelleira fell within the 95% confidence interval of the changes in frequency under the null model of genetic drift for the estimated N_e.

In order to test whether disassortative mating evolved as a response to inbreeding, molecular coancestry coefficients (molecular similarity sensu Walsh and Lynch, 2018, p. 693) were calculated between all pairs of individuals with Metapop2 (López-Cortegano et al., 2019). The average coancestry between individuals involved in mating pairs was compared with the average pairwise coancestry between all six individuals collected within the microarea (mating pair + 4 unmated individuals, Figure 2), and the significance of the difference between these was evaluated by a randomization one-way ANOVA (9999 permutations; Peres-Neto and Olden, 2001). If the mating pair coancestry was significantly lower than the microarea estimate, this would be an indication that disassortative mating evolved to avoid inbreeding.

Estimation of Assortative Mating

Assortative mating was calculated using the I_PSI estimator, which has been shown to be robust in cases of populations with variable frequencies of the traits or datasets with small sample sizes (Rolán-Alvarez and Caballero, 2000; Pérez-Figueroa et al., 2005). The I_PSI was calculated for shell color (L-light, D-dark) with the number of mating pairs from the sampled microareas as

$${\text{I}}_{\text{PSI}}=\frac{({\text{PSI}}_{\text{DD}}+{\text{PSI}}_{\text{LL}}-{\text{PSI}}_{\text{DL}}-{\text{PSI}}_{\text{LD}})}{({\text{PSI}}_{\text{DD}}+{\text{PSI}}_{\text{LL}}+{\text{PSI}}_{\text{DL}}+{\text{PSI}}_{\text{LD}})}(2)$$

where each PSI is the number of observed mating pairs divided by the expected number under random mating for each possible combination. I_PSI takes values from −1 (maximum disassortative mating) to 1 (maximum positive assortative mating) with 0 equaling absolute panmixia. The estimation of assortative mating in the wild may be biased due to the scale of choice effect (SCE). This effect biases the value of assortative mating (typically toward positive values) due to a mismatch between the sampling scale and the real scale of mate choice in nature (Rolán-Alvarez et al., 2015b). Therefore, to minimize SCE bias, the estimates of I_PSI were obtained by pooling mating pairs into homogeneous color groups based on the light color frequency within their microareas: 1) Light group (≥ 66%), 2) Mixed group (33–66%), 3) Dark group (≤ 33%). The I_PSI values from these groups were weight-averaged ($\overline{I_{PSI}}$) each year to reduce the effect of samples’ size differences. Significance for I_PSI values was assessed through bootstrapping (10000 iterations) using JMATING V. 1.0.8 (Carvajal-Rodríguez and Rolán-Alvarez, 2006), while significance for $\overline{I_{PSI}}$ was estimated using a one sample t-test (null hypothesis I_PSI = 0).

Estimation of Sexual Fitness

For each year, sexual fitness was estimated in Abelleira and both locations in the White Sea by comparing shell colors before (unmated individuals) and after (mated individuals) selection (Arnold and Wade, 1984; Rolán-Alvarez and Caballero, 2000; López-Cortegano et al., 2020). The individuals belonged exclusively to one of these two groups. Thus, mated and unmated individuals were statistically independent. We carried out two different statistical methods. First, we used logistic regression to study how variation in shell color (independent variable) could determine fitness (dependent variable). In this case, fitness was considered a binary variable with a value of 0 (unmated individuals) or 1 (mated individuals). This method is commonly employed when fitness is binary and the character under study is a quantitative trait (Brodie et al., 1995; Janzen and Stern, 1998). Significance was evaluated through the Wald test (Sokal and Rohlf, 1995, p. 770). Second, we used the cross-product fitness estimator (W), within sex and sample. This is the most suitable estimator for qualitative traits (Knoppien, 1985; Rolán-Alvarez and Caballero, 2000; Martínez-Rodríguez et al., 2019), and has been shown to estimate long-term fitness better than other methods (see Brommer et al., 2004). One of its major advantages is its frequency independence, which makes it very suitable for studying frequency dependent mechanisms (Martínez-Rodríguez et al., 2019). As an example, the W estimator for dark individuals (D) was calculated as

$$W_{\text{D}}=\frac{{\text{AS}}_{\text{D}}\times {\text{BS}}_{\text{L}}}{{\text{BS}}_{\text{D}}\times {\text{AS}}_{\text{L}}}(3)$$

where AS and BS are the frequency of individuals after and before selection (mated and unmated, respectively) for D (dark) and L (light) colors. The W-values and their significance were assessed through bootstrapping (after 10000 iterations) using the software JMATING, by resampling the observed color frequencies before (BS) and after (AS) selection independently.

Linear regression of the frequency of a color on sexual fitness (given as the logarithm of W in base 10) was carried out to find potential patterns of NFDS (i.e., that fitness decreases with frequency). A single linear regression ANOVA was used to assess significance. Parametric statistical tests were implemented in SPSS Statistics v.22.0 (IBM Corp.).

Contribution of Mate Choice and Mate Competition by Multimodel Comparison

We used multimodel inference (see Burnham et al., 2011) to estimate the contribution of mate choice and mate competition to the observed patterns of shell color polymorphism. These analyses were carried out with the software InfoMating v 0.2 (Carvajal-Rodríguez, 2020), which estimates the mutual mating propensities or MMPs (i.e., the tendency to mate of each particular mating combination: LL, LD, DL and DD) expected from the frequency of mating pairs and unmated individuals under each specific model tested (Carvajal-Rodríguez, 2018). The Akaike information criterion corrected for low sample size (AICc) was used to determine the weight of each model to select the best one (Burnham et al., 2011). Since the MMPs of each model are defined by mate choice and/or mate competition parameters (see Figures 2–4 in Carvajal-Rodríguez, 2020), this allowed us to identify which type of parameter is contributing the most to the observed frequencies. We analyzed Abelleira and the White Sea pooling data from the different years and populations, respectively. Additionally, a second analysis was made for pools of homogeneous groups in each population. Given the number and complexity of the obtained models (see Supplementary Table 2), we classified them into mate choice, mate competition and mixed models, attending to the parameters included in each. For example, mate competition parameters were assumed when the parameter is present in either the whole column or row from Supplementary Table 2, otherwise the parameter was assumed to be mate choice.

FIGURE 3

Figure 3. Expected changes in frequency for the color polymorphism, assuming genetic drift. Red dots represent the allele frequency observed in the year at the x axis while black dots represent the observed frequency of alleles the following year (e.g., the black dot on 2011 corresponds to the frequency in 2012). The bars indicate the 95% confidence interval of the expected frequencies under genetic drift, using the estimated effective population size (N_e) on equation 1 (see section “Materials and Methods”). Thus, black dots outside confidence intervals indicate that the change in frequency from 1 year to the next cannot be explained by genetic drift. Plots indicate results for N_e = 600 (A), 394 (B), and 1076 (C).

FIGURE 4

Figure 4. Relationship between fitness (Log10 W) and the percentage of light individuals for each sex of L. fabalis. Each dot represents a year of study in the Abelleira population (2011–2017) while triangles indicate each of the two populations sampled in the White Sea in 2015. Regression lines are shown for females (solid; R² = 0.667) and males (dotted; R² = 0.004). Notice that, since the W estimator for a given morph is relative to the other morph, using it for one or other morph will yield the same results.

Results

Role of Genetic Drift and Inbreeding

The estimated mean effective population size (N_e) in Abelleira was 600, with 95% confidence interval of 394 and 1076. These values were used to calculate the expected variance in the change of allele frequency over one generation (from 1 year to the next) assuming only genetic drift (Figure 3). Using the mean N_e, the allele frequency change between some years could not be explained by drift alone (2011, 2013, and 2014). Similar results were observed with the 95% confidence interval values of N_e (Figure 3).

The average molecular coancestry coefficient for all mating pairs was 0.345 ± 0.078 (mean ± SD) and the overall average for the pairwise comparisons within microarea (including mated and unmated individuals) was 0.346 ± 0.044. Thus, mated individuals did not show significant differences in relatedness compared to surrounding individuals (ANOVA, F₁,₁₈₈ = 0.03; p = 0.859).

Assortative Mating

All I_PSI estimates obtained from the studied populations were negative (except the dark homogeneous group in location 2 of the White Sea; Table 1), indicating a higher frequency of mates between individuals with different shell color (heterotypic) than with the same color (homotypic) in both populations. Estimates of I_PSI for each of the three homogeneous groups were variable, with values ranging from 0.00 to −0.59. Only 9 (6 in Abelleira and 3 in Russia) out of 27 I_PSI estimates were high in absolute value and significant under bootstrapping (p ≤ 0.05). This is likely due to the substantial reduction in the number of individuals when they are grouped in homogeneous units. After this grouping, the estimates in Abelleira varied from −0.18 to −0.40 with an overall mean of −0.29, which was significant (t₂ = 7.07; p = 0.019), pointing to consistent disassortative mating in this population. Average values in the White Sea were not significant but were very similar to those obtained in Abelleira, thus revealing the same pattern for unrelated populations. This suggests that disassortative mating may be widespread across the two areas.

Sexual Fitness Estimation

Estimates of sexual fitness using the cross-product fitness estimator (W) and the Wald test for each sex are shown in Table 2. Sexual fitness for light males was always lower than that for dark males (with the exception of White Sea 1), but the difference was only significant for Abelleira in 2011, 2012, and 2014. In the remaining years, light males’ fitness was only slightly lower than that of dark males. On the contrary, light females’ sexual fitness tended to be larger compared to that of dark ones, with the exception of White Sea 2 and Abelleira in 2012 (P = 0.002). The cross-product estimator had more statistical power and showed more significant cases than the Wald test under logistic regression. The results were, however, similar from both approaches.

TABLE 2

Table 2. Sexual fitness estimates for females and males using the cross-product estimator (W ± sd) and logistic regression slope (B ± sd) for light individuals.

Figure 4 shows the regression of sexual fitness (as log₁₀ of W) of light color on the percentage of light color individuals. Notice that the result for the dark color is the inverse but would yield the same observations and conclusions. Importantly, a significant pattern of NFDS was observed for females (F₁,₇ = 14.00; P = 0.007) but not in males (F₁,₇ = 0.03; P = 0.872; Figure 4).

Mechanism Inference by Multimodel Comparison

The relative weight of the tested models, and the mutual mating propensities for the best-fit model for each population are presented in Figure 5 (Abelleira and White Sea). In Abelleira, models including a parameter of mate choice and mate competition accounted for 66% of the total weight, while models including only mate choice accounted for 29%. Within homogeneous groups these results were very similar except in the light group where models with only mate choice parameters accounted for 69% of the total weight. In the White Sea, both types of models (mate choice and mixed) presented similar weights. While dark and mixed homogeneous groups showed mate choice models to account for a major part of the total weight (62 and 77%, respectively), mixed models accounted for a 66% of the total weight in the light group.

FIGURE 5

Figure 5. Mutual mating propensities obtained for the best models in Abelleira and the White Sea in both total numbers and homogeneous groups, pooled across all years (Abelleira) and populations (White Sea). Male in columns and female in rows. Values within boxes are the best model mutual mating propensities for the four different mating combinations. Pie charts represent the sum of weights for mate choice, mate competition, mixed and random mating models (see Supplementary Table 2).

The best-fit model in Abelleira included both mate choice and mate competition parameters. The estimate of the mate choice parameter was 0.85 ± 0.069 (Figure 5). Thus, mating between homotypic pairs (DD or LL) had 15% less propensity than mating between heterotypic pairs (DL or LD). This preference decreased to 0.55 ± 0.01 within homogeneous groups (a 45% reduction; Figure 5). The estimate of the mate competition parameter was 1.31 ± 0.106 (Figure 5), which implies a sexual advantage of 31% for dark males over light ones, and 1.45 ± 0.23 within homogeneous groups (Figure 5). Results for the light homogeneous group differed somewhat as no mate competition parameters were defined in the best-fit model, but rather a choice component was found in favor of heterotypic mating pairs and against LL ones.

The best-fit model for the White Sea was defined by mate choice only. The estimated mate choice parameter for this model was 0.77 ± 0.06 (a reduction of 23% in the mating propensity of homotypic pairs). The best-fit models within homogeneous groups were somewhat different. The dark group showed the same pattern as in Abelleira, with a mate choice parameter of 0.42 (a reduction of 58% in the mating propensity of homotypic pairs). The mixed group showed two different mate choice parameters. The first parameter reduces the mating propensity of DD mating pairs, while the second increases it in mating pairs of the DL kind (first letter corresponds to the male color). The light group showed that dark males had 51% less propensity to mate and there was a choice component against LL mating pairs.

Discussion

Littorina fabalis color polymorphism has been shown to be stable across years in the current and previous studies (Figure 1; Rolán-Alvarez and Ekendahl, 1996; Rolán-Alvarez et al., 2012). This is in spite of frequent bottlenecks during winter, at least in NW Spain (Rolán-Alvarez et al., 1995), which could aggravate the effects of genetic drift causing allelic fixation. Here we have presented several examples that support the contribution of a mechanism of NFDS for the maintenance of this polymorphism across the species distribution. Our results indicate that such a mechanism depends on precopulatory mate choice. We have also found some results pointing to directional sexual selection in dark male color. However, such directional selection would also cause a fixation of alleles in the long run and thus its contribution to polymorphism maintenance is not likely. Possible alternative scenarios and implications of the present study are discussed below.

The Stochastic Scenario

Genetic variation is expected to be quickly eroded under stochastic forces when the populations have low effective sizes (Cook, 1992; Hoffman et al., 2006). L. fabalis is a species which suffers strong bottlenecks during winter, possibly due to the detachments of Fucus thalluses during storms (Williams, 1996). This could cause severe reductions in effective size (Rolán-Alvarez et al., 1995) and therefore also presumably in color frequencies. However, in the 7 years of study in Abelleira, frequency fluctuations do not indicate such severe reductions (Figure 1). Therefore, we checked whether the observed changes in allele frequencies could be just explained by genetic drift without the need to claim for any selective mechanism.

Our analysis of genetic drift could be questioned by the facts that the genetic basis of color in Littorinids, including L. fabalis, is not yet completely resolved (Reimchen, 1989; Ekendahl, 1995; McQuaid, 1996; Kozminsky, 2014; Johannesson and Butlin, 2017), and because we have combined all the observed colors in just two categories. Our decision for this latter was based on phenotypic similarity (see Rolán-Alvarez et al., 2015b) and because orange and olive colors had a low frequency in all years of study. However, our assumption that shell color is coded by just one locus is conservative, because the variance in gene frequency across years would be further diminished considering two or more loci under different inheritance models (see simulations by Estévez, 2018, Chapter 3). Our results show that, although genetic drift cannot be discarded when fluctuations in color frequencies are small, genetic drift alone cannot explain the maintenance of the polymorphism over time in our study population. This is in agreement with findings in L. saxatilis where genetic drift cannot explain the persistence of rare colors (Johannesson and Butlin, 2017). In contrast, previous results in Littoraria pallescens showed genetic drift to be the main driver of polymorphism maintenance in the absence of selection (Cook, 1992). However, in other mangrove littorinids, heterogeneous natural selection has gained support as a more likely explanation for such a maintenance (Parsonage and Hughes, 2002). Our results also suggest that there is a frequency threshold over which selection acts to prevent allele fixation. That is, when a morph’s frequency surpasses a boundary, selective mechanisms occur that outrun the relative strength of drift. This has also been proposed in the candy-stripe spider (Oxford, 2005) and points to the alternating roles of selection and drift in maintaining polymorphisms (see O’Hara, 2005).

Overall, this suggests that genetic drift is less important than previously thought in this population and species. There may be two reasons for this. First, these populations are dense and stable with a limited population reduction during winter (personal observation). Second, in the related species L. saxatilis, quick recovery of the populations’ genetic diversity has been observed after bottleneck episodes (Johannesson and Johannesson, 1995; Piñeira et al., 2008). Surprisingly, the rarest colors in L. saxatilis persisted over time in spite of strong population reductions (Johannesson and Butlin, 2017). The fact that pregnant females in this species carry offspring from multiple males may prevent allele loss by genetic drift (Panova et al., 2010). We have no direct evidence of such a process in L. fabalis but its evolutionary proximity to L. saxatilis, and the fact that females mate multiple times, suggests that they may spawn multiple sired offspring. This has in fact been studied in the sibling species L. obtusata where between 4 and 6 sires were found to contribute to females’ broods (Paterson et al., 2001). Therefore, it is likely that each L. fabalis female spawns a genetically diverse brood, and consequently, this can contribute to maintaining a larger effective size despite adult population fluctuations.

Balancing Mechanism Maintaining the L. fabalis Color Polymorphism

Even though there are some studies showing evidence of disassortative mating (Takahashi and Hori, 2008; Hedrick et al., 2016, 2018), this is a rare occurrence in natural populations and has been associated with type I errors (Jiang et al., 2013). Previous studies hypothesized that such a lack of cases could be due to the absence of spatial and temporal dynamics incorporated in the analyses (Rolán-Alvarez et al., 2015b; Fargevieille et al., 2017). Here we have accounted for the scale bias effect (SCE; see Rolán-Alvarez et al., 2015b) and have shown that disassortative mating was observed in all the studied years and in populations at both extremes of the species distribution. Preliminary laboratory evidence also indicates disassortative mating in this species and for the same colors from Abelleira (I_PSI = −0.29 ± 0.102; p = 0.0054; unpublished results). In addition, two other NW populations of L. fabalis from Spain show a high and negative assortative mating for shell color (unpublished results). Thus, L. fabalis provides a good example of consistent disassortative mating for shell color.

Negative mate choice for color in L. fabalis (see below) could be the cause of disassortative mating, as observed by experimental trials in various study systems. For instance, in the butterfly Heliconius numata, experimental trials showed a tendency for disassortative mating which was hypothesized to maintain wing color polymorphisms in spite of positive frequency dependent selection by predators (Chouteau et al., 2017). Another well-known example is the white-throated sparrow, where disassortative mating is almost absolute in the wild (Hedrick et al., 2018) and experimental lab evidence supports negative mate choice (Houtman and Falls, 1994). However, as has been observed in several studies (see for example López-Cortegano et al., 2020, and references therein) results from laboratory trials cannot be readily extrapolated to natural populations. Therefore, one advantage of our study is that it shows what is going on directly in the field, and so it has an evolutionary relevance. On the contrary, a problem of field studies is that they can be affected by several source of bias (Alatalo et al., 1986; Crespi, 1989; Taborsky et al., 2009). Nevertheless, we provided estimates corrected for the SCE bias and we also estimated the strength of mate choice and sexual selection effects by independent methods, and both confirm the present interpretation.

First, our best fit models (Figure 5) indicated negative mate choice in almost all cases (the diagonal for heterotypic pairs showed greater values than that of homotypic pairs). Second, the trend observed in Figure 4 for female sexual fitness strongly showed a pattern of NFDS, which is an expected outcome under male mate choice with preference for the dissimilar color (see Kokko et al., 2007; Wellenreuther et al., 2014). This is in accordance with existing evidence which shows that mate choice in Littorina is indeed preferentially exerted by males (Erlandsson and Kostylev, 1995; Ng et al., 2013, 2019; Saltin et al., 2013). Thus, for L. fabalis there seems to be a pattern of NFDS in which the most frequent female color mate less often. Taken together, these pieces of evidence indicate the occurrence of mate choice for dissimilar color in L. fabalis. Mate competition could be affecting our results but our multimodel inference accounted for its contribution on the observed mating propensities. The multimodel inference results did point to a greater mating advantage in dark males, at least in Abelleira (up to a 31%), suggesting that dark males outcompete light males. Such occurrence, however, only explains why there is a greater frequency of mating pairs where a dark male is involved but not why there is a greater frequency of heterotopic pairs. In other words, that dark males have greater sexual fitness does not explain the evidence that points to negative mate choice.

Potential Causes for Observed Patterns

A common reason to mate preferentially with dissimilar morphs is to avoid inbreeding, which can cause fitness disadvantages in the offspring (Keller and Waller, 2002; Charlesworth and Willis, 2009; Duthie and Reid, 2016). Thus, patterns of disassortative mating are not rare in cases where homozygotes are at selective disadvantage (Tregenza and Wedell, 2000; Mays and Hill, 2004). Here, we tested if relatedness of mating pairs were significantly lower than that of unmated individuals from the same microarea. We did not find significant differences and our results do not support that overall preference for dissimilar morphs have evolved to avoid inbreeding, at least in the studied population. Previous studies in L. saxatilis yielded the same observation (Ng and Johannesson, 2016). However, the level of relatedness at some microsatellite loci does not necessarily reflect the relatedness for the whole genome. Mating can therefore respond to differences in certain genomic regions. For example, in the white throated sparrow, homozygotes for the chromosome allele 2^m suffer low fitness because of the exposition of deleterious alleles in that chromosome (Tuttle et al., 2016). Additionally, in this species, disassortative mating between morphs is almost absolute (Hedrick et al., 2018) and explains the polymorphism maintenance. Both pieces together indicate that negative choice may have evolved to reduce the number of less fit homozygous 2^m types in this bird species. In an analogous way to the white-throat sparrows, the observation that light mating pairs is less frequent than the remaining mating combinations could suggest that homozygosity for light color may have deleterious effects in L. fabalis. A thorough analysis of the genetic basis of shell color in L. fabalis is fundamental for studying this hypothesis.

Shell color might not be the trait under direct choice but rather a correlated trait. Association between color and other traits is typical in nature (McKinnon and Pierotti, 2010). Sacchi (1969, 1972) proposed several such associations including salinity and desiccation tolerance. Association between these traits can be due to genetic inversions which can originate linkage disequilibrium between locally adapted loci. Indeed, L. fabalis studies supported the existence of an inversion which affected size, the allozyme Ark Ginase, and a RAPD locus genotype (Johannesson and Mikhailova, 2004; Kemppainen et al., 2011). This inversion would create a locally adapted supergene codifying for several traits. In this scenario, the occurrence of disassortative mating would benefit the occupation of a wider variety of habitats (van Nes, 1998). Of the possible traits associated with color, variation at loci determining the immune system is a promising candidate. In invertebrates in general, color is usually linked to immune responses (see Scheil et al., 2013 and references therein). In turn, mating and immune characteristics have been found to be greatly linked (Lawniczak et al., 2007). In L. fabalis, male mate choice for different color could have evolved as a side effect of mate choice for individuals with dissimilar immune characteristics, in the same fashion as the MHC complex in vertebrates (e.g., Consuegra and Garcia de Leaniz, 2008; Huchard et al., 2013). In this case, the benefit would be to produce more fit heterozygous offspring for immunity genes. It is still uncertain whether male choice can evolve to improve offspring quality. While this has been suggested in house mice, Mus domesticus (Gowaty et al., 2003), other studies testing this hypothesis have found either a negative (Swierk and Langkilde, 2019) or absence of correlation (Bertram et al., 2017).

Male mate choice can also evolve to exploit other fitness benefits. In many cases, males prefer those females that denote greater fecundity (e.g., Edward and Chapman, 2012; Jaworski et al., 2018). In Littorina, an honest indicator of fecundity is female size and males tend to mate with larger females (Saltin et al., 2013, and references therein). Color can also be an indicator of female condition and thus inform males to increase fitness (Amundsen and Forsgren, 2001; Weiss, 2006). We have no evidence that male mate choice for dissimilar colored females is related to increase fecundity in L. fabalis. To explore this, an experimental approach could be performed where numbers of eggs resulting from homotypic and heterotypic mating pairs are compared. A greater number of eggs from heterotypic mating pairs would indicate male mate choice to increase fecundity in this species.

The observed greater sexual fitness in dark males may be driven by temperature and risks to desiccation, given the known relationship between shell color and sunlight incidence (see Heath, 1975; Cook and Freeman, 1986). Shell color in Littorina could be in fact associated to thermal physiology (reviewed by Rolán-Alvarez et al., 2015a). For example, evidence that sunlight is acting differently on dark and light shells has been claimed in the sibling species L. obtusata (Sergievsky, 1992; Phifer-Rixey et al., 2008) and other gastropods (Cepaea nemoralis, Heath, 1975; Littoraria pallescens, Cook and Freeman, 1986; Nucella lapillus, Etter, 1988). An increment in temperature in dark males may increase their sexual fitness. Elevations in sexual selection intensity due to temperature have in fact been shown in other organisms (Kvarnemo, 1996; Allen et al., 2012). Alternatively, an increase in temperature may trigger behaviors to avoid desiccation (see Miller and Denny, 2011). In our study system, it may cause dark males to move more actively for shelter and thus increase the rate of encounter with females. Interestingly, such sexual fitness advantage is observed in dark males in Abelleira (best models in total results in Figure 5), which represents the southern limit of L. fabalis distribution and where the temperatures are higher, but not in the White Sea, in agreement with the temperature hypothesis.

Conclusion

Shell color polymorphism in L. fabalis is persistent despite expected strong genetic drift. We have otherwise found that color polymorphism in this species follows a NFDS pattern when studying the sexual component of fitness. In particular, this pattern was only present in females. This, together with the occurrence of disassortative mating, is indicative of male mate choice, which is common in species of the Littorina genus. Since these patterns result from observations in the wild, we used multimodel inference to assess the possible role of mate choice and/or mate competition. The inference strongly supports mate choice but also indicates mate competition in males. Several hypotheses have been presented here that can account for either mate choice or mate competition. However, the connection, if any, between them is still to be elucidated and needs further research. In particular, understanding the underlying genetics of color and what is the correlation between different traits, could cast light into these mechanisms and explain the underlying biological causes. Experimental trials in the lab for both mate choice and mate competition can provide further information concerning these biological causes. This study has confirmed a rare pattern, disassortative mating, to occur in very distant populations of L. fabalis, and opens the door for further research into understanding the biological mechanisms behind polymorphism maintenance and mate choice in marine systems.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: doi: 10.6084/m9.figshare.13050494.v1 and doi: 10.6084/m9.figshare.13050509.v1.

Author Contributions

DE performed samplings in Galicia, analyzed the data, designed the tables and figures, wrote and reviewed the content of the manuscript, and designed it for publication. EK performed samplings in the White Sea, provided the data from that location, and reviewed the content of the manuscript. AC-R provided software resources for the analyses of data and reviewed the content of the manuscript. AC contributed to the analysis of genetic drift and reviewed the manuscript. RF reviewed the manuscript, added valuable criticism toward methods, and contributed to improve the writing. JG and ER-A made equal contributions to the manuscript, performed samplings in Galicia, contributed to the writing of the manuscript, and reviewed it for publication. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by Xunta de Galicia (ED431C 2016-037), FONDOS FEDER (unha maneira de facer Europa) and Ministerio de Economía, Industria y Competitividad, Agencia Estatal de Investigación (CGL2016-75904-C2-1-P). JG was funded by a JIN project (Jóvenes Investigadores sin vinculación o con vinculación temporal, Ministerio de Ciencia, Innovación y Universidades, RTI2018-101274-J-I00). RF was financed by the FEDER Funds through the Operational Competitiveness Factors Program – COMPETE and by National Funds through FCT – Foundation for Science and Technology though a research contract within the scope of the project “Hybrabbid” (PTDC/BIA-EVL/30628/2017 and POCI-01-0145-FEDER-030628). UVigo Marine research centre is one of the five research centers from all fields of knowledge recognized as “Galician Singular Research Centres” by the Galician Regional Government through the “Excellence in Research (INUGA)” Program from the Regional Council of Culture, Education and Universities. The Centre’s activity is co-funded by the European Union through the ERDF Operational Program Galicia 2014–2020. The research work of EK was carried out as part of the State Task RF (no. AAAAA19-119022690122-5).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Mary Riádigos for administrative help. We thank José Bella, Roger Butlin, and Mauro Santos for improving later versions of the manuscript and general discussion. We also thank Lisa Germany for improving the language and fluency of the manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.614237/full#supplementary-material

References

Alatalo, R. V., Lundberg, A., and Glynn, C. (1986). Female pied flycatchers choose territory quality and not male characteristics. Nature 323, 152–153. doi: 10.1038/323152a0

CrossRef Full Text | Google Scholar

Allen, B. J., Rodgers, B., Tuan, Y., and Levinton, J. S. (2012). Size-dependent temperature and desiccation constraints on performance capacity: implications for sexual selection in a fiddler crab. J. Exp. Mar. Biol. Ecol. 438, 93–99. doi: 10.1016/j.jembe.2012.09.009

CrossRef Full Text | Google Scholar

Amundsen, T., and Forsgren, E. (2001). Male mate choice selects for female coloration in a fish. Proc. Natl. Acad. Sci. U.S.A. 98, 13155–13160. doi: 10.1073/pnas.211439298

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, M. B. (1994). Sexual Selection. Princeton, NJ: Princeton University Press.

Google Scholar

Arnold, S. J., and Wade, M. J. (1984). On the measurement of natural and sexual selection: applications. Evolution 38, 720–734. doi: 10.2307/2408384

CrossRef Full Text | Google Scholar

Ayala, F. J., and Campbell, C. A. (1974). Frequency-dependent selection. Annu. Rev. Ecol. Syst. 5, 115–138. doi: 10.1146/annurev.es.05.110174.000555

CrossRef Full Text | Google Scholar

Bell, G. (2008). Selection, the Mechanism of Evolution. Oxford: Oxford University Press.

Google Scholar

Bertram, S. M., Harrison, S. J., Ferguson, G. L., Thomson, I. R., Loranger, M. J., Reifer, M. L., et al. (2017). What is driving male mate preference evolution in Jamaican field crickets? Ethology 123, 793–799. doi: 10.1111/eth.12649

CrossRef Full Text | Google Scholar

Bleay, C., Comendant, T., and Sinervo, B. (2007). An experimental test of frequency-dependent selection on male mating strategy in the field. Proc. R. Soc. B 274, 2019–2025. doi: 10.1098/rspb.2007.0361

PubMed Abstract | CrossRef Full Text | Google Scholar

Brodie, E. D., Moore, A. J., and Janzen, F. J. (1995). Visualizing and quantifying natural selection. Trends Ecol. Evol. 10, 313–318. doi: 10.1016/S0169-5347(00)89117-X

CrossRef Full Text | Google Scholar

Brommer, J. E., Gustafsson, L., Pietiäinen, H., and Merilä, J. (2004). Single-generation estimates of individual fitness as proxies for long-term genetic contribution. Am. Nat. 163, 505–517. doi: 10.1086/382547

PubMed Abstract | CrossRef Full Text | Google Scholar

Burnham, K. P., Anderson, D. R., and Huyvaert, K. P. (2011). AIC model selection and multimodel inference in behavioral ecology: some background, observations, and comparisons. Behav. Ecol. Sociobiol. 65, 23–35. doi: 10.1007/s00265-010-1029-6

CrossRef Full Text | Google Scholar

Carvajal-Rodríguez, A. (2018). Non-random mating and information theory. Theor. Popul. Biol. 120, 103–113. doi: 10.1016/j.tpb.2018.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvajal-Rodríguez, A. (2020). Multi-model inference of non-random mating from an information theoretic approach. Theor. Popul. Biol. 131, 38–53. doi: 10.1016/j.tpb.2019.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvajal-Rodríguez, A., and Rolán-Alvarez, E. (2006). JMATING: a software for the analysis of sexual selection and sexual isolation effects from mating frequency data. BMC Evol. Biol. 6:40. doi: 10.1186/1471-2148-6-40

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, J., Pereira, C., Sotelo, G., Costa, D., Galindo, J., and Faria, R. (2015). De novo isolation of 17 microsatellite loci for flat periwinkles (Littorina fabalis and L. obtusata) and their application for species discrimination and hybridization studies. J. Molluscan Stud. 81, 421–425. doi: 10.1093/mollus/eyv014

CrossRef Full Text | Google Scholar

Carvalho, J., Sotelo, G., Galindo, J., and Faria, R. (2016). Genetic characterization of flat periwinkles (Littorinidae) from the Iberian Peninsula reveals interspecific hybridization and different degrees of differentiation. Biol. J. Linn. Soc. 118, 503–519. doi: 10.1111/bij.12762

CrossRef Full Text | Google Scholar

Charlesworth, D., and Willis, J. H. (2009). The genetics of inbreeding depression. Nat. Rev. Genet. 10, 783–796. doi: 10.1038/nrg2664

PubMed Abstract | CrossRef Full Text | Google Scholar

Chouteau, M., Llaurens, V., Piron-Prunier, F., and Joron, M. (2017). Polymorphism at a mimicry supergene maintained by opposing frequency-dependent selection pressures. Proc. Natl. Acad. Sci. U.S.A. 114, 8325–8329. doi: 10.1073/pnas.1702482114

PubMed Abstract | CrossRef Full Text | Google Scholar

Consuegra, S., and Garcia de Leaniz, C. (2008). MHC-mediated mate choice increases parasite resistance in salmon. Proc. R. Soc. B 275, 1397–1403. doi: 10.1098/rspb.2008.0066

PubMed Abstract | CrossRef Full Text | Google Scholar

Cook, L. M. (1992). The neutral assumption and maintenance of color morph frequency in mangrove snails. Heredity 69, 184–189. doi: 10.1038/hdy.1992.112

CrossRef Full Text | Google Scholar

Cook, L. M., and Freeman, P. M. (1986). Heating properties of morphs of the mangrove snail, Littoraria pallescens. Biol. J. Linn. Soc. 29, 295–300. doi: 10.1111/j.1095-8312.1986.tb00281.x

CrossRef Full Text | Google Scholar

Crespi, B. J. (1989). Causes of assortative mating in arthropods. Anim. Behav. 38, 980–1000. doi: 10.1016/S0003-3472(89)80138-1

CrossRef Full Text | Google Scholar

Delph, L. F., and Kelly, J. K. (2014). On the importance of balancing selection in plants. New Phytol. 201, 45–56. doi: 10.1111/nph.12441

PubMed Abstract | CrossRef Full Text | Google Scholar

Do, C., Waples, R. S., Peel, D., Macbeth, G. M., Tillett, B. J., and Ovenden, J. R. (2014). NeEstimator v2: Re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol. Ecol. Resour. 14, 209–214. doi: 10.1111/1755-0998.12157

PubMed Abstract | CrossRef Full Text | Google Scholar

Duthie, A. B., and Reid, J. M. (2016). Evolution of inbreeding avoidance and inbreeding preference through mate choice among interacting relatives. Am. Nat. 188, 651–667. doi: 10.1086/688919

PubMed Abstract | CrossRef Full Text | Google Scholar

Edward, D. A. (2015). The description of mate choice. Behav. Ecol. 26, 301–310. doi: 10.1093/beheco/aru142

PubMed Abstract | CrossRef Full Text | Google Scholar

Edward, D. A., and Chapman, T. (2012). Measuring the fitness benefits of male mate choice in Drosophila melanogaster. Evolution 66, 2646–2653. doi: 10.1111/j.1558-5646.2012.01648.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ekendahl, A. (1995). Microdistribution in the field and habitat choice in aquaria by color morphs of Littorina mariae sacchi & rastelli. J. Molluscan Stud. 61, 249–256. doi: 10.1093/mollus/61.2.249

CrossRef Full Text | Google Scholar

Ekendahl, A. (1998). Colour polymorphic prey (Littorina saxatilis Olivi) and predatory effects of a crab population (Carcinus maenas L.). J. Exp. Mar. Biol. Ecol. 222, 239–246. doi: 10.1016/S0022-0981(97)00147-0

CrossRef Full Text | Google Scholar

Erlandsson, J., and Kostylev, V. (1995). Trail following, speed and fractal dimension of movement in a marine prosobranch, Littorina littorea, during a mating and a non-mating season. Mar. Biol. 122, 87–94. doi: 10.1007/BF00349281

CrossRef Full Text | Google Scholar

Estévez, D. (2018). Causas del Polimorfismo de Color en Poblaciones Naturales de Littorina fabalis. Ph.D. thesis, Universidade de Vigo, Vigo.

Google Scholar

Etter, R. J. (1988). Physiological stress and color polymorphism in the intertidal snail Nucella lapillus. Evolution 42, 660–680. doi: 10.2307/2408858

CrossRef Full Text | Google Scholar

Fargevieille, A., Grégoire, A., Charmantier, A., del Rey Granado, M., and Doutrelant, C. (2017). Assortative mating by colored ornaments in blue tits: space and time matter. Ecol. Evol. 7, 2069–2078. doi: 10.1002/ece3.2822

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzpatrick, B. M., Shook, K., and Izally, R. (2009). Frequency-dependent selection by wild birds promotes polymorphism in model salamanders. BMC Ecol. 9:12. doi: 10.1186/1472-6785-9-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzpatrick, M. J., Feder, E., Rowe, L., and Sokolowski, M. B. (2007). Maintaining a behavior polymorphism by frequency-dependent selection on a single gene. Nature 447, 210–212. doi: 10.1038/nature05764

PubMed Abstract | CrossRef Full Text | Google Scholar

Galindo, J., Martínez-Fernández, M., Rodríguez-Ramilo, S. T., and Rolán-Alvarez, E. (2013). The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail. J. Evol. Biol. 26, 1472–1487. doi: 10.1111/jeb.12152

PubMed Abstract | CrossRef Full Text | Google Scholar

Gosden, T. P., and Svensson, E. I. (2007). Female sexual polymorphism and fecundity consequences of male mating harassment in the wild. PLoS One 2:e580. doi: 10.1371/journal.pone.0000580

PubMed Abstract | CrossRef Full Text | Google Scholar

Gowaty, P. A., Drickamer, L. C., and Schmid-Holmes, S. (2003). Male house mice produce fewer offspring with lower viability and poorer performance when mated with females they do not prefer. Anim. Behav. 65, 95–103. doi: 10.1006/anbe.2002.2026

CrossRef Full Text | Google Scholar

Heath, F. J. (1975). Colour, sunlight and internal temperatures in the land-snail Cepaea nemoralis (L.). Oecologia 19, 29–38. doi: 10.1007/BF00377587

PubMed Abstract | CrossRef Full Text | Google Scholar

Hedrick, P. W. (1983). Genetics of Population. Boston, MA: Jones and Bartlett publishers.

Google Scholar

Hedrick, P. W., Smith, D. W., and Stahler, D. R. (2016). Negative-assortative mating for color in wolves. Evolution 70, 757–766. doi: 10.1111/evo.12906

PubMed Abstract | CrossRef Full Text | Google Scholar

Hedrick, P. W., Tuttle, E. M., and Gonser, R. A. (2018). Negative-assortative mating in the white-throated sparrow. J. Hered. 109, 223–231. doi: 10.1093/jhered/esx086

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffman, E. A., Schueler, F. W., Jones, A. G., and Blouin, M. S. (2006). An analysis of selection on a color polymorphism in the northern leopard frog. Mol. Ecol. 15, 2627–2641. doi: 10.1111/j.1365-294X.2006.02934.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Houtman, A. M., and Falls, J. B. (1994). Negative assortative mating in the white-throated sparrow, Zonotrichia albicollis: the role of mate choice and intra-sexual competition. Anim. Behav. 48, 377–383. doi: 10.1006/anbe.1994.1251

CrossRef Full Text | Google Scholar

Huchard, E., Baniel, A., Schliehe-Diecks, S., and Kappeler, P. M. (2013). MHC-disassortative mate choice and inbreeding avoidance in a solitary primate. Mol. Ecol. 22, 4071–4086. doi: 10.1111/mec.12349

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes, K. A., Houde, A. E., Price, A. C., and Rodd, F. H. (2013). Mating advantage for rare males in wild guppy populations. Nature 503, 108–110. doi: 10.1038/nature12717

PubMed Abstract | CrossRef Full Text | Google Scholar

Janzen, F. J., and Stern, H. S. (1998). Logistic regression for empirical studies of multivariate selection. Evolution 52, 1564–1571. doi: 10.2307/2411330

CrossRef Full Text | Google Scholar

Jaworski, K. E., Lattanzio, M. S., Hickerson, C.-A. M., and Anthony, C. D. (2018). Male mate preference as an agent of fecundity selection in a polymorphic salamander. Ecol. Evol. 8, 8750–8760. doi: 10.1002/ece3.4298

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Bolnick, D. I., and Kirkpatrick, M. (2013). Assortative mating in animals. Am. Nat. 181, 125–138. doi: 10.1086/670160

PubMed Abstract | CrossRef Full Text | Google Scholar

Johannesson, K., and Butlin, R. K. (2017). What explains rare and conspicuous colors in a snail? A test of time-series data against models of drift, migration or selection. Heredity 118, 21–30. doi: 10.1038/hdy.2016.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Johannesson, K., and Ekendahl, A. (2002). Selective predation favouring cryptic individuals of marine snails (Littorina). Biol. J. Linn. Soc. 76, 137–144. doi: 10.1111/j.1095-8312.2002.tb01720.x

CrossRef Full Text | Google Scholar

Johannesson, K., and Johannesson, B. (1995). Dispersal and population expansion in a direct developing marine snail (Littorina saxatilis) following a severe population bottleneck. Hydrobiologia 309, 173–180. doi: 10.1007/BF00014485

CrossRef Full Text | Google Scholar

Johannesson, K., and Mikhailova, N. (2004). Habitat-related genetic substructuring in a marine snail (Littorina fabalis) involving a tight link between an allozyme and a DNA locus. Biol. J. Linn. Soc. 81, 301–306. doi: 10.1111/j.1095-8312.2003.00288.x

CrossRef Full Text | Google Scholar

Johannesson, K., Saltin, S. H., Duranovic, I., Havenhand, J. N., and Jonsson, P. R. (2010). Indiscriminate males: Mating behavior of a marine snail compromised by a sexual conflict? PLoS One 5:e1205. doi: 10.1371/journal.pone.0012005

PubMed Abstract | CrossRef Full Text | Google Scholar

Keller, L. F., and Waller, D. M. (2002). Inbreeding effects in wild populations. Trends Ecol. Evol. 17, 230–241. doi: 10.1016/S0169-5347(02)02489-8

CrossRef Full Text | Google Scholar

Kemppainen, P., Lindskog, T., Butlin, R., and Johannesson, K. (2011). Intron sequences of arginine kinase in an intertidal snail suggest an ecotype-specific selective sweep and a gene duplication. Heredity 106, 808–816. doi: 10.1038/hdy.2010.123

PubMed Abstract | CrossRef Full Text | Google Scholar

Knoppien, P. (1985). Rare male mating advantage. Nature 60, 81–117. doi: 10.1038/272189a0

CrossRef Full Text | Google Scholar

Kokko, H., Brooks, R., Jennions, M. D., and Morley, J. (2003). The evolution of mate choice and mating biases. Proc. R. Soc. B 270, 653–664. doi: 10.1098/rspb.2002.2235

PubMed Abstract | CrossRef Full Text | Google Scholar

Kokko, H., Jennions, M. D., and Houde, A. (2007). Evolution of frequency-dependent mate choice: keeping up with fashion trends. Proc. R. Soc. B 274, 1317–1324. doi: 10.1098/rspb.2007.0043

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozminsky, E. V. (2014). Inheritance of the background shell color in the snails Littorina obtusata (Gastropoda, Littorinidae). Russ. J. Genet. 50, 1038–1047. doi: 10.1134/S1022795414100044

CrossRef Full Text | Google Scholar

Kvarnemo, C. (1996). Temperature affects operational sex ratio and intensity of male-male competition: an experimental study of sand gobies, Pomatoschistus minutus. Behav. Ecol. 7, 208–212. doi: 10.1093/beheco/7.2.208

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawniczak, M. K. N., Barnes, A. I., Linklater, J. R., Boone, J. M., Wigby, S., and Chapman, T. (2007). Mating and immunity in invertebrates. Trends Ecol. Evol. 22, 48–55. doi: 10.1016/j.tree.2006.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

López-Cortegano, E., Pérez-Figueroa, A., and Caballero, A. (2019). METAPOP2: Re-implementation of software for the analysis and management of subdivided populations using gene and allelic diversity. Mol. Ecol. Res. 19, 1095–1100. doi: 10.1111/1755-0998.13015

PubMed Abstract | CrossRef Full Text | Google Scholar

López-Cortegano, E., Carpena-Catoira, C., Carvajal-Rodríguez, A., and Rolán-Alvarez, E. (2020). Mate choice based on body size similarity in sexually dimorphic populations causes strong sexual selection. Anim. Behav. 160, 69–78. doi: 10.1016/j.anbehav.2019.12.005

CrossRef Full Text | Google Scholar

Martínez-Rodríguez, P., Rolán-Alvarez, E., del Mar Pérez-Ruiz, M., Arroyo-Yebras, F., Carpena-Catoira, C., Carvajal-Rodríguez, A., et al. (2019). Geographic and temporal variation of distinct intracellular endosymbiont strains of Wolbachia sp. in the grasshopper Chorthippus parallelus: a frequency-dependent mechanism? Microb. Ecol. 77, 1036–1047. doi: 10.1007/s00248-019-01338-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Mays, H. L., and Hill, G. E. (2004). Choosing mates: good genes versus genes that are a good fit. Trends Ecol. Evol. 19, 554–559. doi: 10.1016/j.tree.2004.07.018

PubMed Abstract | CrossRef Full Text | Google Scholar

McKillup, S. C., and McKillup, R. V. (2008). Apostasy and selection for crypsis in the marine snail Littoraria filosa: an explanation for a balanced colour polymorphism. Biol. J. Linn. Soc. 95, 62–71. doi: 10.1111/j.1095-8312.2008.01032.x

CrossRef Full Text | Google Scholar

McKinnon, J. S., and Pierotti, M. E. R. (2010). Color polymorphism and correlated characters: genetic mechanisms and evolution. Mol. Ecol. 19, 5101–5125. doi: 10.1111/j.1365-294X.2010.04846.x

PubMed Abstract | CrossRef Full Text | Google Scholar

McQuaid, C. (1996). Biology of the gastropod family Littorinidae: 1, Evolutionary aspects. Oceanogr. Mar. Biol. Annu. Rev. 34, 233–262.

Google Scholar

Miller, L. P., and Denny, M. W. (2011). Importance of behavior and morphological traits for controlling body temperature in littorinid snails. Biol. Bull. 220, 209–223. doi: 10.1086/BBLv220n3p209

PubMed Abstract | CrossRef Full Text | Google Scholar

Ng, T. P. T., and Johannesson, K. (2016). No precopulatory inbreeding avoidance in the intertidal snail Littorina saxatilis. J. Molluscan Stud. 82, 213–215. doi: 10.1093/mollus/eyv035

CrossRef Full Text | Google Scholar

Ng, T. P. T., Rolán-Alvarez, E., Dahlén, S. S., Davies, M. S., Estévez, D., Stafford, R., et al. (2019). The causal relationship between sexual selection and sexual size dimorphism in marine gastropods. Anim. Behav. 148, 53–62. doi: 10.1016/j.anbehav.2018.12.005

CrossRef Full Text | Google Scholar

Ng, T. P. T., Saltin, S. H., Davies, M. S., Johannesson, K., Stafford, R., and Williams, G. A. (2013). Snails and their trails: the multiple functions of trail-following in gastropods. Biol. Rev. Camb. Philos. Soc. 88, 683–700. doi: 10.1111/brv.12023

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Hara, R. B. (2005). Comparing the effects of genetic drift and fluctuating selection on genotype frequency changes in the scarlet tiger moth. Proc. R. Soc. B Biol. Sci. 272, 211–217. doi: 10.1098/rspb.2004.2929

PubMed Abstract | CrossRef Full Text | Google Scholar

Olendorf, R., Rodd, F. H., Punzalan, D., Houde, A. E., Hurt, C., Reznick, D. N., et al. (2006). Frequency-dependent survival in natural guppy populations. Nature 441, 633–636. doi: 10.1038/nature04646

PubMed Abstract | CrossRef Full Text | Google Scholar

Olsen, K. C., Ryan, W. H., Winn, A. A., Kosman, E. T., Moscoso, J. A., Krueger-Hadfield, S. A., et al. (2020). Inbreeding shapes the evolution of marine invertebrates. Evolution 74, 871–882. doi: 10.1111/evo.13951

PubMed Abstract | CrossRef Full Text | Google Scholar

Oxford, G. S. (2005). Genetic drift within a protected polymorphism: enigmatic variation in color-morph frequencies in the candy-stripe spider, Enoplognatha ovata. Evolution 59, 2170. doi: 10.1554/05-046.1

CrossRef Full Text | Google Scholar

Panova, M., Boström, J., Hofving, T., Areskoug, T., Eriksson, A., Mehlig, B., et al. (2010). Extreme female promiscuity in a non-social invertebrate species. PLoS One 5:e0009640. doi: 10.1371/journal.pone.0009640

PubMed Abstract | CrossRef Full Text | Google Scholar

Parsonage, S., and Hughes, J. (2002). Natural selection and the distribution of shell color morphs in three species of Littoraria (Gastropoda: Littorinidae) in Moreton Bay, Queensland. Biol. J. Linn. Soc. 75, 219–232. doi: 10.1046/j.1095-8312.2002.00015.x

CrossRef Full Text | Google Scholar

Paterson, I. G., Partridge, V., and Buckland-Nicks, J. (2001). Multiple paternity in Littorina obtusata (Gastropoda, Littorinidae) revealed by microsatellite analysis. Biol. Bull. 200, 261–267. doi: 10.2307/1543508

PubMed Abstract | CrossRef Full Text | Google Scholar

Peres-Neto, P. R., and Olden, J. D. (2001). Assessing the robustness of randomization tests: examples from behavioral studies. Anim. Behav. 61, 79–86. doi: 10.1006/anbe.2000.1576

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez-Figueroa, A., Caballero, A., and Rolán-Alvarez, E. (2005). Comparing the estimation properties of different statistics for measuring sexual isolation from mating frequencies. Biol. J. Linn. Soc. 85, 307–318. doi: 10.1111/j.1095-8312.2005.00491.x

CrossRef Full Text | Google Scholar

Phifer-Rixey, M., Heckman, M., Trussell, G. C., and Schmidt, P. S. (2008). Maintenance of clinal variation for shell color phenotype in the flat periwinkle Littorina obtusata. J. Evol. Biol. 21, 966–978. doi: 10.1111/j.1420-9101.2008.01549.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Piñeira, J., Quesada, H., Rolán-Alvarez, E., and Caballero, A. (2008). Genetic impact of the Prestige oil spill in wild populations of a poor dispersal marine snail from intertidal rocky shores. Mar. Pollut. Bull. 56, 270–281. doi: 10.1016/j.marpolbul.2007.10.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Pusey, A. (1987). Sex-biased dispersal and inbreeding avoidance in birds and mammals. Trends Ecol. Evol. 2, 295–299. doi: 10.1016/0169-5347(87)90081-4

CrossRef Full Text | Google Scholar

Reid, D. G. (1996). Systematics and Evolution of Littorina. London: The Ray Society.

Google Scholar

Reimchen, T. E. (1979). Substratum heterogeneity, crypsis, and color polymorphism in an intertidal snail (Littorina mariae). Can. J. Zool. 57, 1070–1085. doi: 10.1139/z79-135

CrossRef Full Text | Google Scholar

Reimchen, T. E. (1989). Shell color ontogeny and tubeworm mimicry in a marine gastropod Littorina mariae. Biol. J. Linn. Soc. 36, 97–109. doi: 10.1111/j.1095-8312.1989.tb00485.x

CrossRef Full Text | Google Scholar

Rios-Cardenas, O., Bono, L., and Morris, M. R. (2018). Frequency-dependent selection and fluctuations around an equilibrium for alternative reproductive tactics in a swordtail. Anim. Behav. 140, 19–28. doi: 10.1016/j.anbehav.2018.03.018

CrossRef Full Text | Google Scholar

Rolán-Alvarez, E., Austin, C. J., and Boulding, E. G. (2015a). The contribution of the genus Littorina to the field of evolutionary ecology. Oceanogr. Mar. Biol. Annu. Rev. 53, 157–214. doi: 10.1023/a:1023038817912

CrossRef Full Text | Google Scholar

Rolán-Alvarez, E., Carvajal-Rodríguez, A., de Coo, A., Cortés, B., Estévez, D., Ferreira, M., et al. (2015b). The scale-of-choice effect and how estimates of assortative mating in the wild can be biased due to heterogeneous samples. Evolution 69, 1845–1857. doi: 10.1111/evo.12691

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolán-Alvarez, E., and Caballero, A. (2000). Estimating sexual selection and sexual isolation effects from mating frequencies. Evolution 54, 30–36. doi: 10.1111/j.0014-3820.2000.tb00004.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolán-Alvarez, E., and Ekendahl, A. (1996). Sexual selection and non-random mating for shell color in a natural population of the marine snail Littorina mariae (Gastropoda: Prosobranchia). Genetica 97, 39–46. doi: 10.1007/BF00132579

CrossRef Full Text | Google Scholar

Rolán-Alvarez, E., Saura, M., Diz, A. P., Rivas, M. J., Alvarez, M., CortéS, B., et al. (2012). Can sexual selection and disassortative mating contribute to the maintenance of a shell color polymorphism in an intertidal marine snail? Curr. Zool. 58, 463–474. doi: 10.1093/czoolo/58.3.463

CrossRef Full Text | Google Scholar

Rolán-Alvarez, E., Zapata, C., and Alvarez, G. (1995). Distinct genetic subdivision in sympatric and sibling species of the genus Littorina (Gastropoda: Littorinidae). Heredity 74, 1–9. doi: 10.1038/hdy.1995.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Sacchi, C. F. (1969). Recherches sur I’tcologie comparée de: Littorina obtusata (L.) et de Littorina mariae Sacchi and Rast. (Gastropoda Prosobranchia) en Galice et en Bretagne. Sci. Mar. 33, 381–413.

Google Scholar

Sacchi, C. F. (1972). “Recherches sur la valence thermique du couple d’espèces intertidales Littorina obtusata (L.) et L. mariae [=fabalis] Sacchi & Rast. (Gasteropoda, Prosobranchia),” in Proceedings of the Fifth European Marine Biology Symposium, ed. B. Battaglia (Padova: Piccin), 348.

Google Scholar

Sacchi, C. F., and Rastelli, M. (1966). Littornia mariae nov. sp. les différences morphologiques et écologiques entre “nains” et “normaux” chez l’ “espèce” L. obtusata (L.) (Gastr. Prosobr.) et leur signification adaptative et évolutive. Atti della Società Italiana de Scienze naturale e del Museo Civico di Storia naturale di Milano 105, 351–369. doi: 10.15468/39omei

CrossRef Full Text | Google Scholar

Saltin, S. H., Schade, H., and Johannesson, K. (2013). Preference of males for large females causes a partial mating barrier between a large and a small ecotype of Littorina fabalis (W. Turton, 1825). J. Molluscan Stud. 79, 128–132. doi: 10.1093/mollus/eyt003

CrossRef Full Text | Google Scholar

Scheil, A. E., Hilsmann, S., Triebskorn, R., and Köhler, H. R. (2013). Shell color polymorphism, injuries and immune defense in three helicid snail species, Cepaea hortensis, Theba pisana and Cornu aspersummaximum. Results Immunol. 3, 73–78. doi: 10.1016/j.rinim.2013.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Sergievsky, S. O. (1992). “A review of ecophysiological studies of the color polymorphism of Littorina obtusata (L.) and L. saxatilis (Olivi) in the White Sea,” in Proceedings of the Third International Symposium on Littorinid Biology, eds J. Grahame, P. J. Mill, and D. G. Reid (London: The Malacological Society of London), 235–245.

Google Scholar

Sinervo, B., and Lively, C. M. (1996). The rock–paper–scissors game and the evolution of alternative male strategies. Nature 380, 240–243. doi: 10.1038/380240a0

CrossRef Full Text | Google Scholar

Sokal, R. R., and Rohlf, F. J. (1995). Biometry. New York, NY: W. H. Freeman.

Google Scholar

Svensson, E. I., Abbott, J. K., and Härdling, R. (2005). Female polymorphism, frequency dependence, and rapid evolutionary dynamics in natural populations. Am. Nat. 165, 567–576. doi: 10.1086/429278

PubMed Abstract | CrossRef Full Text | Google Scholar

Svensson, E. I., and Connallon, T. (2019). How frequency-dependent selection affects population fitness, maladaptation and evolutionary rescue. Evol. Appl. 12, 1243–1258. doi: 10.1111/eva.12714

PubMed Abstract | CrossRef Full Text | Google Scholar

Swierk, L., and Langkilde, T. (2019). Fitness costs of mating with preferred females in a scramble mating system. Behav. Ecol. 30, 658–665. doi: 10.1093/beheco/arz001

PubMed Abstract | CrossRef Full Text | Google Scholar

Taborsky, B., Guyer, L., and Taborsky, M. (2009). Size-assortative mating in the absence of mate choice. Anim. Behav. 77, 439–448. doi: 10.1016/j.anbehav.2008.10.020

CrossRef Full Text | Google Scholar

Takahashi, T., and Hori, M. (2008). Evidence of disassortative mating in a Tanganyikan cichlid fish and its role in the maintenance of intrapopulation dimorphism. Biol. Lett. 4, 497–499. doi: 10.1098/rsbl.2008.0244

PubMed Abstract | CrossRef Full Text | Google Scholar

Tregenza, T., and Wedell, N. (2000). Genetic compatibility, mate choice and patterns of parentage: invited review. Mol. Ecol. 9, 1013–1027. doi: 10.1046/j.1365-294X.2000.00964.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Trotter, M. V., and Spencer, H. G. (2008). The generation and maintenance of genetic variation by frequency-dependent selection: constructing polymorphisms under the pairwise interaction model. Genetics 180, 1547–1557. doi: 10.1534/genetics.108.088880

PubMed Abstract | CrossRef Full Text | Google Scholar

Tucker, G. M. (1991). Apostatic selection by song thrushes (Turdus philomelos) feeding on the snail Cepaea hortensis. Biol. J. Linn. Soc. 43, 149–156. doi: 10.1111/j.1095-8312.1991.tb00590.x

CrossRef Full Text | Google Scholar

Tuttle, E. M., Bergland, A. O., Korody, M. L., Brewer, M. S., Newhouse, D. J., Minx, P., et al. (2016). Divergence and functional degradation of a sex chromosome-like supergene. Curr. Biol. 26, 344–350. doi: 10.1016/j.cub.2015.11.069

PubMed Abstract | CrossRef Full Text | Google Scholar

Valvo, J. J., Rodd, F. H., and Hughes, K. A. (2019). Consistent female preference for rare and unfamiliar male color patterns in wild guppy populations. Behav. Ecol. 30, 1672–1681. doi: 10.1093/beheco/arz134

PubMed Abstract | CrossRef Full Text | Google Scholar

van Nes, S. (1998). Why Different Morphs of Littorina fabalis Prefer Different Microhabitats. Master’s thesis, University of Gothenburg, Gothenburg.

Google Scholar

Walsh, B., and Lynch, M. (2018). Evolution and Selection of Quantitative Traits. New York, NY: Oxford University Press.

Google Scholar

Waples, R. S., and Do, C. (2008). LDNE: a program for estimating effective population size from data on linkage disequilibrium. Mol. Ecol. Resour. 8, 753–756. doi: 10.1111/j.1755-0998.2007.02061.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Weir, B. S. (1979). Inferences about linkage disequilibrium. Biometrics 35, 235–254. doi: 10.2307/2529947

CrossRef Full Text | Google Scholar

Weir, J. C. (2018). The evolution of color polymorphism in British winter-active Lepidoptera in response to search image use by avian predators. J. Evol. Biol. 31, 1109–1126. doi: 10.1111/jeb.13290

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiss, S. L. (2006). Female-specific color is a signal of quality in the striped plateau lizard (Sceloporus virgatus). Behav. Ecol. 17, 726–732. doi: 10.1093/beheco/arl001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wellenreuther, M., Svensson, E. I., and Hansson, B. (2014). Sexual selection and genetic color polymorphisms in animals. Mol. Ecol. 23, 5398–5414. doi: 10.1111/mec.12935

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, G. A. (1990). The comparative ecology of the flat periwinkles, Littorina obtusata (L.) and L. mariae Sacchi et Rastelli. Field Stud. 7, 469–482.

Google Scholar

Williams, G. A. (1992). The effect of predation on the life histories of Littorina obtusata and Littorina mariae. J. Mar. Biol. Assoc. U.K. 72, 403–416. doi: 10.1017/S0025315400037784

CrossRef Full Text | Google Scholar

Williams, G. A. (1996). “Seasonal variation in a low shore Fucus serratus (Fucales, Phaeophyta) population and its epiphytic fauna,” in Proceedings of the Fifteenth International Seaweed Symposium, eds S. C. Lindstrom D. J. Chapman (Dordrecht: Springer), 191–197. doi: 10.1007/978-94-009-1659-3_26

CrossRef Full Text | Google Scholar

Wright, S. (1922). Coefficients of inbreeding and relationship. Am. Nat. 56, 330–338. doi: 10.1086/279872

CrossRef Full Text | Google Scholar

Wright, S. (1939). The distribution of self-sterility alleles in populations. Genetics 24, 538–552. doi: 10.2307/2406213

CrossRef Full Text | Google Scholar