Frequency-dependent Selection Research Articles

Evolution repeats itself in replicate long-term studies in the wild.

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

Evolution and genetic architecture of sex-limited polymorphism in cuckoos.

Sex-limited polymorphism has evolved in many species including our own. Yet, we lack a detailed understanding of the underlying genetic variation and evolutionary processes at work. The brood parasitic common cuckoo (Cuculus canorus) is a prime example of female-limited color polymorphism, where adult males are monochromatic gray and females exhibit either gray or rufous plumage. This polymorphism has been hypothesized to be governed by negative frequency-dependent selection whereby the rarer female morph is protected against harassment by males or from mobbing by parasitized host species. Here, we show that female plumage dichromatism maps to the female-restricted genome. We further demonstrate that, consistent with balancing selection, ancestry of the rufous phenotype is shared with the likewise female dichromatic sister species, the oriental cuckoo (Cuculus optatus). This study shows that sex-specific polymorphism in trait variation can be resolved by genetic variation residing on a sex-limited chromosome and be maintained across species boundaries.

The narrow footprint of ancient balancing selection revealed by heterokaryon incompatibility genes in Aspergillus fumigatus.

In fungi, fusion between individuals leads to localized cell death, a phenomenon termed heterokaryon incompatibility. Generally, the genes responsible for this incompatibility are observed to be under balancing selection resulting from negative frequency-dependent selection. Here, we assess this phenomenon in Aspergillus fumigatus, a human pathogenic fungus with a very low level of linkage disequilibrium as well as an extremely high crossover rate. Using complementation of auxotrophic mutations as an assay for hyphal compatibility, we screened sexual progeny for compatibility to identify genes involved in this process, called het genes. In total, 5/148 (3.4%) offspring were compatible with a parent and 166/2142 (7.7%) sibling pairs were compatible, consistent with several segregating incompatibility loci. Genetic mapping identified five loci, four of which could be fine mapped to individual genes, of which we tested three through heterologous expression, confirming their causal relationship. Consistent with long-term balancing selection, trans-species polymorphisms were apparent across several sister species, as well as equal allele frequencies within A. fumigatus. Surprisingly, a sliding window genome-wide population-level analysis of an independent dataset did not show increased Tajima's D near these loci, in contrast to what is often found surrounding loci under balancing selection. Using available de novo assemblies, we show that these balanced polymorphisms are restricted to several hundred base pairs flanking the coding sequence. In addition to identifying the first het genes in an Aspergillus species, this work highlights the interaction of long-term balancing selection with rapid linkage disequilibrium decay.

Absence of long-term balancing selection on variation in EuMYB3, an R2R3-MYB gene responsible for the anther-color polymorphism in Erythronium umbilicatum

Balancing selection has been shown to be common in plants for several different types of traits, such as self-incompatibility and heterostyly. Generally, for these traits balancing selection is generated by interactions among individuals or between individuals and other species (e.g., pathogens or pollinators). However, there are phenotypic polymorphisms in plants that do not obviously involve types of interactions that generate balancing selection. Little is known about the extent to which balancing selection also acts to preserve these polymorphisms. Here we ask whether balancing selection preserves an anther-color polymorphism in Erythronium umbilicatum (Liliaceae). We identified a major gene underlying this polymorphism. We then attempted to detect signatures of balancing selection on that gene by developing a new coalescence test for balancing selection. We found that variation in anther color is in large part caused by variation in a paralog of EuMYB3, an anthocyanin-regulating R2R3-MYB transcription factor. However, we found little evidence for balancing selection having acted historically on EuMYB3. Our results thus suggest that plant polymorphisms, especially those not involved in interactions that are likely to generate negative frequency-dependent selection, may reflect a transient state in which one morph will eventually be fixed by either genetic drift or directional selection. Our results also suggest that regulation of the anthocyanin pathway is more evolutionarily labile than is generally believed.

Abstract A027: Detecting and modeling frequency-dependent evolution in serially passaged primary triple negative breast cancer patient derived xenografts

