Frequency-dependent Selection Research Articles

An updated phylogeography and population dynamics of porcine circovirus 2 genotypes: are they reaching an equilibrium?

Porcine circovirus 2 (PCV2) emerged more than three decades ago as one of the most impactful pathogens in the swine industry. Despite being a DNA virus, one of the hallmarks of PCV2 is its high evolutionary rate, which has led to the emergence of different genotypes, each exhibiting varying degrees of evolutionary success. Current knowledge suggests the occurrence of three main waves of genotype dominance, alternating over time (i.e., PCV2a, PCV2b, and PCV2d), alongside less prevalent genotypes. However, although PCV2d is currently the most common genotype nowadays, the others continue being circulating in the pig population. The present study reconsidered the epidemiological and evolutionary patterns of PCV2 genotypes using phylodynamic analyses, benefiting from an almost 10-fold increase in ORF2 sequence availability compared to previous studies. Additionally, a phylogeographic analysis was performed to investigate viral dispersal patterns and frequency, and the selective pressures acting on the capsid protein were estimated and compared among genotypes. While successive emergence of major genotypes was confirmed, this study extends previous findings by revealing subsequent prevalence fluctuations of PCV2a and PCV2b after their initial decline. This evolutionary process may represent an example of balancing selection, specifically negative frequency-dependent selection, where a genotype fitness diminishes as it becomes more common, and vice versa. Variations in genotype- or clade-specific immunity-affected by the local prevalence of viral groups-combined with the periodic introduction of strains that have independently evolved in different regions, may have led to fluctuations in the population dynamics of major genotypes over time. These fluctuations were associated with ongoing evolution and variations in the capsid amino acid profile. These findings have profound implications for future control strategies. Although PCV2d remains the most prevalent and widespread genotype, other genotypes should not be neglected. Control strategies should thus target the entire PCV2 population, with a focus on fostering broader and more cross-protective immunity.

Prediction of post-PCV13 pneumococcal evolution using invasive disease data enhanced by inverse-invasiveness weighting.

After introducing pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in Streptococcus pneumoniae. Predicting which pneumococcal strains will become common in carriage after vaccination can enhance vaccine design, public health interventions, and understanding of pneumococcal evolution. Invasive pneumococcal isolates were collected during 1998-2018 by the Active Bacterial Core surveillance (ABCs). Carriage data from Massachusetts (MA) and Southwest United States were used to calculate weights. Using pre-vaccine data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in carriage to the proportion in invasive data. Genomic data were processed under bioinformatic pipelines to define genetically similar sequence clusters (i.e., strains), and accessory genes (COGs) present in 5-95% of isolates. Weights were applied to adjust observed strain proportions and COG frequencies. The negative frequency-dependent selection (NFDS) model predicted strain proportions by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine COG frequencies. Inverse-invasiveness weighting increased the correlation of COG frequencies between invasive and carriage data in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13; and between different epochs in the invasive data. Weighting the invasive data significantly improved the NFDS model's accuracy in predicting strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R2 increasing from 0.254 before weighting to 0.545 after weighting. The weighting system adjusted invasive disease data to better represent the pneumococcal carriage population, allowing the NFDS mechanism to predict strain proportions in carriage in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance.IMPORTANCEStreptococcus pneumoniae, a common colonizer in the human nasopharynx, can cause invasive diseases including pneumonia, bacteremia, and meningitis mostly in children under 5 years or older adults. The PCV7 was introduced in 2000 in the United States within the pediatric population to prevent disease and reduce deaths, followed by PCV13 in 2010, PCV15 in 2022, and PCV20 in 2023. After the removal of vaccine serotypes, the prevalence of carriage remained stable as the vacated pediatric ecological niche was filled with certain non-vaccine serotypes. Predicting which pneumococcal clones, and which serotypes, will be most successful in colonization after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. While carriage data, which are collected from the pneumococcal population that is competing to colonize and transmit, are most directly relevant to evolutionary studies, invasive disease data are often more plentiful. Previously, evolutionary models based on negative frequency-dependent selection (NFDS) on the accessory genome were shown to predict which non-vaccine strains and serotypes were most successful in colonization following the introduction of PCV7. Here, we show that an inverse-invasiveness weighting system applied to invasive disease surveillance data allows the NFDS model to predict strain proportions in the projected carriage population in the post-PCV13/pre-PCV15 and pre-PCV20 epoch. The significance of our research lies in using a sample of invasive disease surveillance data to extend the use of NFDS as an evolutionary mechanism to predict post-PCV13 population dynamics. This has shown that we can correct for biased sampling that arises from differences in virulence and can enrich the value of genomic data from disease surveillance and advance our understanding of how NFDS impacts carriage population dynamics after both PCV7 and PCV13 vaccination.

