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