Abstract Tumor evolution is a contributing factor to most failed cancer treatments. The intratumor heterogeneity of cellular phenotypes found in the tumor microenvironment acts as an evolutionary substate for cancer cells during tumorigenesis, subsequent growth, and after the patient is in the clinic. Selective forces imposed by certain medications can cause the emergence of a tumor with a resistant phenotype. Chemotherapies that target rapidly dividing cells will select for slower dividing ones which are less sensitive to such agents. However, many selective pressures acting on the tumor microenvironment are yet to be elucidated. For example, the effect one cancer cell has on another’s growth through competitive (i.e., the struggle for resources) or cooperative (e.g., exchange of growth factors) interactions remains poorly understood. Our central hypothesis is that cancer cells with distinct karyotypes will have a varying growth rate influenced by the karyotype landscape of the entire tumor. That is, we expect the fitness of a cancer cell to be largely determined by its karyotype and the karyotypes of its neighbors. Thus, we predict that frequency-dependent selection has a significant impact on clonal evolution in cancer. To test this, we have developed two novel methods for respectively identifying likely instances of frequency-dependent selection and modeling the long-term evolutionary trajectory in such cases. To detect frequency-dependent effects, we investigate the relationship between the frequency of one clone and the growth rate of all others (including itself). If a statistically significant correlation can be found, we proceed to modeling clonal evolution. To do this, we use a penalized least squares method and find the best-fitting parameterization of the replicator equation to the clonally resolved frequency data. This provides us with a payoff matrix describing the intensity of competition between each pair of clones. Furthermore, it allows us to model the long-term evolutionary dynamics of the tumor and to make predictions regarding co-existence of clones or clonal sweeps. Recently, investigators serially passaged primary triple negative breast cancer (TNBC) patient derived xenograft (PDX) tumors in NRG mice for over a year, inferring clonal frequencies at multiple timepoints using single cell whole genome (scWGS) and RNA sequencing (scRNA-seq). From this data, we were able to detect statistically significant frequency-dependent effects in the majority of serially passaged PDX tumors. Of these, most are predicted to result in a clonal sweep. Three of the PDX tumors were treated with cisplatin which revealed clone-specific drug response. This presents an opportunity for adaptive therapy scheduling which takes advantage of fitness tradeoffs under different contexts. A greater understanding of the effects frequency-dependent selection has on tumor evolution will contribute to the design of therapies intended to preclude or delay the emergence of resistance to treatment. Citation Format: Thomas A. Veith, Richard Beck, Noemi Andor. Detecting and modeling frequency-dependent evolution in serially passaged primary triple negative breast cancer patient derived xenografts [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Translating Cancer Evolution and Data Science: The Next Frontier; 2023 Dec 3-6; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_2):Abstract nr A027.

Floral morph variation mediated by clonal growth and pollinator functional groups of Limonium otolepis in a heterostylous fragmented population.

Abstract. Heterostyly, a genetic style polymorphism, is linked to symmetric pollen transfer, vital for its maintenance. Clonal growth typically impacts sexual reproduction by influencing pollen transfer. However, the floral morph variation remains poorly understood under the combined effects of pollinators and clonal growth in heterostyly characterized by negative frequency-dependent selection and disassortative mating. We estimated morph ratios, ramets per genet and heterostylous syndrome and quantified legitimate pollen transfer via clonal growth, pollinators and reciprocal herkogamy between floral morphs in Limonium otolepis, a fragmented population composed of five subpopulations in the desert environment of northwestern China, with small flower and large floral morph variation. All subpopulations but one exhibited pollen-stigma morphology dimorphism. The compatibility between mating types with different pollen-stigma morphologies remained consistent regardless of reciprocal herkogamy. Biased ratios and ramets per genet of the two mating types with distinct pollen-stigma morphologies caused asymmetric pollen flow and varying fruit sets in all subpopulations. Short-tongued insects were the primary pollinators due to small flower sizes. However, pollen-feeding Syrphidae sp. triggered asymmetry in pollen flow between high and low sex organs, with short-styled morphs having lower stigma pollen depositions and greater variation. Clonal growth amplified this variation by reducing intermorph pollen transfer. All in all,pollinators and clonal growth jointly drive floral morph variation. H-morphs with the same stigma-anther position and self-incompatibility, which mitigate the disadvantages of sunken low sex organs with differing from the classical homostyly, might arise from long- and short-styled morphs through a 'relaxed selection'. This study is the first to uncover the occurrence of the H-morph and its associated influencing factors in a distylous plant featuring clonal growth, small flowers and a fragmented population.

Rediversification following ecotype isolation reveals hidden adaptive potential

