Chromosomal inversions have long fascinated evolutionary geneticists because of what they can do structurally. When a stretch of chromosome breaks and reattaches in the opposite orientation, heterozygotes, individuals carrying one inverted copy and one standard copy, produce fewer viable recombinant gametes. The result is that genes trapped inside an inversion tend to be inherited together as a block, shielded from the shuffling effects of recombination. This makes inversions powerful vehicles for adaptation: co-adapted gene combinations that suit a particular environment can be preserved intact rather than broken apart each generation. Yet the same property raises a puzzle. If an inversion carries a combination of genes that is strongly beneficial in one setting, theory suggests it should sweep to fixation, eliminating the standard arrangement. Why, then, do so many inversions persist as polymorphisms, with both arrangements circulating in the same population for long periods?

The Mank lab’s study tackles this question in a system ideally suited to the task. Trinidadian guppies live in rivers where waterfalls create natural barriers, and above and below these barriers, populations experience radically different regimes of predation. Downstream, guppies coexist with predatory fish that exert intense natural selection. Upstream, predators are absent or rare, and the dominant selective force becomes sexual selection, as females choose among males on the basis of their colour ornaments. These paired high-predation and low-predation populations differ across a wide suite of morphological, behavioural and life-history traits, and they have replicated this divergence independently in multiple rivers, a natural experiment repeated by the landscape itself.

Rather than hunting for inversions at known locations, the researchers applied an unbiased method designed to detect them anywhere in the genome. They sampled guppies from three Trinidadian rivers, each containing a high-predation and a low-predation population. The scan revealed 22 inversions scattered across the genome. Using linked-read sequencing technology, which provides long-range information about how genetic variants are physically connected along chromosomes, the team was able to verify the precise breakpoints of all but one of these inversions, a level of validation that lends considerable confidence to the catalogue.

The ages of the inversions proved startling. Population-genetic analysis indicated that they range from roughly 1 to 6 million years old, and that they predate the colonisation of each of the three rivers by the guppies now living there. In other words, these are not recent mutations that arose locally in response to local conditions. They are ancient structural variants that arrived in Trinidad’s waters already established, and have been carried along as guppies founded new river systems. The inversions are also widespread, distributed across the genome rather than clustered on a single chromosome, suggesting that no part of the guppy’s genetic architecture is immune to this kind of structural change.

When the researchers tested whether inversion frequencies differed between high- and low-predation populations, three inversions emerged as significantly associated with the local adaptation syndromes that distinguish the two environments. Guppies from high-predation localities tend to be plainer, leaner and faster; those from low-predation localities are more colourful and carry traits favoured by mate choice. The association between specific inversions and these contrasting life histories suggests that the gene blocks locked inside them contribute to the coordinated shifts that predation pressure produces. Notably, however, none of the three inversions was reciprocally fixed across all three replicate rivers. The same inversion might be common below the waterfalls in one river and far less differentiated in another, a pattern that hints at the contingent, population-specific nature of adaptation even when the underlying selective pressures are similar.

Alongside these three adaptive candidates, the team observed seven additional inversions maintained in all three rivers without any evidence of association with the predation gradient. Their persistence demanded an explanation. Simulations modelling the expected level of inversion polymorphism under a neutral model, in which variants drift by chance alone, revealed that the observed diversity far exceeds what neutrality can accommodate. Something must be actively maintaining both arrangements of these chromosomes, preventing either from drifting to loss or sweeping to fixation.

The most striking clue came from a comparison with earlier work on guppy colouration. The researchers found a significant overlap between the polymorphic inversions and genomic loci previously implicated in variation in male ornament patterns, the pigmented spots and splashes that females scrutinise when choosing mates. This overlap suggests a mechanism long proposed in theory but rarely demonstrated with genomic data in the wild: negative frequency-dependent selection arising from female preference for novelty. If females favour males with rare or unusual colour patterns, then whichever arrangement of an inversion is currently common becomes disadvantageous, and the rarer arrangement gains a mating advantage. The two chromosomal arrangements are thereby trapped in a perpetual seesaw, each rising when the other falls, and neither ever winning outright. Such balancing selection can preserve structural variation for millions of years, exactly the timescales the team inferred for the guppy inversions.

