Dynamic genomes uncover opposite sex determination in the invasive quagga and zebra mussels

This study reveals that the invasive zebra and quagga mussels possess strikingly different sex determination architectures, with zebra mussels exhibiting a polygenic ZZ/ZW system and quagga mussels utilizing a localized XX/XY system driven by a novel male-determining *FoxL2-Y* locus, highlighting unexpected reproductive divergence and rapid evolutionary mechanisms in closely related invasive species.

The Gist

Imagine two very successful, invasive cousins who have taken over freshwater lakes across the globe: the Zebra Mussel and the Quagga Mussel. They are notorious for clogging pipes, damaging boats, and outcompeting native wildlife. Scientists have long wanted to control them, perhaps by using "genetic biocontrol"—a fancy way of saying, "Can we tweak their genes to stop them from reproducing?"

To do that, you first need to know how they decide whether a baby mussel becomes a boy or a girl. For decades, this was a mystery. This new study acts like a detective story, using high-tech genomic tools to crack the case. Here is what they found, explained simply.

The Big Discovery: Cousins, But Different

The researchers expected these two mussels to have the same "blueprint" for making boys and girls because they are close relatives. Instead, they found that their blueprints are completely different. It's like finding out that two identical twins have completely different ways of deciding who gets to be the captain of the ship.

1. The Zebra Mussel: The "Committee" System (Polygenic)

Think of the Zebra Mussel's sex determination system as a large committee meeting.

• How it works: There isn't just one "boss" gene telling the mussel to be male or female. Instead, there are multiple genes scattered across different chromosomes (think of them as different departments in a company) that all have to vote.
• The Evidence: The scientists found signals for sex determination on three different chromosomes. It's a "polygenic" system, meaning it relies on the combined effect of many small factors.
• The Key Player: One of the most important "voters" on this committee is a gene called FoxL2. In almost all animals, FoxL2 is the "Queen Bee" that tells cells to become female. In Zebra Mussels, it seems to be a crucial part of the female vote.
• The Analogy: Imagine trying to stop a committee from making a decision. You'd have to change the minds of many different people in different rooms. This makes it very hard to engineer a genetic "fix" to control their population.

2. The Quagga Mussel: The "Switch" System (Localized XY)

The Quagga Mussel, on the other hand, uses a much simpler, more traditional system, like a light switch.

• How it works: They have a specific, localized region on one chromosome (Chromosome 13) that acts as the master switch. If you have this switch, you are male; if you don't, you are female. This is a classic XX/XY system (like humans).
• The Surprise: The scientists found a "rogue" gene in this switch area. It's a copy of the FoxL2 gene (the female gene), but it's been mutated and moved to the male side. Let's call it FoxL2-Y.
• The Twist: This "FoxL2-Y" gene is broken. It has a "stop sign" (a premature stop codon) in the middle of it, so it can't make a normal protein. However, the scientists think it might act like a radio broadcast rather than a physical worker. It might produce a signal (RNA) that shuts down the normal female FoxL2 gene during the very early stages of development, forcing the embryo to become male.
• The Analogy: Imagine a factory where the "Make Female" machine is usually running. The Quagga Mussel has a broken copy of that machine's manual (FoxL2-Y) that, when read, somehow jams the original machine, forcing the factory to build a "Male" product instead.
• The Implication: Because this system relies on a single "switch" (or a small cluster of genes), it is much easier to target with genetic tools like CRISPR. If you can flip that switch, you could theoretically control the population.

Why Does This Matter?

The study reveals that even closely related species can evolve completely different ways of solving the same problem (making boys and girls) in a very short time (about 10 million years).

• For Science: It shows that genomes are incredibly dynamic. Genes can duplicate, move, break, and get repurposed to create entirely new systems. It's like evolution playing with LEGO bricks, building different structures from the same set of pieces.
• For the Environment: This is huge news for fighting invasive species.
  • Zebra Mussels: Because they use a "committee" system, trying to genetically engineer them to stop reproducing is going to be like trying to stop a committee by changing just one person's mind. It's difficult and complex.
  • Quagga Mussels: Because they use a single "switch" system, they might be a much better candidate for genetic biocontrol. If scientists can figure out how to flip that FoxL2-Y switch, they might be able to crash the population.