Microbial communities play a critical role in ecological processes, and their diversity is key to their functioning. However, little is known about whether communities can regenerate ecological diversity following ecotype removal or extinction and how the rediversified communities would compare to the original ones. Here, we show that simple two-ecotype communities from the E.coli long-term evolution experiment (LTEE) consistently rediversified into two ecotypes following the isolation of one of the ecotypes, coexisting via negative frequency-dependent selection. Communities separated by more than 30,000 generations of evolutionary time rediversify in similar ways. The rediversified ecotype appears to share a number of growth traits with the ecotype it replaces. However, the rediversified community is also different from the original community in ways relevant to the mechanism of ecotype coexistence-for example, in stationary phase response and survival. We found substantial variation in the transcriptional states between the two original ecotypes, whereas the differences within the rediversified community were comparatively smaller, although the rediversified community showed unique patterns of differential expression. Our results suggest that evolution may leave room for alternative diversification processes even in a maximally reduced community of only two strains. We hypothesize that the presence of alternative evolutionary pathways may be even more pronounced in communities of many species where there are even more potential niches, highlighting an important role for perturbations, such as species removal, in evolving ecological communities.

Why the Shaw-Mohler equation works and when it doesn't.

Fitness gain curves were introduced into the framework of the Shaw-Mohler equation, the foundation of sex allocation theory. I return to the Shaw-Mohler equation to consider how it embodies the rare-sex advantage underlying frequency-dependent selection on the sex ratio. The Shaw-Mohler formulation is based on the numbers of males and females randomly mating in a panmictic population. Gain curves are meant to describe reproductive success through male and female functions in hermaphrodites and have been inserted in place of male and female numbers in the Shaw-Mohler equation. In doing so, gain curves bypass consideration of the implicit mating process in the Shaw-Mohler argument and can lead to anomalies like unequal total male and female fitness in a population. If gain curves truly represent fitness gain, equality of total male and female fitness requires a constant sex allocation of equal resource investment into male and female functions. The blurring of input with fitness outcome has led to misinterpretation of what gain curves mean in reproductive ecology. They can describe a particular reproductive ecology, such as diminishing fitness returns on resource investment, but lack causal efficacy with respect to sex allocation.

Bernoulli factories and duality in Wright–Fisher and Allen–Cahn models of population genetics

Mathematical models of genetic evolution often come in pairs, connected by a so-called duality relation. The most seminal example are the Wright–Fisher diffusion and the Kingman coalescent, where the former describes the stochastic evolution of neutral allele frequencies in a large population forwards in time, and the latter describes the genetic ancestry of randomly sampled individuals from the population backwards in time. As well as providing a richer description than either model in isolation, duality often yields equations satisfied by quantities of interest. We employ the so-called Bernoulli factory – a celebrated tool in simulation-based computing – to derive duality relations for broad classes of genetics models. As concrete examples, we present Wright–Fisher diffusions with general drift functions, and Allen–Cahn equations with general, nonlinear forcing terms. The drift and forcing functions can be interpreted as the action of frequency-dependent selection. To our knowledge, this work is the first time a connection has been drawn between Bernoulli factories and duality in models of population genetics.

Embracing the diversity in diverse warning signals

Positive frequency-dependent selection should theoretically lead to monomorphic warning coloration. Instead, numerous examples of polymorphic warning signals exist. Biases – for example, in human perception – hinder our appreciation and research of understanding warning signal diversity. We propose strategies to counter such biases and objectively move our field forward.

Genetic basis of Arabidopsis thaliana responses to infection by naïve and adapted isolates of turnip mosaic virus.

Plant viruses account for enormous agricultural losses worldwide, and the most effective way to combat them is to identify genetic material conferring plant resistance to these pathogens. Aiming to identify genetic associations with responses to infection, we screened a large panel of Arabidopsis thaliana natural inbred lines for four disease-related traits caused by infection by A. thaliana-naïve and -adapted isolates of the natural pathogen turnip mosaic virus (TuMV). We detected a strong, replicable association in a 1.5 Mb region on chromosome 2 with a 10-fold increase in relative risk of systemic necrosis. The region contains several plausible causal genes as well as abundant structural variation, including an insertion of a Copia transposon into a Toll/interleukin receptor (TIR-NBS-LRR) coding for a gene involved in defense, that could be either a driver or a consequence of the disease-resistance locus. When inoculated with TuMV, loss-of-function mutant plants of this gene exhibited different symptoms than wild-type plants. The direction and severity of symptom differences depended on the adaptation history of the virus. This increase in symptom severity was specific for infections with the adapted isolate. Necrosis-associated alleles are found worldwide, and their distribution is consistent with a trade-off between resistance during viral outbreaks and a cost of resistance otherwise, leading to negative frequency-dependent selection.

