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.

Punctuated virus-driven succession generates dynamical alternations in CRISPR-mediated microbe-virus coevolution.

The coevolutionary dynamics of lytic viruses and microbes with CRISPR-Cas immunity exhibit alternations between sustained host control of viral proliferation and major viral epidemics in previous computational models. These alternating dynamics have yet to be observed in other host-pathogen systems. Here, we address the breakdown of control and transition to large outbreaks with a stochastic eco-evolutionary model. We establish the role of host density-dependent competition in punctuated virus-driven succession and associated diversity trends that concentrate escape pathways during control phases. Using infection and escape networks, we derive the viral emergence probability whose fluctuations of increasing size and frequency characterize the approach to large outbreaks. We explore alternation probabilities as a function of non-dimensional parameters related to the probability of viral escape and host competition. Our results demonstrate how emergent feedbacks between host competition and viral diversification render the host immune structure fragile, potentiating a dynamical transition to large epidemics.

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.

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.

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.

Resource Variation Within and Between Patches: Where Exploitation Competition, Local Adaptation, and Kin Selection Meet.

AbstractIn patch- or habitat-structured populations, different processes can favor adaptive polymorphism at different scales. While spatial heterogeneity can generate spatially disruptive selection favoring variation between patches, local competition can lead to locally disruptive selection promoting variation within patches. So far, almost all theory has studied these two processes in isolation. Here, we use mathematical modeling to investigate how resource variation within and between habitats influences the evolution of variation in a consumer population where individuals compete in finite patches connected by dispersal. We find that locally and spatially disruptive selection typically act in concert, favoring polymorphism under a wider range of conditions than when in isolation. But when patches are small and dispersal between them is low, kin competition inhibits the emergence of polymorphism, especially when the latter is driven by local competition for resources. We further use our model to clarify what comparisons between trait and neutral genetic differentiation (QST/FST comparisons) can tell about the nature of selection. Overall, our results help us understand the interaction between two major drivers of polymorphism: locally and spatially disruptive selection, and how this interaction is modulated by the unavoidable effects of kin selection under limited dispersal.

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.

Jeans and language: kin networks and reproductive success are associated with the adoption of outgroup norms.

Traditional norms of human societies in rural China may have changed owing to population expansion, rapid development of the tourism economy and globalization since the 1990s; people from different ethnic groups might adopt cultural traits from outside their group or lose their own culture at different rates. Human behavioural ecology can help to explain adoption of outgroup cultural values. We compared the adoption of four cultural values, specifically speaking outgroup languages/mother tongue and wearing jeans, in two co-residing ethnic groups, the Mosuo and Han. Both groups are learning outgroup traits, including each other's languages through contact in economic activities, education and kin networks, but only the Mosuo are starting to lose their own language. Males are more likely to adopt outgroup values than females in both groups. Females of the two groups are no different in speaking Mandarin and wearing jeans, whereas males do differ, with Mosuo males being keener to adopt them than Han males. The reason might be that Mosuo men experience more reproductive competition over mates, as Mosuo men have larger reproductive skew than others. Moreover, Mosuo men but not others gain fitness benefits from the adoption of Mandarin (they start reproducing earlier than non-speakers). This article is part of the theme issue 'Social norm change: drivers and consequences'.

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.

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.

Emergence of novel non-aggregative variants under negative frequency-dependent selection in Klebsiella variicola.

Klebsiella variicola is an emergent human pathogen causing diverse infections, some of which in the urinary tract. However, little is known about the evolution and maintenance of genetic diversity in this species, the molecular mechanisms and their population dynamics. Here, we characterized the emergence of a novel rdar-like (rough and dry) morphotype which is contingent both on the genetic background and the environment. We show that mutations in either the nitrogen assimilation control gene (nac) or the type III fimbriae regulator, mrkH, suffice to generate rdar-like colonies. These morphotypes are primarily selected for the reduced inter-cellular aggregation as a result of MrkH loss-of-function which reduces type 3 fimbriae expression. Additionally, these clones also display increased growth rate and reduced biofilm formation. Direct competitions between rdar and wild type clones show that mutations in mrkH provide large fitness advantages. In artificial urine, the morphotype is under strong negative frequency-dependent selection and can socially exploit wild type strains. An exhaustive search for mrkH mutants in public databases revealed that ca 8% of natural isolates analysed had a truncated mrkH gene many of which were due to insertions of IS elements, including a reported clinical isolate with rdar morphology. These strains were rarely hypermucoid and often isolated from human, mostly from urine and blood. The decreased aggregation of these mutants could have important clinical implications as we hypothesize that such clones could better disperse within the host allowing colonisation of other body sites and potentially leading to systemic infections.

