The biogeography of the stripped Venus clam, Chamelea gallina, in the Mediterranean Sea indicates limited gene flow and shows evidence of local adaptation.

This study reveals that while the striped Venus clam (*Chamelea gallina*) exhibits high gene flow across the Mediterranean and Northeast Atlantic, local adaptation to environmental gradients drives significant genetic differentiation, highlighting the necessity of incorporating adaptive genomic data to define robust management units for sustainable fisheries.

The Gist

The Big Picture: The "Invisible Wall" in the Ocean

Imagine the Mediterranean Sea as a giant, bustling swimming pool. In this pool lives the Striped Venus Clam (Chamelea gallina), a tasty shellfish that people love to eat. For a long time, scientists thought these clams were like a single, giant family where everyone could easily swim from one end of the pool to the other. They thought the clams were all mixed up, with no real "neighborhoods."

But this new study says: "Not so fast!"

While the clams can swim far (thanks to their babies floating in the currents), they aren't actually all the same. It turns out that the ocean has invisible walls and different "neighborhoods" that make the clams adapt to their specific local environments, even if they can physically travel between them.

The Story in Three Acts

Act 1: The "Passport" vs. The "DNA Test"

To figure this out, the scientists took a DNA "passport" from 226 clams from six different spots: the Gulf of Cadiz (Atlantic side), the Alboran Sea, the Balearic Islands, the Ebro Delta, the Tyrrhenian Sea, and the Adriatic Sea.

They looked at two types of genetic markers:

1. The Neutral Markers (The Passport): These are like a generic ID card that doesn't change much. When they looked at these, the clams from all over looked almost identical. It was like saying, "Hey, everyone in this pool has the same passport!" This suggested that the clams were mixing freely.
2. The Adaptive Markers (The DNA Test): These are the parts of the DNA that actually help the clam survive its specific neighborhood. When they looked at these, the story changed completely. Suddenly, the clams from the Adriatic Sea looked very different from the clams in the Tyrrhenian Sea, and the Atlantic clams were distinct from the Mediterranean ones.

The Analogy: Imagine a group of people who all speak the same language (Neutral DNA). But if you ask them about their local recipes, the person from Italy makes pasta, the person from Spain makes paella, and the person from Greece makes moussaka (Adaptive DNA). Even though they can travel and talk to each other, their "local culture" keeps them distinct.

Act 2: The Ocean's "Traffic Jams"

Why do these differences exist if the clams can swim everywhere? The ocean isn't a smooth, flat surface. It has currents, temperature changes, and salinity shifts that act like traffic jams or borders.

• The Strait of Gibraltar: This is the door between the Atlantic Ocean and the Mediterranean. It's a tough spot to cross, acting like a bouncer at a club.
• The Almeria-Oran Front: This is a giant underwater wall of water where cold, fresh water meets warm, salty water. It's like a river running through the ocean.
• The Adriatic Sea: This is a shallow, enclosed bay with lots of river water pouring in. It's a totally different "room" than the deep, open Tyrrhenian Sea.

The clams that live in the cold, fresh Adriatic have evolved to handle that specific "room." If a clam from the warm Tyrrhenian Sea tries to move there, it might not survive as well. So, over time, the clams in each "room" have tweaked their DNA to fit their specific environment.

Act 3: Why This Matters for Your Dinner Plate

Here is the most important part. These clams are a huge food source, but many fisheries are collapsing because they are being overfished.

For years, managers have treated the whole Mediterranean as one big stock. They thought, "If we catch too many in Italy, the clams from Spain will just swim over and fill the gap."

The study says this is a dangerous mistake.

Because the clams are locally adapted, the "Spanish clams" might not be able to save the "Italian clams" if the Italian population crashes. It's like thinking that if your local bakery burns down, a bakery from a different city with a different recipe will just move in and fix it. They can't; they are too different.

The Takeaway

The ocean looks like one big, connected blue space, but for the Striped Venus Clam, it's actually a patchwork quilt of distinct neighborhoods.

• The Old View: "Everything is mixed up; we can manage it all as one big group."
• The New View: "There are invisible walls and local adaptations. We need to manage each neighborhood separately to keep the clams (and our seafood supply) safe."

By understanding these hidden genetic differences, scientists hope to create better rules to protect the clams from overfishing and climate change, ensuring we can still enjoy them in the future.

Technical Summary

1. Problem Statement

The striped Venus clam (Chamelea gallina) is a commercially vital bivalve in the Mediterranean Sea and Northeast Atlantic, yet its fisheries are facing collapse due to overexploitation and environmental stressors. Effective fisheries management requires defining biologically meaningful stock units. However, previous studies using neutral genetic markers (microsatellites, mtDNA) have yielded conflicting results, suggesting either panmixia (high connectivity) or regional differentiation.

The core problem addressed is the decoupling between neutral genetic connectivity and adaptive divergence. While high dispersal potential (via planktonic larvae) suggests homogenized populations, the Mediterranean Sea is characterized by complex biogeographic boundaries (e.g., the Almería-Oran Front, the Strait of Gibraltar) and environmental gradients (temperature, salinity, nutrients) that may drive local adaptation. The study aims to determine if these environmental barriers create distinct management units that are invisible to neutral markers but critical for the species' evolutionary resilience.

2. Methodology