Genetic basis of Arabidopsis thaliana responses to infection by naïve and adapted isolates of turnip mosaic virus

Plant viruses account for enormous agricultural losses worldwide, and the most effective way to combat them is to identify genetic material conferring plant resistance to these pathogens. Aiming to identify genetic associations with responses to infection, we screened a large panel of Arabidopsis thaliana natural inbred lines for four disease-related traits caused by infection by A. thaliana-naïve and -adapted isolates of the natural pathogen turnip mosaic virus (TuMV). We detected a strong, replicable association in a 1.5 Mb region on chromosome 2 with a 10-fold increase in relative risk of systemic necrosis. The region contains several plausible causal genes as well as abundant structural variation, including an insertion of a Copia transposon into a Toll/interleukin receptor (TIR-NBS-LRR) coding for a gene involved in defense, that could be either a driver or a consequence of the disease-resistance locus. When inoculated with TuMV, loss-of-function mutant plants of this gene exhibited different symptoms than wild-type plants. The direction and severity of symptom differences depended on the adaptation history of the virus. This increase in symptom severity was specific for infections with the adapted isolate. Necrosis-associated alleles are found worldwide, and their distribution is consistent with a trade-off between resistance during viral outbreaks and a cost of resistance otherwise, leading to negative frequency-dependent selection.

Habitat heterogeneity limits prey colour polymorphism maintained via negative frequency-dependent selection.

The persistence of non-neutral trait polymorphism is enigmatic because stabilizing selection is expected to deplete variation. In cryptically coloured prey, negative frequency-dependent selection due to search image formation by predators has been proposed to favour rare variants, promoting polymorphism. However, in a heterogeneous environment, locally varying disruptive selection favours patch type-specific optima, resulting in spatial segregation of colour variants. Here, we address whether negative frequency-dependent selection can overcome selection posed by habitat heterogeneity to promote local polymorphism using an individual-based model. In addition, we compare how prey and predator mobility may modify the outcome. Our model revealed that frequency-dependent predation could strongly promote local prey polymorphism, but only when differences between morphs in patch-specific fitness were small. The effect of frequency-dependent predation depended on the predator adjustment of search image and was hampered by the prey population structure. Gene flow due to prey movement counteracted local selection, promoted local polymorphism to some extent, and relaxed the conditions for polymorphism due to frequency-dependent predation. Importantly, abrupt spatial changes in morph frequencies decreased the probability that mobile frequency-dependent predators could maintain local prey polymorphism. Overall, our study suggests that in a spatially heterogeneous environment, negative frequency-dependent selection may help maintain local polymorphism but only under a limited range of conditions.

How should we distinguish between selectable and circumstantial traits?

There is surprisingly little philosophical work on conceptually spelling out the difference between the traits on which natural selection may be said to act (e.g. "having a high running speed") and mere circumstantial traits (e.g. "happening to be in the path of a forest fire"). I label this issue the "selectable traits problem" and, in this paper, I propose a solution for it. I first show that, contrary to our first intuition, simply equating selectable traits with heritable ones is not an adequate solution. I then go on to argue that two recent philosophical solutions to this problem-due to Peter Godfrey-Smith and Pierrick Bourrat-are unconvincing because they cannot accommodate frequency-dependent selection. The way out of this difficulty is, I argue, to accept that extrinsic properties dependent on relations between intrinsic properties of the population members should also count as selectable traits. I then show that my proposal is legitimized by more than the simple accommodation of frequency-dependent selection.

Evidence for Negative Frequency Dependent Sexual Selection on Eye Color in Europeans

Evidence for Negative Frequency Dependent Sexual Selection on Eye Color in Europeans

The effect of divergent and parallel selection on the genomic landscape of divergence.

While the role of selection in divergence along the speciation continuum is theoretically well understood, defining specific signatures of selection in the genomic landscape of divergence is empirically challenging. Modelling approaches can provide insight into the potential role of selection on the emergence of a heterogenous genomic landscape of divergence. Here, we extend and apply an individual-based approach that simulates the phenotypic and genotypic distributions of two populations under a variety of selection regimes, genotype-phenotype maps, modes of migration, and genotype-environment interactions. We show that genomic islands of high differentiation and genomic valleys of similarity may respectively form under divergent and parallel selection between populations. For both types of between-population selection, negative and positive frequency-dependent selection within populations generated genomic islands of higher magnitude and genomic valleys of similarity, respectively. Divergence rates decreased under strong dominance with divergent selection, as well as in models including genotype-environment interactions under parallel selection. For both divergent and parallel selection models, divergence rate was higher under an intermittent migration regime between populations, in contrast to a constant level of migration across generations, despite an equal number of total migrants. We highlight that interpreting a particular evolutionary history from an observed genomic pattern must be done cautiously, as similar patterns may be obtained from different combinations of evolutionary processes. Modelling approaches such as ours provide an opportunity to narrow the potential routes that generate the genomic patterns of specific evolutionary histories.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

