Genetic population structure and demographic history of Pacific cod in Japanese waters: Implications for stock identification using SNP markers

This study utilizes genome-wide SNP markers to reveal three distinct genetic groups and unique demographic histories of Pacific cod in Japanese waters, demonstrating that carefully selected outlier SNPs significantly enhance stock identification accuracy for sustainable fisheries management.

The Gist

Imagine the Pacific Ocean as a massive, bustling city. In this city, the Pacific cod is a popular resident, swimming everywhere from the cold north to the warmer south. For a long time, fishery managers (the city planners) have tried to manage these fish by drawing lines on a map, assuming that all the cod in one area are the same family, and all the cod in another area are a different family.

But here's the problem: Fish don't carry ID cards, and they don't always stay within the lines humans draw. They mix, mingle, and migrate. So, how do we know who belongs to which "family" to make sure we don't overfish them?

This paper is like a genetic detective story that uses high-tech DNA tools to figure out exactly how these cod families are related in Japanese waters.

The Big Discovery: Three Distinct "Tribes"

The researchers took DNA samples from nearly 500 cod caught at 33 different spots around Japan. They didn't just look at a few genes; they looked at thousands of tiny genetic markers (SNPs) across the entire genome. Think of this as checking thousands of pages in a family history book rather than just looking at a single photo.

They discovered that the cod aren't just one big mixed-up crowd. Instead, they belong to three distinct genetic tribes:

1. The "Broad Range" Tribe (JBR): These are the wanderers. They are found all over the place, from the Pacific coast to the Sea of Japan. They are the most numerous and diverse.
2. The "Northernmost" Tribe (NHH): These are the locals. They stick tightly to the very top of Honshu (near a place called Mutsu Bay) and the southern coast of Hokkaido. They are famous for "natal homing," which means they are like homing pigeons that return to their exact birthplace to have babies.
3. The "Western Sea" Tribe (WSJ): These are the westerners, found exclusively in the Western Sea of Japan. They are genetically quite different from the other two, almost like a separate branch of the family tree that went its own way a long time ago.

The Family History: Surviving the Ice Age

The researchers also acted like time travelers. By analyzing the DNA of a few representative fish, they reconstructed the population's history going back to the last Ice Age (about 11,000 to 100,000 years ago).

Imagine the Ice Age as a massive freeze that forced everyone into small, crowded shelters. When the ice melted, the families expanded again, but they did it differently:

• The Broad Range tribe bounced back the fastest and grew the biggest.
• The Western Sea tribe shrank the most during the ice age and has been growing slowly since.
• The Northernmost tribe was somewhere in the middle.

This proves these aren't just random groups; they have been on separate life journeys for thousands of years.

The "Homing Pigeon" Mystery

One of the coolest findings involves the Northernmost (NHH) tribe. Scientists have long suspected that cod in Mutsu Bay return to that specific bay to spawn. This study confirmed it genetically.

However, there was a twist. Old fishing records (tagging studies from 30 years ago) suggested these fish mostly stayed near Mutsu Bay. But the new DNA data showed they are also found further north along the coast. It's like realizing that while your family reunion is in your hometown, your cousins are actually living in the next town over, too. This suggests that fish populations can change over time, and we need to keep watching them closely.

The Solution: A Better "ID Card" System

The ultimate goal of this research is to help fishery managers. If you catch a fish in a mixed school, how do you know which tribe it belongs to? This is called Genetic Stock Identification (GSI).

• The Old Way: To get a 90% accurate answer, you used to need to check 500+ genetic markers. It was like trying to identify a person by asking them 500 questions. It was slow and expensive.
• The New Way: The researchers found that if you pick the right 8 to 10 specific markers (called "outlier" markers), you get the same 90% accuracy. These special markers are like the person's unique fingerprint or a distinctive scar that instantly tells you who they are.

Why Does This Matter?

Think of fish stocks like a bank account. If you don't know how much money is in the account (how many fish there are) or which specific account you are withdrawing from, you might accidentally bankrupt the family.

By understanding that there are three distinct tribes with different histories and behaviors, managers can:

1. Stop guessing: They can stop treating all cod in Japan as one big group.
2. Protect the vulnerable: If the "Western Sea" tribe is shrinking, they can protect that specific area without worrying about the "Broad Range" tribe.
3. Use better tools: They can now use a small, cheap, and fast DNA test (the 8-marker panel) to monitor the fish, ensuring the fishery remains sustainable for the future.

In short, this paper gives us a genetic map of the Pacific cod, proving that even in the vast ocean, fish have their own neighborhoods, family histories, and unique identities that we must respect to keep the ocean healthy.

Technical Summary

1. Problem Statement

Pacific cod (Gadus macrocephalus) is a critical fishery resource in the North Pacific. While stock management units are defined by administrative and geographical boundaries (e.g., seven sea areas around Japan), these often lack biological validation.

• The Challenge: Marine fish populations often exhibit subtle genetic differentiation due to high gene flow, making stock identification difficult. Traditional methods (microsatellites, mtDNA) lack the resolution to detect fine-scale structures or dynamic spatiotemporal mixing.
• Specific Gaps:
  • The genetic distinctness of local spawning aggregations, such as the Mutsu Bay group in northern Honshu, remains unclear.
  • There is a need to optimize Single Nucleotide Polymorphism (SNP) marker panels for Genetic Stock Identification (GSI) to reduce costs while maintaining high assignment accuracy.
  • The demographic history and divergence timing of these populations relative to the Last Glacial Period (LGP) are not fully understood.

2. Methodology