The evolution and maintenance of trioecy with cytoplasmic male sterility.

Trioecy, the co-existence of females, males and hermaphrodites, is a rare sexual system in plants that may be an intermediate state in transitions between hermaphroditism and dioecy. Previous models have identified pollen limitation as a necessary condition for the evolution of trioecy from hermaphroditism. In these models, the seed-production and pollen production of females and males relative to those of hermaphrodites, respectively, are compromised by self-fertilization by hermaphrodites under pollen- limitation. Here, we investigate the evolution of trioecy via the invasion of cytoplasmic male sterility (CMS) into androdioecious populations in which hermaphrodites co-occur with males and where the male determiner is linked to a (partial) fertility restorer. We show that the presence of males in a population renders invasion by CMS more difficult. However, the presence of males also facilitates the maintenance of trioecy even in the absence of pollen limitation by negative frequency-dependent selection, because males reduce the transmission of CMS by females by siring sons (which cannot transmit CMS). We discuss our results in light of empirical observations of trioecy in plants and its potential role in the evolution of dioecy.

The extraordinary Atlantic silverside and me

Abstract My career is inextricably linked to the remarkable characteristics and astonishing insights provided by a simple yet elegant species, the Atlantic silverside, Menidia menidia. My initial field and laboratory studies provided the first demonstration of temperature-dependent sex determination (TSD) and its adaptive significance in fish, as well as the first empirical confirmation of frequency-dependent sex ratio selection. This early work then led to the discovery of countergradient variation in growth, evolutionary trade-offs associated with growth, pioneered the study of fisheries-induced evolution, and furthered our understanding of the spatial scale of local adaptation in marine species. Most of these topics are now vigorous areas of research encompassing numerous species and systems. Here I describe how these unforeseen discoveries were made and led to other findings on bluefish (Pomatomus saltatrix) and striped bass (Morone saxatilis). An interest in leadership during my early years was later rekindled as I took on administrative roles. It all began with an intriguing observation stemming from my MS thesis that I was compelled to follow. My career illustrates the value of exploring deeply the natural history of lesser-known yet easy-to-study species and to follow unexpected observations even if peripheral to a preconceived plan of research.

The role of mate competition in speciation and divergence: a systematic review.

Competition for mates can play a critical role in determining reproductive success, shaping phenotypic variation within populations, and influencing divergence. Yet, studies of the role of sexual selection in divergence and speciation have focused disproportionately on mate choice. Here, we synthesize the literature on how mate competition may contribute to speciation and integrate concepts from work on sexual selection within populations - mating systems, ecology, and mate choice. Using this synthesis, we generate testable predictions for how mate competition may contribute to divergence. Then, we identify the extent of existing support for these predictions in the literature with a systematic review of the consequences of mate competition for population divergence across a range of evolutionary, ecological, and geographic contexts. We broadly evaluate current evidence, identify gaps in available data and hypotheses that need testing, and outline promising directions for future work. A major finding is that mate competition may commonly facilitate further divergence after initial divergence has occurred, e.g., upon secondary contact and between allopatric populations. Importantly, current hypotheses for how mate competition contributes to divergence do not fully explain observed patterns. While results from many studies fit predictions of negative frequency dependent selection, agonistic character displacement, and ecological selection, results from ~30% studies did not fit existing conceptual models. This review identifies future research aims for scenarios in which mate competition is likely important but has been understudied, including how ecological context and interactions between mate choice and mate competition can facilitate or hinder divergence and speciation.

Female preference counteracts negative frequency dependent selection for a stable polymorphism in a livebearing fish

The maintenance of variation within natural populations is key for natural selection to operate. Polymorphism is an intriguing form of variation that involves the persistence of multiple discrete phenotypes called morphs. Polymorphism is often explained by negative frequency dependent selection (NFDS), under which rare morphs have an advantage, such that no one morph takes over. We tested whether mating polymorphism is maintained by NFDS generated by female choice, in the livebearing poeciliid fish, Girardinus metallicus, whose males are either common, plain morphs or rare, black morphs. Females were treated with one morph for several weeks and tested with both, in mate choice and direct competition assays. Contrary to NFDS via female choice, females preferred the morph with which they were treated. This may disrupt the polymorphism given the rarity of black morphs in the wild, unless black morphs have other advantages: we found that black morphs tended to exhibit higher mating activity, and other studies have demonstrated that they are more aggressive. Interestingly, only black morphs display to females prior to copulation; however, there is little evidence for female preference for this morph or the mating display. These results suggest functions for the mating display of black morphs beyond courtship, including aggressive mate guarding, and prompt a discussion of what constitutes courtship behavior.