Experimental evidence of mate choice as the driving mechanism behind negative assortative mating for shell colour in a marine snail

Abstract How colour polymorphisms are maintained in natural populations constitutes a key subject of study for evolutionary ecologists. One of the mechanisms that has been put forward to explain this phenomenon is negative frequency-dependent selection, which can be caused by different mechanisms. Among these mechanisms is mate choice, or the selection as a potential mate of the rare or dissimilar colour morph by the choosy sex of a given species. In the context of colour polymorphism, mate choice predicts the occurrence of a negative assortative mating pattern, whereby different colour morphs mate more often than expected by chance alone. However, negative assortative mating can also result from other mechanisms, making it, at best, a fallible indicator of a mate choice-driven negative frequency-dependent selection. For this reason, experimental evidence explicitly connecting such a mating pattern to mate choice is necessary to claim that a particular colour polymorphism is being maintained by negative frequency-dependent selection. Here we present experimental evidence backing this relationship for the colour polymorphic marine gastropod Littorina saxatilis from the Ría de Vigo (NW Iberian Peninsula), although how specifically these organisms choose their mates remains unclear. This calls for further experimental efforts to clarify this issue and its relationship to the maintenance of colour polymorphism in L. saxatilis.

Non-additive effects between genotypes: Implications for competitive fitness assays.

Competitive fitness assays are widely used in evolutionary biology and typically rely on a reference strain to compare different focal genotypes. This approach implicitly relies on the absence of interaction between the competing genotypes. In other words, the performance of the reference strain must not depend on the competitor. This report scrutinized this assumption by competing diverged Drosophila simulans populations against a common reference strain. We detected strong evidence for interaction between the competing genotypes: (1) Frequency-dependent selection was common with opposite effects in genetically diverged populations. (2) Temporal heterogeneity of fitness estimates, which can be partially attributed to a competitor-specific delay in the eclosion of the reference strain. We propose that this inconsistent behavior of the reference strain can be considered a specific case of a genotype × environment interaction. Focal populations could modify the environment of the reference strain, either indirectly by altering the microbiome composition and food availability or directly by genotype-specific cannibalism. Our results provide new insights into the interaction of diverged genotypes and have important implications for the interpretation of competitive fitness assays.

ANTI-APOSTATIC FEEDING BEHAVIOR IN TWO SPECIES OF ASSASSIN BUG, THE WHITE-SPOTTED ASSASSIN BUG PLATYMERIS BIGUTTATUS AND RED-SPOTTED ASSASSIN BUG PLATYMERIS RHADAMANTHUS (HEMIPTERA: REDUVIIDAE)

Abstract Predator-prey relationships are a driving force when examining community ecology. One question in predator-prey interactions is how a predator may respond to novel prey phenotypes, a form of frequency dependent prey selection. A predator may selectively target or avoid the novel prey, a behavior referred to as anti-apostatic or apostatic selection respectively. We examined apostatic prey selection by observing feeding behavior in two species of assassin bugs (Hemiptera: Reduviidae), the white-spotted assassin bug (Platymeris biguttatus) and the red-spotted assassin bug (Platymeris rhadamanthus) which were fed unmodified and modified (novel) phenotypes of Turkistan cockroaches, Blatta lateralis (Blattodea: Blattidae). The two species of assassin bugs represent either generalist (P. biguttatus) or specific feeders (P. rhadamanthus). We hypothesized that the generalist feeder would engage in anti-apostatic feeding while the specific feeder would engage in apostatic feeding. Our results indicated that there was no difference in feeding behavior between the two species, however an overall trend of apostatic feeding was observed in response to two of the four novel prey phenotypes. Factors beyond prey phenotype may also be influencing feeding choice in these insects. While the two predator species did not exhibit different feeding behaviors, the observation of apostatic feeding in this genus suggest an overall trend of utilizing specific prey images in feeding.