The study employed a multi-faceted genomic approach combining reduced-representation sequencing, whole-genome sequencing (WGS), and advanced statistical modeling.

• Sampling and Genotyping:

  • Samples: 496 Pacific cod individuals collected from 33 sites around the Japanese archipelago (plus 2 samples from the Yellow Sea).
  • Technique: Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di).
  • Data Processing: Initial 517 samples yielded 3.8 billion reads. After filtering for duplicates, contamination, and low genotyping rates, 496 high-quality samples were retained.
  • SNP Datasets:
    • Dataset 1: 55,702 SNPs (including singletons) for diversity metrics.
    • Dataset 2: 6,035 LD-pruned SNPs (MAF > 0.05, $r^2 < 0.1$) for population structure.
    • Dataset 3: Integrated Dataset 2 with 2 Yellow Sea WGS samples.

• Population Structure Analysis:

  • PCA & DAPC: Principal Component Analysis (PCA) and Discriminant Analysis of Principal Components (DAPC) were used to identify genetic clusters.
  • Phylogenetics: Maximum-likelihood trees (IQ-TREE) and pairwise $F_{ST}$ calculations (Weir & Cockerham) assessed differentiation.
  • Isolation-by-Distance (IBD): Mantel tests using least-cost path distances (avoiding land/shallow water) to correlate genetic and geographic distances.

• Demographic History:

  • WGS: Whole-genome sequencing of three representative individuals (one per genetic group).
  • PSMC: Pairwise Sequentially Markovian Coalescent analysis to reconstruct effective population size ($N_e$) changes since the LGP.

• Marker Selection for GSI:

  • Framework: A resampling-based cross-validation (1,000 iterations) using DAPC.
  • Strategy: Compared "background" SNPs (neutral + outlier) against "outlier" SNPs detected via OutFLANK and pcadapt.
  • Metric: Estimated Correct Assignment Rate (ECAR) to determine the minimum number of markers required for >90% accuracy.

3. Key Results

A. Genetic Population Structure

The analysis revealed three distinct genetic groups with significant differentiation (Global $F_{ST} = 0.056$):

1. Japanese Broad Range (JBR): Widely distributed across Japanese waters.
2. Northernmost Honshu–Hokkaido (NHH): Restricted to the Pacific side of Hokkaido, northern Honshu, and Mutsu Bay. This group was genetically identified for the first time.
3. Western Sea of Japan (WSJ): Exclusive to the Western Sea of Japan coast.

• Differentiation: WSJ showed high differentiation from JBR ($F_{ST} = 0.061$) and NHH ($F_{ST} = 0.074$). JBR and NHH were less differentiated ($F_{ST} = 0.016$).
• Diversity: Nucleotide diversity ($\pi$) was highest in JBR, followed by NHH and WSJ.
• Linkage Disequilibrium (LD): WSJ showed slower LD decay and higher long-range LD, suggesting potential sub-structuring within this group.

B. Demographic History

• Divergence: All three groups diverged during the Last Glacial Period (LGP).
• Trajectories: Post-divergence, all groups experienced a decline followed by recent rapid expansion.
  • JBR: Smallest decline, most pronounced recent expansion (Highest recent $N_e$).
  • WSJ: Greatest historical shrinkage, least recent expansion (Lowest recent $N_e$).
  • NHH: Intermediate trajectory.
• Tajima's D: All groups showed negative skew, consistent with population expansion, with the strongest signal in JBR.

C. Isolation-by-Distance (IBD)

• Pattern: Significant IBD was observed across Japanese waters and within the JBR group.
• Exception: No IBD was detected within the NHH group. This suggests that gene flow within the NHH range is not restricted by geographic distance, likely due to strong natal homing to the Mutsu Bay spawning aggregation, which homogenizes the genetic composition of the feeding grounds.

D. Genetic Stock Identification (GSI) Optimization

• Background SNPs: Required >500 markers to achieve >90% assignment accuracy between JBR and NHH.
• Outlier SNPs: Only 8 or more outlier SNPs were required to achieve comparable (>90%) accuracy.
• Validation: 18 outlier SNPs were consistently detected across 70% of resampling trials. Functional annotation placed many near protein-coding genes (e.g., cytochrome P450 2J4-like), suggesting adaptive divergence.

4. Significance and Implications

• Refined Stock Management: The study provides the first genomic evidence for a distinct NHH stock linked to the Mutsu Bay spawning aggregation. This validates the biological reality of local management units previously based on tagging data alone.
• Cost-Effective Monitoring: The finding that a small panel of outlier SNPs (8 loci) can achieve high assignment accuracy offers a practical, low-cost strategy for routine GSI in fisheries management, replacing the need for hundreds of neutral markers.
• Demographic Insights: The distinct demographic histories (divergence during LGP, varying expansion rates) confirm that these are independent evolutionary units, reinforcing the need for separate stock assessments.
• Biological Realism: The lack of IBD in the NHH group highlights the power of natal homing in maintaining genetic cohesion over large geographic feeding ranges, a critical factor for understanding mixed-stock fisheries.
• Future Directions: The study underscores the need for broader sampling (e.g., Russian waters, Korean Peninsula) to fully resolve the North Pacific population structure and warns that unsampled peripheral populations may obscure global connectivity patterns.

Conclusion

This research successfully integrates population genomics and demographic modeling to redefine Pacific cod management units in Japanese waters. By identifying three distinct genetic groups and demonstrating the efficacy of outlier SNP panels for stock identification, the study provides a robust scientific foundation for sustainable fisheries management and biological stock assessment.