Rapid and repeated evolution of myosin copy number in threespine stickleback

This study demonstrates that threespine sticklebacks have rapidly and repeatedly evolved adaptive expansions of the *MYH3C* gene family on their sex chromosome through tandem duplications, resulting in increased muscle gene expression and distinct ecological divergence between freshwater and marine populations.

The Gist

The Big Picture: Fish, Muscles, and a Genetic "Copy-Paste" Glitch

Imagine the three-spined stickleback fish as a species that is constantly reinventing itself. About 20,000 years ago, when the ice age ended, these fish moved from the ocean into thousands of new freshwater lakes and rivers.

In the ocean, they are like marathon runners: sleek, armored, and built for long-distance swimming. In freshwater, they are like sprinters: they have lost their heavy armor, have deeper bodies, and need to be able to dart away quickly from predators in tight spaces.

This paper discovers how the freshwater fish became such good sprinters. It turns out they didn't just tweak their engines; they found a way to copy and paste a specific part of their genetic manual to build more powerful muscles.

---

The Main Characters: The "Myosin" Factory

Inside every fish, there are genes that act as blueprints for building muscles. One specific family of genes is called Myosin Heavy Chain. Think of these genes as the instruction manuals for building the pistons in a car engine.

• Marine Fish (Ocean): They have a standard set of these manuals. Let's say they have 3 copies of a specific instruction (Gene C3). This is enough for their steady, long-distance swimming.
• Freshwater Fish (Lakes): When they moved to the lakes, they started showing up with 4, 5, or even 6 copies of that same instruction manual.

The Analogy: Imagine you are building a house.

• The ocean fish have one blueprint for the kitchen. They build one kitchen, which is fine for a small family.
• The freshwater fish suddenly have five blueprints for the kitchen. They don't build five separate houses; instead, they use those extra blueprints to build a massive, high-performance kitchen with five ovens. This allows them to cook (swim) much faster.

The Discovery: It Happened Over and Over Again

The researchers looked at fish from all over the world (from Alaska to Iceland). They found that every single time marine fish moved into a freshwater lake, they independently developed this "copy-paste" mutation.

• The "Hotspot" Effect: It's like if you walked into a bakery in New York, London, and Tokyo, and every single baker decided to add a second oven to their shop on the exact same day. This suggests that adding these extra gene copies is a huge advantage for living in freshwater.
• Speed of Change: In some lakes where scientists introduced marine fish just a few years ago, the fish already started copying their genes within six to eight years. Evolution usually takes thousands of years, but here, it happened in the blink of an eye.

How Did They Do It? (The Glitch)

How do you get extra copies of a gene? Nature has a few ways to mess up the copying process, and this fish used two specific tricks:

1. The "Stutter" (MMBIR): Imagine a photocopier that jams. When it tries to copy a page, it gets stuck, then restarts, but it accidentally copies the same paragraph twice before moving on. This happened once to turn the 3-copy version into a 4-copy version.
2. The "Mix-Up" (NAHR): Once the fish had 4 copies, the extra copies were so similar to each other that the cell's machinery got confused during reproduction. It thought, "Oh, I have two of these pages, let me swap them!" This swapping accidentally created a 5th or 6th copy.

It's like having a deck of cards where, every time you shuffle, you accidentally duplicate the Ace of Spades. Eventually, you have a deck with six Aces of Spades instead of one.

Why Does Having More Copies Matter?

The researchers found that having more copies of the gene didn't just make the fish "stronger" in a vague way. It specifically boosted the production of fast-twitch muscle fibers.

• The Engine Analogy: The extra gene copies act like a turbocharger. They don't change the type of engine (it's still a muscle), but they crank up the volume.
• The Result: Freshwater fish can explode with speed to escape predators. Ocean fish, with fewer copies, are built for endurance, not speed.
• The "Volume Knob": The study showed that the extra copies act like a volume knob turned up high. The more copies the fish has, the louder the "muscle building" signal gets, specifically in the muscles used for sprinting.

A Twist: The Boys vs. The Girls

There is one final interesting twist. These genes are located on the X chromosome (the sex chromosome).