Aetiological Naturalism in the Philosophy of Medicine: A Shaky Project

Griffiths and Matthewson (2018) employ the selected effects theory to contend that disease involves the impairment of the normal functioning of biological items. Since the selected effects theory focuses on the past effects of those items, I refer to their proposal as “aetiological naturalism”. In this paper, I argue that aetiological naturalism cannot constitute an adequate theory of disease. This is due to the fact that the selected effects theory, which lies at the heart of aetiological naturalism, is flawed. One promise of the selected effects theory is indeed that it is able to account for our normative intuitions about dysfunctional biological items by grounding them on the concept of natural selection incorporated in the selected effects theory itself, where this promise rests upon its claim that appealing to the fitness-enhancing effect of biological items can always explain why they persisted in a population. However, I contend that the naturalisation of normativity cannot be cashed out in terms of biological items’ past effects by discussing two biological phenomena: phenotypic plasticity and negative frequency-dependent selection. I illustrate that in both cases the selected effects theory cannot recognise adaptive items as having a selected effects function, therefore preventing the possibility to assign them any dysfunction. The normative force of the selected effects theory is consequently much diminished, leading in turn to a weakening of aetiological naturalism, given that such project actually relies on the assumed normative force of the selected effects theory itself.

MSphere of Influence: Predicting the evolution of pathogen populations.

Tatum D. Mortimer works in the field of pathogen population genomics and evolution. In this mSphere of Influence article, she reflects on how "Frequency-dependent selection can forecast evolution in Streptococcus pneumoniae" by Azarian et al. and "Contingency, repeatability, and predictability in the evolution of a prokaryotic pangenome" by Beavan et al. made an impact on her by highlighting the ways in which genomic data can be used to predict pathogen evolution.

Meta-analytical evidence for frequency-dependent selection across the tree of life.

Explaining the maintenance of genetic variation in fitness-related traits within populations is a fundamental challenge in ecology and evolutionary biology. Frequency-dependent selection (FDS) is one mechanism that can maintain such variation, especially when selection favours rare variants (negative FDS). However, our general knowledge about the occurrence of FDS, its strength and direction remain fragmented, limiting general inferences about this important evolutionary process. We systematically reviewed the published literature on FDS and assembled a database of 747 effect sizes from 101 studies to analyse the occurrence, strength, and direction of FDS, and the factors that could explain heterogeneity in FDS. Using a meta-analysis, we found that overall, FDS is more commonly negative, although not significantly when accounting for phylogeny. An analysis of absolute values of effect sizes, however, revealed the widespread occurrence of modest FDS. However, negative FDS was only significant in laboratory experiments and non-significant in mesocosms and field-based studies. Moreover, negative FDS was stronger in studies measuring fecundity and involving resource competition over studies using other fitness components or focused on other ecological interactions. Our study unveils key general patterns of FDS and points in future promising research directions that can help us understand a long-standing fundamental problem in evolutionary biology and its consequences for demography and ecological dynamics.

Mean-field interacting multi-type birth–death processes with a view to applications in phylodynamics

Multi-type birth–death processes underlie approaches for inferring evolutionary dynamics from phylogenetic trees across biological scales, ranging from deep-time species macroevolution to rapid viral evolution and somatic cellular proliferation. A limitation of current phylogenetic birth–death models is that they require restrictive linearity assumptions that yield tractable message-passing likelihoods, but that also preclude interactions between individuals. Many fundamental evolutionary processes – such as environmental carrying capacity or frequency-dependent selection – entail interactions, and may strongly influence the dynamics in some systems. Here, we introduce a multi-type birth–death process in mean-field interaction with an ensemble of replicas of the focal process. We prove that, under quite general conditions, the ensemble’s stochastically evolving interaction field converges to a deterministic trajectory in the limit of an infinite ensemble. In this limit, the replicas effectively decouple, and self-consistent interactions appear as nonlinearities in the infinitesimal generator of the focal process. We investigate a special case that is rich enough to model both carrying capacity and frequency-dependent selection while yielding tractable message-passing likelihoods in the context of a phylogenetic birth–death model.

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion.

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.