The Bottom Line

This paper is a masterclass in how nature reinvents itself. The Zebra and Quagga mussels are neighbors, but they speak different genetic languages when it comes to sex. By decoding these languages, scientists are one step closer to understanding how to manage these invasive pests, turning a complex biological mystery into a potential roadmap for saving our waterways.

Technical Summary

1. Problem Statement

The zebra mussel (Dreissena polymorpha) and the quagga mussel (Dreissena rostriformis bugensis) are two of the most ecologically and economically damaging invasive freshwater bivalves globally. Despite their impact, their sex determination (SD) systems remain entirely uncharacterized. Understanding the genetic architecture of SD is a critical prerequisite for developing genetic biocontrol strategies (e.g., gene drives to skew sex ratios) to manage these populations. Previous studies in bivalves have been limited by a lack of chromosome-scale reference genomes and the prevalence of homomorphic (morphologically indistinguishable) sex chromosomes, making the identification of sex-linked loci difficult. Furthermore, it was unknown whether these two closely related species (diverged ~10–12 million years ago) share a conserved SD mechanism or have evolved distinct architectures.

2. Methodology

The authors employed an integrative multi-omics approach combining chromosome-scale genomics, population resequencing, and transcriptomics:

• Genome Assembly & Karyotyping:
  • Generated a high-quality, chromosome-scale reference genome for the quagga mussel using PacBio HiFi reads (308 Gbp) and Illumina paired-end data from a single male.
  • Scaffolding was performed using Hi-C data, resulting in 16 pseudomolecules.
  • Validated the assembly with karyotyping of metaphase spreads from embryos, confirming 16 pairs of chromosomes.
• Population Resequencing:
  • Conducted whole-genome resequencing (WGS) of 40 individuals (20 males, 20 females) for each species (80 total).
  • Aligned reads to species-specific reference genomes (zebra: existing; quagga: newly assembled).
  • Calculated genome-wide diversity statistics and male-female differentiation ($F_{ST}$).
• Sex-Linked Locus Identification:
  • $F_{ST}$ Scans: Identified regions with high male-female differentiation.
  • Allele Frequency Analysis: Screened for ZW (female heterogametic) and XY (male heterogametic) SNP patterns.
  • Model-Based Testing: Used SDpop to calculate Bayesian Information Criterion (BIC) differences between NULL, XY, and ZW models across 10 kb windows.
  • K-mer Analysis: Performed reference-free k-mer analysis (31-mers) to identify sex-specific sequences independent of mapping bias.
• Candidate Gene Characterization:
  • Analyzed gene content within sex-linked regions, focusing on conserved sex-determination gene families (Fox, DMRT, Sox).
  • Identified a specific FoxL2 paralog in the quagga SD region and compared its structure to canonical FoxL2.
• Transcriptomics:
  • Screened 23 embryonic/developmental RNA-seq libraries and newly generated adult tissue-specific RNA-seq data.
  • Used diagnostic k-mers (spanning specific indels and SNVs) to distinguish transcripts from the highly similar FoxL2 and FoxL2-Y loci, overcoming mapping ambiguity.

3. Key Contributions

• First Chromosome-Scale Genome for Quagga Mussel: Provided a high-quality reference (1.61 Gbp, 16 chromosomes) enabling fine-scale genomic analysis previously impossible due to the fragmented nature of the previous assembly.
• Discovery of Divergent SD Architectures: Revealed that two closely related invasive species utilize fundamentally different sex determination systems:
  • Zebra Mussel: A polygenic ZZ-ZW system with multiple sex-linked regions distributed across at least three linkage groups.
  • Quagga Mussel: A localized XX-XY system controlled by a single ~800 kb region on Linkage Group 13.