The study’s broader significance lies in how it reframes the interplay between the two great engines of evolution. Natural selection, driven here by predation, clearly acts on some inversions, pushing their frequencies up or down in different environments. But sexual selection, driven by female choice, appears to be the force that keeps other inversions polymorphic over geological time, replenishing the raw material of male ornament variation on which mate choice feeds. The two processes do not simply act in opposition or in parallel; they interact in shaping the frequency dynamics of the same structural variants. For evolutionary biologists, the guppy results provide a vivid demonstration that the maintenance of genetic diversity is not merely a story of heterozygote advantage or environmental heterogeneity, but can hinge on the ever-shifting tastes of choosy females.

There are also practical implications for how researchers read the genome. Because inversions suppress recombination across large genomic neighbourhoods, association studies that ignore them can misattribute the effects of genes inside inversions, and conservation geneticists assessing the adaptive potential of fragmented populations may underestimate the hidden structural diversity those populations carry. The guppy work, conducted with an unbiased detection method and validated breakpoints, offers a template for surveying inversions in other natural systems where ancient polymorphisms may be quietly shaping traits, from plumage in birds to life histories in insects. For now, the humble Trinidadian guppy has added another chapter to its scientific biography: a fish whose colours have long charmed biologists turns out to owe part of that charm to chromosomal architecture laid down millions of years before any guppy swam a Trinidadian river.

Subject of Research: Chromosomal inversion polymorphisms and sexually selected traits in natural guppy populations

Article Title: Ancient inversion polymorphisms associate with sexually selected traits across natural guppy populations

Article References: Lin, Y., van der Bijl, W., & Mank, J. E. (2026). Ancient inversion polymorphisms associate with sexually selected traits across natural guppy populations. PLOS Biology, 24(9), e3004024. https://doi.org/10.1371/journal.pbio.3004024

Image Credits: AI Generated

DOI: 10.1371/journal.pbio.3004024

Keywords: guppies, chromosomal inversions, sexual selection, natural selection, local adaptation, balancing selection, frequency-dependent selection, Trinidad, population genetics, male colouration, PLOS Biology, structural variants

The team’s central contribution is a matrix representation of genotypic selection. Instead of writing out separate equations for every possible genotype, the authors encode the fitness consequences of every allele pairing in a single symmetric matrix, which they call the fitness-effect matrix. Each element of this matrix quantifies how much a particular genotype departs from a neutral baseline: a value of zero means the genotype enjoys no selective advantage or disadvantage, while positive or negative values signal fitness gains or losses. Because the matrix is symmetric, with the fitness of the heterozygote carrying allele i and allele j identical to that of the reverse pairing, the entire selection regime at a locus with n alleles is compressed into n(n+1)/2 numbers. This compact encoding exploits a deep mathematical symmetry: the same matrix structures that govern deterministic selection also describe the variance of random genetic drift, allowing the two great engines of evolutionary change to be handled within one coherent formalism.

One of the framework’s most appealing features is its flexibility. The authors show that the matrix representation accommodates an impressive range of selection regimes without modification. Additive selection, where each allele contributes independently to fitness; multiplicative selection, where effects combine like compound interest; frequency-dependent selection, where the advantage of an allele shifts as its prevalence changes; temporally varying selection, where the environment flips fitness rankings across generations; and heterozygote advantage, the famous mechanism that maintains sickle-cell anaemia at high frequencies in malaria-endemic regions, all fall naturally within the same mathematical structure. This generality matters because heterozygote advantage and other multiallelic interactions are precisely the cases where two-allele approximations break down most dramatically, producing stable polymorphisms and oscillatory dynamics that biallelic models misinterpret or miss entirely.