Frequency-dependent seed selection by rodents: Response to seed tannins and sizes

Seed predation by rodents relies on the abundance of available seeds and seed traits like sizes and tannin-content. Here, by tracking the fates of 2400 tagged seeds of two sympatric tree species Quercus wutaishanica and Corylus heterophylla with varying size and tannin content, we assessed rodents’ preference for specific seed traits. Meanwhile, we set the ratio of large to small seeds within species as 9:1, 7:3, 5:5, 3:7, and 1:9, respectively, to examine rodent behavior in response to seed relative frequencies. We found that: 1) rodents preferentially hoarded C. heterophylla seeds with lower-tannins after dispersal while Q. wutaishanica seeds with higher-tannins are more often eaten in situ or after dispersal. Peeling seed coat is an adaptive strategy for rodents to avoid high tannins, particularly for Q. wutaishanica. 2) Rodents prioritize small seeds over large ones for Q. wutaishanica whereas significant differences were not observed for C. heterophylla. Rodents commonly eat small seeds after successful dispersal and prefer to hoard larger ones over long distances. 3) Positive frequency-dependent selection was observed for both seed types. Our results suggest that the selection of seeds by rodents is determined by seed tannins and sizes, and changes in seed size frequencies could alter the interaction between trees and rodents from predation towards mutualism.

Beyond ungulate density: Prey switching and selection by the wolf in a recolonised area

The analysis of mechanisms underlying prey use and selection is crucial for understanding the potential for apex predators to shape ecological processes across trophic levels. Significant implications for the conservation and management of ecosystems can ensue, especially in human-modified landscapes characterised by the recent recovery of apex predators, such as those of Europe. Generalist predators may show frequency-dependent prey selection, leading to prey switching; however, there is no ubiquitous support for positive relationships between prey abundance and prey selection. In a protected area recently recolonised by the wolf Canis lupus and hosting abundant wild prey (>20 individuals/km2), we assessed prey use and selection (>3000 wolf scats analysed) throughout 7 years since wolf settlement (April 2016-to-March 2023), during which prey density was relatively stable. Wolf diet was dominated by wild boar Sus scrofa and fallow deer Dama dama (>70 % of occurrences and volume); the roe deer Capreolus capreolus and coypu Myocastor coypus represented main secondary prey. The wolf showed a clear prey switching throughout the study. In the first two years, the fallow deer was the main prey and was positively selected. Then, wild boar replaced fallow deer as main prey, and was positively selected, with a progressive decrease in the use of fallow deer throughout the years. Selection indices of wild boar were negatively associated with those of fallow deer. The relationship between selection indices and population density was inconsistent across species, being negative for the wild boar and weakly positive for the fallow deer, although the latter showed no sharp decline in the area. Concurrent studies showed increased diurnal activity of fallow deer and a consequent reduced temporal overlap with the wolf throughout years since the settlement of this predator, suggesting temporal avoidance. By supporting a minor role of prey density in mediating prey switching, our results suggest a major role of other factors (e.g., antipredator behavioural responses). Alternative substantial prey are important to buffer apex predators not only towards changes in the abundance of main prey, but also towards the effects of other factors potentially reducing prey detectability, e.g., antipredator behaviour.

Frequency-dependent selection of neoantigens fosters tumor immune escape and predicts immunotherapy response

Cancer is an evolutionary process shaped by selective pressure from the microenvironments. However, recent studies reveal that certain tumors undergo neutral evolution where there is no detectable fitness difference amongst the cells following malignant transformation. Here, through computational modeling, we demonstrate that negative frequency-dependent selection (or NFDS), where the immune response against cancer cells depends on the clonality of neoantigens, can lead to an immunogenic landscape that is highly similar to neutral evolution. Crucially, NFDS promotes high antigenic heterogeneity and early immune evasion in hypermutable tumors, leading to poor responses to immune checkpoint blockade (ICB) therapy. Our model also reveals that NFDS is characterized by a negative association between average clonality and total burden of neoantigens. Indeed, this unique feature of NFDS is common in the whole-exome sequencing (WES) datasets (357 tumor samples from 275 patients) from four melanoma cohorts with ICB therapy and a non-small cell lung cancer (NSCLC) WES dataset (327 tumor samples from 100 patients). Altogether, our study provides quantitative evidence supporting the theory of NFDS in cancer, explaining the high prevalence of neutral-looking tumors. These findings also highlight the critical role of frequency-dependent selection in devising more efficient and predictive immunotherapies.

