🧬 biology

Beavers on the selection-drift seesaw: imprints of the past and present-day processes on the MHC genes variability in the Eurasian beaver

This study demonstrates that the Eurasian beaver's MHC gene variability was shaped by historical positive and balancing selection, severely reduced by a 20th-century population bottleneck, and only partially restored through recent expansion without evidence of current selection.

Original authors: Jan Náhlovský, Marco Heurich, Steven K. Windels, Frank N. Rosell, Christian A. Robstad, Alexander P. Saveljev, Aleš Vorel, Pavel Munclinger

Published 2026-08-31

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Jan Náhlovský, Marco Heurich, Steven K. Windels, Frank N. Rosell, Christian A. Robstad, Alexander P. Saveljev, Aleš Vorel, Pavel Munclinger

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

In the living world, every animal carries a biological toolkit designed to recognize and fight off the invisible army of viruses, bacteria, and parasites that constantly threaten its health. A crucial part of this toolkit is a set of genes known as the major histocompatibility complex, or MHC. Think of these genes as the instruction manual for building a sophisticated security system on the surface of cells. This system scans for foreign invaders and sounds the alarm if it finds something that does not belong. Because the threats in nature are so diverse and constantly changing, the instructions for this security system must be incredibly varied. If a population has many different versions of these genes, it is more likely that at least some individuals will be able to recognize a new disease and survive. However, when a population shrinks drastically, this genetic variety can vanish, leaving the survivors with a much weaker defense.

Scientists have long been interested in how these immune genes change over time, especially in species that have faced near-extinction. The Eurasian beaver offers a perfect, albeit tragic, story for this kind of study. Once roaming freely across Europe and Asia, these industrious builders were hunted to the brink of disappearance by the early 1900s. Only about 1,200 individuals remained, scattered in tiny, isolated pockets of habitat. Conservationists worked tirelessly to save them, moving survivors from these small groups to new areas to rebuild the population. Today, there are roughly 1.5 million beavers again, but they are a mix of descendants from those original, separated groups. Researchers wanted to know what happened to the beavers' immune system during this dramatic crash and recovery. Did the population lose its ability to fight disease? Did the mixing of different groups restore that lost strength? And is the immune system still evolving today, or has it settled into a new normal?

A team of researchers from universities and conservation organizations across Europe and Russia set out to answer these questions by looking directly at the beavers' DNA. They focused on two specific parts of the immune system's instruction manual, known as the DRB and DQA genes. These sections are particularly important because they are the parts of the protein that actually grab onto the invading pathogens. The team collected blood and tissue samples from beavers in the few remaining original "relict" populations that survived the bottleneck, as well as from the newly formed populations that grew from translocations. They also gathered samples from the North American beaver, a close relative, to see how the two species compared. By sequencing the DNA, they could count exactly how many different versions of these immune genes existed in each group and look for signs of how nature had shaped them.

The results painted a clear picture of a species that was battered by history but is now recovering. In the tiny, isolated groups of beavers that survived the 20th-century bottleneck, the genetic variety of the immune system was shockingly low. In some of these groups, every single beaver carried the exact same version of the immune gene, a situation that would be considered a genetic dead end for most other animals. This confirmed that the severe population crash had acted like a powerful filter, wiping out most of the ancient variety through a process known as genetic drift, where chance events determine which genes survive. However, the story did not end there. When the researchers looked at the newly formed populations, where beavers from different original groups had been mixed together, they found that the genetic diversity had bounced back. The new populations carried a wider array of immune gene versions, essentially reassembling the toolkit from the fragments left behind in the different relict groups.

The study also looked back in time to understand how these genes behaved before the population crash. By analyzing the specific changes in the DNA sequences, the researchers found strong evidence that the beaver immune system had been under intense pressure from diseases in the distant past. The genes showed signs of "positive selection," meaning that nature had actively favored new, different versions of the genes to keep up with evolving pathogens. They also found evidence of "balancing selection," a mechanism that helps keep multiple versions of a gene alive in a population over long periods. This suggests that for most of the beaver's history, a diverse immune system was essential for survival. Yet, despite this deep history of adaptation, the researchers found no evidence that the beavers are currently undergoing similar evolutionary changes in their new, large populations. The genes appear to be stable, and the mixing of populations has restored the necessary variety without requiring new mutations to arise right now.

Interestingly, the two immune genes they studied did not behave exactly the same way. One of them, the DRB gene, showed a rich history of adaptation and a good recovery of diversity. The other, the DQA gene, was much more uniform and showed almost no variation, even between the two beaver species. This difference suggests that while the beaver immune system as a whole has a complex history, different parts of it respond to pressure in different ways. The fact that the beavers have managed to rebuild their population numbers so successfully, despite the loss of so much genetic variety during the bottleneck, suggests that the current mix of genes is sufficient for them to thrive in their current environment. However, the researchers note that this restored diversity is still a patchwork of what was left after the crash. While the beavers have made a remarkable comeback, the loss of ancient genetic variety means they might be more vulnerable if they encounter a completely new type of disease in the future, a risk that remains as they continue to expand into new territories.

Technical Summary: MHC Gene Variability in the Eurasian Beaver