The biological motivation for the work is easy to appreciate. The major histocompatibility complex, the gene cluster that underpins the vertebrate immune system’s ability to recognise pathogens, is among the most polymorphic regions of the human genome, with thousands of alleles described. The ABO blood group system, with its three principal alleles and their combinatorial genotypes, shapes transfusion medicine and has been linked to susceptibility to various diseases. Genes underlying monogenic disorders likewise often harbour multiple disease variants. At each of these loci, an individual carries only two alleles, but the population as a whole can carry many, generating a rich spectrum of genotypes whose relative fitnesses determine how variation is maintained or lost. Understanding selection in these settings requires tools that treat the full genotype matrix, not a collapsed two-allele shadow of it.

Describing selection is only half the problem; the other half is measuring it. Here the study makes a second, arguably more practical, contribution. The authors derive a procedure for estimating the fitness-effect matrix from allele frequency trajectories, the kind of data now routinely generated by experimental evolution and by ancient DNA and ecological time-series sampling. The logic is elegant: given a candidate matrix, one can predict how allele frequencies should change from one generation to the next, and then compare those predictions with the observed trajectory. The researchers formalise this comparison in a cost function, weighted by the Shahshahani metric, a geometric measure of distance between allele frequency vectors that arises naturally from the variance structure of genetic drift. When the candidate matrix equals the true one, the cost vanishes; otherwise, discrepancies between predicted and observed frequencies inflate the cost, and a numerical optimisation procedure adjusts the matrix until the best fit is found.

Several technical refinements make the estimation procedure robust in practice. Because relative fitnesses are only defined up to a multiplicative constant, the authors fix one genotype as a reference with a vanishing fitness-effect, removing a redundancy that would otherwise make the problem ill-posed. They also impose symmetry on the matrix throughout the optimisation, ensuring that the heterozygote fitnesses remain consistent regardless of allele ordering. To handle the practical realities of noisy, sparse data, the cost function includes a small constant that prevents infinite weightings when an allele’s frequency dips toward zero, and a ridge penalty that shrinks estimated fitness effects toward zero, suppressing the wild values that small samples can otherwise produce. The authors report that their results are largely insensitive to the first of these tuning parameters, while the optimal strength of the ridge penalty decreases with population size, a pattern they verified systematically across simulations.

The framework also handles an awkward gap that plagues real datasets: sampling. Experimental populations are rarely censused every single generation, and natural populations even less so. The authors extend their method to trajectories in which allele frequencies are known only at irregular intervals. In this case, the prediction is propagated forward through all the unsampled generations, with each intermediate step feeding into the next, until the model reaches the next sampling point. The cost is then computed only at the sampled generations. This means the method can extract fitness information from patchy time-series, provided the underlying dynamics remain governed by the same selection regime across the gaps, an assumption the authors are careful to state.

To demonstrate the approach on real data, the team applied their estimation procedure to time-series data from an experimental yeast evolution study. Yeast populations, with their rapid generations and tractable haplotype tracking, provide an ideal proving ground for microevolutionary inference. The analysis illustrated that multiallelic fitness interactions, including signatures consistent with heterozygote advantage, could be characterised directly from haplotype frequency data, without the need to reduce the system to pairwise biallelic comparisons. This demonstration is significant because it shows the method working end to end: from raw frequency trajectories, through the optimisation machinery, to a fitted matrix of genotype-specific fitness effects that tells a coherent evolutionary story.

The broader implications reach across several fields. For experimental evolution, the framework offers a principled way to quantify fitness landscapes at highly polymorphic loci, where epistatic interactions between alleles within a gene can be mapped and compared across environments. For conservation and evolutionary biology, it provides tools to detect balancing selection at immune genes and other multiallelic systems in natural populations, using the increasingly affordable time-series data generated by ecological genomics programmes. For medical genetics, a better handle on multiallelic selection could sharpen understanding of how disease-associated variation persists in populations despite its costs. The work, funded in part by the European Research Council under the Horizon 2020 programme, is published open access, and the authors have released the full mathematical derivations, including the stochastic dynamics under a Wright-Fisher model and a diffusion approximation, in the article and its supplementary material.