Host adaptation drives genetic diversity in a vector-borne disease system.

The range of hosts a pathogen can infect is a key trait, influencing human disease risk and reservoir host infection dynamics. Borrelia burgdorferi sensu stricto (Bb), an emerging zoonotic pathogen, causes Lyme disease and is widely considered a host generalist, commonly infecting mammals and birds. Yet the extent of intraspecific variation in Bb host breadth, its role in determining host competence, and potential implications for human infection remain unclear. We conducted a long-term study of Bb diversity, defined by the polymorphic ospC locus, across white-footed mice, passerine birds, and tick vectors, leveraging long-read amplicon sequencing. Our results reveal strong variation in host breadth across Bb genotypes, exposing a spectrum of genotype-specific host-adapted phenotypes. We found support for multiple niche polymorphism, maintaining Bb diversity in nature and little evidence of temporal shifts in genotype dominance, as would be expected under negative frequency-dependent selection. Passerine birds support the circulation of several human-invasive strains (HISs) in the local tick population and harbor greater Bb genotypic diversity compared with white-footed mice. Mouse-adapted Bb genotypes exhibited longer persistence in individual mice compared with nonadapted genotypes. Genotype communities infecting individual mice preferentially became dominated by mouse-adapted genotypes over time. We posit that intraspecific variation in Bb host breadth and adaptation helps maintain overall species fitness in response to transmission by a generalist vector.

The maintenance of genetic diversity under host-parasite coevolution in finite, structured populations.

As a corollary to the Red Queen hypothesis, host-parasite coevolution has been hypothesized to maintain genetic variation in both species. Recent theoretical work, however, suggests that reciprocal natural selection alone is insufficient to maintain variation at individual loci. As highlighted by our brief review of the theoretical literature, models of host-parasite coevolution often vary along multiple axes (e.g. inclusion of ecological feedbacks or abiotic selection mosaics), complicating a comprehensive understanding of the effects of interacting evolutionary processes on diversity. Here we develop a series of comparable models to explore the effect of interactions between spatial structures and antagonistic coevolution on genetic diversity. Using a matching alleles model in finite populations connected by migration, we find that, in contrast to panmictic populations, coevolution in a spatially structured environment can maintain genetic variation relative to neutral expectations with migration alone. These results demonstrate that geographic structure is essential for understanding the effect of coevolution on biological diversity.

Ecotype formation and prophage domestication during gut bacterial evolution.

How much bacterial evolution occurs in our intestines and which factors control it are currently burning questions. The formation of new ecotypes, some of which capable of coexisting for long periods of time, is highly likely in our guts. Horizontal gene transfer driven by temperate phages that can perform lysogeny is also widespread in mammalian intestines. Yet, the roles of mutation and especially lysogeny as key drivers of gut bacterial adaptation remain poorly understood. The mammalian gut contains hundreds of bacterial species, each with many strains and ecotypes, whose abundance varies along the lifetime of a host. A continuous high input of mutations and horizontal gene transfer events mediated by temperate phages drives that diversity. Future experiments to study the interaction between mutations that cause adaptation in microbiomes and lysogenic events with different costs and benefits will be key to understand the dynamic microbiomes of mammals. Also see the video abstract here: https://youtu.be/Zjqsiyb5Pk0.

Shell color polymorphism in marine gastropods.

Marine gastropods are characterized by an incredible variation in shell color. In this review, we aim to introduce researchers to previous studies of shell color polymorphism in this group of animals, trying to provide an overview of the topic and highlighting some potential avenues for future research. For this, we tackle the different aspects of shell color polymorphism in marine gastropods: its biochemical and genetic basis, its patterns of spatial and temporal distribution, as well as its potential evolutionary causes. In particular, we put special emphasis on the evolutionary studies that have been conducted so far to reveal the evolutionary mechanisms responsible for the maintenance of shell color polymorphism in this group of animals, as it constitutes the least addressed aspect in existing literature reviews. Several general conclusions can be drawn from our review: First, natural selection is commonly involved in the maintenance of gastropod color polymorphism; second, although the contribution of neutral forces (gene flow-genetic drift equilibrium) to shell color polymorphism maintenance do not seem to be particularly important, it has rarely been studied systematically; third, a relationship between shell color polymorphism and mode of larval development (related to dispersal capability) may exist. As for future studies, we suggest that a combination of both classical laboratory crossing experiments and -Omics approaches may yield interesting results on the molecular basis of color polymorphism. We believe that understanding the various causes of shell color polymorphism in marine gastropods is of great importance not only to understand how biodiversity works, but also for protecting such biodiversity, as knowledge of its evolutionary causes may help implement conservation measures in those species or ecosystems that are threatened.