Problem Statement
Major Histocompatibility Complex (MHC) genes are critical for adaptive immunity and are typically characterized by extreme polymorphism maintained by balancing and positive selection. However, studying MHC in wild populations is methodologically challenging due to frequent locus duplications, high allele counts, and complex recombination patterns. The Eurasian beaver (Castor fiber) presents a unique case study: the species underwent a severe demographic bottleneck in the early 20th century, reducing the global population to approximately 1,200 individuals across isolated relict populations. This event fragmented the species' range and was followed by a rapid recovery supported by translocations from these relict populations. While mitochondrial DNA and microsatellite markers have documented the genetic consequences of this bottleneck and subsequent mixing, the specific dynamics of MHC variability—particularly regarding the interplay of historical selection, the bottleneck effect, and contemporary selection in newly formed mixed populations—remained unclear. Furthermore, previous studies on beaver MHC were limited by partial sequencing and incomplete sampling of eastern relict populations.

Methodology
The study utilized a comprehensive sampling strategy involving 91 individuals from relict populations (including previously unsampled Belarus and Voronezh refugia) and newly formed mixed populations in Central Europe and Russia. Additionally, 30 North American beavers (Castor canadensis) were sampled for comparative analysis.

Genetic Targets: The researchers focused on the second exon of the MHC II DRB and DQA loci.
Sequencing Strategy: To avoid the amplification of pseudogenes reported in previous studies using older primers, new primers were designed based on the North American beaver genome. The study employed Sanger sequencing for most samples and Illumina amplicon sequencing to validate heterozygosity and resolve complex genotypes. RNA was also sequenced to confirm gene expression.
Analytical Framework:
- Recombination: Detected using RDP v.4.101 and GARD.
- Selection in the Past: Positive selection was assessed via Z-tests (Nei-Gojobori and Kumar methods) on Antigen Binding Site (ABS) vs. non-ABS regions. The McDonald-Kreitman (MK) test and its extension by Soni et al. (2022) were used to estimate the proportion of nonsynonymous polymorphisms maintained by balancing selection. Site-specific selection was analyzed using FEL, FUBAR, and MEME methods.
- Population Genetics: Allelic richness, heterozygosity, and $F_{ST}$ were calculated using FSTAT.
- Contemporary Selection: The study tested for Divergent Allele Advantage (DAA) to determine if heterozygotes with highly divergent alleles were overrepresented, a signature of current balancing selection.

Key Results

Locus Structure and Diversity: The study confirmed that the Eurasian beaver possesses a simplified MHC structure with no locus duplications for DRB and DQA, with a maximum of two alleles per individual.
- Relict Populations: Only 12 DRB alleles were found across all Eurasian beaver relict populations. Many alleles were private to specific relict groups, though some sharing occurred between geographically distant populations (e.g., Germany and Belarus).
- North American Comparison: The sister species, C. canadensis, showed higher DRB diversity with 18 alleles found in the total North American sample (including 12 alleles specifically in the Minnesota population), but was monomorphic at the DQA locus.
- Newly Formed Populations: Mixed populations in Central Europe and Russia exhibited increased diversity (4–6 DRB alleles per population) compared to relict populations (1–3 alleles), resulting in significantly higher allelic richness, observed heterozygosity, and gene diversity.
Historical Selection Signals:
- Positive Selection: Strong signals of positive selection were detected in the DRB locus, particularly within the Antigen Binding Site (ABS). The ratio of non-synonymous to synonymous substitutions was high, and the MK test extension indicated that approximately 92% of non-synonymous mutations were maintained by balancing selection.
- Recombination: Clear evidence of intralocus recombination was found in the DRB sequences, complicating phylogenetic reconstruction but not negating the selection signals.
- Locus Divergence: The DQA locus showed significantly lower variability than DRB and appeared to be under purifying selection, contrasting with the positive/balancing selection observed in DRB.
Bottleneck and Drift Effects: The severe bottleneck resulted in a drastic reduction of MHC variability, with some relict populations fixed for single alleles at both loci. This reduction was more pronounced than typically observed in neutral markers, consistent with theoretical expectations that MHC diversity declines faster during bottlenecks.
Contemporary Selection: Despite the restoration of genetic diversity in newly formed mixed populations, the study found no evidence of ongoing positive or balancing selection acting on MHC loci in current populations. Tests for Divergent Allele Advantage (DAA) yielded no significant results, suggesting that the current allele combinations are not being actively shaped by pathogen-driven selection in the sampled timeframe.

Significance and Claims
The paper claims to provide a refined understanding of how MHC variability is shaped by the oscillating forces of selection and drift in a species with a well-documented demographic history.

Simplified Model: The study highlights the Eurasian beaver as a valuable model for immunogenetics due to its non-duplicated, low-variability MHC system, which allows for clearer dissection of evolutionary processes compared to complex, duplicated systems.
Historical vs. Current Dynamics: The authors demonstrate that while historical positive and balancing selection created a diverse ancestral gene pool, the recent bottleneck caused a severe erosion of this diversity through genetic drift.
Restoration via Translocation: The mixing of individuals from different relict populations during reintroduction successfully restored MHC diversity in newly established populations, a phenomenon the authors note as an exceptional case compared to other reintroduced vertebrates where diversity is often lost due to founder effects.
Absence of Current Selection: Crucially, the study modestly concludes that while diversity has been restored, there is currently no detectable signature of natural selection acting on these MHC loci in the contemporary populations. The authors suggest this may be due to the complexity of demographic history masking selection signals, insufficient statistical power, or the possibility that the bottleneck itself eliminated certain pathogens, temporarily altering selective pressures.

The paper does not claim that the restored diversity guarantees long-term resilience against novel pathogens, noting that the retaining low diversity of MHC loci could constitute a potential risk for future disease challenges.

Technical Summary: MHC Gene Variability in the Eurasian Beaver

More like this