• Females (XX): Have two X chromosomes, so they can have a lot of these extra gene copies (up to 12 copies total if both chromosomes are expanded).
• Males (XY): Have only one X chromosome. Even if their X chromosome has the extra copies, they can't have as many total copies as the females.

This means the "sprinting boost" might be stronger in female freshwater fish than in males, potentially leading to differences in how fast the boys and girls can swim.

The Bottom Line

This paper tells the story of a fish that solved a survival problem by accidentally copying its own instruction manual.

Instead of waiting millions of years to slowly evolve a new muscle type, the stickleback fish found a "cheat code" in their DNA. By duplicating a specific gene cluster, they instantly upgraded their engines to become the sprinters of the freshwater world. It's a perfect example of how nature can use simple mistakes (copying errors) to create complex, life-saving adaptations very quickly.

Technical Summary

1. Problem Statement

Copy number variations (CNVs) are a major source of genomic diversity and are linked to both disease and adaptive evolution. However, identifying adaptive CNVs in wild populations is challenging because:

• CNVs are difficult to detect and assemble accurately due to high sequence similarity between duplicated regions.
• It is often unclear whether CNV changes are driven by natural selection or genetic drift.
• The molecular mechanisms driving recurrent CNV expansion in natural populations are not well understood.

The study focuses on the threespine stickleback (Gasterosteus aculeatus), a model organism for studying rapid adaptive evolution following the retreat of glaciers (~20,000 years ago). Marine stickleback repeatedly colonized freshwater habitats, evolving distinct phenotypes (e.g., reduced armor, deeper bodies, burst swimming). The authors sought to identify genomic loci where copy number changes are repeatedly selected during this transition and to elucidate the molecular mechanisms and functional consequences of these changes.

2. Methodology

The researchers employed a multi-faceted approach combining population genomics, long-read sequencing, transcriptomics, and functional assays:

• Population Genomics & Read Depth Analysis:

  • Analyzed whole-genome resequencing data from 96 populations (marine and freshwater) across the Pacific and Atlantic basins.
  • Calculated normalized read depth to estimate copy number of the MYH3C (Myosin Heavy Chain 3 Cluster C) gene cluster on Chromosome XIX.
  • Examined hybrid zones in five river systems to test if copy number differences persist despite gene flow.
  • Analyzed "transplant" populations (marine fish introduced to freshwater lakes in the 1980s–2010s) to observe evolutionary changes in real-time.

• Long-Read Sequencing & De Novo Assembly:

  • Generated high-quality de novo assemblies using PacBio HiFi and Oxford Nanopore long reads for specific marine and freshwater individuals.
  • Sequenced Bacterial Artificial Chromosome (BAC) clones to validate the structure of the MYH3C cluster.
  • Used these assemblies to resolve the complex tandem duplication structures that short-read sequencing could not.

• Phylogenetics & Molecular Evolution:

  • Constructed phylogenetic trees using neutral SNPs to determine if high-copy haplotypes evolved once or repeatedly.
  • Aligned protein sequences to identify amino acid divergence between marine and freshwater alleles and among duplicated copies.

• Expression Analysis:

  • Performed allele-specific expression (ASE) analysis using RNA-seq on F1 hybrids (crossing 3-copy marine and 5-copy freshwater haplotypes) to distinguish expression from different alleles and copies.
  • Utilized PacBio Kinnex (Iso-Seq) long-read sequencing to accurately map transcripts to specific gene copies (C3A, C3B, C3L) which are highly similar in sequence.

• Functional Assays:

  • Created transgenic stickleback by injecting constructs containing the intergenic region between MYH3C2 and MYH3C3 upstream of a GFP reporter to test for enhancer activity.
  • Conducted QTL mapping in a genetic cross to confirm the locus of copy number variation.

3. Key Contributions & Results

A. Recurrent and Rapid Evolution of Copy Number