The study employed a seascape genomics framework to analyze genome-wide Single Nucleotide Polymorphisms (SNPs) across six sampling regions: Gulf of Cádiz (Atlantic), Alboran Sea, Balearic Sea (Gandia and Ebro Delta), Tyrrhenian Sea, and Adriatic Sea ($N=226$ individuals).

• Genotyping & Filtering:

  • Data generated via the DArTseq™ platform.
  • Initial dataset: 20,626 SNPs.
  • Filtering Pipeline: Applied strict quality controls including Minimum Allele Frequency (MAF > 0.05), missing data thresholds (<35%), linkage disequilibrium pruning ($R^2 < 0.8$), and removal of related individuals (siblings).
  • Dataset Partitioning:
    • Neutral Dataset: 1,080 SNPs (loci not under selection, filtered for Hardy-Weinberg Equilibrium).
    • Adaptive Dataset: 213 SNPs (putatively under selection).

• Outlier Detection:

  • Two complementary methods were used to identify loci under selection:
    1. pcadapt: PCA-based method identifying loci associated with genetic axes.
    2. OutFLANK: FST-based method modeling the neutral distribution of differentiation.
  • Result: Only pcadapt identified outliers (213 SNPs); OutFLANK detected none, resulting in no overlap between the two methods.

• Population Structure Analysis:

  • k-means clustering (adegenet) and sparse non-negative matrix factorization (sNMF) (LEA) were applied to both neutral and adaptive datasets to infer the number of genetic clusters ($K$).

• Genotype-Environment Association (GEA):

  • Redundancy Analysis (RDA) was used to link SNP allele frequencies with environmental variables (Sea Surface Temperature, Salinity, Nutrients, etc.).
  • Partial RDA was performed to control for demographic structure (using neutral PC scores as conditioning variables) to reduce false positives.
  • Environmental data sourced from Bio-ORACLE (2000–2014).

3. Key Contributions

• Integration of Adaptive Genomics: The study moves beyond neutral markers to explicitly characterize adaptive divergence in a commercially exploited bivalve, demonstrating that local adaptation can persist despite high gene flow.
• Biogeographic Validation: It validates the role of major Mediterranean oceanographic fronts (Almería-Oran Front, Balearic Front) and basin boundaries (Adriatic vs. Tyrrhenian) as effective barriers to adaptive gene flow.
• Management Implications: It provides a genomic basis for redefining fisheries management units, arguing that current broad-scale management strategies may overlook critical, locally adapted populations.

4. Key Results

A. Population Structure: The Neutral vs. Adaptive Decoupling

• Neutral Loci: Showed weak genetic differentiation.
  • k-means and sNMF suggested a largely panmictic population ($K=1$ or $K=2$ with very subtle structure).
  • Differentiation was primarily observed between the Adriatic Sea and all other sites, and the Gulf of Cádiz and the rest.
  • This confirms high connectivity via larval dispersal across most of the Mediterranean.
• Adaptive Loci: Revealed pronounced genetic structure.
  • k-means identified 5 distinct genetic clusters ($K=5$): Alboran Sea, Adriatic Sea, Gulf of Cádiz, Balearic sites, and Tyrrhenian Sea.
  • sNMF confirmed 4 clusters, grouping Balearic and Alboran sites but separating the Adriatic and Gulf of Cádiz.
  • Conclusion: While neutral markers suggest a single stock, adaptive markers reveal distinct, regionally coherent units.

B. Genetic Diversity

• Heterozygosity: Observed heterozygosity ($H_O$) was generally low but consistent across most sites (mean ~0.10).
• Outliers: The Gandia (Balearic) site showed significantly higher genetic diversity compared to other sites.
• Resilience: Despite historical overfishing and population collapses, neutral genetic diversity remains high, suggesting that demographic collapse has not yet led to significant genetic erosion (a common lag effect in marine species).

C. Genotype-Environment Associations

• Drivers: Three environmental variables were retained as significant predictors of adaptive variation: Mean Sea Surface Temperature (SST), Mean Sea Surface Salinity (SSS), and Mean Nitrate.
• Variance: The RDA explained a small proportion of total variance (<1%), which is typical for polygenic adaptation in large SNP datasets, but the axes were statistically significant ($p=0.001$).
• Candidate SNPs: 19 SNPs on the first axis and 9 on the second were identified as candidate adaptive markers. These loci are likely involved in physiological responses to thermal stress, osmoregulation, and nutrient availability.

5. Significance and Implications

• Fisheries Management: The study argues that management units for C. gallina should be defined by adaptive genetic boundaries rather than neutral connectivity. Treating the Mediterranean as a single stock or broad regional stock risks the loss of locally adapted populations, which are essential for resilience against climate change.
• Climate Change Adaptation: The identification of temperature and salinity as key selective axes highlights that populations are already adapted to specific local conditions. As the Mediterranean warms, these locally adapted units may have different capacities to survive, necessitating region-specific conservation strategies.
• Biogeography: The findings reinforce that the Mediterranean is not a homogenous "sea" but a mosaic of adaptive landscapes shaped by oceanographic fronts. This supports the integration of seascape genomics into biogeographic regionalization schemes.

In summary, the paper demonstrates that local adaptation generates biologically meaningful population structure in Chamelea gallina that is invisible to neutral markers, urging a shift toward adaptive management strategies to ensure the long-term sustainability of this critical fishery.