• Identification of FoxL2-Y: Discovered a novel male-determining candidate locus in quagga mussels (FoxL2-Y), a paralog of the conserved female-associated gene FoxL2. This represents a rare case of a female-associated gene being co-opted for male determination via duplication and divergence.
• Mechanistic Insight: Proposed a model where FoxL2-Y acts as a regulatory RNA or truncated protein during early development to bias sex determination toward the male pathway, potentially driven by haploid selection on gamete-interaction genes (C-lectins) within the SD region.

4. Key Results

Genomic Architecture & Diversity

• Synteny: Despite similar genome sizes and chromosome numbers (16), extensive rearrangements exist between the species (e.g., fission/fusion events involving LG3, LG11, LG12, LG13, LG14).
• Diversity: Quagga mussels exhibit nearly double the nucleotide diversity ($\pi \approx 0.0046$) compared to zebra mussels ($\pi \approx 0.0027$), suggesting higher standing genetic variation.
• Differentiation: Genome-wide $F_{ST}$ between sexes is near zero for both species, confirming the absence of heteromorphic (differentiated) sex chromosomes.

Zebra Mussel: Polygenic ZZ-ZW System

• Signatures: Multiple regions across LG5, LG9, and LG13 showed $F_{ST} \approx 0.5$ and strong support for a ZW model in SDpop.
• Candidate Genes:
  • LG5: Contains three candidate regions, one of which harbors FoxL2 (a known female-determining gene in bivalves).
  • LG9 & LG13: Contain transcription factors, Tudor domain proteins, and long non-coding RNAs (lncRNAs).
• Conclusion: Sex is likely determined by the additive or epistatic effects of multiple loci (polygenic), rather than a single master switch.

Quagga Mussel: Localized XX-XY System

• Signatures: All sex-linked signals converged on a single ~800 kb region on LG13 (Q.13.1).
  • High male-female $F_{ST}$ and XY SNP clustering.
  • Coverage Bias: Males showed higher sequencing coverage than females in this region, indicating the presence of Y-linked sequence absent in females.
  • K-mer Analysis: Confirmed a sharp enrichment of male-specific k-mers and near-zero female-specific k-mers in this region.
• Candidate Gene (FoxL2-Y):
  • Located within the SD region; derived from a duplication of the canonical FoxL2.
  • Structural Divergence: Contains a large intronic repeat insertion, a 6-amino-acid deletion, and a premature stop codon within the forkhead (FH) DNA-binding domain.
  • Expression: Detected at extremely low levels specifically in trochophore and early veliger stages (early development) but absent in adults and later stages.
• Conclusion: Suggests a single master locus system where FoxL2-Y likely functions as a non-canonical regulator (possibly via RNA interference or truncated protein) to suppress the female pathway.

5. Significance

• Evolutionary Biology: Demonstrates the extreme lability of sex determination systems. Closely related species can evolve from a polygenic system to a localized master locus system (or vice versa) over short evolutionary timescales (~10 Myr) without the formation of heteromorphic chromosomes.
• Mechanism of Turnover: Supports the hypothesis that haploid selection (acting on gamete interaction genes like C-lectins enriched in the quagga SD region) can facilitate the rapid turnover of SD systems by linking new sex-determining alleles to loci under strong selection.
• Biocontrol Implications:
  • Quagga Mussel: The discovery of a single, localized XY region with a specific candidate gene (FoxL2-Y) suggests that genetic biocontrol (e.g., CRISPR-based gene drives to skew sex ratios) is theoretically more feasible and tractable in this species.
  • Zebra Mussel: The polygenic nature of their SD system implies that targeting a single gene may be insufficient, as sex determination relies on the cumulative dosage of multiple loci, making biocontrol strategies significantly more complex.
• Genomic Resource: The new chromosome-scale quagga genome serves as a critical resource for future studies on invasive species genomics, adaptation, and reproductive biology.