The impact of habitat loss on molecular signatures of coevolution between an iconic butterfly (Alcon blue) and its host plant (Marsh gentian).

Habitat loss is threatening natural communities worldwide. Small and isolated populations suffer from inbreeding and genetic drift, which jeopardize their long-term survival and adaptive capacities. However, the consequences of habitat loss for reciprocal coevolutionary interactions remain poorly studied. In this study, we investigated the effects of decreasing habitat patch size and connectivity associated with habitat loss on molecular signatures of coevolution in the Alcon blue butterfly (Phengaris alcon) and its most limited host, the marsh gentian (Gentiana pneumonanthe). Because reciprocal coevolution is characterized by negative frequency-dependent selection as a particular type of balancing selection, we investigated how signatures of balancing selection vary along a gradient of patch size and connectivity, using single nucleotide polymorphisms (SNPs). We found that signatures of coevolution were unaffected by patch characteristics in the host plants. On the other hand, more pronounced signatures of coevolution were observed in both spatially isolated and in large Alcon populations, together with pronounced spatial variation in SNPs that are putatively involved in coevolution. These findings suggest that habitat loss can facilitate coevolution in large butterfly populations through limiting swamping of locally beneficial alleles by maladaptive ones. We also found that allelic richness (Ar) of the coevolutionary SNPs is decoupled from neutral Ar in the butterfly, indicating that habitat loss has different effects on coevolutionary as compared with neutral processes. We conclude that this specialized coevolutionary system requires particular conservation interventions aiming at generating a spatial mosaic of both connected and of isolated habitat to maintain coevolutionary dynamics.

Prevalence of Viral Frequency-Dependent Infection in Coastal Marine Prokaryotes Revealed Using Monthly Time Series Virome Analysis.

Viruses infecting marine prokaryotes have a large impact on the diversity and dynamics of their hosts. Model systems suggest that viral infection is frequency dependent and constrained by the virus-host encounter rate. However, it is unclear whether frequency-dependent infection is pervasive among the abundant prokaryotic populations with different temporal dynamics. To address this question, we performed a comparison of prokaryotic and viral communities using 16S rRNA amplicon and virome sequencing based on samples collected monthly for 2 years at a Japanese coastal site, Osaka Bay. Concurrent seasonal shifts observed in prokaryotic and viral community dynamics indicated that the abundance of viruses correlated with that of their predicted host phyla (or classes). Cooccurrence network analysis between abundant prokaryotes and viruses revealed 6,423 cooccurring pairs, suggesting a tight coupling of host and viral abundances and their "one-to-many" correspondence. Although stable dominant species, such as SAR11, showed few cooccurring viruses, a fast succession of their viruses suggests that viruses infecting these populations changed continuously. Our results suggest that frequency-dependent viral infection prevails in coastal marine prokaryotes regardless of host taxa and temporal dynamics. IMPORTANCE There is little room for doubt that viral infection is prevalent among abundant marine prokaryotes regardless of their taxa or growth strategy. However, comprehensive evaluations of viral infections in natural prokaryotic communities are still technically difficult. In this study, we examined viral infection in abundant prokaryotes by monitoring the monthly dynamics of prokaryotic and viral communities at a eutrophic coastal site, Osaka Bay. We compared the community dynamics of viruses with those of their putative hosts based on genome-based in silico host prediction. We observed frequent cooccurrence among the predicted virus-host pairs, suggesting that viral infection is prevalent in abundant prokaryotes regardless of their taxa or temporal dynamics. This likely indicates that frequent lysis of the abundant prokaryotes via viral infection has a considerable contribution to the biogeochemical cycling and maintenance of prokaryotic community diversity.

Balancing selection and the crossing of fitness valleys in structured populations: diversification in the gametophytic self-incompatibility system.

The self-incompatibility locus (S-locus) of flowering plants displays a striking allelic diversity. How such a diversity has emerged remains unclear. In this article, we performed numerical simulations in a finite island population genetics model to investigate how population subdivision affects the diversification process at a S-locus, given that the two-gene architecture typical of S-loci involves the crossing of a fitness valley. We show that population structure slightly reduces the parameter range allowing for the diversification of self-incompatibility haplotypes (S-haplotypes), but at the same time also increases the number of these haplotypes maintained in the whole metapopulation. This increase is partly due to a higher rate of diversification and replacement of S-haplotypes within and among demes. We also show that the two-gene architecture leads to a higher diversity in structured populations compared with a simpler genetic architecture, where new S-haplotypes appear in a single mutation step. Overall, our results suggest that population subdivision can act in two opposite directions: it renders S-haplotypes diversification easier, although it also increases the risk that the self-incompatibility system is lost.