Strain-specific structural variant landscapes shape mutation retention following mutagenesis in Caenorhabditis elegans

This study reveals that in *Caenorhabditis elegans*, strain-specific structural variant architectures significantly influence mutation retention dynamics following mutagenesis, with higher outcrossing rates paradoxically leading to greater structural variant burdens and reduced purging of deleterious mutations.

The Gist

Imagine the genome of a worm (Caenorhabditis elegans) as a massive, intricate instruction manual for building a living organism. Usually, when this manual gets damaged, the organism tries to fix it or, if the damage is too bad, the organism dies off. This is nature's way of "purging" bad errors.

For a long time, scientists thought this cleaning process worked the same way for all types of errors. They believed that if worms mated with each other (outcrossing), it would help shuffle the deck and easily throw out the bad cards (mutations). But this study suggests that the story is much more complicated, especially when the "errors" are huge structural changes rather than just a single typo.

Here is the breakdown of what the researchers found, using some everyday analogies:

1. The Experiment: Stress-Testing the Worms

The researchers took three different strains of worms (let's call them Team N2, Team AB1, and Team CB4856). Think of these teams as different families with different personalities:

• Team N2: The "lab rat." They are very domesticated, rarely mate with others, and mostly reproduce by themselves (selfing).
• Team CB4856: The "wild card." They are from Hawaii, very active, and mate with each other frequently (high outcrossing).
• Team AB1: The "middle child." They are somewhere in between.

The scientists subjected all three teams to a chemical "stress test" (mutagenesis) for five generations. This was like shaking the instruction manual violently to introduce typos and tear out pages. Then, they let the worms recover for three generations without the stress to see what happened.

2. The Big Surprise: The "Wild Card" Kept the Most Damage

You might expect that the team that mates the most (Team CB4856) would be the best at shuffling out the bad mutations. After all, mixing genes usually helps fix errors.

But the opposite happened.

• Team CB4856 (The High-Outcrossers): Even though they mated a lot, they ended up keeping the most damage. They retained huge chunks of torn pages (Structural Variants) and many typos (SNPs).
• Team N2 (The Selfers): Even though they didn't mate much, they managed to keep their instruction manual much cleaner. They purged the damage more effectively.

The Analogy: Imagine you have a messy room.

• Team N2 is like someone who lives alone. They don't invite anyone in, so they can quickly spot the trash, pick it up, and throw it out immediately.
• Team CB4856 is like a person who throws a huge party every night. They think, "If we mix everything around, we'll find the trash!" But instead, the party creates more mess. The guests (new genetic combinations) accidentally hide the trash under piles of clothes, making it harder to get rid of.

3. The "Linkage Block" Problem

Why did the party-goers keep the trash? The study found that the damage wasn't just tiny typos; it was massive structural changes—like whole paragraphs or pages being deleted, duplicated, or flipped upside down.

When these big structural changes happen, they act like glue. They stick nearby sections of the manual together.

• If a "bad page" is glued to a "good page," you can't throw away the bad page without destroying the good one.
• In the highly mating worms (CB4856), these "glued blocks" were so large that the natural selection process couldn't separate the bad mutations from the good ones. The bad mutations were "hitchhiking" on the structural damage.

4. The "Hitchhiking" Effect

The researchers found that in the wild worms (CB4856), the tiny typos (SNPs) were often stuck inside the big structural damage.

• Analogy: Imagine a bad apple (a mutation) is glued inside a giant, heavy box (a structural variant). If you try to throw away the bad apple, you have to throw away the whole box. But if the box is too heavy or valuable, you might decide to keep the whole thing, even with the bad apple inside.
• Because the mating worms kept these "heavy boxes," they ended up keeping a lot of bad apples too.

5. The Transposons: The "Glitter Bombs"

The study also looked at "Transposable Elements" (TEs). Think of these as glitter bombs or sticky notes that can jump around the manual and stick to new places.

• The wild worms (CB4856) had the most active glitter bombs.
• These jumping elements caused a lot of the structural damage.
• Even though the wild worms were good at mixing genes, the sheer volume of these jumping elements created so much chaos that the "cleaning crew" (natural selection) couldn't keep up.

6. The Simulation: A Computer Proof

To make sure this wasn't just a fluke, the scientists built a computer model. They simulated worm populations with different mating rates.

• Result: The computer confirmed that when worms mate a lot, they actually retain more of these big structural errors, especially if the errors are slightly harmful but not deadly. The "glue" effect of the structural variants prevents the bad mutations from being purged.

The Takeaway

This study changes how we think about evolution and cleaning up genetic damage.

• Old Idea: Mating is always good for cleaning up genetic errors.
• New Idea: Mating is great for small errors (typos), but if the damage is huge (structural variants), mating can actually make it harder to get rid of the bad stuff. The "glue" of large structural changes traps the bad mutations, allowing them to survive even in populations that are mixing their genes frequently.

In short: Sometimes, keeping your circle small and stable (like the lab worms) is actually better at cleaning up a massive mess than throwing a big, chaotic party (like the wild worms). The "party" just ends up hiding the mess in the corners.

Technical Summary

1. Problem Statement

Classical evolutionary theory posits that outcrossing (sexual reproduction) facilitates the purging of deleterious mutations through recombination, which separates harmful alleles from beneficial ones. However, this theory has largely focused on Single Nucleotide Polymorphisms (SNPs). Structural Variants (SVs)—including large insertions, deletions, duplications, and inversions—pose a unique challenge because they can span large genomic regions and suppress local recombination, potentially creating "linkage blocks" that hinder the removal of deleterious mutations.

The study addresses a critical gap: How do mating systems (specifically outcrossing rates) and genomic architecture influence the retention of SVs versus SNPs following acute mutagenic stress? It is unclear whether high outcrossing rates effectively purge large-scale structural mutations or if SVs persist due to their ability to suppress recombination, and how this varies across different genetic backgrounds.

2. Methodology

The researchers employed an integrated approach combining experimental evolution, whole-genome sequencing (long- and short-read), and population genetic simulations.

• Experimental Design:

  • Strains: Three genetically distinct C. elegans isolates were used: N2 (Bristol, lab-adapted, low outcrossing), AB1 (intermediate), and CB4856 (Hawaiian, wild isolate, high outcrossing).
  • Mutagenesis: Populations were exposed to two mutagens for five consecutive generations:
    • EMS (Ethyl methanesulfonate): Primarily induces SNPs.
    • Formaldehyde: Induces complex DNA lesions and SVs.
  • Recovery: Following mutagenesis, populations were allowed to recover for three generations without mutagen exposure.
  • Replication: Four replicate populations per strain per treatment (12 total experimental populations).

• Phenotypic Assays:

  • Male Frequency & Outcrossing: Quantified to assess mating behavior changes. Outcrossing was defined by male progeny production thresholds specific to each strain.
  • Relative Fitness: Measured via competitive assays against a GFP-marked reference strain (ST2) to determine if mutagenesis caused irreversible fitness decline.

• Genomic Analysis:

  • Sequencing:
    • PacBio HiFi (Long-read): Used for parental strains and one replicate per strain/treatment to detect large SVs and complex rearrangements.
    • Illumina (Short-read): Used for all replicates to quantify SNPs and validate SVs.
  • Bioinformatics:
    • De novo mutation calling (SNPs via Freebayes; SVs via Sniffles, Manta, and VG toolkit).
    • Transposable Element (TE) Analysis: Used the TransposonUltimate pipeline to distinguish TE-mediated SVs from other structural changes.
    • Breakpoint Analysis: Classified SV breakpoints by microhomology signatures (NHEJ, MMEJ, SSA, HMR) to infer DNA repair mechanisms.

• Computational Modeling:

  • SLiM 4.3 Simulations: Developed a stochastic population model to test how varying outcrossing rates affect the retention of SVs under different fitness landscapes (synergistic epistasis vs. linear).

3. Key Contributions

1. Strain-Specific Mutation Landscapes: Demonstrated that genetic background significantly dictates the accumulation and retention of mutations, independent of the mutagen type.
2. SV Retention Paradox: Challenged the assumption that high outcrossing always purges deleterious variation. Found that the strain with the highest outcrossing rate (CB4856) retained the highest burden of SVs.
3. Coupling of SNPs and SVs: Revealed that SNPs are non-randomly enriched within SV intervals, suggesting that SVs create genomic "islands" where nucleotide mutations accumulate and are shielded from recombination-based purging.
4. TE Dynamics: Identified that while TEs contribute to SVs, they account for a minority of the total SV burden, and their activity is strain-specific (highest in CB4856).

4. Key Results

• Phenotypic Responses:

  • CB4856 exhibited the highest baseline male frequency and outcrossing rates. Mutagenesis reduced male frequency in CB4856 but had no significant effect on N2 or AB1.
  • Fitness Recovery: All strains recovered relative fitness to ancestral levels within three generations, indicating that short-term fitness metrics do not reflect the underlying genomic load.

• Mutation Retention (SNPs vs. SVs):

  • CB4856 (High Outcrossing): Retained the highest number of de novo SNPs and SVs. Notably, SVs in CB4856 spanned ~40–43% of the genome, including megabase-scale inversions and duplications.
  • N2 & AB1 (Low/Moderate Outcrossing): Retained fewer SVs. N2 showed a high enrichment of SNPs within SV regions but a lower total SV burden.
  • Exonic Impact: CB4856 had the highest number of mutations overlapping coding exons, increasing the potential for functional disruption despite fitness recovery.

• Genomic Architecture & Recombination:

  • SNP Enrichment: SNPs were significantly enriched within SV intervals across all strains (2.6–4.9x expectation in N2/AB1; ~1.9x in CB4856).
  • Large SVs: Mutagenesis induced massive structural changes, including a 2.49 Mb inversion in N2 and inversions spanning nearly half the genome in CB4856.
  • Repair Mechanisms: Breakpoint analysis showed that while repair signatures (NHEJ, MMEJ, SSA) varied by strain, they did not correlate directly with SV abundance, suggesting SV burden is driven by genome architecture rather than just repair pathway efficiency.

• Transposable Elements (TEs):

  • TEs accounted for <15% of de novo SVs.
  • Zator elements were the most active transposons across all strains.
  • CB4856 showed ~4x more transposition events than N2 or AB1, correlating with its high outcrossing rate and genomic instability.

• Simulation Findings:

  • Population simulations confirmed that outcrossing rates >30% can lead to the retention of SVs if the mutations have small deleterious effects and exhibit synergistic epistasis. In fully selfing populations, SVs were rapidly purged. This supports the experimental finding that high outcrossing does not guarantee the removal of structural variants.

5. Significance

• Redefining Mutation Purging: The study suggests that the "purging" effect of outcrossing is mutation-size dependent. While outcrossing may effectively purge SNPs, it may inadvertently retain SVs by preventing the recombination necessary to break up large linkage blocks created by the SVs themselves.
• Genomic Architecture Matters: The physical organization of the genome (SV landscape) dictates how mutations are distributed and retained. Strains with high outcrossing rates may accumulate a "hidden" load of structural variants that are not immediately visible in fitness assays.
• Implications for Evolutionary Theory: These findings refine models of mutation accumulation, highlighting that mating systems interact with mutation size to shape evolutionary trajectories. It suggests that in populations with high recombination, large structural variants may persist longer than previously predicted, potentially driving rapid genomic divergence or adaptation.
• Model Organism Context: The results highlight the limitations of relying solely on the N2 strain (low outcrossing) for evolutionary studies, as wild isolates like CB4856 exhibit fundamentally different genomic responses to stress.

In conclusion, the paper establishes that strain-specific genomic architecture and mating behavior jointly determine the retention of structural variants, challenging the notion that outcrossing universally enhances the purging of deleterious genetic load.