• Ecotypic Divergence: Freshwater stickleback consistently possess higher copy numbers of the MYH3C cluster (ranging from 4 to 6 copies) compared to marine stickleback (typically 3 copies).
• Repetition: Phylogenetic analysis of neutral SNPs showed that high-copy haplotypes are interspersed with low-copy haplotypes across the tree, indicating that copy number expansion has evolved independently and repeatedly in different geographic lineages.
• Selection in Real-Time: In lakes where marine fish were introduced to freshwater within the last few decades, MYH3C copy number increased significantly within 6–8 years. Selection coefficients ($s$) were estimated between 0.18 and 0.57, indicating strong positive selection.
• Maintenance Despite Gene Flow: In river systems with hybrid zones, upstream freshwater residents maintained high copy numbers while downstream anadromous (marine) fish retained low copy numbers, suggesting strong selection maintains the divergence despite gene flow.
• Latitudinal Gradient: Copy number in freshwater populations positively correlates with latitude, suggesting adaptation to northern, colder environments.

B. Molecular Structure and Mechanisms of Expansion

• Tandem Duplication: The expansion involves the tandem duplication of a ~17.6 kb region containing a full MYH3C3 gene and partial flanking regions (MYH3C2 and SYT19).
• Mechanism 1 (Initial Expansion): The transition from a 3-copy (ancestral) to a 4-copy haplotype likely occurred via Microhomology-Mediated Break-Induced Replication (MMBIR). Breakpoints showed 4–8 bp of microhomology, characteristic of this mechanism.
• Mechanism 2 (Subsequent Expansion): Further expansion from 4 to 5 or 6 copies likely occurred via Non-Allelic Homologous Recombination (NAHR). The long stretches of high sequence identity between the initial duplicates (C3A and C3L) provided the substrate for unequal crossing over.
• Sex Chromosome Linkage: The cluster is located on the X chromosome (ChrXIX). The Y chromosome retains a 3-copy structure but with a pseudogenized C1. This creates sexual dimorphism in gene dosage, with females (XX) potentially having higher muscle myosin dosage than males (XY) in freshwater populations.

C. Functional Consequences

• Expression Dosage: Higher genomic copy number correlates with increased mRNA expression in skeletal muscle. However, the increase is not strictly proportional to copy number (e.g., 5 copies do not yield 5x expression), likely due to regulatory constraints or epigenetic silencing of some copies.
• Regulatory Enhancer: The duplicated intergenic region between C2 and C3 contains a muscle-specific enhancer. Transgenic assays confirmed this region drives GFP expression specifically in skeletal muscle fibers.
• Protein Sequence: While the duplicated C3 copies (C3A, C3B) are nearly identical, the C3L (last) copy often contains frameshifts or premature stop codons, suggesting some copies may be non-functional or evolving new functions. Marine-freshwater divergence in amino acid sequences is enriched in the SH3 domain, potentially affecting myosin motor kinetics.

4. Significance

• CNVs as Evolutionary Hotspots: This study identifies MYH3C as a "hotspot" for adaptive evolution, demonstrating that CNVs can be repeatedly selected in natural populations just as frequently as single nucleotide polymorphisms (SNPs).
• Mechanistic Insight: It provides a clear model for how complex structural variations arise: an initial MMBIR event creates a low-copy repeat, which then facilitates rapid, recurrent NAHR-mediated expansions.
• Phenotypic Adaptation: The findings link gene dosage changes to specific ecological adaptations. The increased expression of fast-twitch myosin isoforms in freshwater fish likely supports the burst swimming required for predator evasion and foraging in complex freshwater environments, contrasting with the endurance swimming of marine fish.
• Sexual Dimorphism: The study highlights how CNVs on sex chromosomes can drive sexual dimorphism in physiological traits (muscle performance) between males and females.
• Methodological Advance: The work underscores the necessity of long-read sequencing and BAC validation to accurately resolve complex, repetitive genomic regions that are often collapsed or misassembled in standard short-read reference genomes.

In conclusion, the paper demonstrates that the rapid, repeated expansion of the MYH3C gene cluster is a key adaptive mechanism in stickleback, driven by specific molecular mechanisms (MMBIR and NAHR) and resulting in functional changes in muscle physiology that facilitate survival in freshwater environments.