None of this means the two-allele models that dominate the literature are wrong; they remain excellent approximations when a locus genuinely harbours few variants, and their simplicity has powered a century of insight. But as genomic data grows richer and time-series sampling becomes routine, the loci that matter most, from the immune genes that fight our infections to the variants that shape our blood and our diseases, increasingly demand a mathematics that matches their diversity. The fitness-effect matrix of Vellnow, Gossmann and Waxman provides exactly that: a compact, flexible and estimable description of selection that finally looks beyond two alleles and sees the full, tangled richness of genetic variation as evolution actually experiences it.

Subject of Research: A matrix framework for modelling and estimating multiallelic genotypic selection from allele frequency time-series data

Article Title: Beyond two alleles: multiallelic genotypic selection and its estimation from time-series data

Article References: Vellnow, N., Gossmann, T. I., & Waxman, D. (2026). Beyond two alleles: multiallelic genotypic selection and its estimation from time-series data. BMC Genomics. https://doi.org/10.1186/s12864-026-13411-5

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13411-5

Keywords: population genetics, multiallelic loci, genotypic selection, time-series data, heterozygote advantage, genetic drift, fitness estimation, Wright-Fisher model, experimental evolution, yeast, BMC Genomics, balancing selection

The expectations were not arbitrary. Genes involved in allorecognition, the ability to distinguish self from non-self among members of the same species, are classically assumed to be driven by balancing selection, the evolutionary force that maintains many different versions of a gene within a population. This pattern is abundantly documented in the major histocompatibility complex of vertebrates, in fungal heterokaryon incompatibility loci such as het-c, in plant self-incompatibility systems, and in the allorecognition genes of hydroids and sponges. Under this model, rare alleles enjoy a selective advantage, frequency-dependent dynamics keep many variants circulating, and the protein products of the gene diversify into ever more distinctive forms, each tuned to recognize a different repertoire of non-self partners. When the researchers set out to characterize BHF across natural and laboratory populations of B. schlosseri, they anticipated finding precisely these signatures.

Instead, they found something stranger. The team analyzed full-length BHF alleles from nineteen colonies, drawn both from laboratory-born stocks and wild populations, and uncovered thirty-three distinct alleles. At the DNA level, the locus is a riot of variation: nucleotide diversity is pronounced, alleles vary extensively in length, and the genealogy of the sequences shows clear evidence of intragenic recombination, with different segments of the gene telling contradictory evolutionary stories. Yet when the researchers translated the alleles into their protein products, the picture collapsed into remarkable uniformity. Those thirty-three alleles encode only seventeen distinct protein variants. The explosive divergence seen in the nucleotide sequence is largely silent at the level of the protein, a pattern that immediately suggested strong functional constraints keeping the encoded molecule stable even as its underlying blueprint drifts.

The geographic distribution of these alleles added another layer of intrigue. Highly divergent alleles, sequences so different from one another that they appear to have separated in the distant past, tend to coexist within the same individual colony, a hallmark of long-lived polymorphism. At the same time, identical alleles were found to be shared across continental-scale distances, a pattern difficult to reconcile with purely local evolution. In a species whose larvae disperse only over short ranges, such global allele sharing points to movement mediated by humans: B. schlosseri is a notorious fouling organism that hitches rides on ship hulls and in ballast water, and previous work by the same group and others has shown that historical navigation routes have left detectable footprints on the species’ seascape genetics. The new findings suggest that this human-mediated gene flow has effectively homogenized BHF alleles across ocean basins, mixing variant repertoires that might otherwise have diverged in isolation.

To probe the evolutionary forces at work, the researchers subjected the allele alignments to a battery of statistical tests designed to detect natural selection at individual codons. Contrary to the predictions of the diversifying-allele model, no evidence of balancing selection or directional positive selection emerged from the analyses. What the tests did reveal were signatures of purifying selection, the process by which natural selection removes harmful mutations and conserves protein function. In other words, the vast majority of the nucleotide variation in BHF is either silently synonymous or otherwise constrained from changing the protein it encodes, and the few amino-acid substitutions that do occur appear tolerated only within narrow limits. For a gene whose job is, by definition, to discriminate among an enormous diversity of non-self surfaces, this conservation of protein structure is genuinely counterintuitive.