Inversions Can Accumulate Balanced Sexual Antagonism: Evidence from Simulations and Drosophila Experiments.

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with In(3R)K favoring survival and In(3L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, In(3L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas In(3R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

Frequency-Dependent Ecological Interactions Increase the Prevalence, and Shape the Distribution, of Preexisting Drug Resistance

The evolution of resistance remains one of the primary challenges for modern medicine, from infectious diseases to cancers. Many of these resistance-conferring mutations often carry a substantial fitness cost in the absence of treatment. As a result, we would expect these mutants to undergo purifying selection and be rapidly driven to extinction. Nevertheless, preexisting resistance is frequently observed from drug-resistant malaria to targeted cancer therapies in non-small-cell lung cancer (NSCLC) and melanoma. Solutions to this apparent paradox have taken several forms, from spatial rescue to simple mutation supply arguments. Recently, in an evolved resistant NSCLC cell line, we found that frequency-dependent ecological interactions between ancestor and resistant mutant ameliorate the cost of resistance in the absence of treatment. Here, we hypothesize that frequency-dependent ecological interactions in general play a major role in the prevalence of preexisting resistance. We combine numerical simulations with robust analytical approximations to provide a rigorous mathematical framework for studying the effects of frequency-dependent ecological interactions on the evolutionary dynamics of preexisting resistance. First, we find that ecological interactions significantly expand the parameter regime under which we expect to observe preexisting resistance. Next, even when positive ecological interactions between mutants and ancestors are rare, these resistant clones provide the primary mode of evolved resistance because even weak positive interaction leads to significantly longer extinction times. We then find that even in the case where mutation supply alone is sufficient to predict preexisting resistance, frequency-dependent ecological forces still contribute a strong evolutionary pressure that selects for increasingly positive ecological effects (negative frequency-dependent selection). Finally, we genetically engineer several of the most common clinically observed resistance mechanisms to targeted therapies in NSCLC, a treatment notorious for preexisting resistance. We find that each engineered mutant displays a positive ecological interaction with their ancestor. As a whole, these results suggest that frequency-dependent ecological effects can play a crucial role in shaping the evolutionary dynamics of preexisting resistance. Published by the American Physical Society 2024

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/2,142 (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.

When does antimicrobial resistance increase bacterial fitness? Effects of dosing, social interactions, and frequency dependence on the benefits of AmpC β-lactamases in broth, biofilms, and a gut infection model.

One of the longstanding puzzles of antimicrobial resistance is why the frequency of resistance persists at intermediate levels. Theoretical explanations for the lack of fixation of resistance include cryptic costs of resistance or negative frequency-dependence but are seldom explored experimentally. β-lactamases, which detoxify penicillin-related antibiotics, have well-characterized frequency-dependent dynamics driven by cheating and cooperation. However, bacterial physiology determines whether β-lactamases are cooperative, and we know little about the sociality or fitness of β-lactamase producers in infections. Moreover, media-based experiments constrain how we measure fitness and ignore important parameters such as infectivity and transmission among hosts. Here, we investigated the fitness effects of broad-spectrum AmpC β-lactamases in Enterobacter cloacae in broth, biofilms, and gut infections in a model insect. We quantified frequency- and dose-dependent fitness using cefotaxime, a third-generation cephalosporin. We predicted that infection dynamics would be similar to those observed in biofilms, with social protection extending over a wide dose range. We found evidence for the sociality of β-lactamases in all contexts with negative frequency-dependent selection, ensuring the persistence of wild-type bacteria, although cooperation was less prevalent in biofilms, contrary to predictions. While competitive fitness in gut infections and broth had similar dynamics, incorporating infectivity into measurements of fitness in infections significantly affected conclusions. Resistant bacteria had reduced infectivity, which limited the fitness benefits of resistance to infections challenged with low antibiotic doses and low initial frequencies of resistance. The fitness of resistant bacteria in more physiologically tolerant states (in biofilms, in infections) could be constrained by the presence of wild-type bacteria, high antibiotic doses, and limited availability of β-lactamases. One conclusion is that increased tolerance of β-lactams does not necessarily increase selection pressure for resistance. Overall, both cryptic fitness costs and frequency dependence curtailed the fitness benefits of resistance in this study.