The authors propose that the peculiar contrast between nucleotide and protein diversity at BHF arises from the combined effects of several processes acting simultaneously. Intragenic recombination continually shuffles DNA segments between alleles, generating new nucleotide combinations and inflating sequence diversity without necessarily altering protein sequences, particularly when crossovers occur in regions where the variants are functionally equivalent or where synonymous sites dominate. Human-mediated transport spreads the resulting allele repertoire globally, ensuring that any given colony can harbor partners for its most divergent alleles. Meanwhile, linkage to nearby loci that genuinely are under balancing selection could drag neutral and nearly neutral variation at BHF along for the ride, inflating apparent diversity at the locus without any direct selective diversification of the gene itself. This hitchhiking scenario echoes mechanisms documented in other systems, where polymorphisms maintained at one gene cast a long shadow of diversity across adjacent regions of the genome.

The technical groundwork for the study rested on methods refined over decades of research into the B. schlosseri fusion-rejection system, first described genetically by Abe Sabbadin in 1962 and later linked to a MHC-like gene system by Scofield and colleagues in 1982. The chromosomal location of the histocompatibility locus was resolved in 2005 by De Tomaso and coworkers, and the identity of BHF itself was established by Voskoboynik and colleagues in a landmark 2013 Science paper. Recombination analysis in the new study employed modern tools including RDP5, which screens sequence data for exchange events using multiple independent detection algorithms, while tests of codon-specific selection drew on maximum-likelihood frameworks that estimate synonymous and non-synonymous substitution rates across the gene. Exon- and intron-based phylogenies were compared to expose topological incongruence, the fingerprint of recombination, and predicted pairwise interaction outcomes were modeled on the assumption that colonies sharing at least one protein allele would fuse.

The broader significance of the finding extends beyond a single marine invertebrate. If extensive protein diversification is not a universal driver of allorecognition gene evolution, then the textbook equation of self/non-self recognition systems with runaway diversifying selection may need qualification. It raises the possibility that some recognition systems achieve functional diversity through means other than amino-acid innovation, for instance through the combinatorial pairing of a modest set of protein variants, or through differences in expression and regulation that are invisible to sequence-based analyses. It also cautions researchers studying other invertebrate allorecognition loci, from the Alr genes of hydroids to the aggregation factors of sponges, against assuming that polymorphism at the nucleotide level necessarily translates into adaptive protein diversification, since recombination, demographic mixing and linkage can each generate misleading signals.

For the field of comparative immunology, the BHF locus remains an indispensable touchstone precisely because it sits at the boundary between innate and adaptive immunity, in an organism whose blood-forming system has recently been shown to resemble that of mammals in surprising complexity. The new work does not diminish the gene’s predictive power over fusion and rejection; it simply reframes how that power is maintained. Purifying selection, recombination, and a globe-trotting lifestyle, it turns out, can conspire to keep a recognition gene both endlessly variable in its DNA and stubbornly conservative in its protein, a combination that ensures B. schlosseri will continue to puzzle and instruct immunologists for years to come.

Subject of Research: Evolutionary genetics of the Botryllus histocompatibility factor allorecognition locus in the colonial tunicate Botryllus schlosseri

Article Title: Intra-colony divergence and global allele sharing reflect purifying selection and recombination at the Botryllus histocompatibility factor locus

Article References: Intra-colony divergence and global allele sharing reflect purifying selection and recombination at the Botryllus histocompatibility factor locus. (n.d.). https://doi.org/10.1007/s00251-026-01413-2

Image Credits: AI Generated

DOI: 10.1007/s00251-026-01413-2

Keywords: Botryllus schlosseri, allorecognition, BHF, purifying selection, intragenic recombination, histocompatibility, tunicate immunity, balancing selection, polymorphism, self/non-self recognition, immunogenetics, gene flow